Models of post-translational protein translocation

Chaperone-assisted translocation of flexible polymers in three dimensions

Polymer translocation through a nanometer-scale pore assisted by chaperones binding to the polymer is a process encountered in vivo for proteins. Studying the relevant models by computer simulations is computationally demanding. Accordingly, previous studies are either for stiff polymers in three dimensions or flexible polymers in two dimensions. Here, we study chaperone-assisted translocation of flexible polymers in three dimensions using Langevin dynamics. We show that differences in binding mechanisms, more specifically, whether a chaperone can bind to a single site or multiple sites on the polymer, lead to substantial differences in translocation dynamics in three dimensions. We show that the single-binding mode leads to dynamics that is very much like that in the constant-force driven translocation and accordingly mainly determined by tension propagation on the cis side. We obtain β≈1.26 for the exponent for the scaling of the translocation time with polymer length. This fairly low value can be explained by the additional friction due to binding particles. The multiple-site binding leads to translocation the dynamics of which is mainly determined by the trans side. For this process we obtain β≈1.36. This value can be explained by our derivation of β=4/3 for constant-bias translocation, where translocated polymer segments form a globule on the trans side. Our results pave the way for understanding and utilizing chaperone-assisted translocation where variations in microscopic details lead to rich variations in the emerging dynamics.

Quantifying the role of chaperones in protein translocation by computational modeling

The molecular chaperone Hsp70 plays a central role in the import of cytoplasmic proteins into organelles, driving their translocation by binding them from the organellar interior. Starting from the experimentally-determined structure of the E. coli Hsp70, we computed, by means of molecular simulations, the effective free-energy profile for substrate translocation upon chaperone binding. We then used the resulting free energy to quantitatively characterize the kinetics of the import process, whose comparison with unassisted translocation highlights the essential role played by Hsp70 in importing cytoplasmic proteins.

Polymer translocation: the first two decades and the recent diversification

Probably no other field of statistical physics at the borderline of soft matter and biological physics has caused such a flurry of papers as polymer translocation since the 1994 landmark paper by Bezrukov, Vodyanoy, and Parsegian and the study of Kasianowicz in 1996. Experiments, simulations, and theoretical approaches are still contributing novel insights to date, while no universal consensus on the statistical understanding of polymer translocation has been reached. We here collect the published results, in particular, the famous-infamous debate on the scaling exponents governing the translocation process. We put these results into perspective and discuss where the field is going. In particular, we argue that the phenomenon of polymer translocation is non-universal and highly sensitive to the exact specifications of the models and experiments used towards its analysis.

Polymer translocation through a nanopore driven by binding particles: influence of chain rigidity

We investigate the influence of chain rigidity on the dynamics of polymer translocation in the presence of binding particles (BPs) through a nanopore using two-dimensional Langevin dynamics simulations. With increasing chain rigidity κ, the mean translocation time 〈τ〉 increases monotonically due to an increase in the radius of gyration and a decrease in the center of mass velocity. Particularly for weak binding, we further find that 〈τ〉 shows a power-law behavior with the persistence length lp. Analysis indicates a scaling relation between the average velocity of the center of mass of a chain 〈vc.m.〉 and lp. As the chain becomes stiffer, the distribution of the translocation time τ approximates the Gaussian distribution and gets broader with the peak position being shifted towards longer translocation time. The corresponding translocation coordinate smax of the maximum waiting time gets smaller with increasing chain rigidity. Finally, under an extremely low BP concentration, 〈τ〉 shows a minimum for small κ, while it decreases monotonically for large κ with increasing binding energy. Our results suggest a nontrivial effect of the intrinsic property of chains on the dynamics of polymer translocation driven by BPs.

BiP clustering facilitates protein folding in the endoplasmic reticulum

The chaperone BiP participates in several regulatory processes within the endoplasmic reticulum (ER): translocation, protein folding, and ER-associated degradation. To facilitate protein folding, a cooperative mechanism known as entropic pulling has been proposed to demonstrate the molecular-level understanding of how multiple BiP molecules bind to nascent and unfolded proteins. Recently, experimental evidence revealed the spatial heterogeneity of BiP within the nuclear and peripheral ER of S. cerevisiae (commonly referred to as 'clusters'). Here, we developed a model to evaluate the potential advantages of accounting for multiple BiP molecules binding to peptides, while proposing that BiP's spatial heterogeneity may enhance protein folding and maturation. Scenarios were simulated to gauge the effectiveness of binding multiple chaperone molecules to peptides. Using two metrics: folding efficiency and chaperone cost, we determined that the single binding site model achieves a higher efficiency than models characterized by multiple binding sites, in the absence of cooperativity. Due to entropic pulling, however, multiple chaperones perform in concert to facilitate the resolubilization and ultimate yield of folded proteins. As a result of cooperativity, multiple binding site models used fewer BiP molecules and maintained a higher folding efficiency than the single binding site model. These insilico investigations reveal that clusters of BiP molecules bound to unfolded proteins may enhance folding efficiency through cooperative action via entropic pulling.

Simulating the pulling of stalled elongated peptide from the ribosome by the translocon

The nature of the coupling between the stalling of the elongated nascent peptide chain in the ribosome and its insertion through the translocon is analyzed, focusing on the recently discovered biphasic force that overcomes the stalling barrier. The origin of this long-range coupling is explored by coarse-grained simulations that combine the translocon (TR) insertion profile and the effective chemical barrier for the extension of the nascent chain in the ribosome. Our simulation determined that the inserted H segment is unlikely to climb the TR barrier in parallel with the peptide synthesis chemical step and that the nascent chain should first overcome the chemical barriers and move into the ribosome-TR gap region before the insertion into the TR tunnel. Furthermore, the simulations indicate that the coupled TR-chemistry free energy profile accounts for the biphasic force. Apparently, although the overall elongation/insertion process can be depicted as a tug-of-war between the forces of the TR and the ribosome, it is actually a reflection of the combined free-energy landscape. Most importantly, the present study helps to relate the experimental observation of the biphasic force to crucial information about the elusive path and barriers of the TR insertion process.

Long-timescale dynamics and regulation of Sec-facilitated protein translocation

We present a coarse-grained modeling approach that spans the nanosecond- to minute-timescale dynamics of cotranslational protein translocation. The method enables direct simulation of both integral membrane protein topogenesis and transmembrane domain (TM) stop-transfer efficiency. Simulations reveal multiple kinetic pathways for protein integration, including a mechanism in which the nascent protein undergoes slow-timescale reorientation, or flipping, in the confined environment of the translocon channel. Competition among these pathways gives rise to the experimentally observed dependence of protein topology on ribosomal translation rate and protein length. We further demonstrate that sigmoidal dependence of stop-transfer efficiency on TM hydrophobicity arises from local equilibration of the TM across the translocon lateral gate, and it is predicted that slowing ribosomal translation yields decreased stop-transfer efficiency in long proteins. This work reveals the balance between equilibrium and nonequilibrium processes in protein targeting, and it provides insight into the molecular regulation of the Sec translocon.

Pushing, pulling and trapping--modes of motor protein supported protein translocation

Protein translocation across the cellular membranes is an ubiquitous and crucial activity of cells. This process is mediated by translocases that consist of a protein conducting channel and an associated motor protein. Motor proteins interact with protein substrates and utilize the free energy of ATP binding and hydrolysis for protein unfolding, translocation and unbinding. Since motor proteins are found either at the cis- or trans-side of the membrane, different mechanisms for translocation have been proposed. In the Power stroke model, cis-acting motors are thought to push, while trans-motors pull on the substrate protein during translocation. In the Brownian ratchet model, translocation occurs by diffusion of the unfolded polypeptide through the translocation pore while directionality is achieved by trapping and refolding. Recent insights in the structure and function of the molecular motors suggest that different mechanisms can be employed simultaneously.

Directed motion emerging from two coupled random processes: translocation of a chain through a membrane nanopore driven by binding proteins

We investigate the translocation of a stiff polymer consisting of M monomers through a nanopore in a membrane, in the presence of binding particles (chaperones) that bind onto the polymer, and partially prevent backsliding of the polymer through the pore. The process is characterized by the rates: k for the polymer to make a diffusive jump through the pore, q for unbinding of a chaperone, and the rate qκ for binding (with a binding strength κ); except for the case of no binding κ = 0 the presence of the chaperones gives rise to an effective force that drives the translocation process. In more detail, we develop a dynamical description of the process in terms of a (2+1)-variable master equation for the probability of having m monomers on the target side of the membrane with n bound chaperones at time t. Emphasis is put on the calculation of the mean first passage time [Formula: see text] as a function of total chain length M. The transfer coefficients in the master equation are determined through detailed balance, and depend on the relative chaperone size λ and binding strength κ, as well as the two rate constants k and q. The ratio γ = q/k between the two rates determines, together with κ and λ, three limiting cases, for which analytic results are derived: (i) for the case of slow binding ([Formula: see text]), the motion is purely diffusive, and [Formula: see text] for large M; (ii) for fast binding ([Formula: see text]) but slow unbinding ([Formula: see text]), the motion is, for small chaperones λ = 1, ratchet-like, and [Formula: see text]; (iii) for the case of fast binding and unbinding dynamics ([Formula: see text] and [Formula: see text]), we perform the adiabatic elimination of the fast variable n, and find that for a very long polymer [Formula: see text], but with a smaller prefactor than for ratchet-like dynamics. We solve the general case numerically as a function of the dimensionless parameters λ, κ and γ, and compare to the three limiting cases.

Computer simulation of polypeptide translocation through a nanopore

A simplified model of polypeptide chains was designed and studied by means of computer simulations. Chains were represented by a sequence of united atoms located at the positions of the alpha-carbons. A further assumption was the lattice approximation for the chains. We used a (310) lattice, which was found useful for studying properties of proteins. The force field used consisted of a long-range contact potential between amino-acid residues and a local preference for forming alpha-helical states. The chain consisted of two kinds of residues: hydrophilic (P) and hydrophobic (H) ones forming model helical septets--HHPPHPP--in a sequence. The chains were placed near an impenetrable surface with a square hole in it. The size of the hole was comparable or smaller than the size of a chain. The properties of these model chains were determined using the Monte-Carlo simulation method. During the simulations, translocation of the chain through the hole in the wall was observed. The influence of the chain length, the temperature differences on both sides of the wall and the force field on the chain properties were investigated. It was shown that the translocation time scales as N(2.2) and it was found that the presence of the local helical potential significantly slows down the process of translocation. [Figure: see text]. The snapshots of typical chain's conformation obtained during the simulation for chain consisted of N = 60. The values of the local potential epsilon(loc) = -8.

Chaperone-assisted translocation

We investigate the translocation of a stiff polymer through a nanopore in a membrane, in the presence of binding particles (chaperones) that bind reversibly to the polymer on both sides of the membrane. A bound chaperone covers one (univalent binding) or many (multivalent binding) binding sites. Assuming that the diffusion of the chaperones is fast compared to the rate of translocation we describe the process by a one-dimensional master equation. We expand previous models by a detailed study of the effective force in the master equation, which is obtained by the appropriate statistical mechanical average over the chaperone states. The dependence of the force on the degree of valency (the number of binding sites occupied by a chaperone) is studied in detail. We obtain finite size corrections (to the thermodynamical expression for the force), which, for univalent binding, can be expressed analytically. We finally investigate the mean velocity for translocation as a function of chaperone binding strength and size. For both univalent and multivalent binding simple results are obtained for the case of a sufficiently long translocating polymer.

Precursors with altered affinity for Hsp70 in their transit peptides are efficiently imported into chloroplasts

Protein import into chloroplasts is postulated to occur with the involvement of molecular chaperones. We have determined that the transit peptide of ferredoxin-NADP(H) reductase precursor binds preferentially to an Hsp70 from chloroplast stroma. To investigate the role of Hsp70 molecular chaperones in chloroplast protein import, we analyzed the import into pea chloroplasts of preproteins with decreased Hsp70 binding affinity in their transit peptides. Our results indicate that the precursor with the lowest affinity for Hsp70 molecular chaperones in its transit peptide was imported to chloroplasts with similar apparent Km as the wild type precursor and a 2-fold increase in Vmax. Thus, a strong interaction between chloroplast stromal Hsp70 and the transit peptide seems not to be essential for protein import. These results indicate that in chloroplasts the main unfolding force during protein import may be applied by molecular chaperones other than Hsp70s. Although stromal Hsp70s undoubtedly participate in chloroplast biogenesis, the role of these molecular chaperones in chloroplast protein translocation differs from the one proposed in the mechanisms postulated up to date.

Molt cycle–dependent molecular chaperone and polyubiquitin gene expression in lobster

Lobster claw muscle undergoes atrophy in correlation with increasing ecdysteroid (steroid molting hormone) titers during premolt. In vivo molecular chaperone (constitutive heat shock protein 70 [Hsc70], heat shock protein 70 [Hsp70], and Hsp90) and polyubiquitin messenger ribonucleic acid (mRNA) levels were examined in claw and abdominal muscles from individual premolt or intermolt lobsters. Polyubiquitin gene expression was assayed as a marker for muscle atrophy. Both Hsc70 and Hsp90 mRNA levels were significantly induced in premolt relative to intermolt lobster claw muscle, whereas Hsp70 mRNA levels were not. Hsp90 gene expression was significantly higher in premolt claw muscle when compared with abdominal muscle. Polyubiquitin mRNA levels were elevated in premolt when compared with intermolt claw muscle and significantly elevated relative to premolt abdominal muscle.

The brownian ratchet and power stroke models for posttranslational protein translocation into the endoplasmic reticulum

A quantitative analysis of experimental data for posttranslational translocation into the endoplasmic reticulum is performed. This analysis reveals that translocation involves a single rate-limiting step, which is postulated to be the release of the signal sequence from the translocation channel. Next, the Brownian ratchet and power stroke models of translocation are compared against the data. The data sets are simultaneously fit using a least-squares criterion, and both models are found to accurately reproduce the experimental results. A likelihood-ratio test reveals that the optimal fit of the Brownian ratchet model, which contains one fewer free parameter, does not differ significantly from that of the power stroke model. Therefore, the data considered here cannot be used to reject this import mechanism. The models are further analyzed using the estimated parameters to make experimentally testable predictions.