A systems approach to model natural variation in reactive properties of bacterial ribosomes

Background Natural variation in protein output from translation in bacteria and archaea may be an organism-specific property of the ribosome. This paper adopts a systems approach to model the protein output as a measure of specific ribosome reactive properties in a ribosome-mediated translation apparatus. We use the steady-state assumption to define a transition state complex for the ribosome, coupled with mRNA, tRNA, amino acids and reaction factors, as a subsystem that allows a focus on the completed translational output as a measure of specific properties of the ribosome. Results In analogy to the steady-state reaction of an enzyme complex, we propose a steady-state translation complex for mRNA from any gene, and derive a maximum specific translation activity, Ta(max), as a property of the ribosomal reaction complex. Ta(max) has units of a-protein output per time per a-specific mRNA. A related property of the ribosome, T˜a(max⁡), has units of a-protein per time per total RNA with the relationship T˜a(max⁡) = ρa Ta(max), where ρa represents the fraction of total RNA committed to translation output of Pa from gene a message. Ta(max) as a ribosome property is analogous to kcat for a purified enzyme, and T˜a(max⁡) is analogous to enzyme specific activity in a crude extract. Conclusion Analogy to an enzyme reaction complex led us to a ribosome reaction model for measuring specific translation activity of a bacterial ribosome. We propose to use this model to design experimental tests of our hypothesis that specific translation activity is a ribosomal property that is subject to natural variation and natural selection much like Vmax and Km for any specific enzyme.


Results:
In analogy to the steady-state reaction of an enzyme complex, we propose a steady-state translation complex for mRNA from any gene, and derive a maximum specific translation activity, T a(max) , as a property of the ribosomal reaction complex. T a(max) has units of a-protein output per time per a-specific mRNA. A related property of the ribosome, , has units of a-protein per time per total RNA with the relationship = ρ a T a(max) , where ρ a represents the fraction of total RNA committed to translation output of P a from gene a message. T a(max) as a ribosome property is analogous to k cat for a purified enzyme, and is analogous to enzyme specific activity in a crude extract.
Conclusion: Analogy to an enzyme reaction complex led us to a ribosome reaction model for measuring specific translation activity of a bacterial ribosome. We propose to use this model to design experimental tests of our hypothesis that specific translation activity is a ribosomal property that is subject to natural variation and natural selection much like V max and K m for any specific enzyme.

Background
Measures of primary sequence variation in ribosomal RNA, ribosomal protein, and non-ribosomal enzyme protein show that the corresponding gene sequences provide a record of natural variation associated with phylogeny in bacteria and archaea [1,2]. Recent analysis of a select subset of 32 protein-coding genes, mostly associated with transcription and translation machinery, established phy-logenetic relationships among the host genomes of 69 bacterial and archaeal organisms, and the basis for a molecular tree of variation calibrated to geological time [1]. The phylogenetic tree constructed in that study largely resembles trees based upon analyses of 16S rDNA sequences [3].
In addition to the phylogenetic perspective based on sequence variation in ribosomal RNA and ribosomal protein encoding genes, we suggest that this sequence variation underlies differences in the performance of the translational machinery. Modeling variation in the kinetics of translation, the primary energy-requiring process of microbes, should lead to an enhanced ecological perspective of microbial life, much like enzyme kinetics has provided the basis for physiological models of metabolism.
Variation in kinetic properties of enzymes has been associated with phylogenetic variation in closely related bacteria but data are less readily available for the variation in enzyme families. In some few studies, K m and V max have been used as functional parameters of enzymes to infer neutral or adaptive properties associated with allelic variation in natural populations [4,5]. Hereafter, we refer to the kinetic parameters of enzymes as specific reactive properties. These studies suggest that we should expect phylogenetic variation in measures of specific reactive properties, e.g. K m , V max , and k cat , for a family of functionally conserved enzymes. Though a recent, systematic study of the natural variation in specific reactive properties of enzyme families does not appear to be available, a brief survey of literature supports the expectation for natural variation of K m , V max in an enzyme family. For example, the published values for K m of biosynthetic threonine deaminase (EC 4.2.1.16) vary from 0.25 mM in Thiobacillus acidophilus and 3.2 mM in Proteus morganii, to 8.0 mM in Escherichia coli, 13.9 mM in Geobacillus stearothermophilus, and 14.0 mM in Thermus sp. [6][7][8][9][10]. Similarly, the V max value for purified threonine deaminase varies over a range from 2.1 μmol/min/mg in Pseudomonas putida to 210 μmol/min/mg in E. coli and even higher in other organisms. Overall, the reactive properties of the biosynthetic threonine deaminase (4.2.1.16) family of enzymes appear to reflect a natural variation over an approximate range greater than two orders of magnitude. Simply stated, enzymes and enzyme complexes that evolve to carry out the same reactions in a set of organisms, should be expected to catalyze these reactions with sometimes substantial differences in steady-state rates and reaction efficiencies, owing to natural variation in sequence associated with reaction properties. Variation in reaction properties may confer a fitness advantage on organisms and be fixed rapidly in populations, or be selectively neutral and accumulate in populations as a result of genetic drift.
In analogy to the natural variation of conserved reactive properties of proteins, ribosomal RNA sequences undergo natural variation that forms the basis for phylogenetic grouping of all organisms and reflects the physiological adaptations associated with the sequence variation [9].
Regions of rRNA that are most highly conserved are speculated to be associated with the ribosomal reaction properties moreso than the regions which are more free to vary. Ribosomal proteins also reflect sufficiently conserved variation that allows for their use in phylogenetic grouping [1]. Since phylogenetic variation in enzyme protein sequence is associated with variation in measurable, reactive properties of the enzymes, the suggestion seems reasonable that measurable, reactive properties of ribosomes are associated with regions of rRNA and ribosomal protein sequences that undergo natural variation. Thus, reactive properties of ribosomes should also be expected to vary in ways similar to the natural variation observed for enzymes [11].
A small, phylogenetically varied subset of bacteria, grouped by fast or slow growth rates in response to nutrient availability, demonstrated variation in a measurable property of ribosomal activity termed translational power by Dethlefsen and Schmidt [12]. This translational power displayed natural variation over approximately a five-fold range [11]. By inference from those studies, we postulate that specific reactive properties of ribosomes are measurably constant for any specific bacterial phylotype or species in a constant environment, and that each reactive property is subject to natural variation under selective constraints to maintain function.
In the present study we consider the ribosome and associated translational machinery as a co-evolving complex of functionally conserved molecules with measurably varied reactive properties in distinct bacterial populations. Consequently, we expect to find that ribosomes from different bacteria vary in their steady-state output of proteins as a result of selective pressures that have fine-tuned the translational machinery for optimal activity in a given environment. In analogy to a multi-enzyme complex, we model ribosomal catalyzed translation activity as a steady-state process involving the entire translational complex as a transition-state intermediate such that the measurably important output is a functional protein. The translational apparatus complex is, then, a subsystem of the organism and produces the mass ratios of proteins required for maintenance and growth. Treatment of the translational complex as a subsystem operating in a steady-state is the key factor that makes measurement of protein output sufficient to characterize subsystem properties associated with a ribosome. Our model suggests that variation in the reactive properties of ribosomes can be detected by measuring the steady-state output of specific, functional proteins. This model substitutes simpler, systems measures to replace the more complex, step-wise, mechanistic measures used elsewhere, [e.g. [13,14]]. Details of the design strategy to develop this measure follow.

Strategy to measure reactive properties of ribosomes
In order to measure specific reactive properties of ribosomes, the design of a model of the translational apparatus as a bacterial subsystem must take into account all factors that can influence the rate of protein output that establishes gene expression levels and the mass ratio of proteins for steady-state growth. When all such factors can be held constant, then measurement of translational output records a ribosome-specific property. Under any steady state growth condition the properties of transcription and translation machinery are constant. A few definitions, here, are useful to build a model. For this model we define expression level as the steady state mass ratio of protein, P a , from a specific gene compared to the total cell protein, P, i.e. P a /P.
For this study, transcription processes in steps 1 and 2 are held constant, and steps 3-11 form the basis to model the steady-state output activity of the translation apparatus. Fig. 1 and the accompanying description refer to these steps 3-11.

Rationale for the model
In the exponential state of bacterial growth in a closed system, the rate of production of cell mass, dx/dt, is proportional to the initial cell mass, x, with the proportionality represented by the instantaneous growth rate constant, k, as In this steady state, the rate of production of any component of cell mass is commonly and reasonably assumed to be a constant fraction of the total cell mass, and k would be the expected proportionality constant for production of any component of the cell. This is a reasonable assumption and subject to a simple proof that need not be shown here. From this inference, as suggested by Neidhardt et al. [19], k defines the steady-state rate of protein production in Eq. (2).
In the steady state translation output of a genome represented by Eq. (2), the protein output per unit total RNA is where T = the protein translation output of a genome, measured in units of total protein, P, per unit total RNA, R, in a bacterium in steady state growth with instantaneous growth rate, k [12,20]. This relationship approximates a property of average translation activity per unit total RNA for an organism in a specific state of growth. We might also refer to T as the translation completion rate in recognition that not all starts of polypeptide synthesis reach completion [18]. Considering the protein product from a single gene as a regulated, fractional component of cell mass, Eq. (2) implies that the change in the mass of protein, P a , from any gene a in genome A, i.e. a ∈ A, is It follows that a translation output T a , from any gene a per unit total RNA results from translation activity by the subset of the ribosomes translating the message from gene a. Ribosomal RNA constitutes approximately 85% of total RNA in bacteria [21], so total RNA is mostly ribosomal.
Ribosome reaction model, analogy to an enzyme complex Figure 1 Ribosome reaction model, analogy to an enzyme complex. Abbreviations: R 30 , 30s ribosomal subunit; m a , mRNA from gene a; R 50 , 50s ribosomal subunit; aa, amino acid; tRNA aa , tRNA isoaccepting species for amino acid aa; aa~tRNA aa , aminoacylated tRNA; Rx, reaction; and P a , protein from translation of mRNA from gene a.
Let ρ a represent the fraction of total RNA committed to translation output of P a from gene a message, then the amount of R committed to producing P a is ρ a R. Thus ρ a R approximates the fraction of ribosomes translating the mRNA from gene a m a . Since Eq. (3) defines the general relationship of the overall translational output to total RNA in a growing organism or culture, i.e. the translational power, we propose Eq. (5) as a special case of Eq.
(3) to represent the specific translational output from any one specific gene, a.

Ribosome analogy to an enzyme complex
Through simplifying assumptions, treat the 70S ribosome as a multi-subunit enzyme that catalyzes reactions to polymerize amino acids to make a functional protein product, P a . A summary representation of this process is diagrammed as a chemical reaction sequence ( Fig. 1). For any gene, treat the initiation process and assembly of the 70S ribosome as analogous to activation of the ribosome as an enzyme. The formation of the initiation complex involves binding of 30S ribosome (R 30 ) with mRNA from gene a, (m a ), at a Shine-Dalgarno (SD) or other ribosome binding site; f-met~tRNA; the initiation factors IF1, IF2 and IF3; GTP; and subsequent assembly with 50S ribosome (R 50 ) to form the activated ribosome, R 70 m a . In this analogy, reaction cofactors are not rate limiting, m a is an activator, and aminoacylated tRNAs (aa~tRNA aa ) are substrate molecules. The amino acid polymerization reactions that follow utilize the pools of aa~tRNA aa to produce the polypeptides for the functional P a protein. Treat all steps of polypeptide chain elongation reactions as a single transition state intermediate. The overall reaction product is P a , formed in a step involving release from the ribosome-mRNA complex and assembly of the active protein.
We consider the last step as rate-determining. Measurement of the rate of increase of P a reactivity is the measure of the rate of production of functional proteins from translation of gene a message.
Continuing the analogy, the active ribosome complex, R 70 m a , would be expected to vary the rate of polypeptide output in response to the variables of input, especially the pool of aa~tRNA aa , and the codon composition of the gene transcript, m a . The rate of the translation output would be expected to vary as a function of the degree of match of the aa~tRNA aa composition and abundance with the relative codon frequency of the transcript. We posit that a major source of the variation of specific translation activity is the degree of mismatch between the relative codon frequency and the tRNA abundance and composition. The greater the degree of mismatch, the higher the frequency of codon-specific pauses and the greater is the probability of premature termination of polypeptide synthesis, i.e. a higher drop-off frequency. Thus, within a genome, A, the translation activity of its ribosomes would be expected to vary over a definable range up to a maximum translation activity that is a property of the ribosome.
The rate of protein output in this Ribosome Reaction Model is a function of the catalytic capacity of the activated ribosome complex acting upon the pool of aa~tRNA aa as shown by Eq. (6).
We can see that the rate of production of P a is the initial reaction velocity of output from the gene a message, thus From Eqs. (6-7) we see that Eqs. (4)(5)(6) imply that , and from our defi- Activated ribosomes within an organism will therefore be expected to display a property of maximum specific translation activity, T a(max) , under optimal conditions of reactant availability. By inference, it may be possible to measure the catalytic efficiency for use of any specific aa~tRNA aa in translation of m a by measure of k cat /K m , i.e. T a(max) /K m . We might expect, then, that the translation properties of ribosomes in a R 70 m a complex are constant, are organism specific, and therefore are genome-specific. k cat for the complex would be expected to vary according to m a properties that influence the expression level.

Single gene translation activity, factors and determinants
Modifying Eq. (5) to introduce a term for the expression level of gene a, define the specific translation activity for gene a as where λ a = P a /P is the expression level. The λ a term shows the fractional contribution of functional protein from any gene, a, to the total protein content of the organism in any steady state of a growth cycle. Thus, λ a is the mass ratio of protein P a compared to the total cellular protein. The dP a / dt varies among the genes of a genome such that a gene or gene set that is highly expressed has a higher production rate for P a than genes expressed less frequently. Therefore, multiple expression states must be expected to exist according to the functional role of a gene or set of genes in the life cycle of an organism. The consequence of these multiple expression states is that the steady state translation of transcripts from any low expression gene contributes fewer proteins than a high expression gene to the total cellular protein.
The half-life of proteins that comprise P would be expected to vary over some range and that variation will constitute some frequency distribution. Considering the variation in stability of any protein P a , a more accurate representation of λ a = P a /P may be In theory, the measurable total protein mass, P, will be experimentally indistinguishable from ∑P a even if the t 1/2 of one or several P a is extremely short. Therefore, we expect that P ≈ ∑P a . However, in the case of an extremely short half-life for any P a , the expression of P a would appear to be low and lead to an erroneous conclusion that the mass ratio λ a is low. Though a condition of high instability of P a may be rare, calculation of the initial mass ratio from Eq. (11) can be accomplished by substituting the initial protein output P a' for P a , where c a is a protein stability constant, and t is the time increment between completion of protein synthesis and sampling. These relationships derive from treating the decay of the functional protein P a as P a = P a' exp-(c a t). Note that initially, i.e. when t = 0, the measured value and initial value for P a are equal, and when c a is small, i.e. stable protein, P a ≈ P a' . Thus P a' is an estimate of the original protein mass when t, c a , and P a are known.
Since not all translation starts, , reach completion for a gene message owing to premature termination, ψ a in Eq.
(13) represents a translation completion efficiency for each gene.
Then δ, a so-called "drop-off frequency" sometimes used to describe translation efficiency, may be defined as We may now also define specific translation activity for a single gene a as and when considering the instability of proteins from any gene a,

Ribosome loading and transcript stability
We will next treat, specifically, the initiation step, also called the ribosome loading step. Thus far we have considered factors affecting the specific translation activity, T a . The overall translation rate, however, also depends upon T kP a s a the frequency of ribosome loading onto mRNA, the stability of that mRNA, and the probability and frequency of the gene transcription. The first step of the reaction sequence ( Fig. 1) represents ribosomal loading onto transcript m a from gene a ∈ A, described by rate equation Eq. This association constant can be thought of as a loading factor that is a ratio of bound to free mRNA from gene a, and can be related to the probability of message binding to ribosome. The amount of m a undergoing translation compared to the total m a is the relative frequency of loading and approximates the probability of ribosomal loading, i.e.
where f l is the relative loading frequency and π a is the probability of message from gene a binding to ribosomes for translation. The specific translation rate, T a , weighted by the fraction of total message loaded, π a ([m a ] + [R 70 m a ]) = π a m a(total) , gives an overall expectation for translation rate of gene a message to make the functional protein, P a (Eq. (21)).
T rate = π a m a(total) T a (21) Factors and growth conditions that affect transcription initiation are too numerous and variable to treat in the context of overall translation rate. Therefore, the last step we consider in this treatment of the overall specific trans-lation rate, T rate , is transcript stability. The higher the stability of m a , the greater is the output of P a per time, and conversely, the lower the stability of m a , the lower the overall translation rate. Analogous to treatment of protein stability in Eq. (12), m a represents the amount of mRNA from gene a that is available to bind a ribosome. It is easy to see that instability of m a reduces the amount of transcript loaded for translation and reduces the overall translation rate. If the stability of mRNA, i.e. half-life, for gene a is known, it should be possible to use Eq. (22) to estimate the actual m a available for loading, where ω a is a constant of mRNA stability for gene a, and m a' is the initial amount of message.

Maximum translation rate
Ideal conditions that include optimal mRNA stability and protein stability would be expected to maximize production of proteins from a single, protein-coding gene. Every translation start results in a completed protein under optimal conditions when ψ a = 1.0. Thus, an expected maximum, steady state specific translation activity by any ribosome to make any P a may be defined as a maximum specific translation activity To approximate a T max for any gene a under non-ideal laboratory or natural conditions we may use Eq. (18) to correct for protein instabilities and translation efficiency to establish that is the maximum speed at which ribosomes can translate any gene in one organism at a steady state of growth. Such optimal conditions imply existence of a K a' , a maximum efficiency of ribosome binding to mRNA such that where and this also implies a corresponding π a' for the K a' . It follows, then, from Eqs. (24) and (25) that a maximum overall translation output rate is given by T rate (max) = π a' m a(total) T max .
( 2 6  It is now possible to account for the multiple effect(s) of growth conditions on T a , including factors that affect translation efficiency, e.g. the drop-off rate and completion rate.
As defined in Eq. (24), the maximum specific translation output activity for a ribosome, T max , defines a genomespecific property as suggested by Dethlefsen and Schmidt [11]. As a general property, T max defines a capacity for production of protein, and that capacity will variably express in relation to the transcript composition from a particular gene. A convenient and practical way to look at specific translation output capacity for a single gene tied to the steady state growth rate of the organism is to combine the factors that cause variability from gene to gene within a genome into the expression for maximum specific translation output activity per unit total RNA in translation of m a . From Eq. (23) for expression of protein from a specific gene, a in genome A, we can see that and in this form it is not necessary to know or measure protein starts from gene a, i.e. , as long as the expres- The maximum specific translation output rate of functional, gene a protein can be measured in multiple organisms as a means to compare this ribosome-specific property among the set of organisms by using Eq. (27). As long as the functional proteins measured derive from an orthologous set of genes, the actual measurement of T a(max) should reveal the variation in the specific translation rate property of ribosomes from one organism with respect to that property in every other organism. Recognizing that the approximation is always less than but fractionally related to T a(max) by ρ a to a specific gene a, we suggest that this approximation will provide a basis to detect, quantify and compare variation in the specific translation rate property of a ribosome. We suggest that measures the specific translation activity of ribosomes analogous to enzyme specific activity measurements in crude extracts, e.g. V max /P, and that T a(max) is analogous to the k cat of a purified enzyme, e.g. V max /P a . We propose that Eq. (27) is an idealized mathematical model of the organism-specific translation output rate as a specific, functional property of the ribosome and that Eq. (28) is a readily measurable, fractional approximation.

Discussion and conclusion
Analogy to the reaction of an enzyme complex formed the basis to derive Eq. (27) as a simple model for specific translation activity as a specific reaction property of a ribosome. The measure of T a(max) is made by the measure of kP/ρ a R and λ a /ψ a for the species-specific reaction property of a ribosome. Factors that make the measure speciesspecific are the expression level of the protein measured, captured in λ a as the relative ratio of P a to the total protein P, and the translation completion efficiency, ψ a . This model of translation activity of a ribosome provides a way to compare ribosome function in one bacterium to another, in bacteria that differ in their growth responses to changing environments or growth conditions. A key postulate in deriving this model is the steady-state assumption that for the time interval of constant growth conditions, the level of biosynthetic intermediates is constant and sufficient to support growth at a constant rate. The steady-state assumption, in turn, permits treatment of the ribosomal reaction complex (Fig. 1) as a cellular reaction subsystem, and to measure its output as an indicator of the subsystem function and therefore as a property of the subsystem. These two simplifications make it possible to treat ribosomal protein synthesis as a "black box" subsystem comprised by many component activities, and to use only the output of functional protein to measure this subsystem property. Any gene product can be used to measure a functional property of the protein synthesis subsystem. For example, the specific activity of a biosynthetic pathway enzyme can be used to measure the reaction property of specific protein output per unit of translating RNA associated with a single gene in a crude codon usage previously reported could be used. As an example, the codon adaptation index of Sharp and Li [21] has been shown to correlate directly with the expressionlevel of protein-coding genes in Escherichia coli and may be used as a value for ψ a in this context [16]. We also suggest that a Pearson cross correlation measure of similarity between the relative codon frequency vector of a test gene, a, and a reference vector, b, of a high-expression gene set, i.e. r γ (a|b), may be well suited and sufficient for this measure.
Measurement of the activity of the functional protein output from a gene permits use of any protein-coding gene or set of genes to calibrate specific translation activity as the measurable reactive property or parameter of ribosomes, using T a(max) when ρ a is known or using when ρ a is unknown. It is reasonable to expect that measurement of as the specific translation activity for m a is fractionally related to T a(max) just as the specific activity of an enzyme measured in a crude cell extract, V max /P, is fractionally related to k cat = V max /P a for the purified enzyme.
So, it is appropriate to treat the measurement of the specific translation activity, , for bacterial ribosomes under defined growth conditions as one would treat measuring the specific activity, V max /P, for an enzyme under defined reaction conditions in a crude cell extract.
The ribosome reaction models for this study expect values for T a(max) and to vary according to the gene-specific mRNA undergoing translation and the specific reactive properties of the ribosomes within an organism.
Owing to the ρ a relationship between T a(max) and , what is true for is also true for T a(max) . Since it is likely that a value for ρ a will frequently be unknown, this discussion will treat as a proxy for T a(max) . A basic assumption of our model is that the reactive properties of ribosomes are constant within an organism under a defined set of growth conditions. The model pioneered by Kurland [20,18] and uniquely applied by Dethlefsen and Schmidt [11] is based upon treating the generalized protein production per unit time per total RNA (kP/R) as constant for a specific growth condition of an organism. If the apparent difference in translational power between ribosomes of a bacterium A compared to B can be attributed to differences in the ribosomes [11], we suggest that may improve the accuracy of measurement of the ribosomal differences and offer a way to study the natural variation in the properties associated with ribosomal reactivity. Our proposed improvement in accuracy is based upon the measurement of translational output, or power, from the mRNA from an individual gene in our ribosome reaction model. This model predicts that the natural variation in specific translation activity for mRNA from any genes within a genome should show a central tendency that is characteristic for that organism, such that for specific genes a 1 through a n in the genome of an organism A. This means, of course, that while the specific translation activity for different genes will be expected to vary within a genome, an average specific translation activity will represent the ribosome property for the host organism. According to Eq.

… …
izes the natural variation in this ribosomal reactive property.
In addition to using the reaction rate measurement T a (max) and to characterize ribosomes, it should be possible to measure a K m for the ribosome reaction model by using in vitro translation with rate-limiting concentrations of an aminoacylated tRNA as conceptually depicted in Fig. 1. Measurement of a K m for a charged tRNA would then allow characterization of a ribosome reaction efficiency as T a(max) /K m or /K m . A further advantage of this model is that it corrects for various specific errors, e.g. translational frameshift and other output errors, by basing the ribosome reactivity upon production of functional enzymes and measuring those functions directly.
A useful model should reflect the behavior of a system under naturally occurring as well as experimental conditions. One observation is that ribosomes are rarely saturated under a specific growth condition [11]. Achievement of a steady-state proposed by this model requires only a constant balance of reactants to achieve a stable output over time. This steady-state model applies for all degrees of saturation. This Ribosome Reaction model, , behaves such that high speed of translation is associated with high values for the growth rate, k, and expression level λ a , and low values for the translation efficiency ψ a . Conversely, slower translation speed is associated with slower growth, k, lower expression levels, λ a , and higher translation efficiency, ψ a . In other words, organisms that are rapid responders to change of environmental conditions will tend to have sloppier translation than slow responders [12], and our Ribosome Reaction model supports these translation relationships. Development of this model provides a potential boost to successful inference of physiology from genomic sequence. If the relationships of Eq. (27) and the proportionate measurement of Eq. (28) hold in analysis of codon usage and gene expression level, it may become possible to associate specific ribosomal characteristics to bacterial lifestyle and ecological niche, i.e. to establish that ribosomes are customized for their host environment. For example, using a sufficient number of measurements of threonine deaminase specific activity in steadystate growth of a phylogenetically varied set of bacteria, along with codon usage patterns of these genes, the calcu-lated value of T a(max) may distribute into expected ranges that are consistent with specific lifestyles from mostly dormant to rapid responders. This model provides an accessible tool and potential to measure a physiological bias of ribosomes to phylogeny and the environment in which they operate.

Abbreviations
Except in Fig. 1 and noted in the Fig. 1 legend, only standard abbreviations were used in this paper.