A Chemical Kinetic Basis for Measuring Translation Initiation and Elongation Rates from Ribosome Profiling data

Analysis methods based on simulations and optimization have been previously developed to estimate relative translation rates from next-generation sequencing data. Translation involves molecules and chemical reactions; hence, bioinformatics methods consistent with the laws of chemistry and physics are more likely to produce accurate results. Here, we derive simple equations based on chemical kinetic principles to measure the translation-initiation rate, transcriptome-wide elongation rate, and individual codon translation rates from ribosome profiling experiments. Our methods reproduce the known rates from ribosome profiles generated from detailed simulations of translation. Applying our methods to data from S. cerevisiae and mouse embryonic stem cells we find that the extracted rates reproduce previously reported correlations with various molecular quantities. A codon can exhibit up to 26-fold variability in its translation rate depending upon its context with in a transcript. This broad distribution means that the average translation rate of a codon is not representative of the rate at which most instances of that codon are translated. We also find that mouse embryonic stem cells have a global translation speed that is almost two-fold faster than previously reported. This large variability in translation rates suggests that translational regulation might be used by cells to a greater degree than previously thought. RNA-Seq data corresponding to Ribosome profiling data (GSE75322) of two biological replicates of S288C yeast cells