Regulation of multispanning membrane protein topology via post-translational annealing

  1. Reid C Van Lehn
  2. Bin Zhang
  3. Thomas F Miller III  Is a corresponding author
  1. California Institute of Technology, United States

Abstract

The canonical mechanism for multispanning membrane protein topogenesis suggests that protein topology is established during cotranslational membrane integration. However, this mechanism is inconsistent with the behavior of EmrE, a dual-topology protein for which the mutation of positively charged loop residues, even close to the C-terminus, leads to dramatic shifts in its topology. We use coarse-grained simulations to investigate the Sec-facilitated membrane integration of EmrE and its mutants on realistic biological timescales. This work reveals a mechanism for regulating membrane-protein topogenesis, in which initially misintegrated configurations of the proteins undergo post-translational annealing to reach fully integrated multispanning topologies. The energetic barriers associated with this post-translational annealing process enforce kinetic pathways that dictate the topology of the fully integrated proteins. The proposed mechanism agrees well with the experimentally observed features of EmrE topogenesis and provides a range of experimentally testable predictions regarding the effect of translocon mutations on membrane protein topogenesis.

https://doi.org/10.7554/eLife.08697.001

eLife digest

Proteins are long chains of smaller molecules called amino acids, and are built inside cells by a molecular machine called the ribosome. Many important proteins must be inserted into the membrane that surrounds each cell in order to carry out their role. As these proteins are being built by the ribosome, they thread their way into a membrane-spanning channel (called the translocon) from the inner side of the membrane. Short segments of these integral membrane proteins (called transmembrane domains) then become embedded in the membrane, while other parts of the protein remain on either side of the membrane.

For a membrane protein to work properly, the end of each of its transmembrane domains must be on the correct side of the membrane (i.e., the protein must obtain the correct ‘topology’). The conventional model for this process suggests that topology is fixed when the first transmembrane domain of a protein is initially integrated into the membrane, while the ribosome is still building the protein. This model can explain most integral membrane proteins, which only have a single topology. However, it cannot explain the family of membrane proteins that have an almost equal chance of adopting one of two different topologies (so-called ‘dual-topology proteins’).

Van Lehn et al. have now used computer modeling to simulate how a bacterial protein called EmrE (which is a dual-topology protein) integrates into the membrane via the translocon. The results reveal that a few transmembrane domains in EmrE do not fully integrate into the membrane while the ribosome is building the protein. Instead, these transmembrane domains slowly integrate after the ribosome has finished its job.

These findings contradict the conventional model and suggest that some membrane proteins only become fully integrated after the protein-building process is complete. The next step in this work is to experimentally test predictions from the computer simulations.

https://doi.org/10.7554/eLife.08697.002

Introduction

Integral membrane proteins (IMPs) are central to cellular functions that include signal transduction, transport across the cell membrane, and energy conversion. Performing these roles requires integration of the IMPs into the membrane with the correct topology (i.e., the correct orientation of the fully integrated IMP relative to the membrane). In most cases, membrane integration proceeds via the Sec translocon, a conserved protein-conducting channel located in the endoplasmic reticulum membrane in eukaryotes or in the plasma membrane in bacteria (White and von Heijne, 2004). During this process, the ribosome or other molecular motor docks to the cytoplasmic opening of the translocon, feeding the nascent protein into the translocon channel (Shao and Hegde, 2011); conformational changes in the lateral gate (LG) helices of the translocon then allow sufficiently hydrophobic segments of the nascent protein to integrate as transmembrane domains (TMD) (Hessa et al., 2005; Egea and Stroud, 2010; Zhang and Miller, 2010; Gogala et al., 2014). The orientation of a single TMD relative to the membrane is determined by factors that include the hydrophobicity of the TMD and the charge and length of the soluble loops that flank the TMD (Goder and Spiess, 2001, 2003; Devaraneni et al., 2011). However, the extent to which these factors influence the topology of multispanning IMPs is less clear.

The conventional model of multispanning IMP topogenesis assumes that a single dominant topology is established via the successive integration of TMDs that thread back-and-forth across the membrane in alternating orientations (Blobel, 1980; Wessels and Spiess, 1988; Sadlish et al., 2005). In this cotranslational model, the dominant IMP topology is determined by the orientation of the N-terminal TMD and is primarily dictated by the features of that leading TMD (Hartmann et al., 1989; Borel and Simon, 1996; Dale et al., 2000). However, the cotranslational model is challenged by dual-topology proteins, which exhibit both possible orientations of the fully integrated IMP with respect to the membrane in approximately 1:1 stoichiometry (Rapp et al., 2006, 2007). The most thoroughly studied dual-topology protein is the bacterial multidrug transporter EmrE (Chen et al., 2007), which can be biased in favor of a single dominant topology by introducing positive charges to any of its soluble loops (Rapp et al., 2006, 2007; Seppälä et al., 2010). The dominant topology of each EmrE mutant retains the loop with the additional positive charges in the cytoplasm (Seppälä et al., 2010), apparently satisfying the empirical trend known as the ‘positive-inside’ rule which notes that the combined charges of the cytoplasmic loops (i.e., K+R bias) of an IMP correlates with its dominant topology (von Heijne, 1986). Surprisingly, adding charges to even C-terminal loops can influence the dominant topology of EmrE, suggesting that such mutations have a long-range effect on the orientation of previously-translated TMDs. This finding is inconsistent with the cotranslational model and raises interesting questions about IMP topogenesis. At what point is IMP topology established with respect to ribosomal translation? Are TMD orientations locked-in during the period in which the nascent IMP is attached to the ribosome (i.e., cotranslationally) or do TMD orientations remain subject to change even upon completion of ribosomal translation (i.e., post-translationally)?

In this work, we simulate the topogenesis of EmrE and its mutants to address limitations in the cotranslational model of IMP topogenesis by understanding when IMP topology is established (co- or post-translationally) and how topology is regulated. We use a coarse-grained (CG) model that enables access to a timescale of minutes while retaining sufficient chemical accuracy to capture the forces that drive membrane integration (Zhang and Miller, 2012a). The distribution of topologies predicted by the simulations are in good agreement with previous experimental findings (Rapp et al., 2007; Seppälä et al., 2010). The simulation results show that TMDs in the dual-topology mutants do not completely integrate by the end of translation; instead, the slow post-translational flipping of loops across the membrane allows misintegrated TMDs to reorient and insert into the membrane. The fully integrated topology is determined by the position of the loop that undergoes flipping most slowly. This work elucidates the mechanism by which dual-topology protein topology is established, reconciles dominant protein topologies with the positive-inside rule, and predicts the role that the translocon plays in mediating multispanning IMP topogenesis. Other examples of post-translational topological changes in diverse multispanning IMP systems suggest that this mechanism may have generality beyond EmrE (Lu et al., 2000; Lambert and Prange, 2001; Kanki et al., 2002; Skach, 2009; Öjemalm et al., 2012; Bowie, 2013; Virkki et al., 2014).

Coarse-grained model

The cotranslational integration and topogenesis of EmrE and its mutants is simulated using a recently developed CG model (Zhang and Miller, 2012a), which we employ essentially unchanged from its original introduction. Figure 1 illustrates the CG representation of a nascent protein and the protocol for simulating membrane integration. The ribosome, translocon, and nascent protein are all composed of CG beads. Each bead has a diameter of σ = 0.8 nm to represent approximately three amino-acid residues. This bead diameter is similar to the Kuhn length of polypeptides (Staple et al., 2008) so that the nascent protein can be treated as a freely jointed chain. The surrounding solvent and lipid bilayer are included implicitly, a technique that is used in other CG models of the translocon (Rychkova and Warshel, 2013). The time-evolution of nascent protein configurations is calculated using Brownian dynamics with a 100 ns timestep. The kinetics of the LG are modeled as stochastic transitions between a closed conformation, which prevents the nascent protein from exiting from the channel interior to the membrane, and an open conformation, which removes the barrier to membrane entry. All bead positions are projected onto the plane that passes along the translocon channel axis between the helices forming the LG. This off-lattice 2D approximation reflects the cylindrical geometry of the channel and is inspired by previous models of biopolymer translocation through nanopores (Huopaniemi et al., 2006). Beads representing the ribosome enclosure and translocon are placed to approximate their structures (Van den Berg et al., 2004; Frauenfeld et al., 2011). Two negative charges are placed on a bead at the cytosolic end of the translocon LG, whereas two positive charges are placed on a bead at the periplasmic end of the translocon LG. This charge distribution reflects the position of conserved charged residues (White and von Heijne, 2004) near the translocon LG that have been previously shown to affect single-spanning protein topogenesis (Goder et al., 2004). Full details of the model are provided in Appendix 1.

Figure 1 with 2 supplements see all
Schematic of Sec-mediated cotranslational integration of EmrE and corresponding simulation representation.

(A) At top, an illustration of the structural motifs in EmrE, including indication of the charged residues in the soluble loops with black circles and the transmembrane domain (TMD)/loop numbering scheme that is employed in the text; below, the corresponding sequence of coarse-grained (CG) beads that represent the EmrE amino-acid sequence. TMDs and loops are assigned based on the hydropathy plot and consensus topology prediction shown in Figure 1—figure supplement 1. (B) At top, a schematic illustration of the sequential integration of TMDs to obtain a multispanning Nperi/Cperi topology, in which both the N- and C-terminal loops are positioned in the periplasm, according to the cotranslational model; below, representative simulation snapshots of EmrE as the nascent protein grows during translation, integrates into the membrane, and exits the channel in the Nperi/Cperi multispanning topology. The nascent protein is colored according to the legend at top, the ribosome is brown, and the translocon is green with translocon charges labeled explicitly.

https://doi.org/10.7554/eLife.08697.003

The CG model is well-suited to simulating the kinetics of cotranslational IMP integration, a process that is challenging for atomistic models (Zhang and Miller, 2010; Gumbart et al., 2011; Zhang and Miller, 2012b; Rychkova and Warshel, 2013) due to the large system size (>100,000 atoms) and the long timescale (minutes) of translation. We note that the model does not include nascent protein secondary/tertiary structure, charged lipids, protein chaperones, or an electrostatic potential across the membrane. However, the model does include explicit LG/translation dynamics, electrostatic interactions with the translocon, water/bilayer transfer free energies, and a direct mapping between the nascent protein sequence and the CG representation. The model thus captures the major physicochemical features of the translocon-membrane system (White and von Heijne, 2004). Moreover, the model has been shown to accurately predict features of single-spanning IMP integration and topogenesis (Zhang and Miller, 2012a), including the sigmoidal dependence of stop-transfer efficiency on TMD hydrophobicity (Hessa et al., 2005), the inversion of signal-anchor orientation during translation (Goder and Spiess, 2003), and the effect of translation rates and sequence features on signal-anchor orientation (Goder and Spiess, 2003). In particular, the model has been shown (Zhang and Miller, 2012a) to correctly describe integration processes that are governed either by thermodynamics (Hessa et al., 2005) or kinetics (Goder and Spiess, 2003), and it has provided a means of understanding the competition between such effects. The model has also been shown to correctly predict the dominant topology for a three-TMD multispanning IMP with a strong positive-inside bias (Zhang and Miller, 2012a). The strong agreement between simulation and experimental results presented in this work further indicates that IMP topological determinants are captured at this CG resolution.

EmrE protein

The EmrE amino-acid sequence includes four hydrophobic domains and five hydrophilic loops, according to both the hydropathy plot and consensus topology prediction shown in Figure 1—figure supplement 1. The hydropathy plot was calculated using the Wimley–White hydrophobicity scale (Wimley et al., 1996). The black line in the hydropathy plot indicates the water–octanol transfer free energy per residue and the overlaid red line shows a moving average using a 7-residue window. The consensus topology prediction was generated by the TOPCONS 1.0 server (Bernsel et al., 2009) and agrees with previous representations of EmrE structural elements (Seppälä et al., 2010). Shaded regions in the hydropathy plot indicate the predicted TMDs and loops.

In the CG model, each TMD is represented by four CG beads and each soluble loop is represented by five CG beads, as seen in Figure 1A. The CG beads assume one of four types as determined by the associated amino-acid residues in the nascent protein; these CG bead-types include V (moderately hydrophobic), L (very hydrophobic), Q (neutral-hydrophilic), and K (positively charged). Among these types, the CG beads vary with respect to their charge and their water/membrane transfer free energies (Appendix table 1). In the hydropathy profile, the N-terminal TMD (TMD1) is less hydrophobic than the other three TMDs, so its beads are assigned the V bead type. All other TMD beads are assigned the L bead type. Beads in each soluble loop are assigned to either the K or Q bead type, depending on the location of positive charges in the amino-acid sequence; positive charges are highlighted in red in the EmrE wild-type amino-acid sequence in Figure 1—figure supplement 1. Each K bead type is assigned a +2 charge, following previous work (Zhang and Miller, 2012a). Negative charges are excluded from the CG representation of EmrE, because EmrE exhibits a small number of such charges (Figure 1—figure supplement 1) and because the experimentally studied EmrE mutations focus only on the addition/removal of positively charged residues (Seppälä et al., 2010). Nonetheless, the effect of negatively charged residues in the CG simulation was explicitly tested in Figure 5—figure supplement 1 and was found to be minor. Similarly, the results of the simulations are robust with respect to changes in the modeling of TMD1 hydrophobicity (Figure 5—figure supplement 1) and loop length (Figure 3—figure supplement 3).

Using the CG model, we consider a series of EmrE mutants from Rapp et al. (2007) and Seppälä et al. (2010). We include EmrE mutants with single charge mutations—K3, T28R, A52K, L85R, and R111—from Seppälä et al. (2010) and EmrE mutants with single dominant topologies—EmrE(Ncyto) and EmrE(Nperi)—from Rapp et al. (2007). We also consider a series of mutants in which the protein has either zero positive charge or positive charges in only a single loop—nEmrE, nK3, nT28R1, nT28R2, nT28R, nA52K, nL85R, and nR111—from Seppälä et al. (2010). This list includes all 16 of the EmrE and nEmrE mutants with single added charges studied experimentally by Seppälä et al. (2010); mutants with added C-terminal His residues or an extra TMD are not considered. Finally, we include a ‘cotranslationally-biased’, or CB, mutant that has elongated, 10-bead hydrophilic loops and two positives charges in the first, third, and fifth loops to create a strong K+R bias that favors a Ncyto/Ccyto topology (i.e., with both the N-terminal and C-terminal loops in the cytoplasm) according to the positive-inside rule (von Heijne, 1986; Rapp et al., 2006); this protein is expected to be strongly biased towards membrane integration via the cotranslational mechanism, providing a useful comparison with the other EmrE mutants. The CG representation of each mutant is listed in Appendix table 2; for each mutant, charge mutations are reflected by changing between Q-type and K-type beads at the appropriate point in the sequence. Despite its simplicity, we emphasize that the CG representation captures the major features of EmrE and its mutants, including the number of TMDs/loops and the distribution of charges.

Simulation protocol

As illustrated in Figure 1B, the dynamics of the ribosome/nascent protein/translocon complex is directly simulated using the CG model. Each CG trajectory is initiated with a short nascent protein attached to the ribosome exit channel; as a function of time, the nascent protein grows in length (while remaining attached to the ribosome) until it completes translation and is released from the ribosome. The dynamics of the nascent protein continue to be simulated until the protein reaches a fully integrated topology.

Simulations are initialized from equilibrated configurations of the nascent protein, initially comprised of 9 CG beads, with the C-terminus attached to the ribosome exit channel (Figure 1—figure supplement 2). Translation is performed by adding a new CG bead to the C-terminus of the nascent protein and attaching it to the ribosome exit channel; the previous C-terminus is released from the exit channel. The simulation is then continued for 125 ms before the next bead is added, a simulation time which corresponds to a translation rate of 24 residues/s (Bilgin et al., 1992). At the end of translation, the C-terminus is released from the ribosome exit channel and simulations are continued until all beads in the TMDs are at least 4.5σ from the origin and integrated with either a Ncyto/Ccyto or Nperi/Cperi topology. The ribosome remains bound to the translocon for the duration of all simulations (Potter and Nicchitta, 2002; Schaltetzky and Rapoport, 2006). The distance threshold ensures that the final configuration of the protein has exited from both the ribosome and translocon channel.

The trajectory termination criteria are designed to examine the effects of the Sec-facilitated membrane integration process on EmrE topogenesis. Specifically, it is assumed that upon reaching configurations in which all of the TMDs are integrated into the membrane, the protein topology remains irreversibly fixed for all subsequent times; physical processes that may lead to this irreversibility include the dimerization of EmrE proteins to form functional channels in the membrane (Lloris-Garcerá et al., 2012) or the degradation of undimerized EmrE proteins prior to topological inversion (Woodall et al., 2015). Given the symmetry of the membrane-protein interactions in the absence of the translocon, if the CG trajectories were allowed to run for infinitely long times to reach full equilibration after diffusing away from the translocon, the relative probability of the Ncyto/Ccyto and Nperi/Cperi topologies would be equal, regardless of the protein sequence. The employed trajectory termination criteria thus isolate the role of the non-equilibrium integration process in determining IMP topology. Demonstration of the robustness of the reported results to the cutoff values employed in the trajectory termination criteria are provided in the Robustness checks for the trajectory termination criteria section of the ‘Materials and methods’.

The integration and orientation of a TMD is interpreted from the positions of hydrophobic beads in each TMD and the third bead in each hydrophilic loop. The coordinate system is defined with the x-axis perpendicular to the bilayer (Figure 1—figure supplement 2). The origin is placed at the center of the channel such that negative x-values indicate cytoplasmic positions. A TMD is considered integrated if −2σx ≤ 2σ for all four hydrophobic beads, corresponding to positions within the implicit bilayer, and if all y-positions are outside of the translocon interior. A loop is considered to be in the cytoplasm if the position of the reference bead satisfies x < −σ and in the periplasm if x > σ. The Ncyto/Ccyto topology is reached if the first, third, and fifth loops are positioned in the cytoplasm and the second and fourth loops are positioned in the periplasm. The Nperi/Cperi topology has the opposite loop positions as shown in Figure 1B.

For each mutant, 250 independent trajectories are performed for a total of 4000 CG trajectories and nearly 6000 min of aggregate simulation time. Error bars measure the standard error between 2 blocks of 125 simulated trajectories. Complete system configurations are saved every 50 ms while loop positions and TMD orientations are saved every 1 ms.

Results

Simulations match experimental observations of topology

For all 16 of the EmrE and nEmrE mutants with single added charges studied by Seppälä et al. (2010), Figure 2 compares the experimentally observed IMP topologies with the prediction from the CG model. Specifically, the figure compares the fraction of fully integrated proteins that adopt the Ncyto/Ccyto topology, with the remainder in the Nperi/Cperi topology. The top and bottom rows show variants of EmrE and nEmrE respectively. Each mutant differs only in the number and location of charges in the hydrophilic loops. A schematic of each mutant drawn with the dominant topology predicted from simulations is included; positive charges are indicated as filled-in circles with additional charges relative to EmrE (top row) or nEmrE (bottom row) highlighted in red. The topologies determined experimentally in Seppälä et al. (2010) are expressed as the fraction of Ncyto/Ccyto topologies by dividing the cell activity of each protein coexpressed with the EmrE(Nperi) mutant by the total growth of the protein coexpressed with either the EmrE(Nperi) or EmrE(Ncyto) mutant (Seppälä et al., 2010), as described in the Experimental interpretation of EmrE topology section of the ‘Materials and methods’.

Figure 2 with 2 supplements see all
Topologies determined from simulations (blue) and compared to the experiments of Seppälä et al. (2010) (red), reporting the fraction of fully integrated integral membrane protein (IMP) configurations in the Ncyto/Ccyto topology.

Error bars indicate the standard error measured from independent blocks of simulations or taken from Seppälä et al. (2010). The dominant topology for each mutant is indicated schematically with additional positive charges relative to EmrE (top) or nEmrE (bottom) drawn as red dots.

https://doi.org/10.7554/eLife.08697.006

It is clear from Figure 2 that the simulations are in excellent qualitative agreement with the experiments by correctly predicting the near 1:1 stoichiometry of wild-type EmrE and identifying the dominant topology for nearly all of the proteins considered. Figure 2—figure supplement 2 illustrates that the distribution of topologies determined experimentally and the distribution of topologies measured from the simulations are linearly correlated (Pearson correlation coefficient, r = 0.92); points lying in the two shaded quadrants of the graph correspond to proteins for which the simulations and experiments predict consistent topologies. All mutants, with the exception of A52K, have the same dominant topology in the simulations as in the experiments within the statistical error. The agreement between simulations and experiments suggests that the CG model correctly reproduces the essential molecular features of topogenesis; in the following, we analyze the ensembles of CG trajectories that give rise to these computed IMP topologies.

Dual-topology proteins exhibit slow post-translational integration

To investigate the molecular processes that govern the establishment of EmrE topology, we first examine the kinetics by which fully integrated topologies are reached. As a function of time, Figure 3A shows the fraction of CG trajectories in which the studied protein reaches a fully integrated topology for several EmrE mutants and the CB mutant. 0 s corresponds to the end of translation and negative values of time correspond to the period that precedes the end of ribosomal translation in which the nascent protein is still attached to the ribosome. Over 90% of the CB mutant trajectories reach the Ncyto/Ccyto topology within 3 s following the completion of translation and thus rapidly integrate as expected for the cotranslational model (Blobel, 1980; Wessels and Spiess, 1988; Sadlish et al., 2005); mechanistic features of individual TMD integration steps are discussed in the Cotranslational integration pathways section of the ‘Materials and methods’. In contrast, all variants of EmrE reach a fully integrated topology much more slowly, requiring hundreds of seconds for some CG trajectories to fully integrate (see also Figure 3—figure supplement 2).

Figure 3 with 3 supplements see all
Kinetics of EmrE topogenesis.

(A) Fraction of CG trajectories in which all TMDs are fully integrated in a multispanning topology, plotted as a function of time for several mutants. (B) Fraction of CG trajectories in which each TMD is integrated, plotted as a function of time for the cotranslationally-biased (CB) mutant (top) and EmrE (bottom). The snapshots show an example of a simulation in which TMD4 of EmrE does not integrate during translation. In both panels, 0 s corresponds to the end of translation and negative values of time correspond to the period that precedes the end of ribosomal translation.

https://doi.org/10.7554/eLife.08697.009

The slow post-translational integration of the dual-topology EmrE mutants is due to the fact that a significant fraction of trajectories exhibit configurations in which some TMDs are not fully integrated at the end of translation. As a function of time, Figure 3B shows the fraction of CG trajectories in which each TMD is integrated for both the CB mutant (top) and EmrE (bottom). TMDs in the CB mutant integrate sequentially with near 100% efficiency during translation, which is consistent with the standard cotranslational model of topogenesis (c.f. Figure 1) and explains the rapid timescale for fully integrating into a multispanning topology shown in Figure 3A. In contrast, the TMDs of EmrE exhibit only partial integration, even at long times following the completion of translation. Snapshots of a typical misintegrated TMD in EmrE are shown in Figure 3B. Various experiments have indicated that such configurations with misintegrated TMDs arise due to frustration from charges placed in consecutive loops (Gafvelin and von Heijne, 1994), the strong orientational preference of a neighboring TMD (Öjemalm et al., 2012), or the weak stop-transfer efficiency of marginally hydrophobic TMDs (Moss et al., 1998). Consistent with these experimental observations, the simulations in Figure 3B find that the weakly hydrophobic TMD1 of EmrE integrates the least efficiently, followed by TMD4 which is flanked by two charged loops.

The proposed mechanism

Kinetic annealing of the end-of-translation ensemble

Analysis of the simulated CG trajectories reveals a straightforward molecular mechanism by which the multispanning topology of EmrE and its mutants is established. This mechanism, which we refer to as kinetic annealing of the end-of-translation (EOT) ensemble, is illustrated in Figure 4 and involves two basic steps. In the first step, the cotranslational integration (or misintegration) of each TMD leads to an ensemble of IMP configurations associated with the time at which ribosomal translation completes; we call this set of configurations the EOT ensemble. In the second step of the proposed mechanism, configurations in the EOT ensemble anneal toward a fully integrated topology as a function of time as loops post-translationally flip across the membrane. The rate at which the soluble loops undergo post-translational flipping is a key determinant of the fully integrated topology. We will show that this mechanism explains the unexpected elements of EmrE topogenesis observed experimentally, including the topogenic effect of C-terminal mutations (Seppälä et al., 2010).

The proposed mechanism of kinetic annealing of the EOT ensemble.

(A) Representative configurations visited by EmrE, T28R, and nR111 mutants within 1 s of simulation time following the end of ribosomal translation; this set of configurations for each mutant is called the EOT ensemble. The most-charged loop is indicated in red to highlight its position relative to the membrane. (B) The average loop positions from the EOT ensemble of each mutant is presented in terms of the fraction of configurations for which each loop occupies the cytoplasm. Adding charge to a loop biases towards a cytoplasmic position. (C) Post-translational annealing of representative EOT configurations for each mutant. Horizontal arrows indicate possible transitions as loops stochastically flip across the membrane to correct misintegrated TMDs. In each case, the position of the slowest-flipping loop (loop in red) determines the fully integrated topology by retaining its initial EOT position.

https://doi.org/10.7554/eLife.08697.013

The first step of the proposed mechanism is presented in Figure 4A and Figure 4B in greater detail. As illustrated in Figure 4A, the EOT ensemble of each mutant is determined cotranslationally as TMDs exit the translocon. Differences in loop charges in the various mutants leads to variation in the corresponding EOT ensembles, because electrostatic interactions between highly-charged loops and the translocon favor their cytoplasmic retention (Goder et al., 2004; Zhang and Miller, 2012a). Figure 4A shows representative members of the EOT ensemble for the EmrE, T28R, and nR111 mutants with the most-charged loop of each mutant highlighted in red. The EOT ensemble is defined as the set of configurations visited by a given nascent protein within 1 s of simulation time following the termination of ribosomal translation. The schematics indicate how various TMDs integrate or misintegrate to give rise to heterogeneity in the EOT ensemble of configurations, while the loops with added charges preferentially obtain cytoplasmic positions. Figure 4B further quantifies the cytoplasmic bias of charged loops by showing the EOT ensemble averaged loop positions with respect to the membrane for all five loops in each mutant, expressed as the fraction of EOT configurations with a given loop in the cytoplasm. The increased cytoplasmic localization exhibited by the L2 and L5 loops in T28R and nR111 respectively highlights the effect of adding positive charges. Similarly, the reduced cytoplasmic retention of L2 and L4 in nR111 relative to EmrE is due to the removal of charges from these loops.

The second step of the proposed mechanism is presented in Figure 4C in greater detail. For each of the three mutants, the figure schematically illustrates the post-translational kinetics of two representative configurations from the EOT ensemble. Black horizontal arrows indicate how the flipping of soluble loops across the membrane leads to transitions between intermediate configurations. The most-charged loop is again highlighted in red. Each configuration post-translationally anneals toward a fully integrated topology as loops stochastically flip across the membrane to correct the misintegrated TMDs. The soluble loops undergo flipping at different rates, with charged loops flipping more slowly. The slowest-flipping loop thus determines the fully integrated topology that is most kinetically-accessible from a given configuration in the EOT ensemble, because the other loops will more rapidly flip. The EmrE examples (Figure 4C, left) demonstrate how the equal distribution of L2 positions with respect to the membrane in the EOT ensemble leads to two different fully integrated topologies, giving rise to the dual-topology behavior. The T28R examples (Figure 4C, right) show that increasing the charge of L2, thereby biasing its cytoplasmic localization in the EOT ensemble (Figure 4B), leads to a dominant Nperi/Cperi topology. Finally, the nR111 examples illustrate how C-terminal charges can have a long-range topogenic effect by biasing the fully integrated proteins towards a dominant Ncyto/Ccyto topology.

The proposed mechanism predicts that the final topological distribution of each EmrE mutant is determined by both the distribution of configurations in the EOT ensemble and the available post-translational kinetic pathways that lead to fully integrated protein topologies. In the following, we provide detailed analysis of the simulated CG trajectories to support these elements of the proposed mechanism.

Charge mutations bias loop positions in the EOT ensemble

To investigate the first step of the proposed mechanism (Figure 4A and Figure 4B), we examine the degree to which changing the number of charges in a given soluble loop shifts the position of that loop in the EOT ensemble. Figure 5 presents the average position with respect to the membrane of each individual loop of EmrE in the EOT ensemble (blue bars); as well as the corresponding average position of each loop in the mutant for which that loop includes an additional positive charge (red bars). Loop positions are expressed as the fraction of CG trajectories in which the loop is positioned in the cytoplasm. In each case, the addition of positive charge to a soluble loop leads to an increase in its degree of cytosolic localization, as is consistent with previous simulations and experiments of single-spanning TMDs (Goder et al., 2004; Zhang and Miller, 2012a). These results support the first step of the proposed mechanism and show that interactions of the nascent protein with its translocon/ribosome/membrane environment lead to significant shifts in the EOT ensemble of configurations.

Figure 5 with 1 supplement see all
In blue, the average loop positions from the EOT ensemble of EmrE is presented in terms of the fraction of configurations for which each loop occupies the cytoplasm.

In red, the corresponding average loop position from the EOT ensemble of EmrE mutants; for each loop, the presented result is for the charge mutation associated with that loop.

https://doi.org/10.7554/eLife.08697.014

Rate of loop-flipping depends on charge mutations

To investigate the second step of the proposed mechanism (Figure 4C), we examine the molecular processes by which configurations in the EOT ensemble reach a fully integrated topology. The energetic cost for flipping a hydrophilic loop across the hydrophobic membrane increases with the hydrophilicity of the loop; the loop-flipping frequency observed during simulations is thus expected to decrease for loops with larger numbers of charges. Figure 6 shows the computed loop-flipping frequencies for each loop in the EmrE mutants. In this analysis, loop-flipping events are determined by comparing loop positions with respect to the membrane in 1-ms time intervals, as described in the Calculation of loop-flipping frequency section of the ‘Materials and methods’. The number of charges in each loop are marked with dots. As expected, highly-charged loops exhibit a decreased loop-flipping frequency. The figure also reveals that the terminal L1 and L5 loops have a lower loop-flipping frequency than the intermediate L2-4 loops. Loop-flipping events are not found to be strongly concerted, as two or more loops were observed to flip concurrently in only 0.015% of all 1-ms time intervals in which at least one loop-flipping event was observed. However, the loop-flipping frequency of a given loop is impacted by the orientation of its neighboring TMDs; on average, a loop with a single misintegrated neighboring TMD flips 1.5 times more frequently than the same loop with zero misintegrated neighboring TMDs, while a loop with two misintegrated TMDs flips 3.7 times more frequently than the same loop with zero misintegrated neighboring TMDs. Additional details on these calculations are presented in the Calculation of loop-flipping frequency section of the ‘Materials and methods’.

Loop-flipping frequencies computed from the CG trajectories.

Darker squares correspond to less frequent loop-flipping events according to the logarithmic color map. Gray dots indicate the number of positively charged residues on each loop. For each mutant, the red box indicates the slowest-flipping loop.

https://doi.org/10.7554/eLife.08697.016

The most important feature in Figure 6 is the identification of a slowest-flipping loop for each mutant (red boxes). The slowest-flipping loop typically exhibits a loop-flipping frequency that is orders of magnitude slower than the other loops, although in four cases (K3, L85R, nEmrE, and nT28R1), two loops have slow loop-flipping frequencies that are within a factor of two. The variation in loop-flipping frequencies explains the difference in kinetics in Figure 3A, such that mutants with multiple slow-flipping loops (K3) reach a fully integrated topology more slowly than mutants with a single slowest-flipping loop (EmrE, T28R) or mutants largely devoid of charge (nK3). These results confirm that the loop-flipping frequency of a given loop depends strongly on the number of charges on that loop, indicating that charge mutations can impact the determination of the slowest-flipping loop.

Position of slowest-flipping loop in EOT ensemble determines fully integrated topology

We now investigate the degree to which the position of the slowest-flipping loop in the EOT ensemble correlates with its position in the fully integrated topology. For the simulated CG trajectories, Figure 7A demonstrates strong correlation (R2 = 0.85) between the position of the slowest-flipping loop in the EOT ensemble and the corresponding position in the fully integrated configuration reached at the end of the CG trajectory. The results in Figure 7A indicate that the complexity of post-translational kinetics can be distilled to a much simpler picture in which the key parameter is the location of the slowest-flipping loop at the end of ribosomal translation. The fully integrated topology is almost completely determined at the time at which ribosomal translation ends, despite the fact that the kinetics of loop-flipping takes hundreds of seconds to complete.

Figure 7 with 1 supplement see all
Relating configurations in the EOT ensemble to the fully integrated IMP topology.

(A) Comparison of the average position of the slowest-flipping loop in the EOT ensemble to the average position of that same loop in the ensemble of fully integrated configurations. The average loop positions are presented in terms of the fraction of configurations for which each loop occupies the cytoplasm. The black dashed line indicates perfect correlation. (B) Fraction of Ncyto/Ccyto topologies for the CG trajectories that have reached full integration, excluding all trajectories for which at least one TMD is misintegrated, plotted as a function of time following the end of translation. The dashed lines indicate the fraction of Ncyto/Ccyto topologies for each mutant after all CG trajectories reach full integration, corresponding to the results from Figure 2.

https://doi.org/10.7554/eLife.08697.017

In Figure 7A, the K3 and L85R mutants deviate most significantly from the plotted correlation between the EOT ensemble and the final topology; as seen in Figure 6, these two mutants exhibit a pair of slow loop-flipping frequencies rather than a single, well-separated slowest loop-flipping frequency. For a more detailed analysis of these special cases that involve a pair of slow loop-flipping frequencies, we direct the reader to the Alternative definition of the slowest-flipping loop position for mutants with two slow-flipping loops section of the ‘Materials and methods’ and the corresponding results in Figure 7—figure supplement 1. However, we emphasize that the close agreement between the results in Figure 7A and Figure 7—figure supplement 1 indicate that our conclusions regarding the strong correlation between the EOT ensemble and the final topology are robust with respect to the details of the definition of the slowest-flipping loop.

The results in Figures 2, 3 neglect the possibility that misintegrated proteins may be degraded prior to reaching a fully integrated topology. Several bacterial proteases that degrade membrane proteins have been characterized which provides insight into the approximate degradation timescale (Dalbey et al., 2011). For example, FtsH is a membrane-embedded protease that degrades misassembled IMPs over timescales ranging from 2 min (for SecY) to 15 min (for YccA) in Escherichia coli (Ito and Akiyama, 2005), and even longer timescales for degradation have been observed in eukaryotic systems (Buck and Skach, 2005; Feige and Hendershot, 2013); very recently, FtsH was also shown to degrade undimerized EmrE on a sub-30 min timescale (Woodall et al., 2015). In comparison to the simulated trajectories (Figure 3—figure supplement 2), these degradation timescales are relatively slow, supporting the assumption that IMP integration and post-translational annealing reaches completion prior to significant degradation. Nonetheless, if degradation of EmrE occurs on faster timescales, it could potentially impact the reported topologies from the simulations. To investigate this effect, Figure 7B shows the relative fraction of Ncyto/Ccyto and Nperi/Cperi protein topologies for the CG trajectories that have reached fully integrated topologies as a function of time, excluding all trajectories for which at least one TMD is misintegrated. If it is assumed that fully integrated proteins are resistant to degradation (or that rapid dimerization following the full integration of EmrE protects the proteins from degradation [Woodall et al., 2015]), then each point in Figure 7B represents the distribution of topologies that would be observed if all misfolded proteins were uniformly degraded at the corresponding time. Data are shown for degradation times ranging from 5 s to 100 s following the end of translation; the dashed lines indicate the overall fraction of Ncyto/Ccyto topologies for each mutant after all CG trajectories reach fully integrated topologies, corresponding to the results from Figure 2. Figure 7B shows that the distribution of topologies is nearly constant with respect to degradation time, preserving the correlation between the position of the slowest-flipping loop at the end of translation and in the fully integrated topology. These results suggest that the predicted distribution of protein topologies from simulation is relatively robust with respect to possible degradation processes that occur on the same timescale as post-translational annealing.

Discussion

The results of our CG simulations support a mechanism for multispanning IMP topogenesis in which an ensemble of configurations with misintegrated TMDs undergo kinetically-controlled TMD reorientations to reach a fully integrated topology. Introducing charge mutations to the soluble loops of a multispanning IMP leads to shifts in both the distribution of loop positions in the EOT ensemble (Figure 5) and changes in the kinetics of loop-flipping events that lead to the fully integrated topologies (Figure 6). The combination of these effects is found to govern the observed distribution of fully integrated topologies in the CG simulations (Figure 7). This proposed mechanism explains the experimental finding that adding charges to any of the soluble loops of EmrE, even a loop near the C-terminus, affects the observed topology (Seppälä et al., 2010). The proposed mechanism also agrees with recent experiments that find EmrE to undergo partial topological rearrangements that correspond to the loop-flipping events described here (Woodall et al., 2015). Furthermore, the mechanism can explain deviations from the positive-inside rule if the position of the slowest-flipping loop in the EOT ensemble enforces a topology in which the majority of the positive charges are in periplasmic loops, as seen for the K3 mutant (Figure 2).

In addition to explaining existing experimental data for the topogenesis of the EmrE mutants, the proposed mechanism yields a number of new and experimentally testable predictions. A simple overarching prediction of the mechanism is that changes to the ribosome or translocon that affect the EOT ensemble may lead to significant shifts in topology. Figure 8 shows the average position of the slowest-flipping loop in the EOT ensemble after slowing translation from 24 residues/s to 6 residues/s to model the addition of cycloheximide (Goder and Spiess, 2003), removing the periplasmic positive charge from the channel, or removing the cytoplasmic negative charge from the channel (Goder et al., 2004). For single-spanning IMPs, the rate of translation and the removal of translocon charges were previously found to significantly affect TMD orientation in both simulations and experiments (Goder and Spiess, 2003; Goder et al., 2004; Zhang and Miller, 2012a). We find that slowing translation has a minimal effect on the mutants studied here, and Figure 8—figure supplement 1 confirms this finding for other translation rates. Given that these EOT loop positions are unchanged, and given that the post-translational dynamics is unaffected by the ribosomal translation rate, these results suggest that changing translation rate will not affect the final distribution of fully integrated topologies. In contrast, Figure 8 shows that removing either the cytoplasmic or periplasmic charge on the translocon significantly decreases the cytoplasmic retention of the slowest-flipping loops by increasing the periplasmic accessibility of highly charged loops. Most notably, it is found that for two of the EmrE mutants (indicated in dashed boxes) the translocon charge mutations dramatically shift the slowest-flipping loop position in the EOT ensemble from being primarily cytosolic to being primarily periplasmic, suggesting that the dominant topology for these EmrE mutants will be similarly reversed by the translocon charge mutations. These changes in IMP topology due to channel mutations are experimentally testable predictions of the proposed mechanism.

Figure 8 with 1 supplement see all
Predictions from the CG simulations in terms of changing the rate of ribosomal translation and introducing translocon mutations.

In each case, the figure presents the average position of the slowest-flipping loop in the EOT ensemble. The average loop position corresponds to the fraction of configurations for which that loop occupies the cytoplasm. The first column (dark gray) shows loop positions from the control simulation protocol. The remaining three columns show loop positions associated with the fourfold slowing of ribosomal translation (blue), associated with removal of the positive charge on the periplasmic opening of the translocon (red), and associated with removal of the negative charge on the cytoplasmic opening of the translocon (orange). The dashed boxes emphasize EmrE mutants that are predicted to exhibit dramatic inversions of the dominant topology upon translocon mutations.

https://doi.org/10.7554/eLife.08697.019

A notable aspect of the CG model is the absence of asymmetric features in the membrane or environment that favor either the Ncyto/Ccyto or Nperi/Cperi topology under equilibrium thermodynamic conditions, such as the electrostatic potential across the inner membrane of E. coli or an asymmetric distribution of charged lipids (Bogdanov et al., 2008; Vitrac et al., 2013; Bogdanov et al., 2014). In the CG model, neglecting interactions with the Sec translocon, both the Ncyto/Ccyto or Nperi/Cperi topologies are energetically equivalent and would be observed with equal probability if simulations were continued for an infinitely long time. The prediction of a dominant topology by the CG model arises from the initial distribution of configurations in the EOT ensemble (due to interactions of the nascent protein with the translocon complex) and from the available kinetic pathways that allow the configurations in the EOT ensemble to reach fully integrated topologies. We note that the changes in topology predicted in Figure 8 would be unexpected from a model in which the dominant topology of an IMP is determined by thermodynamic equilibration, since the equilibrium distribution of the protein topologies would be unaffected by transient interactions with the translocon or ribosome during initial membrane integration.

We further note that direct comparison of the experimental and simulation timescales for the kinetic annealing of misintegrated TMDs is limited by both the accuracy of the CG model as well as the neglect of external chaperone proteins, such as TRAP, TRAM, or other members of the Sec complex (Sommer et al., 2013; Zhu et al., 2013; Aviram and Schuldiner, 2014; Jung et al., 2014), that may catalyze loop-flipping. However, since the topological predictions of the proposed mechanism are primarily sensitive to which soluble loop flips most slowly—as opposed to the actual timescale of loop-flipping—we expect that the presented conclusions are relatively robust with respect to these effects. This robustness is directly illustrated in Figure 7B, which shows that the relative fraction of proteins that reach each fully integrated topology is nearly constant as a function of time.

Conclusions

In this work, we utilize a recently developed CG computational approach (Zhang and Miller, 2012a) that enables the direct simulation of Sec-facilitated membrane integration of proteins on biological timescales to investigate the topogenesis of the dual-topology EmrE protein and its mutants. In addition to demonstrating excellent agreement with the experimentally observed topologies of EmrE and its mutants (Seppälä et al., 2010), the simulations reveal a novel mechanism for the regulation of topogenesis in multi-spanning membrane proteins, in which initially misintegrated configurations of the proteins undergo post-translational annealing to reach final, fully integrated topologies. The energetic barriers associated with this post-translational annealing process enforce kinetic pathways that dictate the topology of the fully integrated proteins. The inclusion of charged residues on the soluble loops of the IMP can lead to significant changes in the distribution of fully integrated topologies by both altering the ensemble of protein configurations at the end of ribosomal translation, as well as by altering the available kinetic pathways that lead to fully integrated topologies.

This analysis leads to a number of experimentally testable predictions regarding IMP topogenesis. In particular, the results of Figure 8 predict that the mutation of charged residues near the cytoplasmic or periplasmic openings of the translocon channel will lead to significant shifts in the observed topology of several EmrE mutants. More generally, we note that any effect of channel mutations on the fully integrated IMP topology would indicate that kinetic effects during translation influence topogenesis, as suggested by the proposed mechanism. Additionally, we predict that the introduction of IMP mutations that significantly alter the EOT ensemble with respect to the cytosolic localization of the slowest-flipping soluble loop, either by introducing charge mutations or by changing TMD hydrophobicity, will influence the multispanning IMP topology.

Although the current manuscript primarily focuses on the mechanism of topogenesis in the dual-topology EmrE protein, the mechanism and simulation analysis presented here has broader implications for topogenesis in other multispanning IMPs. For EmrE and its mutants, we find that a significant fraction of the IMP configurations are misintegrated upon completion of ribosomal translation and undergo subsequent post-translational annealing to reach fully integrated topologies. In contrast, a CB mutant exhibits an essentially fully integrated ensemble of configurations at the time that ribosomal translation completes. For other IMPs, a combination of these behaviors may well be expected (Lu et al., 2000; Lambert and Prange, 2001; Kanki et al., 2002; Skach, 2009; Öjemalm et al., 2012; Bowie, 2013; Virkki et al., 2014), with some fraction of the nascent protein configurations reaching fully integrated topologies at the completion of ribosomal translation and some fraction reaching misintegrated configurations that subsequently undergo post-translational annealing. Indeed, the importance of chaperone proteins such as YidC or Sec62 that post-translationally rescue misintegrated TMDs (Sommer et al., 2013; Zhu et al., 2013; Aviram and Schuldiner, 2014; Jung et al., 2014) may be connected to this necessary process of annealing initially misintegrated IMP configurations towards fully integrated topologies. The emerging understanding of the role of the Sec translocon in regulating IMP topogenesis, as well as advances in the methodologies for probing and modifying interactions between the nascent protein and the translocon complex, hold intriguing possibilities for the prediction and control of protein folding in cellular environments.

Materials and methods

Calculation of loop-flipping frequency

Request a detailed protocol

Loop-flipping frequencies are calculated by monitoring loop positions with respect to the membrane in 1-ms time intervals. A loop-flipping event is counted if a given loop switches from a cytoplasmic to periplasmic (or vice versa) position according to the definitions in the main text and if the y-position of the reference bead in the loop (defined in the main text) is at least 4.5σ after the flip (to exclude counting the rapid flipping of loops within the translocon channel). The loop-flipping frequencies presented in Figure 6 are obtained by averaging the frequency of loop-flipping events in each trajectory. Sufficiently infrequent loop flips will not occur in every trajectory, but such events are observed in the combined ensemble of trajectories.

To examine the effect of neighboring TMDs on the loop-flipping frequency, we separately calculate the loop-flipping frequency of each loop for configurations in which zero, one, or two of the neighboring TMDs is misintegrated (discussed in the Rate of loop-flipping depends on charge mutations section of the ‘Results’).

Experimental interpretation of EmrE topology

Request a detailed protocol

In Seppälä et al. (2010), the dominant topologies of EmrE mutants are determined by measuring the growth of E. coli cells in the presence of ethidium bromide (EtBr). EtBr is toxic to E. coli, but antiparallel EmrE dimers, in which the two monomers forming the dimer have opposite topologies, confer drug resistance. EmrE dimerization can also be suppressed by including an E14D mutation. The topology of an EmrE variant with the E14D mutation can thus be inferred by coexpressing the mutant with another EmrE variant of known topology, as any resulting drug resistance (and cell growth) can be attributed to the formation of antiparallel dimers. To enable a direct comparison between the topologies measured from simulations and the experimental results, we convert the experimentally-measured cell activities from Seppälä et al. (2010) to the fraction of Ncyto/Ccyto topologies by assuming a linear relationship between cell growth and the population of antiparallel EmrE dimers. The fraction of Ncyto/Ccyto topologies is calculated as

(1) f(Ncyto/Ccyto)=A(Nperi)A(Ncyto)+A(Nperi),

where A(Ncyto) and A(Nperi) are the experimentally-measured cell activities for cells coexpressing the EmrE(Ncyto) and EmrE(Nperi) mutants, respectively. Greater cell growth in the presence of the EmrE(Nperi) mutant, which exhibits a single dominant Nperi/Cperi topology, indicates that the mutant of interest exhibits a larger fraction of the opposite Ncyto/Ccyto topology, and vice versa for growth in the presence of the EmrE(Ncyto) mutant. Experimental values for the activities of the EmrE and nEmrE mutants are taken from Figure 2 and Figure S1 of Seppälä et al. (2010), respectively; these values are used to compute the fraction of Ncyto/Ccyto topologies reported in Figure 2 of the current paper. Values for the activities of the (Nout(E14D) + Nin) and the Nout constructs from Figure 1 of Seppälä et al. (2010) are used to approximate the topology of the EmrE(Nperi) mutant, while the activities of the (Nin(E14D) + Nout) and the Nin constructs from the same figure are used to approximate the topology of the EmrE(Ncyto) mutant. Error bars are approximated via standard error propagation techniques based on Equation 1.

Alternative definition of the slowest-flipping loop position for mutants with two slow-flipping loops

Request a detailed protocol

In Figure 7A, the average position of the slowest-flipping loop relative to the membrane in the EOT ensemble is compared with the average position of that same loop in the ensemble of fully integrated configurations at the end of the CG trajectories. Four mutants (K3, L85R, nEmrE, and nT28R1), however, have two slow-flipping loops with similar loop-flipping frequencies (Figure 6), and two of these mutants (K3 and L85R) deviate most significantly in terms of the correlation in Figure 7A.

To better understand the effect of multiple slow-flipping loops on the correlation between the EOT ensemble and the final topology, the current section provides additional analysis in which a more sophisticated definition of the ‘slowest-flipping loop’ is employed for the four mutants that exhibit a pair of slow-flipping loops. Below, we present this alternative definition, which leads to a slightly better correlation between the EOT ensemble and the ensemble of fully integrated configurations, as plotted in Figure 7—figure supplement 1.

The alternative definition of the slowest-flipping loop for mutants with two slow-flipping loops is given by ϕEOT and ϕFI, which report on the average position of the two slow-flipping loops in the EOT ensemble and in the ensemble of fully integrated configurations, respectively.

The quantity ϕFI reports on the average position of the two slow-flipping loops in the ensemble of fully integrated configurations at the end of the CG trajectories. For the L85R and nEmrE mutants, the two slow-flipping loops (L2/L4 and L1/L5, respectively) reach positions on the same side of the membrane in either fully integrated topology; for these mutants, ϕFI is defined as the fraction of fully integrated configurations for which both slow-flipping loops are positioned in the cytoplasm. For the K3 and nT28R1 mutants, the two slow-flipping loops (L1 and L2) reach positions on opposite sides of the membrane in either fully integrated topology; for these mutants, ϕFI is defined as the fraction of fully integrated configurations for which L1 is positioned in the cytoplasm and L2 is positioned in the periplasm. For the nEmrE, K3, and nT28R1 mutants, ϕFI is equivalent to the fraction of CG trajectories that reach the fully integrated Ncyto/Ccyto topology, whereas for the L85R mutant, ϕFI is equivalent to the fraction of CG trajectories that reach the fully integrated Nperi/Cperi topology.

The quantity ϕEOT reports on the average position of the two slow-flipping loops in the EOT ensemble. For each mutant, ϕEOT is defined as

(2) ϕEOTL85R=0.5[f(cyto)L2+f(cyto)L4]ϕEOTnEmrE=0.5[f(cyto)L1+f(cyto)L5]ϕEOTK3=0.5[1+f(cyto)L1f(cyto)L2]ϕEOTnT28R1=0.5[1+f(cyto)L1f(cyto)L2],

where f (cyto)Li is the fraction of configurations in the EOT ensemble for which loop Li is in the cytoplasm. As for the previous definition of ϕFI, this definition accounts for the fact that the two slow-flipping loops of the L85R and nEmrE mutants reach the same side of the membrane in the fully integrated topologies, while the two slow-flipping loops of the K3 and nT28R1 mutants reach opposite sides of the membrane in the fully integrated topologies. The definition in Equation 2 additionally assumes that the post-translational annealing of misintegrated configurations in the EOT ensemble is equally rate-limited by the two slow-flipping loops.

Using these alternative definitions for the position of the slowest-flipping loop (i.e., ϕFI and ϕEOT), Figure 7—figure supplement 1 compares the average position of the slowest-flipping loop in the EOT ensemble to the average position of that same loop in the ensemble of fully integrated configurations. Having more carefully accounted for the effect of both slow-flipping loops in the K3, L85R, nEmrE, and nT28R1 mutants, this figure reveals a slight improvement in the correlation (R2 = 0.88 vs R2 = 0.85) in comparison to the results in Figure 7A.

Robustness checks for the trajectory termination criteria

Request a detailed protocol

Alternative trajectory termination criteria are tested to ensure the robustness of the simulated distribution of multispanning topologies presented in Figure 2. As a first alternative, the original set of CG trajectories are extended by 50 s, and the distribution of topologies is determined from the position of the slowest-flipping loop at the end of the extended trajectories. As a second alternative, the distribution of topologies is calculated from the subset of original CG trajectories that reach fully integrated topologies in which all beads in the TMDs are at least 20σ, rather than 4.5σ, from the origin. These robustness checks are presented in Figure 2—figure supplement 1 and exhibit excellent correlation with the results obtained using the original protocol.

Additionally, Figure 7B shows results in which the CG trajectories are terminated at a range of fixed times following the end of ribosomal translation. Again, the results using this alternative trajectory termination criterion are in good agreement with the results obtained using the original protocol (indicated in dashed lines in Figure 7B).

Cotranslational integration pathways

Request a detailed protocol

From the ensemble of CG trajectories, it is possible to examine the pathways by which individual TMDs undergo Sec-facilitated cotranslational integration. In particular, following the definitions of Cymer et al. (2014), it is possible to characterize each cotranslational TMD integration event as corresponding to either the ‘channel-sliding’, ‘interface-sliding’, or ‘in-out’ pathways. Simulation snapshots illustrating the three pathways are shown in Figure 3—figure supplement 1.

Each pathway is defined in terms of the series of intermediate states that are visited by the TMD prior to membrane integration. To characterize these intermediate states, the following geometric regions are defined (see Figure 1—figure supplement 2). The channel region is defined as that for which −2σx ≤ 2σ and −2σy ≤ 2σ, the membrane region is defined as that for which −2σx ≤ 2σ and y > 2σ, the ribosome region is defined as that for which −11σx < −2σ and −8.5σy ≤ 4.5σ, and the cytoplasm region is defined as the region outside of the ribosome for which x < −2σ. Finally, a bead is considered to overlap the LG if it is within a distance of σ to any lateral-gate bead.

We now define the four intermediate states. Intermediate state 1 (IS1) is that for which the TMD partially enters the channel; it is defined as the set of configurations for which at least two TMD beads are in the channel region and zero TMD beads are in the membrane region. Intermediate state 2 (IS2) is that for which the TMD fully spans the membrane while in the channel; it is defined as the set of configurations for which all four TMD beads are in the channel region and the two hydrophilic beads that flank the TMD occupy opposite sides of the membrane. Intermediate state 3 (IS3) is that for which the TMD accesses the membrane interior via the LG; it is defined as the set of configurations for which at least one TMD bead occupies the membrane region, the remaining three TMD beads occupy either the channel or membrane regions, and at least one TMD bead overlaps with the LG. Intermediate state 4 (IS4) is that for which the TMD accesses the cytoplasm region without accessing the channel region; it is defined as the set of configurations for which each of the four TMD beads occupies either the ribosome, membrane, or cytoplasm regions and for which at least one of the hydrophilic beads that flank the TMD occupies the cytoplasm region.

In this analysis, cotranslational TMD integration events are defined as those for which the TMD reaches a membrane integrated configuration (for which all four beads of the TMD span the membrane region and the two hydrophilic flanking beads occupy opposite sides of the membrane) before reaching a misintegrated configuration (for which both hydrophilic flanking beads occupy the same side of the membrane and for which all TMD beads and both flanking beads lie outside of the channel and ribosome regions). Using the definitions of intermediate states, the cotranslational integration pathways are defined as follows. In the ‘channel-sliding’ pathway, the TMD partially enters the channel, then crosses the LG, then fully integrates into the membrane; a trajectory thus exhibits this pathway if a TMD visits IS1, IS2, and membrane integration in chronological order and without visiting any other intermediate states. In the ‘interface-sliding’ pathway, the TMD enters the cytoplasm through the gap between the translocon and ribosome, prior to undergoing membrane integration; a trajectory thus exhibits this pathway if a TMD visits IS4 on the way to membrane integration. In the ‘in-out’ pathway, the TMD fully spans the channel prior to membrane integration; a trajectory thus exhibits this pathway if a TMD visits IS3 on the way to membrane integration without visiting IS4.

At right, Figure 3—figure supplement 1 shows the relative fraction of cotranslational TMD integration events that exhibit each of these three pathways. It is clear that the dominant cotranslational integration pathway for all four TMDs in both the EmrE and nEmrE mutants is the ‘channel-sliding’ pathway. This same pathway was also observed in the previous study of single-spanning proteins using the CG model (Zhang and Miller, 2012a) and similar configurations were observed in long-timescale atomistic molecular dynamics simulations of the early stages of cotranslational membrane integration (Zhang and Miller, 2012b). We find that only a small number of CG trajectories exhibit the ‘interface-sliding’ pathway. Finally, we note that the dominant cotranslational integration pathway is likely to depend on the IMP sequence, and the ‘channel-sliding’ behavior may be less dominant in other IMPs with less hydrophobic TMDs.

Appendix 1

The CG simulation model used in this work has been previously used to study single-spanning membrane proteins (Zhang and Miller, 2012a). Here, the key features of the model are summarized, and the minor modifications introduced in the current work are emphasized.

Coordinate system

The CG model projects all system coordinates onto a 2D plane passing through the LG of the translocon channel. The coordinate system is defined with the x-axis passing through the channel perpendicular to the bilayer and the y-axis running parallel to the bilayer as shown in Figure 1—figure supplement 2. The origin is set at the center of the channel such that negative x-values correspond to cytoplasmic locations. The LG is positioned in the positive y-direction and facilitates membrane entry.

Interactions

Interactions in the CG model are defined in terms of the reduced units ϵ and σ, where ϵ = 0.833 kT sets the energy scale and σ = 0.8 nm sets the length scale. σ is equivalent to the diameter of a CG bead. T is equal to 300 K in all simulations and k is the Boltzmann constant.

The nascent protein is treated as a freely jointed chain with connectivity enforced by a finite extension nonlinear elastic potential between consecutive beads,

(3) Ubond(r)=1/2KR02ln(1r2/R02),

where r is the distance between beads, K = 7ϵ/σ2 is the spring constant and R0 = 2σ.

Short-ranged nonbonding interactions are modeled using a cutoff Lennard-Jones (LJ) potential,

(4) ULJ(r)={4ϵij[(σr)12(σr)6]+ϵcrrcl<rrcr0rrcl,r>rcr,

where ϵij defines the strength of the interaction between bead i and j, rcl is a cutoff at low values of r, and rcr is a cutoff at large values of r. ϵcr is a constant to shift the potential to 0 at rcr. The LJ interactions among beads of the nascent protein and between beads of the nascent protein and the ribosome are chosen to be purely repulsive by setting rcr = 21/6σ and ϵij = ϵ. Weak attractive interactions are included between of the nascent protein and beads in the translocon by setting rcr = 2.5σ and ϵij = 1.5ϵ. rcl is set to 0 for all interactions between CG beads with the exception of interactions between the nascent protein and the LG in the open state, in which case rcl = σ to allow beads to exit the channel and enter the membrane.

Electrostatic interactions between charged beads are modeled using the Debye-Hückel potential,

(5) UDH(r)={σqiqj(βr)1exp(r/κ)rσqiqjβ1exp(σ/κ)r<σ,

where κ = 1.4σ approximates the Debye screening length in typical biological media, β = 1/kT, and q is the charge of a CG bead. The electrostatic potential is capped from below to prevent singularities as beads exit the channel and interact with charges on the open LG.

The effect of the implicit bilayer is modeled by calculating the solvation energy of CG beads,

(6) Usolv(x,y)=giS(x;ϕx,ψx)[1S(y;ϕy,ψy)],

where gi is the transfer free energy of bead i and S(x,ϕ,ψ) is a switching function that defines the membrane region,

(7) S(x;ϕ,ψ)=14[1+tanh(xϕ0.25σ)]  [1tanh(xψ0.25σ)].

This form of the solvation energy and switching function defines the implicit bilayer as the region where ϕx < x < ψx and y < ψy or y > ϕy, where ϕx = −2.0σ, ψx = 2.0σ, ϕy = −1.5σ, and ψy = 1.5σ. The transfer free energy, g, for each bead type is approximated from the Wimley–White hydrophobicity scale which measures water–octanol transfer free energies (Wimley et al., 1996). Values of g for the different bead types are summarized in Appendix table 1. Bead types for all EmrE variants studied in this work are listed in Appendix table 2.

Appendix table 1

CG bead charges (q) and water/membrane transfer free energies (g)

https://doi.org/10.7554/eLife.08697.021
QLVK
q0.00.00.02.0
g/ϵ2.0−4.0−2.06.0
Appendix table 2

CG bead sequences for each of the 16 EmrE mutants and the cotranslationally-biased mutant

https://doi.org/10.7554/eLife.08697.022
ProteinL1TMD1L2TMD2L3TMD3L4TMD4L5
EmrEQQQQQVVVVKQKQQLLLLQQQQQLLLLKQQQQLLLLKQQQQ
EmrE(Ncyto)QQQQQVVVVKQQQQLLLLQQQQQLLLLQQQQQLLLLKQQQK
EmrE(Nperi)QQQQQVVVVKKKQQLLLLQQQQQLLLLKKQQQLLLLQQQQQ
K3KQQQQVVVVKQKQQLLLLQQQQQLLLLKQQQQLLLLKQQQQ
T28RQQQQQVVVVKKKQQLLLLQQQQQLLLLKQQQQLLLLKQQQQ
A52KQQQQQVVVVKQKQQLLLLQQKQQLLLLKQQQQLLLLKQQQQ
L85RQQQQQVVVVKQKQQLLLLQQQQQLLLLKKQQQLLLLKQQQQ
R111QQQQQVVVVKQKQQLLLLQQQQQLLLLKQQQQLLLLKQQQK
nEmrEQQQQQVVVVQQQQQLLLLQQQQQLLLLQQQQQLLLLQQQQQ
nK3KQQQQVVVVQQQQQLLLLQQQQQLLLLQQQQQLLLLQQQQQ
nT28R1QQQQQVVVVKQQQQLLLLQQQQQLLLLQQQQQLLLLQQQQQ
nT28R2QQQQQVVVVKQKQQLLLLQQQQQLLLLQQQQQLLLLQQQQQ
nT28RQQQQQVVVVKKKQQLLLLQQQQQLLLLQQQQQLLLLQQQQQ
nA52KQQQQQVVVVQQQQQLLLLQQKQQLLLLQQQQQLLLLQQQQQ
nL85RQQQQQVVVVQQQQQLLLLQQQQQLLLLKKQQQLLLLQQQQQ
nR111QQQQQVVVVQQQQQLLLLQQQQQLLLLQQQQQLLLLKQQQQ
CB mutantL1TMD1L2TMD2L3TMD3L4TMD4L5
KQQQKLLLLQQQQQQQQQQLLLLKQQQQQQQQKLLLLQQQQQQQQQQLLLLKQQQQQQQQK
  1. Moderately hydrophobic (V), very hydrophobic (L), neutral-hydrophilic (Q) and charged (K) beads are colored as light gray, gray, black, and red respectively. Loops and TMDs are labeled.

System dynamics

The time-evolution of the nascent protein is modeled with Brownian dynamics with the equation of motion,

(8) xi(t+Δt)=xi(t)βDV(x(t))xiΔt+2DΔtηi,

where x(t)i is the position (in two dimensions) of bead i at time t, V(x(t)) is the potential energy of the system, D is the isotropic bead diffusion coefficient, Δt is the time step, and ηi is random noise drawn from a Gaussian distribution with zero mean and a variance of one. Only beads in the nascent protein are subject to Brownian dynamics. D is set to a value of 758.7 nm2/s and Δt is set to 100 ns.

LG dynamics are modeled by attempting to stochastically open or close the LG at every simulation timestep with a probability popen = kopenΔt or pclose = kcloseΔt, where

(9) kopen=1τLGexp(βΔGtot)1+exp(βΔGtot),
(10) kclose=1τLG11+exp(βΔGtot).

τLG = 500 ns is the timescale for attempting LG opening/closing and ΔGtot is the change in the free energy for opening the LG. ΔGtot depends on the presence of the nascent protein beads in the channel and is defined as

(11) ΔGtot=i=1Mgi+ΔE+ΔGemptyχempty,

where M is the number of beads occupying the translocon, ΔE is the difference between the total LG/protein LJ interactions in the closed state and total LG/protein LJ interactions in the open state, ΔGempty = 16ϵ is the free energy cost for opening the LG when there is no nascent protein in the channel, and χempty is the fraction of the channel that is empty for a given timestep. The first term promotes LG opening when hydrophobic beads enter the channel, the second term prevents the LG from closing when occluded by a nascent protein, and the third term promotes LG closing once the nascent protein exits the channel. Additional details on the development and numerical testing of the CG model are provided in Zhang and Miller (2012a).

Wall potential

A modification that appears in the current implementation of the CG model is a wall potential that prevents the nascent protein from returning to the translocon once it completes translation and diffuses a given distance away from the channel. The potential has the form

(12) Uwall(y)={1/2Kwall(y10σ)2y<10σ0y10σ,

where the spring constant, Kwall, is set to 10 ϵ/σ2. The potential is only added to the system when all beads of the nascent protein have y-positions greater than 12σ. Inclusion of the wall potential was found to avoid artifacts associated with the nascent protein interacting with the translocon long after exiting the channel. These artifacts were expected to be accentuated in the CG model due to its reduced dimensionality; nonetheless, the results are qualitatively unchanged if the wall potential is not included.

References

    1. Blobel G
    (1980) Intracellular protein topogenesis
    Proceedings of the National Academy of Sciences of USA 77:1496–1500.
    https://doi.org/10.1073/pnas.77.3.1496
    1. Schaltetzky J
    2. Rapoport TA
    (2006)
    Ribosome binding to and dissociation from translocation sites of the endoplasmic reticulum membrane
    Molecular Biology of the Cell 17:3860–3869.
    1. von Heijne G
    (1986)
    The distribution of positively charged residues in bacterial inner membrane proteins correlates with the trans-membrane topology
    The EMBO Journal 5:3021–3027.

Decision letter

  1. Gerhard Hummer
    Reviewing Editor; The Max Planck Institute of Biophysics, Germany

eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see review process). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.

Thank you for submitting your work entitled “Regulation of multispanning membrane protein topology via post-translational annealing” for peer review at eLife. Your submission has been favorably evaluated by Randy Schekman (Senior Editor) and three reviewers, one of whom, Gerhard Hummer, is a member of our Board of Reviewing Editors.

he reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.

Co-translational and post-translational steps in Sec-facilitated insertion of the EmrE integral membrane protein are studied using a coarse-grained simulation model. Consistent with experiment, the simulations produce a mixture of orientations of the protein in the membrane that is sensitive to mutations in the loops connecting the transmembrane helices. Transient topologies are found to anneal by flipping loops across the membrane during post-translational annealing. The resulting orientation distribution thus depends strongly on the electrostatic character of these loops. This work sheds light on the mechanism of membrane protein insertion, its kinetics, and the factors determining the orientation and its variability.

All three referees are overall positive but raise a number of concerns regarding the details of the model, the interpretation of the simulation results, and the presentation in the paper. I find the issues raised by the referees to be relevant. Addressing them in a revision should result in a stronger paper. Following is a brief summary of their reports that emphasizes the main points as seen from my perspective:

Reviewer #1 is concerned that the outcome of the simulations was sensitive to the somewhat arbitrary simulation termination criterion. Other concerns center on the modeling of the peptides, including the assignment of CG sequences, loop length, and the neglect of negative charge (however, the latter may not be particularly relevant for the EmrE sequence).

Reviewer #2 points out a number of relevant connections to the literature, and asks for clarifications/rectifications in the presentation of the model and the simulation results.

Reviewer #3 is concerned that the reported agreement with experiment, while visually strong, may not be quite as convincing in a more quantitative statistical analysis. The reviewer is also concerned about the factors determining the rates of loop flipping in the model and in real systems. The role of degradation (which could affect slowly relaxing EOT structures) should also be discussed in more detail.

In the following, I highlight the points in the reports I consider most relevant.

1) The terms topology, single-topology, and dual-topology should be explained to the general reader. The exact meaning of these technical terms is not immediately obvious and has caused some confusion even among the experts (cf. dual vs. antiparallel topology).

2) The simulations are terminated when a dissociated Ncyto/Ccyto or Nperi/Cperi state with 4 inserted TMDs is reached. The authors should briefly discuss (a) if the results are sensitive to variations in the termination criterion (e.g. by continuing the simulations for an additional set time, by asking for a larger distance from the origin, by running the simulations for a fixed time after initial release etc.) and (b) if there is a justification for this criterion (e.g. that subsequent changes to the topology occur on much longer time scales). It would also be useful to indicate the length of the longest simulation at which 0% of misintegrated proteins is observed.

3) Proteins with mixed orientations in the membrane appear to be the exception. The authors should thus briefly discuss whether by applying the same coarse-graining recipe to other proteins, and using the same simulation protocol, one indeed finds well defined orientations.

4) TMD1 was treated differently from the other 3 TMDs. What happens if one uses the same bead type for all TMDs? (The text states that L-type beads were used, but Supplementary file 1C lists them as V-type beads. The legend of the table also has the assignment of moderate and strong hydrophobicity of L and V beads switched vis-à-vis the text).

5) A discussion of the present results in light of earlier studies should prove insightful. In particular, Cymer et al. in a recent review (J Mol. Biol., 427:999-1022, 2015) raised questions about the simplistic model of release of TM helices in the correct topology. They noted that (a) the assembly process is driven by thermodynamics and (b) at an elongation rate of about 50 ms per amino acid, there is ample time for a chain to sample what could be a very complicated energy landscape. Here the authors indeed find that the insertion pathway and topology of a given EmrE mutant are not immutable. This connection should be discussed.

6) The pathways of insertion should be discussed in more detail. Figure 1 shows the 2D model used by Miller et al. The nascent chain can exit from the translocon either by a gap between translocon and membrane or through the translocon. Although the upper panel of Figure 1B suggests passage through the translocon, the lower panel shows the nascent chain passing from ribosome directly into the membrane. The paper does not address which of these two pathways is preferred. I think the authors should comment on this and, if possible, provide more information on the pathways observed in the simulations. This is important, because Cymer et al. (2015) suggested that the preferred path of the nascent chain might be along the translocon lateral gate rather than through the translocon as suggested by the current cartoons. The lower panel of Figure 1B suggests that the lateral gate pathway may be the preferred path.

7) A more thorough statistical analysis of the data would strengthen the paper. P-values, rank correlations, or other measures of correlation could be calculated for the data presented in Figure 2 so that the reader can understand the correlation and statistical significance of the simulation results compared to the experimental results.

8) What is the role of neighboring protein elements on the flipping frequency of a given loop? How do these correlated motions affect the flipping frequencies presented in Figure 6?

9) A more detailed discussion of the possible role of degradation would be important. Based on experimentally determined timescales of degradation (a brief comment is presented in the last paragraph of the Discussion), and the time scales of insertion and annealing observed here, can one estimate the degraded fraction and the effect on the orientation distribution?

Along the same lines, in the third paragraph of the subsection “Position of slowest-flipping loop in EOT ensemble determines fully integrated topology”, it is not immediately clear how “the fraction of CG trajectories that have reached fully integrated multispanning configurations” (Figure 7B) is related to degradation of sequences as a function of time prior to reaching a fully integrated multispanning topology. This paragraph needs to be edited for greater clarity to reveal the authors' assumptions more clearly.

[Editors' note: further revisions were requested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled “Regulation of multispanning membrane protein topology via post-translational annealing” for further consideration at eLife. Your revised article has been favorably evaluated by Randy Schekman (Senior Editor), a Reviewing Editor, and two reviewers. The manuscript has been improved but there are some remaining issues that need to be addressed before acceptance, as outlined below:

Please discuss the recent study of Woodall et al. (Nat. Comms, 6:8099; DOI: 10.1038/ncomms9099) on EmrE topologies.

Following are the full reports of the reviewers.

Reviewer #2:

The authors have, in my opinion, responded extremely well to the previous reviews. Critical questions raised have been addressed by the performance of additional simulations and analyses. The conclusions of the paper have been strengthened significantly as a result.

A new issue has arisen, however. Woodall et al. (Nat. Comms 6:8099; DOI: 10.1038/ncomms9099) have just reported (26 August) experiments directed at answering the question of how a cell can generate equal populations of two opposite topologies of EmrE. An important conclusion is that the FtsH protease of E. coli cleans up misfolded or unpaired EmrEs to create equal populations of both EmrE topologies. The authors should address these new findings, particularly in the context of their Figure 7B.

Reviewer #3:

The authors have satisfactorily addressed my concerns. I recommend the manuscript's publication.

https://doi.org/10.7554/eLife.08697.023

Author response

1) The terms topology, single-topology, and dual-topology should be explained to the general reader. The exact meaning of these technical terms is not immediately obvious and has caused some confusion even among the experts (cf. dual vs. antiparallel topology).

We agree that clarification and unambiguous usage of these terms is essential.

Throughout the manuscript, we use “topology” to indicate the orientation of a fully integrated, multispanning protein with respect to the cell membrane. A given IMP can thus become fully integrated in one of two topologies (which for EmrE correspond to Nperi/Cperi and Ncyto/Ccyto). To emphasize that “topology” refers to the fully integrated proteins, as opposed to partially integrated configurations of the proteins, we at times refer to the “fully integrated topology” – however, in both cases, we are referring to the same notion and the more verbose term is used to avoid confusion.

“Dual-topology proteins” refers to those which exhibit both possible orientations of the fully integrated IMP with respect to the cell membrane in approximately 1:1 stoichiometry (Rapp et al., Nat. Struct. Biol., 2006; Rapp et al., Science, 2007).

The term “single-topology” has been eliminated from the manuscript, as it is not needed and only creates potential confusion.

The term “antiparallel” is used only in reference to dimers of IMPs with opposite topologies.

To address this comment in the revised manuscript, the first two paragraphs of the Introduction have been modified to clarify the usage of “topology” and “dual-topology”. The text in the subsection “Alternative definition of the slowest-flipping loop position for mutants with two slow-flipping loops” has been modified to clarify the usage of “antiparallel.” We have eliminated usage of the phrase “overall topology” throughout the text to avoid confusion.

2) The simulations are terminated when a dissociated Ncyto/Ccyto or Nperi/Cperi state with 4 inserted TMDs is reached. The authors should briefly discuss (a) if the results are sensitive to variations in the termination criterion (e.g. by continuing the simulations for an additional set time, by asking for a larger distance from the origin, by running the simulations for a fixed time after initial release etc.) and (b) if there is a justification for this criterion (e.g. that subsequent changes to the topology occur on much longer time scales). It would also be useful to indicate the length of the longest simulation at which 0% of misintegrated proteins is observed.

With reference to point (a), we find that the distribution of topologies presented in Figure 2 is not sensitive to variations in the trajectory termination criteria. Per the recommendations of the referees, we tested (i) the effect of extending each original CG trajectory by 50 s and (ii) the effect of increasing the distance from the origin that was used to define the trajectory cutoff in the original protocol. These new robustness checks are presented in Figure 2–figure supplement 1 and show no significant deviations from the original results in Figure 2. Similarly, the results in Figure 7 show that terminating simulations at a set time after the end of translation does not affect the distribution of topologies in Figure 2.

To address this point in the manuscript, details of the new calculations have been added to the new section “Robustness checks for the trajectory termination criteria” in the Materials and methods.

To address point (b) in the manuscript, we provide biophysical justification for the trajectory termination criteria, adding the following explanation in the subsection “Simulation protocol”:

“The trajectory termination criteria are designed to examine the effects of the Sec-facilitated membrane integration process on EmrE topogenesis. […] Demonstration of the robustness of the reported results to the cutoff values employed in the trajectory termination criteria are provided in the Robustness checks for the trajectory termination criteria section of the Materials and methods.”

The choice of trajectory termination criteria is further justified by the robustness of the reported results with respect to the detailed cutoff values, as discussed with regard to point (a).

Finally, in regard to the reviewers’ suggestion to report the longest simulation at which 0% of misintegrated proteins is observed, we have included a table in Figure 3–figure supplement 2 that reports the amount of simulation time required for 50%, 90%, and 95% of the CG trajectories to reach fully integrated topologies for each mutant. This table provides the reader with a useful sense of the range of lengths of the trajectories. (We prefer not to report the “longest trajectory”, since this is not a well-defined property of the converged ensemble and will become arbitrarily long as the number of trajectories in the ensemble of simulated trajectories increases.)

3) Proteins with mixed orientations in the membrane appear to be the exception. The authors should thus briefly discuss whether by applying the same coarse-graining recipe to other proteins, and using the same simulation protocol, one indeed finds well defined orientations.

This is a good point. The coarse-grained model has indeed been tested for the prediction of topology in membrane proteins with a single dominant topology.

In previous work (B. Zhang and T. F. Miller III, Cell Rep., 2012), the topogenesis of single-spanning signal anchor proteins was shown to exhibit similar trends as in the experimental studies of Goder and Spiess (EMBO J., 2003). It was further shown in this previous work that a multispanning (three TMD) membrane protein with a strong positive charge bias obtained the correct, single dominant topology according to the positive-inside rule, in agreement with expectations.

In the current work, the EmrE, EmrE(Ncyto), and EmrE(Nperi) mutants, which have topologies that have been thoroughly characterized experimentally (Rapp et al., Science, 2007; Seppala et al., Science, 2010), exhibit the expected distribution of topologies. The cotranslationally-biased mutant provides another example where incorporating a strong positive charge bias leads to a single dominant topology as expected from the positive-inside rule. These examples underscore the ability of the CG model to produce well-defined topologies.

To emphasize past findings in the revised manuscript, the last paragraph of the subsection “Coarse-grained model” has been expanded to read: “Moreover, the model has been shown to accurately predict features of single-spanning IMP integration and topogenesis […] The strong agreement between simulation and experimental results presented in this work further indicates that IMP topological determinants are captured at this CG resolution.”

4) TMD1 was treated differently from the other 3 TMDs. What happens if one uses the same bead type for all TMDs? (The text states that L-type beads were used, but Supplementary file 1C lists them as V-type beads. The legend of the table also has the assignment of moderate and strong hydrophobicity of L and V beads switched vis-à-vis the text).

The representation for TMD1 used in the model is based on the hydropathy plot presented in Figure 1–figure supplement 1, which shows that the first TMD is significantly less hydrophobic than the following three.

Regardless, we agree that it is worth investigating whether this lower hydrophobicity of TMD1 impacts the reported conclusions; having performed additional simulations to investigate this point, we find that it does not have a significant impact.

To test the robustness of our results with respect to TMD1 hydrophobicity, Figure 5–figure supplement 1 presents the end-of-translation (EOT) ensemble averaged loop positions for the six EmrE mutants. The loop positions are compared between mutants with either a moderately hydrophobic (V-type beads) or very hydrophobic (L-type beads) TMD1. The positions of loops L3-L5 are largely unchanged in all cases. There is a slight increase in L2 cytoplasmic retention for all six mutants; this finding agrees with the increase in type I (Nperi/Ccyto) integration with increasing TMD hydrophobicity that was identified in previous work (B. Zhang and T. F. Miller III, Cell Rep., 2012). Given the strong correlation (Figure 7) between EOT loop positions and the final, fully integrated topology, these results suggest that increasing TMD1 hydrophobicity does not qualitatively affect the reported conclusions.

We have revised the manuscript to include this additional robustness check as follows. We have added Figure 5–figure supplement 1 as described above. Furthermore, in the subsection “EmrE protein”, we have added: “Similarly, the results of the simulations are robust with respect to changes in the modeling of TMD1 hydrophobicity (Figure 5–figure supplement 1) and loop length (Figure 3–figure supplement 3).”

As per the reviewers’ suggestion, we have also rechecked the manuscript and figures to confirm that the reported bead-types for TMD1 are correct.

5) A discussion of the present results in light of earlier studies should prove insightful. In particular, Cymer et al. in a recent review (J Mol. Biol., 427:999-1022, 2015) raised questions about the simplistic model of release of TM helices in the correct topology. They noted that (a) the assembly process is driven by thermodynamics and (b) at an elongation rate of about 50 ms per amino acid, there is ample time for a chain to sample what could be a very complicated energy landscape. Here the authors indeed find that the insertion pathway and topology of a given EmrE mutant are not immutable. This connection should be discussed.

We agree that this review raises interesting questions that should be cited in the manuscript.

As the reviewers point out, the translation process is slow, and equilibrium thermodynamic effects play a major role in determining the insertion pathway of a nascent chain, which is also noted in the review by Cymer et al. Nonetheless, some experiments, such as the work of Goder and Speiss (EMBO J., 2003), also demonstrate a clear role for kinetic effects, in addition to purely equilibrium effects, on nascent chain topogenesis.

The CG model has previously been demonstrated (B. Zhang and T. F. Miller III, Cell Rep., 2012) to correctly describe membrane integration processes that are governed by both thermodynamic and/or kinetic effects. In particular, the CG model correctly reproduces the stop-transfer efficiency of TMDs with varying hydrophobicity (T. Hessa et al., Nature, 2005), a process that is dominated by thermodynamics. The CG model also correctly reproduces the effect of translation rate and nascent chain length on signal anchor topogenesis (Goder and Spiess, EMBO J., 2003), both of which exhibit clear kinetic effects. The CG model is therefore able to model pathways dominated by both kinetics and thermodynamics as discussed in the review by Cymer et al.

To emphasize these points in the revised text, we have expanded the last paragraph of the subsection “Coarse-grained model” to read: “Moreover, the model has been shown to accurately predict features of single-spanning IMP integration and topogenesis […] and it has provided a means of understanding the competition between such effects.” Note that this revision includes explicit citation of the Cymer et al. review (2014).

6) The pathways of insertion should be discussed in more detail. Figure 1 shows the 2D model used by Miller et al. The nascent chain can exit from the translocon either by a gap between translocon and membrane or through the translocon. Although the upper panel of Figure 1B suggests passage through the translocon, the lower panel shows the nascent chain passing from ribosome directly into the membrane. The paper does not address which of these two pathways is preferred. I think the authors should comment on this and, if possible, provide more information on the pathways observed in the simulations. This is important, because Cymer et al. (2015) suggested that the preferred path of the nascent chain might be along the translocon lateral gate rather than through the translocon as suggested by the current cartoons. The lower panel of Figure 1B suggests that the lateral gate pathway may be the preferred path.

As suggested by the reviewers, we have performed additional analysis of the simulation data to characterize the cotranslational integration pathways. We find that the preferred pathway for cotranslational integration is for a TMD to partially enter the channel, then integrate into the membrane via the lateral gate.

Our analysis of the simulated trajectories examines three pathways by which by which individual TMDs undergo Sec-facilitated cotranslational integration. These three pathways are categorized according to the definitions of Cymer et al. In the “channel-sliding” pathway, the TMD partially enters the channel, then crosses the lateral gate, then fully integrates into the membrane. In the “interface-sliding” pathway, the TMD enters the cytoplasm through the gap between the translocon and ribosome, prior to undergoing membrane integration. In the “in-out” pathway, the TMD fully spans the channel prior to membrane integration.

Figure 3–figure supplement 1 in the revised manuscript shows that the dominant cotranslational integration pathway for all four TMDs in both the EmrE and nEmrE mutants is the “channel-sliding” pathway. This same pathway was also observed in our previous study of single-spanning proteins using the CG model (B. Zhang and T. F. Miller III, Cell Rep., 2012) and similar configurations were observed in previous atomistic molecular dynamics simulations (B. Zhang and T. F. Miller III, JACS, 2012). These observations agree in part with the proposal by Cymer et al. that TMDs do not need to fully enter the channel prior to lateral gate opening and membrane integration. However, we observe only a few CG trajectories in which TMDs exhibit the “interface-sliding” pathway. Finally, we note that the “channel-sliding” behavior may be less dominant in other IMPs with less hydrophobic TMDs.

To provide this analysis of the cotranslational integration pathways in the revised manuscript, we have added Figure 3–figure supplement 1 and a new section titled “Cotranslational integration pathways” to the Materials and methods section. We have also updated trajectory snapshots in Figure 1B to depict the dominant integration pathway and added to the subsection “Dual-topology proteins exhibit slow post-translational integration”: “Over 90% of the CB mutant trajectories reach the Ncyto/Ccyto topology within three seconds following the completion of translation and thus rapidly integrate as expected for the cotranslational model (Blobel, 1980; Wessels and Spiess, 1988; Sadlish et al., 2005); mechanistic features of individual TMD integration steps are discussed in the Cotranslational integration pathways section of the Materials and methods.”

7) A more thorough statistical analysis of the data would strengthen the paper. P-values, rank correlations, or other measures of correlation could be calculated for the data presented in Figure 2 so that the reader can understand the correlation and statistical significance of the simulation results compared to the experimental results.

We agree with this comment and have revised the manuscript to include improved statistical analysis.

First, for the data in Figure 2, we calculated the Pearson correlation coefficient to quantify the linear correlation between the simulation and experimental data. This calculation yielded a value of 0.92, confirming a strong linear correlation. We have also graphically illustrated this correlation in Figure 2–figure supplement 2 of the revised text. This plot is divided into quadrants – data points lying in the top right and bottom left quadrants (shaded in light gray) show where simulations and experiments both predict the same dominant topology. All mutants, with the exception of A52K, have the same dominant topology in the simulations as in the experiments within the statistical error. These measures confirm the agreement between simulation/experiment discussed in the text and strengthen the conclusions in the manuscript.

In the revised text, we have added Figure 2–figure supplement 2 and the following passage to the subsection “Simulations match experimental observations of topology”: “Figure 2-figure supplement 2 illustrates that the distribution of topologies determined experimentally and the distribution of topologies measured from the simulations are linearly correlated (Pearson correlation coefficient, r = 0.92); points lying in the two shaded quadrants of the graph correspond to proteins for which the simulations and experiments predict consistent topologies.”

8) What is the role of neighboring protein elements on the flipping frequency of a given loop? How do these correlated motions affect the flipping frequencies presented in Figure 6?

As suggested by the reviewers, we have performed additional analysis of the CG trajectories to address the role of neighboring protein elements and correlated loop motions on the loop-flipping frequencies.

Our new analysis indicates that multiple loops do not flip concurrently within the same time interval, although loops are more likely to flip when neighboring TMDs are misintegrated. These results show that while loop-flipping events are not concerted on the millisecond timescale, the loop-flipping frequency of a given loop is impacted by the orientation of its neighboring TMDs.

To examine the concurrence of loop-flipping events, we calculated the frequency with which multiple loops undergo flipping within the same 1 ms interval. It was found that pairs of loops undergo concurrent flipping with a frequency that is only 0.015% of the frequency of single loop-flipping events. Three or more loops were never observed to concurrently flip. These results indicate that loop-flipping events are not concerted on the millisecond timescale.

To examine the effect of neighboring TMDs on the loop-flipping frequency, we separately calculate the loop-flipping frequency for each loop for configurations in which zero, one, or two of the neighboring TMDs is misintegrated. On average, it was found that a loop with a single misintegrated neighboring TMD flips 1.5 times more frequently than the same loop with zero misintegrated neighboring TMDs, while a loop with two misintegrated neighboring TMDs flips 3.7 times more frequently than the same loop with zero misintegrated neighboring TMDs. These findings show that the orientation of neighboring TMDs does indeed impact the loop-flipping frequency.

To address these points in the revised manuscript, we have added to the subsection “Calculation of loop-flipping frequency”: “Loop-flipping events are not found to be strongly concerted, as two or more loops were observed to flip concurrently in only 0.015% of all 1-ms time intervals in which at least one loop-flipping event was observed. […] Additional details on these calculations are presented in the Calculation of loop-flipping frequency section of the Materials and methods.” We have also added additional details to the “Calculation of flipping frequency” section.

9) A more detailed discussion of the possible role of degradation would be important. Based on experimentally determined timescales of degradation (a brief comment is presented in the last paragraph of the Discussion), and the time scales of insertion and annealing observed here, can one estimate the degraded fraction and the effect on the orientation distribution?

Along the same lines, in the third paragraph of the subsection “Position of slowest-flipping loop in EOT ensemble determines fully integrated topology”, it is not immediately clear how “the fraction of CG trajectories that have reached fully integrated multispanning configurations” (Figure 7B) is related to degradation of sequences as a function of time prior to reaching a fully integrated multispanning topology. This paragraph needs to be edited for greater clarity to reveal the authors' assumptions more clearly.

We agree that this issue of protein degradation in this context is interesting and worth clarifying. While the exact pathway for EmrE degradation is unknown, several bacterial proteases that degrade membrane proteins have been characterized which provides insight into the approximate degradation timescale (R. Dalbey et al., 2011). Notably, FtsH is a membrane-embedded protease that degrades misassembled IMPs over timescales ranging from 2 minutes (for SecY) to 15 minutes (for YccA) in E. coli (K. Ito and Y. Akiyama, Annu. Rev. Microbiol., 2005). In comparison to the simulated trajectories (Figure 3–figure supplement 2), these degradation timescales are relatively slow, supporting the assumption that IMP integration and post-translational annealing reaches completion prior to significant degradation. Nonetheless, if it is assumed that misfolded proteins are uniformly subject to degradation, even on relatively short (i.e. 5-100 s) timescales, we provide analysis in Figure 7B to suggest that degradation still does not significantly affect the distribution of topologies predicted from the simulations.

Our analysis of the likely effects of degradation are presented in Figure 7B. This figure presents the relative fraction of Ncyto/Ccyto and Nperi/Cperi protein topologies for the CG trajectories that have reached fully integrated multispanning topologies as a function of time, excluding all trajectories for which at least one TMD is misintegrated If it is assumed that fully integrated proteins are resistant to degradation, then each point in Figure 7B represents the distribution of topologies that would be observed if all misfolded proteins were uniformly degraded at the corresponding time. The figure shows that the distribution of topologies is nearly constant with respect to degradation time, preserving the correlation between the position of the slowest-flipping loop at the end of translation and in the fully integrated multispanning topology.

To clarify these points in the manuscript, we have expanded the last paragraph of the Results.

[Editors' note: further revisions were requested prior to acceptance, as described below.]

Reviewer #2:

The authors have, in my opinion, responded extremely well to the previous reviews. Critical questions raised have been addressed by the performance of additional simulations and analyses. The conclusions of the paper have been strengthened significantly as a result.

A new issue has arisen, however. Woodall et al. (Nat. Comms 6:8099; DOI: 10.1038/ncomms9099) have just reported (26 August) experiments directed at answering the question of how a cell can generate equal populations of two opposite topologies of EmrE. An important conclusion is that the FtsH protease of E. coli cleans up misfolded or unpaired EmrEs to create equal populations of both EmrE topologies. The authors should address these new findings, particularly in the context of their Figure 7B.

We are pleased that all previous issues raised by this reviewer have been addressed in the revised manuscript. We agree that the new paper referenced by the reviewer merits discussion in the context of our work, and the manuscript has been revised to introduce several points of contact with the study, particularly in regard to the FtsH degradation results discussed in Figure 7B.

We have modified the sentence in the third paragraph of the subsection “Simulation protocol” of the revised manuscript to read: “Specifically, it is assumed that upon reaching configurations in which all of the TMDs are integrated into the membrane, the protein topology remains irreversibly fixed for all subsequent times; physical processes that may lead to this irreversibility include the dimerization of EmrE proteins to form functional channels in the membrane (Lloris-Garcerá et al., 2012) or the degradation of undimerized EmrE proteins prior to topological inversion (Woodall, Yin and Bowie, 2015).”

We have also modified the subsection “Position of slowest-flipping loop in EOT ensemble determines fully integrated topology” to read: “Several bacterial proteases that degrade membrane proteins have been characterized which provides insight into the approximate degradation timescale […] then each point in Figure 7B represents the distribution of topologies that would be observed if all misfolded proteins were ussniformly degraded at the corresponding time.”

Finally, we have added to the first paragraph of the Discussion: “The proposed mechanism also agrees with recent experiments that find EmrE to undergo partial topological rearrangements that correspond to the loop-flipping events described here (Woodall, Yin and Bowie, 2015).”

https://doi.org/10.7554/eLife.08697.024

Article and author information

Author details

  1. Reid C Van Lehn

    Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, United States
    Contribution
    RCVL, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  2. Bin Zhang

    Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, United States
    Contribution
    BZ, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    Competing interests
    The authors declare that no competing interests exist.
  3. Thomas F Miller III

    Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, United States
    Contribution
    TFM, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article
    For correspondence
    tfm@caltech.edu
    Competing interests
    The authors declare that no competing interests exist.

Funding

National Institute of General Medical Sciences (NIGMS) (Ruth Kirschtein Postdoctoral Fellowship 1F32GM113334-01)

  • Reid C Van Lehn

The funder had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

Research reported in this publication was supported by the National Institute of General Medical Sciences of the National Institutes of Health under Ruth L Kirschstein National Research Service Award award number 1F32GM113334-01. Computational resources were provided by the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the US Department of Energy under Contract No. DE-AC02-05CH11231, and by the National Science Foundation under Grant No. CHE-1040558.

Reviewing Editor

  1. Gerhard Hummer, The Max Planck Institute of Biophysics, Germany

Version history

  1. Received: May 14, 2015
  2. Accepted: September 25, 2015
  3. Accepted Manuscript published: September 26, 2015 (version 1)
  4. Version of Record published: November 6, 2015 (version 2)

Copyright

© 2015, Van Lehn et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 3,256
    Page views
  • 760
    Downloads
  • 34
    Citations

Article citation count generated by polling the highest count across the following sources: Scopus, Crossref, PubMed Central.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Reid C Van Lehn
  2. Bin Zhang
  3. Thomas F Miller III
(2015)
Regulation of multispanning membrane protein topology via post-translational annealing
eLife 4:e08697.
https://doi.org/10.7554/eLife.08697

Further reading

    1. Biochemistry and Chemical Biology
    2. Microbiology and Infectious Disease
    Rui-Qiu Yang, Yong-Hong Chen ... Cheng-Gang Zou
    Research Article Updated

    An imbalance of the gut microbiota, termed dysbiosis, has a substantial impact on host physiology. However, the mechanism by which host deals with gut dysbiosis to maintain fitness remains largely unknown. In Caenorhabditis elegans, Escherichia coli, which is its bacterial diet, proliferates in its intestinal lumen during aging. Here, we demonstrate that progressive intestinal proliferation of E. coli activates the transcription factor DAF-16, which is required for maintenance of longevity and organismal fitness in worms with age. DAF-16 up-regulates two lysozymes lys-7 and lys-8, thus limiting the bacterial accumulation in the gut of worms during aging. During dysbiosis, the levels of indole produced by E. coli are increased in worms. Indole is involved in the activation of DAF-16 by TRPA-1 in neurons of worms. Our finding demonstrates that indole functions as a microbial signal of gut dysbiosis to promote fitness of the host.

    1. Biochemistry and Chemical Biology
    2. Structural Biology and Molecular Biophysics
    Daniel Muñoz-Reyes, Levi J McClelland ... Maria Jose Sanchez-Barrena
    Research Article

    The Neuronal Calcium Sensor 1, an EF-hand Ca2+ binding protein, and Ric-8A coregulate synapse number and probability of neurotransmitter release. Recently, the structures of Ric-8A bound to Ga have revealed how Ric-8A phosphorylation promotes Ga recognition and activity as a chaperone and guanine nucleotide exchange factor. However, the molecular mechanism by which NCS-1 regulates Ric-8A activity and its interaction with Ga subunits is not well understood. Given the interest in the NCS-1/Ric-8A complex as a therapeutic target in nervous system disorders, it is necessary to shed light on this molecular mechanism of action at atomic level. We have reconstituted NCS-1/Ric-8A complexes to conduct a multimodal approach and determine the sequence of Ca2+ signals and phosphorylation events that promote the interaction of Ric-8A with Ga. Our data show that the binding of NCS-1 and Ga to Ric-8A are mutually exclusive. Importantly, NCS-1 induces a structural rearrangement in Ric-8A that traps the protein in a conformational state that is inaccessible to Casein Kinase II-mediated phosphorylation, demonstrating one aspect of its negative regulation of Ric-8A-mediated G-protein signaling. Functional experiments indicate a loss of Ric-8A GEF activity towards Ga when complexed with NCS-1, and restoration of nucleotide exchange activity upon increasing Ca2+ concentration. Finally, the high-resolution crystallographic data reported here define the NCS-1/Ric-8A interface and will allow the development of therapeutic synapse function regulators with improved activity and selectivity.