Physical modeling of ribosomes along messenger RNA: Estimating kinetic parameters from ribosome profiling experiments using a ballistic model
Carole Chevalier,Jérôme Dorignac,Yahaya Ibrahim,Armelle Choquet,Alexandre David,Julie Ripoll,Eric Rivals,Frédéric Geniet,Nils-Ole Walliser,John Palmeri,Andrea Parmeggiani,Jean-Charles Walter
DOI: https://doi.org/10.1371/journal.pcbi.1011522
2023-10-21
PLoS Computational Biology
Abstract:Gene expression is the synthesis of proteins from the information encoded on DNA. One of the two main steps of gene expression is the translation of messenger RNA (mRNA) into polypeptide sequences of amino acids. Here, by taking into account mRNA degradation, we model the motion of ribosomes along mRNA with a ballistic model where particles advance along a filament without excluded volume interactions. Unidirectional models of transport have previously been used to fit the average density of ribosomes obtained by the experimental ribo-sequencing (Ribo-seq) technique in order to obtain the kinetic rates. The degradation rate is not, however, accounted for and experimental data from different experiments are needed to have enough parameters for the fit. Here, we propose an entirely novel experimental setup and theoretical framework consisting in splitting the mRNAs into categories depending on the number of ribosomes from one to four. We solve analytically the ballistic model for a fixed number of ribosomes per mRNA, study the different regimes of degradation, and propose a criterion for the quality of the inverse fit. The proposed method provides a high sensitivity to the mRNA degradation rate. The additional equations coming from using the monosome (single ribosome) and polysome (arbitrary number) ribo-seq profiles enable us to determine all the kinetic rates in terms of the experimentally accessible mRNA degradation rate. Genetic translation, involving the transport of ribosomes, is one of the main steps of genetic expression. The molecular mechanisms that regulate ribosome motion along messenger RNA (mRNA) remain to be elucidated. Ribo-sequencing experiments are now capable of providing access to the gene-dependent density of ribosomes along mRNA. The physical modeling of ribosome motion is carried out using 1D transport models, such as TASEP-like ones, which enables one to obtain the kinetic (initiation and elongation) rates from the Ribo-seq densities. Currently used methods, however, do not allow access to mRNA degradation rates and previous experiments lack input to perform a self-consistent quantitative analysis. Here, we propose a new experimental and modeling protocol to obtain the kinetic rates (including mRNA degradation) with better accuracy. Ribo-seq is performed after splitting the mRNA population into groups that depend on the number of translating ribosomes (from 1 to 4). We thus obtain four Ribo-seq profiles per gene properly normalized without requiring independent experiments. We show that these low ribosome number profiles display a higher sensitivity to mRNA degradation than the full polysome profile and discuss the biological translation regimes of several organisms.
biochemical research methods,mathematical & computational biology