Role of cotranslational folding for β-sheet-enriched proteins: A perspective from molecular dynamics simulations

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.

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.

Review of Electrostatic Force Calculation Methods and Their Acceleration in Molecular Dynamics Packages Using Graphics Processors

Molecular dynamics (MD) simulations probe the conformational repertoire of macromolecular systems using Newtonian dynamic equations. The time scales of MD simulations allow the exploration of biologically relevant phenomena and can elucidate spatial and temporal properties of the building blocks of life, such as deoxyribonucleic acid (DNA) and protein, across microsecond (μs) time scales using femtosecond (fs) time steps. A principal bottleneck toward extending MD calculations to larger time scales is the long-range electrostatic force measuring component of the naive nonbonded force computation algorithm, which scales with a complexity of (N, number of atoms). In this review, we present various methods to determine electrostatic interactions in often-used open-source MD packages as well as the implementation details that facilitate acceleration of the electrostatic interaction calculation.