- Research article
- Open Access
A systems approach to model natural variation in reactive properties of bacterial ribosomes
BMC Systems Biology volume 2, Article number: 62 (2008)
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.
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, , has units of a-protein per time per total RNA with the relationship = ρ a Ta(max), where ρ a represents the fraction of total RNA committed to translation output of P a from gene a message. Ta(max)as a ribosome property is analogous to kcat for a purified enzyme, and is analogous to enzyme specific activity in a crude extract.
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.
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 phylogenetic relationships among the host genomes of 69 bacterial and archaeal organisms, and the basis for a molecular tree of variation calibrated to geological time . The phylogenetic tree constructed in that study largely resembles trees based upon analyses of 16S rDNA sequences .
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 Vmax 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 , Vmax, and kcat, 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 , Vmax in an enzyme family. For example, the published values for Km of biosynthetic threonine deaminase (EC 188.8.131.52) 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–10]. Similarly, the Vmax 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 (184.108.40.206) 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 . 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 . 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 .
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 . This translational power displayed natural variation over approximately a five-fold range . 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. Pa/P.
The instantaneous mass ratio or expression level of a protein from any gene in a bacterium under steady state growth conditions is determined by probabilities affecting the following steps: 1, transcription initiation (promoter strength, promoter configuration, transcription factors; and metabolites); 2, transcription termination; 3, mRNA stability; 4, ribosome binding site (e.g. Shine-Dalgarno sequence) affinity for ribosome; 5, ribosome loading at the Shine-Dalgarno (SD) or other sequences; 6, relative codon frequency pattern, i.e. codon usage; 6, translation initiation (start codon; initiation factors); 7, protein elongation (GTP, aminoacylated tRNA, elongation factors), and termination; 8, frequency of polypeptide completion; 9, protein folding; 10, protein stability; and 11, metabolite balance regulating enzyme activity and synthesis [15–18].
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
dx/dt = kx.
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. , 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 . 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
dP a /dt = kP a .
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 , so total RNA is mostly ribosomal. 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 (R30) 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 (R50) to form the activated ribosome, R70m a . In this analogy, reaction cofactors are not rate limiting, m a is an activator, and aminoacylated tRNAs (aa~tRNAaa) are substrate molecules. The amino acid polymerization reactions that follow utilize the pools of aa~tRNAaa 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, R70m a , would be expected to vary the rate of polypeptide output in response to the variables of input, especially the pool of aa~tRNAaa, 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~tRNAaa 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~tRNAaa as shown by Eq. (6).
dP a /dt = k3[R70m a ]
We can see that the rate of production of P a is the initial reaction velocity of output from the gene a message, thus
dP a /dt = v a .
From Eqs. (6–7) we see that
k3 = v a /[R70m a ].
Eqs. (4–6) imply that , and from our definition of ρ a R as the amount of R involved with translation of mRNA from gene a, we infer that ρ a R ≈ [R70m a ] and thus k3 ≈ T a , the specific translation activity. The units for k3 are units of P a formed per unit time per unit R involved in translating message m a from gene a. These units are the same as the units for T a , the specific translation activity in Eq. (5). Under conditions where [R70m a ] is maximum, all mRNA from gene a is loaded onto ribosomes, and the maximum output rate is achieved, then v a becomes Va(max). For conditions producing Va(max), k3 becomes analogous to kcat for an enzyme reaction, as in Eq. (9).
kcat = Va(max)/[R70m a ].
Activated ribosomes within an organism will therefore be expected to display a property of maximum specific translation activity, Ta(max), under optimal conditions of reactant availability. By inference, it may be possible to measure the catalytic efficiency for use of any specific aa~tRNAaa in translation of m a by measure of kcat/Km, i.e. Ta(max)/Km. We might expect, then, that the translation properties of ribosomes in a R70m a complex are constant, are organism specific, and therefore are genome-specific. kcat 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 t1/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 Pa'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 = Pa'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 ≈ Pa'. Thus Pa'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
For unstable P a , i.e. when the protein stability constant c a is large, then it would be appropriate to substitute Pa'for P a from Eq. (13) and to introduce an initial to substitute for to get
The rationale for Eq. (16) is analogous to that for Eq. (12), and is an estimate of the original protein mass from completion of all translation starts when t, c a , and are known. Generally, it would be appropriate to use the efficiency factor of Eq.(13) and reserve Eq. (15) for probably rare instances of a highly unstable P a .
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 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. (19).
In Eq. (19), k R is a third-order rate constant with units of (t-1)(conc-1)(conc-1). The conditions of growth result in [m a ] << [R30] and [m a ] << [R50]. These conditions suggest that the association reaction may display first-order behavior. Let k+1 be a pseudo first-order rate constant, where k+1 = k R [R30] [R50], and substitute into Eq. (19) to obtain
We see then that Eq. (19a) shows the expected pseudo first-order behavior of m a binding to ribosomes. At any steady state of growth, ribosomal activation can be represented by an association constant, K a , for the message m a from gene a, as
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 ] + [R70m a ]) = π a ma(total), gives an overall expectation for translation rate of gene a message to make the functional protein, P a (Eq. (21)).
T rate = π a ma(total)T a
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 translation rate, Trate, 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 ma'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 Tmax 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 Ka', a maximum efficiency of ribosome binding to mRNA such that
where and this also implies a corresponding πa'for the Ka'. It follows, then, from Eqs. (24) and (25) that a maximum overall translation output rate is given by
Trate (max) = πa'ma(total)Tmax.
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, Tmax, defines a genome-specific property as suggested by Dethlefsen and Schmidt . As a general property, Tmax 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 expression level, λ a , and the translation completion efficiency, ψ a , are known. In the section on Discussion and Conclusions we suggest a way to estimate a value for ψ a . When ρ a is known, i.e. the approximate fraction of ribosomes translating m a , the value of Ta(max)will be highest because it is a measure of the protein expression per unit of total RNA translating m a . When ρ a is not known, a fractional approximation of Ta(max)is
which has units of units a-protein per time per unit total RNA. It is clear that <Ta(max)because the approximation of Eq. (28) compares the gene a-protein output rate to the total translational capacity as opposed to the fraction of capacity, ρ a , committed to m a translation in Eq. (27). Thus, ρ a relates to Ta(max), by = ρ a Ta(max).
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 Ta(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 Ta(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. Vmax/P, and that Ta(max)is analogous to the kcat of a purified enzyme, e.g. Vmax/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 Ta(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 species-specific 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 extract of a pure culture of a bacterium. If the kcat of the enzyme is known, i.e. reaction units per mg enzyme, the expression level, λ a , is the mass ratio of the amount of enzyme protein, measured by activity, to the total extract protein. The translation efficiency correlated to ψ a can be calibrated to a scale of codon usage relative to the highest expression set of genes. The definition of ψ a requires values for the codon usage measure to vary over the range 0 ≤ ψ a ≤ 1.0, and several possible calculations of bias in codon usage previously reported could be used. As an example, the codon adaptation index of Sharp and Li  has been shown to correlate directly with the expression-level of protein-coding genes in Escherichia coli and may be used as a value for ψ a in this context . 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 Ta(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 Ta(max)just as the specific activity of an enzyme measured in a crude cell extract, Vmax/P, is fractionally related to kcat = Vmax/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, Vmax/P, for an enzyme under defined reaction conditions in a crude cell extract.
The ribosome reaction models for this study expect values for Ta(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 Ta(max)and , what is true for is also true for Ta(max). Since it is likely that a value for ρ a will frequently be unknown, this discussion will treat as a proxy for Ta(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  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 , 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 a1 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 s(28), the steady state condition under which the relationships of Eq. (29) would hold requires that λ a /ψ a be approximately constant for all a i . This means, then, that if specific translation activity is an approximately constant property of the ribosomes of an organism, the ratio of the expression level to codon usage bias must also be approximately constant. This inference forms the basis of postulates that are experimentally testable and should reveal the natural variation in within an organism.
Measurement of should also provide the basis to compare and assess the natural variation in ribosome reactivity using the same gene or set of genes in all organisms. For that purpose, crude extracts of steady-state growth cultures will yield values for each defined growth condition to compare with for a different organism under a different growth condition. Expect to observe, measuring the same protein for n different bacterial genera, that comprise a skewed frequency distribution that characterizes the natural variation in this ribosomal reactive property.
In addition to using the reaction rate measurement Ta(max)and to characterize ribosomes, it should be possible to measure a Km 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 Km for a charged tRNA would then allow characterization of a ribosome reaction efficiency as Ta(max)/Km or /Km. 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 . 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 , 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 steady-state growth of a phylogenetically varied set of bacteria, along with codon usage patterns of these genes, the calculated value of Ta(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.
- Except in Fig. 1 and noted in the Fig. 1 legend:
only standard abbreviations were used in this paper.
Battistuzzi FU, Feijao A, Hedges SB: A genomic timescale of prokaryote evolution: insights into the origin of methanogenesis, phototrophy, and the colonization of land. BMC Evol Biol. 2004, 4: 44- 10.1186/1471-2148-4-44
Hansmann S, Martin W: Phylogeny of 33 ribosomal and six other proteins encoded in an ancient gene cluster that is conserved across prokaryotic genomes: influence of excluding poorly alignable sites from analysis. Int J Syst Evol Microbiol. 2000, 50 Pt 4: 1655-1663.
Woese CR: Interpreting the universal phylogenetic tree. Proc Natl Acad Sci USA. 2000, 97: 8392-8396. 10.1073/pnas.97.15.8392
Demetrius L: Selective neutrality and enzyme kinetics. J Mol Evol. 1997, 45: 370-377. 10.1007/PL00006242
Wang I-N, Dykhuizen DE: Variation of enzyme activities at a branched pathway involved in the utilization of gluconate in Escherichia coli. Evolution. 2001, 55: 897-908. 10.1554/0014-3820(2001)055[0897:VOEAAA]2.0.CO;2
Eisenstein E: Cloning, expression, purification, and characterization of biosynthetic threonine deaminase from Escherichia coli. J Biol Chem. 1991, 266: 5801-5807.
Higa EH, Ramaley RF: Purification and properties of threonine deaminase from the X-1 isolate of the genus Thermus. J Bacteriol. 1973, 114: 556-562.
Muramatsu N, Nosoh Y: Some catalytic and molecular properties of threonine deaminase from Bacillus stearothermophilus. J Biochem. 1996, 80 (3): 485-490.
Proteau G, Silver M: Effects of amino acids on Thiobacillus acidophilus. II. Threonine deaminase. Can J Microbiol. 1980, 26: 385-388.
Yoshida H, Hanada K, Ohsawa H, Kumagai H, Yamada H: Biosynthetic threonine deaminase from Proteus morganii. Agric Biol Chem. 1982, 46: 1035-1042.
Dethlefsen L, Schmidt TM: Performance of the translational apparatus varies with the ecological strategies of bacteria. J Bacteriol. 2007, 189: 3237-3245. 10.1128/JB.01686-06
Dethlefsen L, Schmidt TM: Differences in codon bias cannot explain differences in translational power among microbes. BMC Bioinformatics. 2005, 6: 3- 10.1186/1471-2105-6-3
Drew DA: A mathematical model for prokaryotic protein synthesis. Bull Math Biol. 2001, 63: 329-351. 10.1006/bulm.2000.0225
Heyd A, Drew DA: A mathematical model for elongation of a peptide chain. Bull Math Biol. 2003, 65: 1095-1109. 10.1016/S0092-8240(03)00076-4
Dong H, Nilsson L, Kurland CG: Co-variation of tRNA abundance and codon usage in Escherichia coli at different growth rates. J Mol Biol. 1996, 260: 649-663. 10.1006/jmbi.1996.0428
Karlin S, Mrazek J, Campbell AM: Codon usages in different gene classes of the Escherichia coli genome. Mol Microbiol. 1998, 29: 1341-1355. 10.1046/j.1365-2958.1998.01008.x
Karlin S, Mrazek J: Predicted highly expressed genes of diverse prokaryotic genomes. J Bacteriol. 2000, 182: 5238-5250. 10.1128/JB.182.18.5238-5250.2000
Kurland CG, Hughes D, Ehrenberg M: Limitations of translational accuracy. Escherichia coli and Salmonella Cellular and Molecular Biology. Edited by: Neidhardt F, Curtiss III R, Ingraham JL, Lin ECC, Low KB, Magasanik B, Reznikoff WS, Riley M, Schaechter M, Umbarger HE. 1996, 1: 979-1004. Washington, D. C.: A SM Press
Neidhardt FC, Ingraham JL, Schaechter M: Physiology of the Bacterial Cell: A Molecular Approach. 1990, Sinauer Associates
Kjeldgaard NO, Kurland CG: The distribution of soluble and ribosomal RNA as a function of growth rate. J Molec Biol. 1963, 6: 341-348.
Karpinets TV, Greenwood DJ, Sams CE, Ammons JT: RNA:protein ratio of the unicellular organism as a characteristic of phosphorus and nitrogen stoichiometry and of the cellular requirement of ribosomes for protein synthesis. BMC Biology. 2006, 4: 30- 10.1186/1741-7007-4-30
Sharp PM, Li W-H: The codon adaptation index – a measure of directional synonymous codon usage bias and its potential applications. Nucleic Acids Res. 1987, 15: 1281-1295. 10.1093/nar/15.3.1281
The authors are indebted to L. Dethlefsen for stimulating discussions and extensive encouragement in early stages of this work. The website databases BRENDA http://www.brenda-enzymes.info/ and BioCyc http://www.biocyc.org/ were used as initial information sources for comparative enzyme properties. This study was supported by NSF grants MCB-0135880 to TMS and JHJ, IBN-0421900 to TMS, and DBI-0139577 to JHJ.
JHJ, TMS and PAH fully shared in all aspects of this study to develop a quantitative model of bacterial translation to measure natural variation in specific protein output by ribosomes. JHJ led in the model development, TMS contributed specific expertise in ribosomal properties, and PAH brought perspectives on enzyme reaction models into the study. All authors shared in all phases of evolution of the concepts, approaches, and the final construct of the model, as well as preparation of the manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
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.
About this article
Cite this article
Jackson, J.H., Schmidt, T.M. & Herring, P.A. A systems approach to model natural variation in reactive properties of bacterial ribosomes. BMC Syst Biol 2, 62 (2008). https://doi.org/10.1186/1752-0509-2-62
- Codon Usage
- Natural Variation
- Protein Output
- Translation Efficiency
- Translation Activity