The ribosome’s interactions with mRNA govern its translation rate and the effects of post-transcriptional regulation. Long, structured 5′ untranslated regions (5′ UTRs) are commonly found in bacterial mRNAs, though the physical mechanisms that determine how the ribosome binds these upstream regions remain poorly defined. Here, we systematically investigate the ribosome’s interactions with structured standby sites, upstream of Shine–Dalgarno sequences, and show that these interactions can modulate translation initiation rates by over 100-fold. We find that an mRNA’s translation initiation rate is controlled by the amount of single-stranded surface area, the partial unfolding of RNA structures to minimize the ribosome’s binding free energy penalty, the absence of cooperative binding and the potential for ribosomal sliding. We develop a biophysical model employing thermodynamic first principles and a four-parameter free energy model to accurately predict the ribosome’s translation initiation rates for 136 synthetic 5′ UTRs with large structures, diverse shapes and multiple standby site modules. The model predicts and experiments confirm that the ribosome can readily bind distant standby site modules that support high translation rates, providing a physical mechanism for observed context effects and long-range post-transcriptional regulation.

Riboswitches are shape-changing regulatory RNAs that bind chemicals and regulate gene expression, directly coupling sensing to cellular actuation. However, it remains unclear how their sequence controls the physics of riboswitch switching and activation, particularly when changing the ligand-binding aptamer domain. We report the development of a statistical thermodynamic model that predicts the sequence-structure-function relationship for translation-regulating riboswitches that activate gene expression, characterized inside cells and within cell-free transcription–translation assays. Using the model, we carried out automated computational design of 62 synthetic riboswitches that used six different RNA aptamers to sense diverse chemicals (theophylline, tetramethylrosamine, fluoride, dopamine, thyroxine, 2,4-dinitrotoluene) and activated gene expression by up to 383-fold. The model explains how aptamer structure, ligand affinity, switching free energy and macromolecular crowding collectively control riboswitch activation. Our model-based approach for engineering riboswitches quantitatively confirms several physical mechanisms governing ligand-induced RNA shape-change and enables the development of cell-free and bacterial sensors for diverse applications.

RNA folding plays an important role in controlling protein synthesis as well as other cellular processes. Existing models have focused on how RNA folding energetics control translation initiation rate under equilibrium conditions but have largely ignored the effects of nonequilibrium RNA folding. We introduce a new mechanism, called "ribosome drafting", that explains how a mRNA's folding kinetics and the ribosome's binding rate collectively control its translation initiation rate. During cycles of translation, ribosome drafting emerges whenever successive ribosomes bind to a mRNA faster than the mRNA can refold, maintaining it in a nonequilibrium state with an acceleration of protein synthesis. Using computational design, time-correlated single photon counting, and expression measurements, we demonstrate that slow-folding and fast-folding RNA structures with equivalent folding energetics can vary protein synthesis rates by 1000-fold. We determine the necessary conditions for ribosome drafting by characterizing mRNAs with rationally designed ribosome binding rates, folding kinetics, and folding energetics, confirming the predictions of a nonequilibrium Markov model of translation. Our results have widespread implications, illustrating how competitive folding and assembly kinetics can shape the gene expression machinery's sequence-structure-function relationship inside cells.

A mRNA's translation rate is controlled by several sequence determinants, including the presence of RNA structures within the N-terminal regions of its coding sequences. However, the physical rules that govern when such mRNA structures will inhibit translation remain unclear. Here, we introduced systematically designed RNA hairpins into the N-terminal coding region of a reporter protein with steadily increasing distances from the start codon, followed by characterization of their mRNA and expression levels in Escherichia coli. We found that the mRNAs’ translation rates were repressed, by up to 530-fold, when mRNA structures overlapped with the ribosome's footprint. In contrast, when the mRNA structure was located outside the ribosome's footprint, translation was repressed by <2-fold. By combining our measurements with biophysical modeling, we determined that the ribosomal footprint extends 13 nucleotides into the N-terminal coding region and, when a mRNA structure overlaps or partially overlaps with the ribosomal footprint, the free energy to unfold only the overlapping structure controlled the extent of translation repression. Overall, our results provide precise quantification of the rules governing translation initiation at N-terminal coding regions, improving the predictive design of post-transcriptional regulatory elements that regulate translation rate.

Centralized facilities for genetic engineering, or "biofoundries", offer the potential to design organisms to address emerging needs in medicine, agriculture, industry, and defense. The field has seen rapid advances in technology, but it is difficult to gauge current capabilities or identify gaps across projects. To this end, our foundry was assessed via a timed "pressure test", in which 3 months were given to build organisms to produce 10 molecules unknown to us in advance. By applying a diversity of new approaches, we produced the desired molecule or a closely related one for six out of 10 targets during the performance period and made advances toward production of the others as well. Specifically, we increased the titers of 1-hexadecanol, pyrrolnitrin, and pacidamycin D, found novel routes to the enediyne warhead underlying powerful antimicrobials, established a cell-free system for monoterpene production, produced an intermediate toward vincristine biosynthesis, and encoded 7802 individually retrievable pathways to 540 bisindoles in a DNA pool. Pathways to tetrahydrofuran and barbamide were designed and constructed, but toxicity or analytical tools inhibited further progress. In sum, we constructed 1.2 Mb DNA, built 215 strains spanning five species ( Saccharomyces cerevisiae, Escherichia coli, Streptomyces albidoflavus, Streptomyces coelicolor, and Streptomyces albovinaceus), established two cell-free systems, and performed 690 assays developed in-house for the molecules.