A systems approach to model natural variation in reactive properties of bacterial ribosomes
© Jackson et al; licensee BioMed Central Ltd. 2008
Received: 10 March 2008
Accepted: 13 July 2008
Published: 13 July 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 18.104.22.168) 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 (22.214.171.124) 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].
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).
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 .
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
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.
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.
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 .
Ribosome loading and transcript stability
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
where ω a is a constant of mRNA stability for gene a, and ma'is the initial amount of message.
Maximum translation rate
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.
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.
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.
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.
- 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-44PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Woese CR: Interpreting the universal phylogenetic tree. Proc Natl Acad Sci USA. 2000, 97: 8392-8396. 10.1073/pnas.97.15.8392PubMed CentralView ArticlePubMedGoogle Scholar
- Demetrius L: Selective neutrality and enzyme kinetics. J Mol Evol. 1997, 45: 370-377. 10.1007/PL00006242View ArticlePubMedGoogle Scholar
- 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;2View ArticlePubMedGoogle Scholar
- Eisenstein E: Cloning, expression, purification, and characterization of biosynthetic threonine deaminase from Escherichia coli. J Biol Chem. 1991, 266: 5801-5807.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- Muramatsu N, Nosoh Y: Some catalytic and molecular properties of threonine deaminase from Bacillus stearothermophilus. J Biochem. 1996, 80 (3): 485-490.Google Scholar
- Proteau G, Silver M: Effects of amino acids on Thiobacillus acidophilus. II. Threonine deaminase. Can J Microbiol. 1980, 26: 385-388.View ArticlePubMedGoogle Scholar
- Yoshida H, Hanada K, Ohsawa H, Kumagai H, Yamada H: Biosynthetic threonine deaminase from Proteus morganii. Agric Biol Chem. 1982, 46: 1035-1042.View ArticleGoogle Scholar
- 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-06PubMed CentralView ArticlePubMedGoogle Scholar
- 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-3PubMed CentralView ArticlePubMedGoogle Scholar
- Drew DA: A mathematical model for prokaryotic protein synthesis. Bull Math Biol. 2001, 63: 329-351. 10.1006/bulm.2000.0225View ArticlePubMedGoogle Scholar
- 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-4View ArticlePubMedGoogle Scholar
- 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.0428View ArticlePubMedGoogle Scholar
- 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.xView ArticlePubMedGoogle Scholar
- 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.2000PubMed CentralView ArticlePubMedGoogle Scholar
- 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 PressGoogle Scholar
- Neidhardt FC, Ingraham JL, Schaechter M: Physiology of the Bacterial Cell: A Molecular Approach. 1990, Sinauer AssociatesGoogle Scholar
- Kjeldgaard NO, Kurland CG: The distribution of soluble and ribosomal RNA as a function of growth rate. J Molec Biol. 1963, 6: 341-348.View ArticleGoogle Scholar
- 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-30PubMed CentralView ArticlePubMedGoogle Scholar
- 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.1281PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.