Downstream mRNA Secondary Structure, Not Codon Elongation Supply, Coordinates Co-Translational Protein Folding Across the Human Ribosomal Exit Tunnel
- 1School of Interdisciplinary Engineering & Sciences (SINES), National University of Sciences & Technology (NUST), Sector H-12, Islamabad, Pakistan
Abstract
How ribosomes pace translation to assist nascent protein folding remains an open question in molecular biology. While synonymous codon selection is widely hypothesized to regulate elongation rates to facilitate domain organization, distinguishing genuine translational kinetics from baseline amino acid preferences has proven technically difficult. Here, we analyze a non-redundant cohort of 1,270 high-resolution human crystal structures (410,151 residues) mapped to their native mRNA transcripts. When we mathematically isolate synonymous codon choices from amino acid identity using orthogonal linear projection, standard codon-supply metrics i.e. the tRNA Adaptation Index (tAI) and the Codon Adaptation Index (CAI), show negligible independent spatial coupling with downstream protein structure. Their uncorrected correlations predominantly reflect local amino acid chemistry rather than physical translation pacing. In contrast, downstream mRNA secondary structure stability (minimum free energy, MFE) displays a subtle but consistent correlation that survives amino acid control. Across 100,000 whole-proteome permutations per pair, this MFE signal centers at an offset of +15 to +16 codons across multiple independent physical properties, including residue packing density (r = -0.0699, Z = -23.64, p < 10-5), crystallographic rigidity (B-factor, r = +0.0614, Z = +13.74, p < 10-5), and solvent burial (SASA, r = +0.0492, Z = +18.98, p < 10-5). This +15 codon offset corresponds directly to the physical dimensions of the eukaryotic 80S ribosome: the path from the peptidyl transferase center to the internal uL4/uL22 constriction neck (~10 amino acids) plus the downstream mRNA helicase entry channel (~5 codons). While the overall effect size is modest, accounting for approximately 0.49% of local packing variance, its spatial specificity and consistency across independent structural metrics suggest that downstream mRNA stability acts as a localized mechanical brake during early chain compaction.
Related
- A common structure in recurrent networks supports neural sequence generation locally and in downstream neurons
- PHATE-Derived Ancestry Coordinates Capture Continuous Population Structure in Human Genomic Data
- Transfer Matrix and Tensor Approaches to Nucleic Acid Secondary Structure Thermodynamics
- Use of FTIR to study secondary structure of texturized plant proteins by high moisture extrusion cooking, a comprehensive review
- Beyond Ramachandran distributions: quantifying secondary structure effects on amino acids and their pairs in the full torsional space
BibTeX
@misc{hassan2026downstream,
title = {{Downstream mRNA Secondary Structure, Not Codon Elongation Supply, Coordinates Co-Translational Protein Folding Across the Human Ribosomal Exit Tunnel}},
author = {M. Saad Ul Hassan and I. Jabeen and Y. S. Kiani},
howpublished = {bioRxiv},
year = {2026},
month = sep,
doi = {10.64898/2026.09.27.753847},
url = {https://www.biorxiv.org/content/10.64898/2026.09.27.753847v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Downstream mRNA Secondary Structure, Not Codon Elongation Supply, Coordinates Co-Translational Protein Folding Across the Human Ribosomal Exit Tunnel
AU - M. Saad Ul Hassan
AU - I. Jabeen
AU - Y. S. Kiani
PY - 2026
DA - 2026/09/27
PB - bioRxiv
DO - 10.64898/2026.09.27.753847
UR - https://www.biorxiv.org/content/10.64898/2026.09.27.753847v1
AB - How ribosomes pace translation to assist nascent protein folding remains an open question in molecular biology. While synonymous codon selection is widely hypothesized to regulate elongation rates to facilitate domain organization, distinguishing genuine translational kinetics from baseline amino acid preferences has proven technically difficult. Here, we analyze a non-redundant cohort of 1,270 high-resolution human crystal structures (410,151 residues) mapped to their native mRNA transcripts. When we mathematically isolate synonymous codon choices from amino acid identity using orthogonal linear projection, standard codon-supply metrics i.e. the tRNA Adaptation Index (tAI) and the Codon Adaptation Index (CAI), show negligible independent spatial coupling with downstream protein structure. Their uncorrected correlations predominantly reflect local amino acid chemistry rather than physical translation pacing. In contrast, downstream mRNA secondary structure stability (minimum free energy, MFE) displays a subtle but consistent correlation that survives amino acid control. Across 100,000 whole-proteome permutations per pair, this MFE signal centers at an offset of +15 to +16 codons across multiple independent physical properties, including residue packing density (r = -0.0699, Z = -23.64, p < 10-5), crystallographic rigidity (B-factor, r = +0.0614, Z = +13.74, p < 10-5), and solvent burial (SASA, r = +0.0492, Z = +18.98, p < 10-5). This +15 codon offset corresponds directly to the physical dimensions of the eukaryotic 80S ribosome: the path from the peptidyl transferase center to the internal uL4/uL22 constriction neck (~10 amino acids) plus the downstream mRNA helicase entry channel (~5 codons). While the overall effect size is modest, accounting for approximately 0.49% of local packing variance, its spatial specificity and consistency across independent structural metrics suggest that downstream mRNA stability acts as a localized mechanical brake during early chain compaction.
ER -