Environmental perturbations lift the degeneracy of the genetic code to regulate protein levels in bacteria

Abstract: The genetic code underlying protein synthesis is a canonical example of a degenerate biological system. Degeneracies in physical and biological systems can be lifted by external perturbations thus allowing degenerate systems to exhibit a wide range of behaviors. Here we show that the degeneracy of the genetic code is lifted by environmental perturbations to regulate protein levels in living cells. By measuring protein synthesis rates from a synthetic reporter library in Escherichia coli, we find that environmental perturbations, such as reduction of cognate amino acid supply, lift the degeneracy of the genetic code by splitting codon families into a hierarchy of robust and sensitive synonymous codons. Rates of protein synthesis associated with robust codons are up to hundred-fold higher than those associated with sensitive codons under these conditions. We find that the observed hierarchy between synonymous codons is not determined by usual rules associated with tRNA abundance and codon usage. Rather, competition among tRNA isoacceptors for aminoacylation underlies the robustness of protein synthesis. Remarkably, the hierarchy established using the synthetic library also explains the measured robustness of synthesis for endogenous proteins in E. coli. We further found that the same hierarchy is reflected in the fitness cost of synonymous mutations in amino acid biosynthesis genes and in the transcriptional control of sigma factor genes. Our study reveals that the degeneracy of the genetic code can be lifted by environmental perturbations, and it suggests that organisms can exploit degeneracy lifting as a general strategy to adapt protein synthesis to their environment.

\section{Results}
\subsection{Degeneracy lifting upon amino acid limitation}

We considered synonymous codons for seven amino acids: Leu, Arg, Ser, Pro, Ile,
Gln, and Phe. This set of seven amino acids is representative of the degeneracy
of the genetic code, in that it includes six-, four-, three- and two-fold
degenerate codon families. We constructed a library of 29 \textit{yellow
fluorescent protein} (\textit{yfp}) gene variants, each of which had between six
and eight synonymous mutations for one of the seven amino acids (Fig.
\ref{fig1}A). In this library, we designed each \textit{yfp} variant to
characterize the effect of one specific codon on protein synthesis. We expressed
the \textit{yfp} variants constitutively at low gene dosage (2 copies /
chromosome, Fig. \ref{fig1}B) in \textit{E. coli} strains that were auxotrophic
for one or more amino acids. We monitored growth and YFP synthesis in these
strains during amino acid-rich growth as well as during limitation for each of
the seven amino acids (\hyperref[methods]{Methods}).

During amino acid-rich growth, our measurements revealed that protein synthesis
rates were highly similar across \textit{yfp} variants, with less than 1.4-fold
variation within all codon families (Fig. \ref{fig1}D, grey bars). Thus, under
rich conditions, the degeneracy of the genetic code remains intact with respect
to protein synthesis. Strikingly, under amino acid-limited growth, codon
families split into a hierarchy of YFP synthesis rates (Fig. \ref{fig1}C,
\ref{fig1}D). We found that some synonymous codons, such as CTA for leucine,
were highly sensitive to environmental perturbation, causing YFP synthesis rates
to be near zero in response to the limitation of these codons’ cognate amino
acids. Conversely, other synonymous codons, such as CTG for leucine, were more
robust to the same perturbation with synthesis rates of YFP up to 100-fold
higher than the sensitive ones\footnote{We define codons as robust when the
synthesis rate from the corresponding \textit{yfp} variant during cognate amino
acid
limitation is higher than the average synthesis rate within that codon family.
Similarly, we define codons as sensitive when the synthesis rate from the
corresponding \textit{yfp} variant during cognate amino acid limitation is lower
than the
average synthesis rate within that codon family.}. In addition to fluorescence,
this difference in
robustness was reflected in protein levels measured with Western blotting (Fig.
\ref{fig_s1}). Notably, even a single substitution to
a perturbation-sensitive codon in the \textit{yfp} coding sequence resulted in
more than a 2-fold difference in YFP synthesis rate during limitation for the
cognate amino acid, without any effect on synthesis rate during amino acid-rich
growth (Fig. \ref{fig_s2}). Only those codons that were cognate to the limiting amino acid
caused splitting of YFP synthesis rates (Fig. \ref{fig_s3}). Interestingly, the splitting
was more acute for codon families with six-fold degeneracy (Leu, Arg, Ser),
while splitting was weaker for codon families with four-, three- and two-fold
degeneracies (Fig. \ref{fig1}D, first row vs. second row). These results support
the idea that greater degeneracy typically allows systems to exhibit wider range
of responses to environmental perturbations \cite{edelman_degeneracy_2001}. In
subsequent experiments, we focused on the two codon families, leucine and
arginine that displayed the largest range of splitting. These two families
constitute 16\% of codons across the genome of \textit{E. coli}.

\subsection{Intracellular determinants of the hierarchy among synonymous codons}
We sought to identify the intracellular parameters that determine the observed
hierarchy of degeneracy splitting during amino acid limitation. To this end, we
quantified the robustness of synthesis rate to amino acid limitation as the
ratio of YFP synthesis rates between amino acid-limited and amino acid-rich
growth phases. Protein synthesis rate is known to be correlated with codon usage
and tRNA abundance during artificial over-expression of proteins
\cite{robinson_codon_1984, gustafsson_codon_2004}. However, we found that
robustness of YFP synthesis to amino acid limitation was not correlated with
either codon usage or tRNA abundance ($r^2 = 0.08$ and $0.00$, squared Spearman
rank-correlation, Fig. \ref{fig_s4}). We then considered determinants of protein synthesis
that might be important specifically during amino acid limitation. tRNA
isoacceptors are uniformly charged (aminoacylated) at about 80\% under amino
acid-rich conditions \cite{kruger_aminoacylation_1998, sorensen_charging_2001}.
However during perturbations such as amino acid limitation, some tRNA
isoacceptors cognate to the amino acid are almost
fully charged while other isoacceptors in the same family have charged fractions
that are close to zero \cite{dittmar_selective_2005, sorensen_over_2005}. A
theoretical model proposed that such selective charging arises from differences
in the relative supply and demand for charged tRNA isoacceptors
\cite{elf_selective_2003}. While it is unclear how this mechanism could solely
control protein levels, charged tRNA play an essential role as substrates for
the elongation of ribosomes across individual codons \cite{ibba_quality_1999}.
Consequently, we hypothesized that selective charging of tRNA isoacceptors also
underlies the observed splitting in synthesis rates among \textit{yfp} variants.
Consistent with this hypothesis, charged fractions of leucine and arginine tRNA
isoacceptors during limitation of cognate amino acid starvation measured in a
previous work \cite{dittmar_selective_2005} were correlated with the robustness
of synthesis rates from \textit{yfp} variants after accounting for codon-tRNA
assignments ($r^2 = 0.78$, Fig. \ref{fig_s5}).

We experimentally tested whether varying the concentration of charged tRNA could
change the hierarchy of protein synthesis rates initially revealed by amino acid
limitation. To this end, we co-expressed each one of the leucine or arginine
tRNA isoacceptors together with each of the six leucine or arginine variants of
\textit{yfp}, respectively (Fig. \ref{fig2}). Previous work
\cite{sorensen_over_2005} showed that overexpression of a single tRNA
isoacceptor cognate to a limiting amino acid enables it to compete better in the
common charging reaction against other isoacceptors. As a result, charged tRNA
concentration of the overexpressed isoacceptor increases, while charged tRNA
concentrations of the remaining isoacceptors for that amino acid decrease or
remain unchanged \cite{sorensen_over_2005}. We found that \textit{yfp} variants
constructed with perturbation-sensitive codons exhibited higher synthesis rates
upon co-expression of tRNA isoacceptors cognate to those perturbation-sensitive
codons (Fig. \ref{fig2}, A-B , bottom three rows, solid black-outlined squares).
Conversely, \textit{yfp}
variants with perturbation-robust codons exhibited lower protein synthesis rates
upon co-expression of non-cognate tRNA isoacceptors (Fig. \ref{fig2}, A-B, top
three rows, non-outlined squares). These two patterns of changes in YFP
synthesis rate mirror previously measured changes in charged tRNA concentration
upon tRNA co-expression \cite{sorensen_over_2005}, thereby suggesting that the
observed hierarchy in synthesis rates of \textit{yfp} variants are tightly
coupled with the concentrations of cognate charged tRNA isoacceptors during
amino acid limitation. By contrast, tRNA co-expression did not affect synthesis
rates from \textit{yfp} variants in the absence of perturbation, i.e., during
amino acid-rich growth (Fig. \ref{fig2}C). We observed several codon-tRNA pairs
with mismatches at the wobble position but that do not satisfy known
wobble-pairing rules (Table S9), and that showed an increase in YFP synthesis
rate upon co-expression of the tRNA isoacceptor during amino acid limitation
(Fig. \ref{fig2}, A-B , dashed black-outlined squares).

\subsection{A codon robustness index for endogenous proteins}
We investigated whether the hierarchy of synthesis rates measured for the
synthetic \textit{yfp} variants also governs the synthesis of endogenous
proteins of \textit{E.
coli}. We first devised a general parameter, hereafter called the codon
robustness index (CRI), to characterize the robustness of any protein’s
synthesis rate to an environmental perturbation associated with limitation of a
specific amino acid (Fig. \ref{fig3}A). We defined CRI as a product of
codon-specific weights wcodon, and we inferred these weights from the synthesis
robustness of \textit{yfp} variants to limitation for their cognate amino acid
(Fig. \ref{fig3}B). Our formulation of CRI is based on the simplifying
assumption that each codon decreases protein synthesis rate by a factor
$w_{codon}$ that is independent of the codon’s intragenic location, the presence
of other codons in the coding sequence, or the specific cellular role of the
encoded protein. By definition, $w_{codon}$ is unity for codons that are not
cognate to the limiting amino acid, and perturbation-robust codons have a
higher $w_{codon}$ value than perturbation-sensitive codons for the limiting
amino acid.

To test the predictive power of CRI, we selected 92 \textit{E. coli} open
reading frames (ORFs) that span a broad range of leucine CRI values and
functional categories (Fig. \ref{fig_s7}, Table \ref{92orfs}). We expressed the corresponding
proteins constitutively as N-terminus fusions with YFP\footnote{The YFP fusion
partner in the 92 ORF-\textit{yfp} fusions used for testing Leu CRI was encoded
by the CTG variant of \textit{yfp} that has the highest, most robust synthesis
rate during leucine limitation. Similarly, the YFP fusion partner in the 56
ORF-\textit{yfp} fusions used for testing Arg CRI was encoded by the AGA variant
of \textit{yfp} that has the highest, most robust synthesis rate during arginine
limitation.} in an \textit{E. coli}
strain auxotrophic for leucine (Fig. \ref{fig3}C, Inset,). Upon leucine
limitation, we found a strong correlation between the robustness of protein
synthesis rates from the 92 ORF-\textit{yfp} fusions and their leucine CRI
values (Fig. \ref{fig3}C, $r^2=0.61$, $P=10^{-23}$, squared Spearman
rank-correlation). Similarly, arginine CRI was also strongly correlated with
robustness of a library of 56 ORF-\textit{yfp} fusions during arginine
limitation ($r^2=0.59$, $P=10^{-12}$, Fig. \ref{fig_s8}, Table \ref{56orfs}). By contrast, standard
measures of translation efficiency under amino acid-rich conditions such as
codon adaptation index \cite{sharp_codon_1987}, tRNA adaptation index
\cite{dos_reis_solving_2004} or folding energy of the mRNA around the start
codon \cite{kudla_coding-sequence_2009} displayed only a weak correlation with
protein synthesis rate from the ORF-\textit{yfp}
fusions during amino acid-rich growth ($r^2 = 0.10$, 0.08, and 0.02 resp., Fig.
\ref{fig_s9}). We further found that changes in Leu CRI calculated from the \textit{yfp}
data could predict both the effect of tRNA co-expression and that of synonymous
mutations on protein synthesis from \textit{E. coli} ORFs during leucine
limitation (Fig. \ref{fig3}D, Fig. \ref{fig_s10}). Importantly, similar to our results using
\textit{yfp} reporters, neither tRNA co-expression nor synonymous mutations for
\textit{E.coli} ORF-\textit{yfp} fusions had a significant effect on the
synthesis rates from these ORFs during leucine-rich growth in absence of
environmental perturbations (Fig. \ref{fig_s11}). Thus the degeneracy of the genetic code
underlies the levels of endogenous protein production only during response to
environmental perturbations.

\subsection{Consequences of degeneracy lifting for fitness and gene regulation}
Degeneracy splitting in physical systems can be exploited to encode information
related to the environmental context \cite{gershenfeld_bulk_1997,
chuang_experimental_1998}. We asked whether bacteria might similarly exploit the
degeneracy splitting of genetic code during response to amino acid limitation.
Hence we tested whether the expression of amino acid biosynthesis genes that
enable bacteria to adapt to amino acid limitation is affected by the hierarchy
between robust and sensitive codons. We found that mutating codons that are
perturbation-robust to those that are perturbation-sensitive in the
leucine-biosynthesis genes \textit{leuA}, \textit{leuC} and \textit{leuD}, and
the arginine biosynthetic gene \textit{carA} decreased their protein synthesis
rate during cognate amino acid limitation, but not during amino acid-rich growth
(Fig. \ref{fig_s12}). Interestingly, in the case of \textit{leuA} and \textit{carA}, the
same synonymous mutations also resulted in a fitness cost for prototrophic
strains upon downshift from amino acid-rich to amino acid-poor conditions (Fig.
\ref{fig4}A). Thus synonymous mutations can have a
significant fitness cost during an environmental perturbation, which is distinct
from that measured under nutrient-rich conditions in the absence of any
perturbation \cite{kudla_coding-sequence_2009, lind_mutational_2010}. However
swapping codons that are perturbation-robust to those that are
perturbation-sensitive in other biosynthesis genes (see \textit{argA} and
\textit{leuC} in Fig. \ref{fig4}A) did not significantly affect fitness,
suggesting that the hierarchy of robust and sensitive codons might be
selectively utilized by bacteria to regulate genes within a single metabolic
pathway.

Perturbations associated with amino acid limitation in \textit{E. coli} can
result in two distinct outcomes, depending on the environmental conditions: On
one hand, when substrates used in amino acid biosynthesis are still abundant in
the environment, the cell up-regulates corresponding biosynthesis genes to
mitigate the limitation of amino acids and resume growth. On the other hand, in
the absence of substrates for amino acid biosynthesis, \textit{E. coli} can
survive a prolonged period in amino acid-poor environments through a cellular
response mediated by sigma factors \cite{traxler_global_2008,
durfee_transcription_2008}. We found that genes encoding several stress-response
sigma factors (\textit{rpoS}, \textit{rpoE} and \textit{rpoH}) are enriched in
TTA and TTG, the leucine codons that ensure robust protein synthesis during
leucine limitation (Fig. \ref{fig4}B, \textit{top panel}). By contrast, genes
for the
housekeeping sigma factor (\textit{rpoD}) and a few minor sigma factors
(\textit{fecI}, \textit{fliA}) are enriched for CTC and CTT, which are sensitive
to leucine limitation. This contrasting pattern is observed for leucine (but
not for arginine), and is further mirrored by the change in transcript abundance
for sigma factor genes in response to leucine limitation (Fig. \ref{fig4}B,
\textit{bottom panel}). Hence degeneracy splitting in the genetic code might be
exploited
in concert with transcriptional control to regulate protein levels.

\section{Discussion}
In summary, we have found that the degeneracy of the genetic code does not have
a role in regulating protein synthesis during amino acid-rich growth. By
contrast, the splitting of this degeneracy upon reduction in amino acid supply
has a potent effect on protein synthesis that results in up to 100-fold
differences in protein synthesis rates between synonymous gene variants. Such a
large role for synonymous codons in protein synthesis is surprising given that
other post-transcriptional mechanisms such as protein degradation are known to
play a significant role upon amino acid limitation \cite{kuroda_role_2001}. We
identified competition between tRNA isoacceptors for aminoacylation as a key
determinant of the hierarchy of protein synthesis rates during amino acid
limitation. Low concentration of a charged tRNA isoacceptor can cause ribosomes
to selectively pause at its cognate codon\footnote{A recent genome-wide study
found increased ribosome pausing at serine codons during serine-limited growth
of \textit{E. coli}. Interestingly, ribosomes paused significantly only at four
out of the six serine codons, and these four codons are precisely the same ones
that caused YFP synthesis rate to be sensitive to serine limitation in our
experiments (Fig. \ref{fig_s13}).}  and trigger ribosome jamming
\cite{tuller_evolutionarily_2010}, translation-recoding
\cite{zaher_primary_2011}, mRNA cleavage \cite{christensen_rele_2003,
li_cleavage_2008, garza-sanchez_amino_2008} or feedback-transcriptional control
\cite{henkin_regulation_2002, proshkin_cooperation_2010}. It will be
interesting to find the relative contribution of these different molecular
processes to the degeneracy lifting uncovered here\footnote{We measured the change in mRNA levels of different \textit{yfp} variants in response to amino acid limitation. Changes in mRNA levels were correlated with corresponding changes in YFP synthesis rates upon amino acid limitation (Fig. \ref{fig_s14}), but were smaller than expected.}. Here, we have
investigated
the effect of a specific environmental perturbation associated with amino acid
limitation in the bacterium \textit{E. coli}. However, this type of perturbation
plays a crucial role in the lifecycle of other bacteria such as
\textit{Myxococcus xanthus} and \textit{Bacillus subtilis} that undergo
differentiation cued by amino acid limitation \cite{dworkin_recent_1996,
ochi_initiation_1981}. Protein synthesis during such differentiation events
might also be regulated by degeneracy lifting of the genetic code. Moreover,
degeneracy lifting could be important during protein synthesis in eukaryotes,
where clinically-important conditions such as neoplastic transformation and drug
treatment are often accompanied by a reduction in amino acid supply
\cite{ye_gcn2-atf4_2010, zhou_high_2009}. Therefore lifting the degeneracy of
the genetic code might emerge as a general strategy for biological systems to
expand their repertoire of responses to environmental
perturbations.

\begin{materials}
\label{methods}
Summary of key methods are given below. Detailed methods for all experiments and
analyses are included in the \hyperref[appendix]{Appendix}.

\subsection{Bacterial strains}
All strains used in this study were obtained from the E.coli Genetic Stock
Center (CGSC), Yale University. Different auxotrophic strains were used
depending on the amino acid that was limiting in the growth medium (Table S5).
Strains were stored as 20\% glycerol stocks at -80$^\circ$C either in 1ml
cryo-vials or in 96-well plates (3799, Costar). For experiments involving over
25 strains, a temporary 20\% glycerol stock was stored at -20$^\circ$C in
96-well PCR plates.

\subsection{Plasmids}
The \textit{pZ} series of plasmids \cite{lutz_independent_1997} were used for
expression of all genes constructed for this study. General features of the
plasmid backbones are described here. Details on individual gene constructs that
were inserted into these plasmid backbones, including DNA sequences and plasmid
maps, are in Appendix. A low-copy plasmid, \textit{pZS*11} [\textit{SC101*} ori
(3-4 copies/cell), AmpR (\textit{bla} gene) and a constitutive
\textit{P\subsc{L}tetO1}
promoter] was used for expression of all fluorescent reporter genes and their
fusions. The synthetic ribosome binding site (RBS) in the original
\textit{pZS*11} backbone was replaced by a modified T7-based RBS that resulted
in efficient expression of most coding sequences. A medium-copy plasmid,
\textit{pZA32} [\textit{p15A} ori (10-12 copies/cell), ChlR (\textit{cat} gene)
and \textit{P\subsc{L}lacO1} promoter] was used for expression of all tRNA
genes. Strains with
\textit{pZA32} plasmids were grown with 1mM IPTG to ensure constitutive
expression of all
tRNA genes. Standard plasmids \textit{pUC18} and \textit{pUC19} were used as
intermediate cloning
vectors for site-directed
mutagenesis.

\subsection{Gene synthesis and cloning}
A single \textit{yfp} sequence was built de novo (synthesis by Genscript, USA).
All subsequent \textit{yfp} variants were constructed using a site-directed
mutagenesis kit (Stratagene). tRNA genes and \textit{E. coli} ORFs were
amplified from the chromosome of a wild-type \textit{E. coli} strain (MG1655) by
PCR (Details on cloning and genes sequences in Appendix).

\subsection{Amino acid limitation experiments}
Overnight cultures were inoculated from glycerol stocks or fresh colonies and
grown in a MOPS-based rich-defined medium with 800$\mu M$ of 19 amino acids and
10mM serine at 30$^\circ$C with shaking. For experiments involving amino acid
limitation, overnight cultures were diluted 1:1000 into a similar rich-defined
medium as the overnight cultures. However the amino acid, whose limitation was
to be induced, was added at a reduced concentration and supplemented with its
methyl-ester analog (see Table S6 for exact concentrations). Amino acid
methyl-esters are inefficiently metabolized analogs of the corresponding amino
acids and have been previously used for steady growth of \textit{E. coli} under
amino acid limiting conditions \cite{yelverton_function_1994,gallant_role_2004}
(see Figs. S15 and S16 for the effect of methyl-ester on growth and measured
robustness during amino acid limited growth). Slight variations in the initial
concentration of either the limiting amino acid or its methyl-ester only results
in shifting of the transition to a higher or lower cell density
without appreciable changes in growth rate (see Notes S1 and S2). Growth and
fluorescence were quantified using a standard 96-well plate reader integrated
with a robotic system. Further details on growth protocols are given in
Appendix.

\subsection{Analysis of cell density and fluorescence time series}
Matlab R2009 (Mathworks) was used for all analyses unless otherwise mentioned.
All correlations and P-values reported in this work were calculated using the
Matlab command ‘corr’ with the ‘Type’ option set to either ‘Spearman’ or
‘Pearson’ as appropriate. Growth and fluorescence time series were fit with
exponential and linear curves in the amino acid rich and amino acid limited
growth regimes, respectively, and the onset time of amino acid limited growth
was automatically inferred from their intersection. Protein synthesis rate, $S$
was calculated as:
\begin{align}
\text{Protein synthesis rate } S= \frac{1}{\text{Absorbance}}  \times  \frac{d
(\text{Fluorescence})} {d (\text{time})}
\end{align}

First, the above formula was evaluated at the onset time of amino acid limited
growth using the exponential fits for absorbance and fluorescence data in the
amino acid rich growth regime. Next, the same formula was evaluated at the onset
time using the linear fits for absorbance and fluorescence data in the amino
acid limited growth regime. These two values correspond to the protein synthesis
rates reported for the amino acid rich and amino acid limited growth regimes
(such as the data in Fig. \ref{fig1}D). Further details of this analysis are
given in
Appendix.

\subsection{Calculation of CRI}
CRI for a protein coding sequence corresponding to a limiting amino acid was
calculated by multiplying the wi values for codons cognate to the limiting amino
acid in that sequence. $w_i$ values shown in Fig. \ref{fig3}B were calculated
using the
robustness of protein synthesis of the corresponding \textit{yfp} variants
during cognate amino acid limitation (Fig. \ref{fig1}D). Based on our non-cognate
limitation experiment (Fig. S2), the $w_i$ values for all codons other than
those
cognate to the limiting amino acid are set to be equal to 1.
For illustration, we demonstrate the calculation of $w_i$ for the six leucine
codons
below. The exact same procedure was followed for other synonymous codon
families. Taking log\subsc{2} $w_i  \equiv W_i$ for each codon, and log\subsc{2}
(robustness
during amino acid limited growth) $  \equiv SR$ for each \textit{yfp} variant,
\begin{align}
7 \times W_{CTA}+15 \times W_{CTG} &=SR_{yfp,CTA}\\
7 \times W_{CTC}+15 \times W_{CTG} &=SR_{yfp,CTC}\\
22 \times W_{CTG}&=SR_{yfp,CTG}\\
7 \times W_{CTT}+15 \times W_{CTG} &=SR_{yfp,CTT}\\
7 \times W_{TTA}+15 \times W_{CTG} &=SR_{yfp,TTA}\\
7 \times W_{TTG}+15 \times W_{CTG} &=SR_{yfp,TTG}
\end{align}
The multiplicative factors on the LHS in front of $W_i$ correspond to the
frequency of the Leu codon $i$ in the corresponding Leu variant of \textit{yfp}
(see Fig. \ref{fig1}A). The RHS is the measured (log\subsc{2}) robustness of protein
synthesis
from the corresponding \textit{yfp} variant during Leu limitation (see Fig. \ref{fig1}D).
These equations were solved simultaneously to determine the $w_i$ value for each
Leu codon. Revised $w_i$ values based on \textit{yfp} measurements in the
presence of \supsc{GAG}Leu2 tRNA (Fig. \ref{fig2}) were used for calculation of
Leu CRI in the
case of GAGLeu2 tRNA co-expression with \textit{E. coli} ORFs (Fig.
\ref{fig3}D).

\end{materials}

\begin{acknowledgments}
We acknowledge J. Shapiro for work on developing the preliminary growth assays,
and J. Gallant and M. Cashel for useful correspondence on optimizing the growth
assays. A.S. is grateful to C. C. Guet for help with cloning and suggesting the
tRNA co-expression experiments, and J. Moffitt for pointing him to the CP78 E.
coli strain. We thank B. Stern, A. Murray, and V. Denic for critical questions.
We thank K. Dave for editorial advice and assistance, and the members of the
Cluzel lab and FAS Center for Systems Biology for experimental support. We
acknowledge M. Aldana, L. David, D. A. Drummond, J. Elf, M. Kim, L. Marshall, M.
Mueller, E. O’Shea, E. Wallace, J. Weissman, K. Wood, and B. Zid for comments on
earlier versions of the manuscript. We are grateful to the anonymous referees
for critical comments on the tRNA co-expression experiment.
\end{acknowledgments}

Genomics (q-bio.GN); Soft Condensed Matter (cond-mat.soft); Biological Physics (physics.bio-ph)