Ribosome exit tunnel electrostatics

Role of Sec61α2 Translocon in Insulin Biosynthesis

Translocational regulation of proinsulin biosynthesis in pancreatic β-cells is unknown, although several studies have reported an important accessory role for the Translocon-Associated Protein complex to assist preproinsulin delivery into the endoplasmic reticulum via the heterotrimeric Sec61 translocon (comprising α, β, and γ subunits). The actual protein-conducting channel is the α-subunit encoded either by Sec61A1 or its paralog Sec61A2. Although the underlying channel selectivity for preproinsulin translocation is unknown, almost all studies of Sec61α to date have focused on Sec61α1. There is currently no evidence to suggest that this gene product plays a major role in proinsulin production, whereas genome-wide association studies indicate linkage of Sec61A2 with diabetes. Here, we report that evolutionary differences in mouse preproinsulin signal peptides affect proinsulin biosynthesis. Moreover, we find that, although some preproinsulin translocation can proceed through Sec61α1, Sec61α2 has a greater impact on proinsulin biosynthesis in pancreatic β-cells. Remarkably, Sec61α2 translocon deficiency exerts a significant inhibitory effect on the biosynthesis of preproinsulin itself, including a disproportionate increase of full-length nascent chain unreleased from ribosomes. This study not only reveals novel translocational regulation of proinsulin biosynthesis but also provides a rationale for genetic evidence suggesting an important role of Sec61α2 in maintaining blood glucose homeostasis.

ARTICLE HIGHLIGHTS

Reversing the relative time courses of the peptide bond reaction with oligopeptides of different lengths and charged amino acid distributions in the ribosome exit tunnel

The kinetics of the protein elongation cycle by the ribosome depends on intertwined factors. One of these factors is the electrostatic interaction of the nascent protein with the ribosome exit tunnel. In this computational biology theoretical study, we focus on the rate of the peptide bond formation and its dependence on the ribosome exit tunnel electrostatic potential profile. We quantitatively predict how oligopeptides of variable lengths can affect the peptide bond formation rate. We applied the Michaelis-Menten model as previously extended to incorporate the mechano-biochemical effects of forces on the rate of reaction at the catalytic site of the ribosome. For a given pair of carboxy-terminal amino acid substrate at the P- and an aminoacyl-tRNA at the A-sites, the relative time courses of the peptide bond formation reaction can be reversed depending on the oligopeptide sequence embedded in the tunnel and their variable lengths from the P-site. The reversal is predicted to occur from a shift in positions of charged amino acids upstream in the oligopeptidyl-tRNA at the P-site. The position shift must be adjusted by clever design of the oligopeptide probes using the electrostatic potential profile along the exit tunnel axial path. These predicted quantitative results bring strong evidence of the importance and relative contribution of the electrostatic interaction of the ribosome exit tunnel with the nascent peptide chain during elongation.

Mechanochemical forces regulate the composition and fate of stalled nascent chains

The ribosome-associated quality control (RQC) pathway resolves stalled ribosomes. As part of RQC, stalled nascent polypeptide chains (NCs) are appended with CArboxy-Terminal amino acids (CAT tails) in an mRNA-free, non-canonical elongation process. CAT tail composition includes Ala, Thr, and potentially other residues. The relationship between CAT tail composition and function has remained unknown. Using biochemical approaches in yeast, we discovered that mechanochemical forces on the NC regulate CAT tailing. We propose CAT tailing initially operates in an "extrusion mode" that increases NC lysine accessibility for on-ribosome ubiquitination. Thr in CAT tails enhances NC extrusion by preventing formation of polyalanine, which can form α-helices that lower extrusion efficiency and disrupt termination of CAT tailing. After NC ubiquitylation, pulling forces on the NC switch CAT tailing to an Ala-only "release mode" which facilitates nascent chain release from large ribosomal subunits and NC degradation. Failure to switch from extrusion to release mode leads to accumulation of NCs on large ribosomal subunits and proteotoxic aggregation of Thr-rich CAT tails.

BEMM-GEN: A Toolkit for Generating a Biomolecular Environment-Mimicking Model for Molecular Dynamics Simulation

Understanding the influence of the cellular environment on protein conformations is crucial for elucidating protein functions within living cells. In studies using molecular dynamics (MD) simulation, carbon nanotubes and hydrophobic cages have been widely used to emulate the cellular environment inside specific large biomolecules such as ribosome tunnels and chaperones. However, recent studies suggest that these uniform hydrophobic models may not adequately capture the environmental effects inside each biomolecule. Based on these facts, it is necessary to generate spherical and cylindrical models with varied chemical properties corresponding to the components within target biomolecules. We developed a biomolecular environment-mimicking model generator (BEMM-GEN) that generates spherical and cylindrical models with user-specified chemical properties and allows the integration of arbitrary protein conformations into the generated models. BEMM-GEN provides model and protein complex structures, along with the corresponding parameter files for MD simulation (AMBER and GROMACS), and users immediately run their MD simulation based on the generated input files. BEMM-GEN can be freely downloaded and installed via a Python package manager (pip install BEMM-gen). The source code files and a user manual for operation are provided on GitHub (https://github.com/y4suda/BEMM-GEN).

Ribosome Tunnel Environment Drives the Formation of α-Helix during Cotranslational Folding

Protein conformations in cells are not solely determined by amino acid sequences; they also depend on cellular environments. For instance, the ribosome tunnel induces its specific α-helix formation during cotranslational folding. Owing to the link between these temporally α-helix and biological functions, the mechanism of α-helix formation inside the ribosome tunnel has been previously explored. Consequently, the conformational restrictions of the tunnel were considered one of the driving forces of α-helix formation. Conversely, the ribosomal tunnel environment, including its chemical properties, appears to influence the α-helix formation. However, a comprehensive analysis of the ribosome tunnel environment's impact on the α-helix formation has not been conducted yet due to challenges in experimentally controlling it. Therefore, as a new computational approach, we proposed a ribosome environment-mimicking model (REMM) based on the radius and components of the experimentally determined ribosome tunnel structures. Using REMM, we assessed the impact of the ribosome tunnel environment on α-helix formation. Herein, we employed carbon nanotubes (CNT) as a reference model alongside REMM because CNT reproduce conformational restrictions rather than the ribosome tunnel environment. Quantitatively, the ability to reproduce the α-helix of nascent peptides in the experimental structure was compared between the CNT and REMM using enhanced all-atom molecular dynamics simulations. Consequently, the REMM more accurately reproduced the α-helix of the nascent peptides than the CNT, highlighting the significance of the ribosome tunnel environment in α-helix formation. Additionally, we analyzed the properties of the peptide inside each model to reveal the mechanism of ribosome tunnel-specific α-helix formation. Consequently, we revealed that the chemical diversities of the tunnel are essential for the formation of backbone-to-backbone hydrogen bonds in the peptides. In conclusion, the ribosome tunnel environment, with the diverse chemical properties, drives its specific α-helix formation. By proposing REMM, we newly provide the technical basis for investigating the protein conformations in various cellular environments.

Single-residue effects on the behavior of a nascent polypeptide chain inside the ribosome exit tunnel

Nascent polypeptide chains (NCs) are extruded from the ribosome through an exit tunnel (ET) traversing the large ribosomal subunit. The ET's irregular and chemically complex wall allows for various NC-ET interactions. Translational arrest peptides (APs) bind in the ET to induce translational arrest, a property that can be exploited to study NC-ET interactions by Force Profile Analysis (FPA). We employed FPA and molecular dynamics (MD) simulations to investigate how individual residues placed in a glycine-serine repeat segment within an AP-stalled NC interact with the ET to exert a pulling force on the AP and release stalling. Our results indicate that large and hydrophobic residues generate a pulling force on the NC when placed ≳10 residues away from the peptidyl transfer center (PTC). Moreover, an asparagine placed 12 residues from the PTC makes a specific stabilizing interaction with the tip of ribosomal protein uL22 that reduces the pulling force on the NC, while a lysine or leucine residue in the same position increases the pulling force. Finally, the MD simulations suggest how the Mannheimia succiniproducens SecM AP interacts with the ET to promote translational stalling.

Comprehensive quantitative modeling of translation efficiency in a genome-reduced bacterium

Translation efficiency has been mainly studied by ribosome profiling, which only provides an incomplete picture of translation kinetics. Here, we integrated the absolute quantifications of tRNAs, mRNAs, RNA half-lives, proteins, and protein half-lives with ribosome densities and derived the initiation and elongation rates for 475 genes (67% of all genes), 73 with high precision, in the bacterium Mycoplasma pneumoniae (Mpn). We found that, although the initiation rate varied over 160-fold among genes, most of the known factors had little impact on translation efficiency. Local codon elongation rates could not be fully explained by the adaptation to tRNA abundances, which varied over 100-fold among tRNA isoacceptors. We provide a comprehensive quantitative view of translation efficiency, which suggests the existence of unidentified mechanisms of translational regulation in Mpn.

A simple geometrical model of the electrostatic environment around the catalytic center of the ribosome and its significance for the elongation cycle kinetics

The central function of the large subunit of the ribosome is to catalyze peptide bond formation. This biochemical reaction is conducted at the peptidyl transferase center (PTC). Experimental evidence shows that the catalytic activity is affected by the electrostatic environment around the peptidyl transferase center. Here, we set up a minimal geometrical model fitting the available x-ray solved structures of the ribonucleic cavity around the catalytic center of the large subunit of the ribosome. The purpose of this phenomenological model is to estimate quantitatively the electrostatic potential and electric field that are experienced during the peptidyl transfer reaction. At least two reasons motivate the need for developing this quantification. First, we inquire whether the electric field in this particular catalytic environment, made only of nucleic acids, is of the same order of magnitude as the one prevailing in catalytic centers of the proteic enzymes counterparts. Second, the protein synthesis rate is dependent on the nature of the amino acid sequentially incorporated in the nascent chain. The activation energy of the catalytic reaction and its detailed kinetics are shown to be dependent on the mechanical work exerted on the amino acids by the electric field, especially when one of the four charged amino acid residues (R, K, E, D) has previously been incorporated at the carboxy-terminal end of the peptidyl-tRNA. Physical values of the electric field provide quantitative knowledge of mechanical work, activation energy and rate of the peptide bond formation catalyzed by the ribosome. We show that our theoretical calculations are consistent with two independent sets of previously published experimental results. Experimental results for E.coli in the minimal case of the dipeptide bond formation when puromycin is used as the final amino acid acceptor strongly support our theoretically derived reaction time courses. Experimental Ribo-Seq results on E. coli and S. cerevisiae comparing the residence time distribution of ribosomes upon specific codons are also well accounted for by our theoretical calculations. The statistical queueing time theory was used to model the ribosome residence time per codon during nascent protein elongation and applied for the interpretation of the Ribo-Seq data. The hypo-exponential distribution fits the residence time observed distribution of the ribosome on a codon. An educated deconvolution of this distribution is used to estimate the rates of each elongation step in a codon specific manner. Our interpretation of all these results sheds light on the functional role of the electrostatic profile around the PTC and its impact on the ribosome elongation cycle.

Translation rate underpins specific targeting of N-terminal transmembrane proteins to mitochondria

Protein biogenesis is a complex process, and complexity is greatly increased in eukaryotic cells through specific targeting of proteins to different organelles. To direct targeting, organellar proteins carry an organelle-specific targeting signal for recognition by organelle-specific import machinery. However, the situation is confusing for transmembrane domain (TMD)-containing signal-anchored (SA) proteins of various organelles because TMDs function as an endoplasmic reticulum (ER) targeting signal. Although ER targeting of SA proteins is well understood, how they are targeted to mitochondria and chloroplasts remains elusive. Here, we investigated how the targeting specificity of SA proteins is determined for specific targeting to mitochondria and chloroplasts. Mitochondrial targeting requires multiple motifs around and within TMDs: a basic residue and an arginine-rich region flanking the N- and C-termini of TMDs, respectively, and an aromatic residue in the C-terminal side of the TMD that specify mitochondrial targeting in an additive manner. These motifs play a role in slowing down the elongation speed during translation, thereby ensuring mitochondrial targeting in a co-translational manner. By contrast, the absence of any of these motifs individually or together causes at varying degrees chloroplast targeting that occurs in a post-translational manner.

Geometric differences in the ribosome exit tunnel impact the escape of small nascent proteins

The exit tunnel is the subcompartment of the ribosome that contains the nascent polypeptide chain and, as such, is involved in various vital functions, including regulation of translation and protein folding. As the geometry of the tunnel shows important differences across species, we focus on key geometrical features of eukaryote and prokaryote tunnels. We used a simple coarse-grained molecular dynamics model to study the role of the tunnel geometry in the post-translational escape of short proteins (short open reading frames [sORFs]) with lengths ranging from 6 to 56 amino acids. We found that the probability of escape for prokaryotes is one for all but the 12-mer chains. Moreover, proteins of this length have an extremely low escape probability in eukaryotes. A detailed examination of the associated single trajectories and energy profiles showed that these variations can be explained by the interplay between the protein configurational space and the confinement effects introduced by the constriction sites of the ribosome exit tunnel. For certain lengths, either one or both of the constriction sites can lead to the trapping of the protein in the "pocket" regions preceding these sites. As the distribution of existing sORFs indicates some bias in length that is consistent with our findings, we finally suggest that the constraints imposed by the tunnel geometry have impacted the evolution of sORFs.