ATP-induced asymmetric pre-protein folding as a driver of protein translocation through the Sec machinery
Abstract
Transport of proteins across membranes is a fundamental process, achieved in every cell by the ‘Sec’ translocon. In prokaryotes, SecYEG associates with the motor ATPase SecA to carry out translocation for pre-protein secretion. Previously, we proposed a Brownian ratchet model for transport, whereby the free energy of ATP-turnover favours the directional diffusion of the polypeptide (Allen et al., 2016). Here, we show that ATP enhances this process by modulating secondary structure formation within the translocating protein. A combination of molecular simulation with hydrogendeuterium-exchange mass spectrometry and electron paramagnetic resonance spectroscopy reveal an asymmetry across the membrane: ATP-induced conformational changes in the cytosolic cavity promote unfolded pre-protein structure, while the exterior cavity favours its formation. This ability to exploit structure within a pre-protein is an unexplored area of protein transport, which may apply to other protein transporters, such as those of the endoplasmic reticulum and mitochondria.
Introduction
The encapsulation and compartmentalisation of cells has necessitated the evolution of machineries that conduct proteins across membranes, including for protein secretion and organellar import. Usually, protein transport occurs before the nascent protein has folded. This paper explores how the protein folding process per se may be exploited to drive protein translocation.
The bulk of protein secretion and membrane protein insertion is conducted by the ubiquitous Sec translocon. In bacteria, this comprises SecY, SecE and usually SecG, with the protein-conducting pore running through the centre of SecY. This complex can associate with either the ribosome for co-translational protein translocation (Blobel and Dobberstein, 1975) – the main pathway for nascent membrane protein insertion in bacteria (Müller et al., 2001; Jungnickel et al., 1994; Ulbrandt et al., 1997) – or with the motor protein SecA for post-translational secretion of pre-proteins (Hartl et al., 1990), which contain an N-terminal cleavable signal sequence. In the latter case, the fully synthesised pre-protein is maintained in an unfolded conformation by chaperones, such as SecB, and SecA itself (Hartl et al., 1990; Arkowitz et al., 1993). The protein must then fold during or after the translocation process.
Post-translational translocation of the unfolded pre-protein occurs through a contiguous channel formed through SecA and SecY (Figure 1A–C; Zimmer et al., 2008; Van den Berg et al., 2004; Li et al., 2016; Bauer and Rapoport, 2009). The first step in this process is the ATP-dependent ‘initiation’ phase, whereby the pre-protein is targeted via its cleavable N-terminal signal sequence to SecA, and subsequently to the SecYEG complex (Fessl et al., 2018). Transfer of the signal sequence to SecY unlocks the protein-channel, enabling the intercalation of the rest of the pre-protein (Corey et al., 2016; Hizlan et al., 2012). Next, the pre-protein is fed through the SecY channel in a process driven by both ATP and the proton-motive-force (PMF) (Brundage et al., 1990), with each mechanism requiring specific interactions with cardiolipin (Gold et al., 2010; Corey et al., 2018; Hendrick and Wickner, 1991). Finally, the signal sequence is cleaved and the pre-protein is either folded or trafficked onwards.
The SecY pore is constricted centrally by a ring of six hydrophobic residues (Van den Berg et al., 2004) (Figure 1A; red mesh), which is flanked on the cytoplasmic and exterior sides by large solvated cavities, resulting in an hourglass-like configuration (Van den Berg et al., 2004) (Figure 1A). This pore ring is capped by the so-called ‘plug’ motif (Van den Berg et al., 2004), which sits on the external face of the pore (Figure 1A; red helix), and helps to seal the channel against undesirable ion flux (Park and Rapoport, 2011). Upon SecA binding, SecY rearranges to form a more open state, involving the relocation of the plug (Zimmer et al., 2008), likely controlled by the ATPase activity of SecA (Fessl et al., 2018; Allen et al., 2016). This state still contains solvated cytoplasmic and exterior cavities, but with a profoundly changed shape of the channel (Figure 1B).
Recently, additional information on the translocation process has come from a pseudo-translocation state captured in a crystal structure of Geobacillus thermodentrificans SecYE bound to Bacillus subtilis SecA (Li et al., 2016); remodelled for Figure 1C. This structure was determined by engineering a region of pre-protein into the SecA two helix finger (2HF), such that it resides within the SecY channel. Amongst many observations, the data reveal that the central pore of SecY tightly clasps the pre-protein (Figure 1—figure supplement 1A).
Previously, we proposed a Brownian ratchet mechanism for ATP-driven protein secretion, where stochastic pre-protein diffusion is biased into directional force by the action of the SecA ATPase (Allen et al., 2016). We demonstrated that ATPase activity was able to regulate the open (ATP) and closed (ADP) states of the SecY channel; thus allowing/preventing passage of specific regions of pre-protein, such as large or charged residues, or small helical regions of secondary structure. The SecA 2HF can sense these regions, and communicate their presence back to SecA to trigger nucleotide exchange, thereby providing the means to bias Brownian motion (see Figure 8 and Video 1 in Allen et al., 2016).
Here, we extend this model to include structural changes within the channel and translocating pre-protein. We employ atomistic molecular dynamics (MD) simulation on a remodelled version of the pseudo-translocation crystal structure, along with electron paramagnetic resonance (EPR) spectroscopy and hydrogen deuterium exchange mass spectrometry (HDX-MS) on E. coli SecA-SecYEG. Together, the data reveal structural changes between the cytoplasmic and exterior cavities of SecY during the SecA ATPase cycle. These changes are instigated by a widening of the cytoplasmic cavity in the ATP-bound state. In the simulation data, this widening results in a reduced degree of pre-protein secondary structure when compared to the exterior cavity. We show this asymmetry is strongly diminished in the ADP-bound state, suggesting pre-protein transport is, in part, driven by ATP-dependent control of secondary structure formation.
Results
Asymmetric secondary structure formation of translocon-engaged pre-protein
To investigate the structural changes in the pre-protein during the translocation process, we remodelled the pseudo-translocation state crystal structure (PDB 5EUL [Li et al., 2016]) into a physiological complex (i.e. a non-fusion protein) containing SecA, SecYE and 76 residues of unfolded pre-protein, hereafter referred to as ‘SecA-SecYE-PP’; see Materials and methods section for full modelling detail (Figure 1—figure supplement 1B–C and Table 1). This structure was built into a solvated lipid bilayer, and simulated over 1 µs with SecA occupied by either ADP or ATP. The simulations were stable within the core SecY region (Figure 2—figure supplement 1) and, crucially, the SecY pore remained tightly formed around the bound pre-protein (Figure 2—figure supplement 2A).
A range of hydrogen-bonded secondary structure motifs are observed within the SecY channel, mainly α-helix in the ATP simulation and 310-helix in the ADP simulation (Figure 2—figure supplement 2B). Strikingly, however, most of this secondary structure appears to form in the exterior cavity, with little on the cytosolic side of the pore ring (e.g. Figure 2A).
Deformation analysis, which reports on local conformational flexibility in structural ensembles (see Materials and methods including Figure 2—figure supplement 3A for details), reveals a high degree of pre-protein perturbation within SecY, particularly in the cytoplasmic cavity (Figure 2B). As the pre-protein is observed to fold primarily in the exterior cavity, this suggests that Sec may be actively keeping the stretch of pre-protein in the cytosolic cavity unfolded.
Quantification of cross-translocon pre-protein secondary structure
To quantify the pre-protein folding more precisely, we carried out a series of additional simulations. Analysis of the initial trajectories (Figure 2—figure supplement 2B) suggests that a time scale of 110 ns is sufficient to kinetically sample helix formation. This matches time scales previously reported for model protein secondary structure formation (Davis et al., 2015), although it should be noted that any potential slower folding events, in the µs range, will not be sampled. We set up 36 independent simulations of ~110 ns (about 4 μs in total) per nucleotide state, varying both lipid composition and pre-protein sequence orientation (see Materials and methods for full details).
The simulations for each nucleotide state were analysed for hydrogen-bonded secondary structure in the 18 residue stretch of pre-protein through SecY. In the ATP-bound state there is >10 fold more pre-protein structure apparent in the exterior cavity compared to the cytoplasmic cavity (Figure 2C; blue data). In contrast, in the ADP-bound state no significant asymmetry with respect to pre-protein secondary structure was observed (Figure 2C; red data). This suggests an ATP-dependent and highly asymmetric influence of the translocon on the folded state of the pre-protein; set up to permit folding on the exterior side of the channel only.
Pre-protein folding asymmetry is enforced by the translocon
Further simulations were conducted to determine if this observed asymmetry was caused by the translocon, or was an intrinsic property of the chosen polypeptide sequence. The coordinates for the unfolded 18 residue stretch were extracted from the simulation and built into separate water boxes, without either SecA-SecYE or the membrane. Mild positional restraints were added to the ends of the peptide to replicate the constraints applied by the translocon and the bound signal sequence. These systems were then treated in the same way as the SecA-SecYE-PP simulations (i.e. simulated for 110 ns and each half analysed for secondary structure). The data show a small difference between the ADP and ATP states, likely a result of the mild positional restraints on the ends of the polypeptide, added to keep the substrate in a relatively transport-like conformation. This does not affect the primary outcome from this experiment, which is that no significant asymmetry in pre-protein formation is observed in either the simulation sets, be they derived from the ADP or ATP simulations (Figure 2D).
Therefore, the asymmetric secondary structure formation must be a consequence of the ATP-bound translocon, rather than the intrinsic properties of the pre-protein per se.
Unfolding of the pre-protein requires the correct signal sequence
Next, we conducted a critical control experiment, by analysing the dependence of asymmetric pre-protein structure on the interaction with a functional signal sequence; required to ‘unlock’ the complex (Corey et al., 2016). For this, we constructed a system in which the engaged pre-protein possesses a defective signal sequence (‘SSΔ4’; missing residues 5–9), but retains its strategic position in contact with the lipid bilayer (Hizlan et al., 2012; Briggs et al., 1986; McKnight et al., 1991). Analysis of this complex reveals that the ATP-induced cross-channel asymmetry is somewhat reduced (Figure 2C and Figure 2—figure supplement 4B). However, more interestingly the degree of pre-protein secondary structure is significantly higher in both cavities; particularly in the cytosolic cavity (Figure 2E). Thus, a productive interaction of the signal sequence with the ATP-associated translocon seems to be required to reduce pre-protein folding in the cytosolic cavity. This suggests that the ability of the translocon to asymmetrically influence the folding propensity of the translocating pre-protein is subject to both activation by ATP and the signal sequence.
No apparent role for specific protein-lipid interactions
Given the known dependency of protein secretion on anionic phospholipids, particularly cardiolipin (Gold et al., 2010; Corey et al., 2018; Hendrick and Wickner, 1991), we decided to look at their effects on the observed translocon-induced asymmetry. The data presented in Figure 2C combine simulations with different lipid compositions, including data in the presence of physiological concentrations of phosphatidylglycerol (PG) lipid and cardiolipin. For these, coarse-grained simulations were run on the post-1 µs ATP and ADP snapshots using the Martini force field (Monticelli et al., 2008; Marrink et al., 2007). Following 1 µs of simulation, multiple protein-lipid specific interactions were identified (Figure 2—figure supplement 4C–D), as per previous data (Corey et al., 2018). These systems were then converted back to an atomistic description (Stansfeld and Sansom, 2011) for further analyses.
Comparison of pre-protein folding in the presence and absence of PG/cardiolipin reveals no significant difference for any of the cavities (Figure 2—figure supplement 4E). Therefore, there appears to be no role for specific protein-lipid interactions in generating a cross-membrane pre-protein folding asymmetry.
Perturbed translational dynamics of water molecules within the SecY translocon
The thermodynamics and kinetics of protein folding can be affected by perturbed water dynamics (Lucent et al., 2007), and water molecules within SecY have previously been shown to exhibit reduced mobility (Capponi et al., 2015). Therefore, to examine the role of water in the pre-protein folding process, 31 structural snapshots were extracted from each of the 1 µs SecA-SecYE-PP ATP and ADP simulations. Short simulations were run to model the water dynamics accurately (see Materials and methods for details).
To measure the translational water dynamics, mean squared displacement (MSD) calculations were employed (as per Capponi et al., 2015; see Materials and methods for details). When applied to the water molecules through SecA-SecYE-PP, there is a clear pattern of perturbation throughout SecY, with the water dynamics at the centre of the channel severely restricted (Figure 3 and Figure 3—figure supplement 1). Comparison of the water molecules on either side of the pore reveal an asymmetry in the ATP-bound system – that is a higher degree of translational diffusion (disorder) in the cytoplasmic cavity than the exterior cavity – but not in the ADP-bound system.
Additional calculations were run for the translocon in a resting state, without pre-protein (Figure 3—figure supplement 2), which resembles the SecA-SecYE-PP in the ADP state, with no asymmetry between cavities. As with the analysis of pre-protein secondary structure formation, these data mark the cytoplasmic cavity of the ATP and pre-protein activated translocon as unusual in its environment. This observation suggests that the arrangement and ordering of waters within the SecY translocon may be contributing to the formation of pre-protein structure, with this effect – and hence the degree of pre-protein structure – lessened in the cytoplasmic cavity when bound to ATP.
However, it must be noted that it is unclear as to whether the perturbed water dynamics are affecting pre-protein folding/unfolding, or are a consequence of the folded state present.
Nucleotide-dependent variation in the geometry of the exterior and interior cavities of the translocon
The conformational entropy of an extended polypeptide is reduced when confined, which should favour more compact, folded states (Zhou and Dill, 2001). This raises the possibility that the ATP-driven opening of the SecY translocon (Allen et al., 2016) might have a direct impact on pre-protein secondary structure, whereby opening of the SecY cavities could promote pre-protein unfolding.
To test this, the cavity sizes of the translocon associated with either ATP or ADP were assessed throughout the 1 µs simulations. The pre-protein was removed from the translocon and the dimensions of the cavities on either side were measured for 31 structural snapshots from each simulation using the HOLE program (Smart et al., 1996). For consistency, the analyses were initiated at a set point between the residues Ile-78 and Ile-275 at the centre of the SecY pore (see Figure 1—figure supplement 1A). The cavity volumes were measured for 6.5 Å on either side of the pore, towards the cytosol (inside) or exterior (Figure 4A; Figure 4—figure supplement 1A). The relative sizes of these defined regions reveal an asymmetry between cytoplasmic and exterior cavities in the ATP-bound simulations (Figure 4B; blue data), but not in the ADP-bound simulations (Figure 4B; red data). The ATP cytoplasmic cavity is the outlier, being at least 10% larger than any of the others. This in turn reduces the degree of pre-protein contact with SecY (Figure 4C), suggesting a role in SecY-pre-protein contact in regulating secondary structure formation.
The pre-protein folding analyses (above; Figure 2) do not reveal the structural basis for asymmetric control of pre-protein folding; specifically, is it a consequence of promoted folding at the exterior, or prevention at the cytosolic cavity? The geometric analyses here suggest it is the latter, caused by a widening of the cytoplasmic cavity in the ATP-bound state.
Cavity size regulation occurs in the absence of pre-protein
To test whether the ATP-driven increase in size of the cytosolic chamber (Figure 4) is dependent on pre-protein, we analysed previously produced simulation data (Allen et al., 2016) of the SecA-SecYEG complex (PDB code 3DIN [Zimmer et al., 2008]) with ATP or ADP bound, without pre-protein. Cavity size analyses reveal the same effect as described above: that is the cytoplasmic cavity of SecY is much larger than either the periplasmic cavity in the ATP-bound state or either cavity in the ADP-bound state (Figure 4—figure supplement 1B–C). Evidently, the cytoplasmic cavity opens up in response to ATP irrespective of the presence of pre-protein.
Validation of MD simulations showing asymmetric SecY cavity size by EPR
Experimental measurement of localised and transient secondary structure changes within a highly dynamic Sec-engaged pre-protein is a challenging prospect, arguably beyond our current capabilities. Instead, we chose to experimentally analyse the observed asymmetry in SecY cavity sizes. To this end, we applied the EPR technique of double electron-electron resonance (DEER; otherwise known as PELDOR) spectroscopy, which allows detailed conformational sampling of the distances between two spin labels attached at specific positions to a molecule (Milov et al., 1981; Martin et al., 1998; Jeschke, 2012).
Firstly, we conducted MD simulations of the T. thermophilus SecYEG resting state structure (PDB 5AWW [Tanaka et al., 2015]) labelled with nitroxide-containing spin label ‘MTS’ (MTS; S-(1-oxyl-2,2,5,5-tetramethyl-2,5-dihydro-1H-pyrrol-3-yl)methyl methanethiosulfonate) (Berliner et al., 1982; Figure 5—figure supplement 1A) covalently attached to two cysteines on either sides of the SecY channel: residues Ile-58 and Ser-101 (Figure 5—figure supplement 1B; Fessl et al., 2018). The resultant simulation provides an ensemble of modelled DEER data for the labelled translocon based on the experimentally determined closed structure.
Next, we recreated the effect of SecA-ATP binding through targeted MD, using SecY from the SecA-SecYEG crystal structure (Zimmer et al., 2008) as a template, and with the targeting force constants being increased every 100 ps over 800 ps (Figure 5—figure supplement 1C). This allowed us to produce an ensemble of conformations of spin-labelled complexes based on the experimental open structure. Comparing the two ensembles, we observe that transitioning to the open structure reduces the inter-label distance substantially (Figure 5A–C). This is mainly caused by resizing of the central cavity, which permits one of the spin labels to flip inside (Figure 5B). Note that the cavity sizes are asymmetric in the open ensemble but not in the resting state (Figure 5D) exactly as for the other simulation sets (Figure 4B).
Experimental DEER spectra were then obtained by attaching spin labels at the equivalent sites on E. coli SecYEG, and measuring the inter-nitroxide distance (Martin et al., 1998; Jeschke, 2012; Haugland et al., 2016; Jeschke et al., 2006; Todd et al., 1989) for SecYEG alone; with SecA and ADP; and with SecA and the non-hydrolysable ATP analogue AMPPNP. The results reveal a strong nucleotide-dependent effect on channel conformation (Figure 5E–F and Figure 5—figure supplement 2). The most probable distance between the two sites in the SecA-SecYEG complex with ADP is ~2.7 nm, with a full width at half height (FWHH) of ~0.4 nm (Figure 5F), matching the MD data extremely well (Figure 5C). This value is identical to the distance for SecYEG alone – supporting earlier findings that the resting SecY and ADP-bound SecA-SecYEG states are very similar with respect to the channel region (Allen et al., 2016). Addition of AMPPNP, meanwhile, causes the two positions to rearrange, with the spin labels coming much closer together to give a modal distance of ~2.3 nm and broader FWHH of ~0.9 nm (Figure 5F) – again, in good agreement with the MD data (Figure 5C). These experimental distances are only an indirect measure of cavity opening, however they do provide direct experimental validation for the MD data, and confirm that ATP does indeed have a potent effect on the conformation of the channel.
HDX-MS reveals asymmetry in SecY cavities
To further explore the environment and conformational dynamics of the protein-channel during the ATPase cycle, we carried out differential HDX-MS experiments on E. coli SecA-SecYEG. HDX-MS reports on the exchange of hydrogen to deuterium in backbone amides of a protein or protein complex (Englander and Kallenbach, 1983; Konermann et al., 2011; Engen, 2009). The exchange rates are dependent on the protein solvent accessibility and hydrogen bonding (Wales and Engen, 2006). The main strength of HDX-MS is it allows non-invasive monitoring of the dynamics of an unmodified protein complex at peptide-level of resolution, enabling comparison between distinct protein states. Here, we directly compare the HDX rates between the ADP- and ATP-bound states, aiming to gain quantitative information on how nucleotide affects global conformational changes in the SecA-SecYEG complex, including networks within the SecY cavities.
We performed differential HDX-MS experiments on SecYEG saturated with SecA and either ADP or AMPPNP (hereafter referred to as ATP, for clarity). The nucleotide-bound conditions require SecA in a slight molar excess over SecY, so any data pertinent to SecA was ignored. By mapping the difference in deuterium uptake between the ATP- and ADP-bound states (ΔATP–ADP) on the crystal structure of SecA-SecYEG (PDB 3DIN), we showed significant changes primarily located in the periplasmic and cytoplasmic regions (Figure 6A and Figure 6—figure supplement 1A–C). Specific peptides were identified with significant ΔHDX (99% confidence interval), as useful reporters for changes in these two regions of interest (Figure 6B and Figure 6—figure supplement 1A).
Peptides lining the cytoplasmic cavity (CC) show deprotection induced by the presence of ATP (4 out of 4 peptides; Figure 6B–C and Figure 6—figure supplement 1B). This is consistent with a destabilisation of this region and may indicate a more ‘open’ cavity. In contrast, from the two peptides identified that line the exterior cavity (PC; not including the ‘plug’, see below), one of them revealed protection and the other a mild deprotection suggesting a more limited effect of ATP on the exterior cavity (Figure 6B,D and Figure 6—figure supplement 1C). This is consistent with the MD simulations which also showed mild dynamic changes in these regions.
Interestingly, detailed inspection of the equivalent regions in the MD simulation data reveals striking similarities with the HDX-MS results. For instance, a peptide in the cytosolic region between TMS 6–7 (CC/1), which lines the protein channel is subject to a massive increase in HDX with ATP (~50% compared to ~10% with ADP; Figure 6C and Figure 6—figure supplement 1B). This region is very highly conserved (Figure 6C and Figure 6—figure supplement 1G) and in the simulation data exhibits a substantial change in its structure between the ATP and ADP associated states (Figure 6—figure supplement 1H). The ATP-dependent mobilisation of this area does indeed appear to increase the size of the cavity and therefore could influence the folded state of the adjacent translocating pre-protein. Together, these analyses provide compelling evidence for a nucleotide-driven asymmetric resizing of the SecY cavities.
Overall the data confirm the presence of long-range conformational changes from the nucleotide binding site on SecA to the protein-channel (Figure 1C) (Fessl et al., 2018; Allen et al., 2016; Robson et al., 2007). Many of these differences occur in regions of SecYEG known to be susceptible to remote nucleotide modulation (Figure 6—figure supplement 1D). For example, the increased HDX (ATP versus ADP) observed in the ‘plug’ (extended plug loop; EPL) of SecY and the amphipathic helix (AH) of SecE are consistent with their known mobilisation in the translocon when activated by ATP (Corey et al., 2016). The contrasting decrease in HDX seen in the two largest cytosolic loops (CL) of SecY (Figure 6—figure supplement 1D) reflects the tighter association with SecA in the ATP bound state (Robson et al., 2007; Robson et al., 2009). Therefore, our HDX-MS experiments clearly support the notion of ATP-induced asymmetry in the SecYEG channel.
Pre-protein secondary structure prevents transit through SecY pore
To determine whether or not the pre-protein secondary structures observed in the simulation data impede passage through the SecY pore, a series of steered MD simulations were carried out. In these, a directional pulling force was applied to a stretch of the pre-protein through SecYE, with SecA and the rest of the pre-protein removed (Figure 7A). The pre-protein contained either an α-helical region – as sampled by equilibrium MD in the ATP- or ADP-bound states (Figure 2), and stabilised here using a distance restraint between positions i and i + 4 – or had the helix abolished using a short steered MD simulation. These substrates were then pulled through the SecY pore, and the passage time of 10 independent repeats recorded. Note that the translocation of α-helical pre-protein is unlikely to occur in the physiological system. These simulations are designed not to determine time to pass a folded region through the SecY pore, rather to report on the ease with which different pre-protein structural states might freely diffuse through.
The results show that the helix-containing substrate takes considerably longer to pass (~5–10 times slower) than the substrate with no helix (Figure 7B) in both the ATP- and ADP-bound states. Additionally, the ADP-bound state is generally more restrictive to transport than the ATP-bound state, supporting the concept of nucleotide regulated pore opening and closure (Allen et al., 2016). Although this setup clearly does not represent the physiological translocation process, the data indicate a significant barrier to the transit of folded regions of pre-protein through the protein channel.
Discussion
While it is widely accepted that pre-proteins need to be unfolded during transport (Arkowitz et al., 1993; Bonardi et al., 2011), it is unclear as to the degree of this unfolding. The general assumption is that the tertiary structure forms after the translocation process, but what of the secondary structure? Does it form prior, during or after protein transport? This study seeks to investigate this problem and explore how the folding process per se might be exploited in the mechanism of protein transport.
A high degree of secondary structure has been reported for extended pre-proteins as they exit the ribosome (Lu and Deutsch, 2005; Hardesty and Kramer, 2001), and many proteins fold further once they emerge into the periplasm (Antonoaea et al., 2008). Furthermore, previous studies have intimated that pre-protein can form secondary structure within the translocon (Gumbart et al., 2011; Zhang et al., 2017). Here, we have combined computational and structural experiments on the SecA-SecY complex, aimed at elucidating the importance of pre-protein folding for the process of protein translocation.
The data shown here reveal a propensity for the pre-protein to form secondary structure in both the cytoplasmic and exterior cavities of the channel, when associated with SecA bound to ADP. Strikingly, when ATP is bound to SecA, this propensity becomes highly asymmetric, with secondary structure strongly disfavoured in the cytoplasmic cavity (Figure 2). In this state, the cavity also exhibits a marked change in environment and increased size (Figures 3B and 4B), confirmed by EPR DEER spectroscopy (Figure 5) and HDX MS (Figure 6). This also causes a decrease in contact between the translocon and pre-protein (Figure 4C). These observations suggest a specific role for ATP (and its hydrolysis to ADP) in regulating the structure of the translocating pre-protein: ATP binding causes the SecY cytoplasmic cavity to expand, reducing the amount of secondary structure, relative to the amount formed in the exterior cavity.
Steered MD simulations confirm that regions of secondary structure are much less likely to pass the narrow SecY pore (Figure 7). Taken together, these effects could contribute significantly to any translocation model that incorporates pre-protein diffusion within the channel, including full diffusion-based models (Allen et al., 2016) and hybrid diffusion/power stroke models (Bauer et al., 2014). Interestingly, the data also argues against full power-stroke based models (Erlandson et al., 2008a), as the pre-protein conformational changes observed here are apparently unrelated to motions of the 2HF.
The results described here allow us to extend our previously proposed ‘Brownian ratchet’ model for protein translocation (Allen et al., 2016), outlined in Figure 8. The initiation process involves ATP dependent activation of SecYEG by SecA and the pre-protein signal sequence (Fessl et al., 2018; Corey et al., 2016; Hizlan et al., 2012; Gold et al., 2013), and results in an intercalated diffusible mature region of pre-protein (Figure 8). When the translocon is bound to ADP, unstructured polypeptide in this region can flow freely within the channel, which explains why ‘backsliding’ can occur in the ADP bound state (Erlandson et al., 2008b). Meanwhile, secondary structure can equally form in the exterior or cytosolic cavities, producing no bias for forward or backward translocation. However, the appearance of secondary structure in the cytosolic cavity promotes nucleotide exchange, ADP for ATP (Allen et al., 2016) (Figure 8, bottom left). As we show in this study, the ATP bound translocon strongly favours secondary structure formation in the exterior cavity, but not so in the cytosolic cavity. Given that the folded domains cannot backtrack, this favours the forward diffusion of polypeptide (Figure 8, upper middle). Successive rounds of hydrolysis and exchange continue until transport is complete and the signal sequence is removed by proteolysis (Figure 8, centre right). Since, during translocation the steady-state complex is dominated by the ATP-bound state (Robson et al., 2009) (Figure 8, green box), the observed asymmetry is maintained throughout (upper states predominate); thus preventing backsliding and favouring productive transport.
This new mechanism for forward transport through the channel, along with additional stimulation from the trans-membrane PMF (Brundage et al., 1990), could achieve the required very fast rates required for efficient protein secretion. This mechanism would also explain observed differences in the rates of transport (Fessl et al., 2018; Liang et al., 2009; Tomkiewicz et al., 2006) and stoichiometry of ATP: pre-protein (Tomkiewicz et al., 2006; van der Wolk et al., 1997), because they would be dependent on the variable propensities of different protein sequences to form secondary structures. Interestingly, it could also explain the higher translocation efficiency of the prl mutants (Bieker et al., 1990), which have been shown to exhibit a more open conformation (Corey et al., 2016), thus favouring the ATP-bound conformation.
Presently, it is only possible to simulate time scales (μs) much shorter than the physiological rates of translocation (40 aa sec−1 [Fessl et al., 2018]). Nevertheless, the kinetically fast event we sample – the formation of pre-protein structure – may well be critical for the transport mechanism.
It should be noted that the identified secondary structures might not reflect the final structure of the transported protein, especially for the β-barrel proteins of the outer membrane. However, they may correspond to folding intermediates proposed to be important for translocation competence through the Sec machinery (Tsirigotaki et al., 2018). The reliance on protein folding to aid the ratcheting process might predict that an inherently disordered protein substrate would fail to transport, as indeed it does in all of the bacterial and eukaryotic counterparts (Sec61; both yeast and mammalian [Gonsberg et al., 2017]). In cases where the propensity for folding is too great, such that it cannot be prevented as described, then other mechanisms could prevail. The E. coli thioredoxin DsbA does just that and is transported across the membrane co-translationally (Huber et al., 2005). Other fast folding proteins could be selected for secretion through the TAT machinery capably of transporting unfolded proteins (Berks, 2015).
Given the universal nature of secondary structure in polypeptides, such conformational control could be central to many other protein translocating machines, including chaperonins, organellar protein importers and the secretion systems of bacterial pathogens.
Materials and methods
Modelling the SecA-SecYEG-PP complex
Request a detailed protocolThe bulk of the data here are collected from a SecA-SecYEG-PP complex in an ATP or ADP state. Models for the simulations were built using chains A, Y and E of the crystal structure 5EUL (Li et al., 2016) as starting coordinates. Several small missing loops were added using the Modeller program (Sali and Blundell, 1993), detailed in Table 1. Some larger loops were added manually based on similar structures of SecA-SecYEG (PDB code 3DIN; Zimmer et al., 2008). These include HWD 1 and 2, two periplasmic SecY loops and a cytoplasmic SecY loop. The details of these loops are given in Table 1, and the loops are shown in Figure 1—figure supplement 1B. We reasoned that this would produce a better model than loops built simply based from sequence.
The original structure has a non-physiological arrangement of the SecA 2HF and bound pre-protein substrate, where these are represented by a single, continuous chain. To correct for this, we remodelled the tip of the 2HF based on 3DIN, to generate a physiological 2HF conformation, which is stable during simulation (Figure 2—figure supplement 1). In the SecA NBS, the ADP-BeFx molecule was replaced with either ADP or ATP (Piggot et al., 2012). The substrate was extended in an unfolded conformation through the SecA ring, via known crosslinking sites (Bauer and Rapoport, 2009). Substrate building was done with PyMOL (Delano, 2002).
Molecular dynamics simulations
Request a detailed protocolAll simulations were run using GROMACS 4.6.4 or 5.1.2 (Berendsen et al., 1995). Simulations were built using the SecA-SecYEG-PP complex described above, or alternatively the M. jannaschii SecYEβ crystal structure (Van den Berg et al., 2004) or the T. thermophilus SecYEG crystal structure (Tanaka et al., 2015).
In most cases, the protein and solvent atoms were described in the OPLS all-atom force field (Jorgensen et al., 1996), with the simulations being run in an OPLS united-atom POPC membrane (Ulmschneider and Ulmschneider, 2009). The protein-membrane structures were built into simulation boxes with periodic boundary conditions in all dimensions and solvated with explicit SPC water and sodium and chloride ions to a neutral charge and concentration of 0.15 M. For the 1 µs SecA-SecYE-PP and M. jannaschii SecYEβ simulations, the systems were energy minimized using the steepest descents method over 2 × 5000 steps, then equilibrated with positional restraints on heavy atoms for 1 ns in the NPT ensemble at 300 K with the Bussi-Donadio-Parrinello thermostat (Bussi et al., 2007) and semi-isotropic Parrinello-Rahman pressure coupling (Parrinello and Rahman, 1981; Nosé and Klein, 1983).
Production simulations were run without positional restraints with 2 fs time steps over 500–1 µs on the UK HPC facility ARCHER, using time awarded by HECBioSim. Bond lengths were constrained using the LINCS method. Non-bonded interaction cut-offs were calculated using the Verlet method, with the neighbour search list updated every 20 steps. Long-range electrostatic interactions were calculated using the particle mesh Ewald method and a cut-off of 1.0 nm was applied for van der Waals and short range electrostatic interactions. All simulations reached a steady state as judged by their root-mean-squared-deviation from the starting structure (Figure 2—figure supplement 1).
For some of the folding data, membranes consisting of 63% POPE, 32% POPG and 5% cardiolipin in a coarse-grained description were built around the protein using the MemProtMD program (Stansfeld et al., 2015). Following 1 µs simulation, snapshots were obtained in which multiple acidic lipids were seen to bind to the complex (Figure 2—figure supplement 4C–D). These were then converted to an atomistic description (Stansfeld and Sansom, 2011), and the systems simulated using the Charmm36 force field (Best et al., 2012).
Deformation analysis on pre-protein
Request a detailed protocolDeformation energies were calculated from the SecA-SecYEG-PP ADP and ATP MD trajectories. Deformation analysis based on normal mode (vibrational) analysis provides a measure for the amount of local flexibility in the protein structure - that is atomic motion relative to neighbouring atoms. We calculated deformation energies and atomic fluctuations of the first three non-trivial modes in Bio3D (Hinsen, 1998; Grant et al., 2006) and visualized the results in PyMOL (Delano, 2002).
Targeted MD
Request a detailed protocolFor the systems built from T. thermophilus SecYEG (Tanaka et al., 2015), simulations were initially run of the resting state of the channel (i.e. SecYEG alone), with MTS spin labels attached to cysteine residues engineered at positions 58 and 101. MTS parameters were produced using ACPYPE (Sousa da Silva and Vranken, 2012). MD simulation of the MTS in water produces the appropriate geometries (Figure 5—figure supplement 1A).
Over 110 ns simulation, distance analyses were carried about between the centre-of-mass of the ON groups of each spin label. After this, targeted MD was carried out on the SecY subunit to sample the conformational arrangement of 3DIN SecY (chain C [Zimmer et al., 2008]). Force constants were applied to the CA atoms, and increased at 100 ps intervals from 10 to 1000 kJ−1 mol−1 nm−2 restraints over 1 ns. Following this, the systems were simulated for a further 70 ns, with the positional restraints eased to 10 kJ−1 mol−1 nm−2. This therefore allowed us to simulate a SecA-bound state of the complex, without needing SecA present. Upon opening, the MTS label was able to flip into the centre of SecY. This process was modelled through a manual switching of rotamer state for the backbone cysteine. Again, distance analyses were performed between the spin labels.
Modelling a defective signal sequence
Request a detailed protocolTo investigate the effect of signal sequence interaction on the system, we modelled in a known defective signal sequence (Hizlan et al., 2012; Emr et al., 1980), where a conserved four hydrophobic residue stretch is removed (here, KKTAIAIAVALAGFATVAS). We modelled this in PyMOL from a 1 µs snapshot of the ATP-bound SecA-SecYE-PP system. To reduce input bias, we allowed the N-terminal KK to remain in contact with the lipid phosphate groups and simply removed the four residues from the protein. We then repositioned the flanking residues slightly to allow bond formation, and relaxed the system using energy minimization. We then simulated the complex out to >400 ns for further analyses.
Secondary structure formation analyses
Request a detailed protocolFormation of secondary structure was determined both using visual analysis and quantified using the dictionary of secondary structure of proteins (DSSP) algorithm (Kabsch and Sander, 1983; Joosten et al., 2011). This was implemented using the Gromacs utility do_dssp, which provides an estimation of secondary structure for each reside based on dihedral angle, including ‘no structure’, ‘α-helix’, ‘310-helix’ and ‘hydrogen-bonded turn’.
Extending simulations to investigate secondary structure formation
Request a detailed protocolTo provide broad conformational sampling and reduce any input bias from the starting configurations, multiple simulations were run from different starting coordinates. For this, new systems were modelled based on snapshots of the ADP- and ATP-bound systems at 1 µs. To abolish all the secondary structure within the pre-protein, the nine residues on each side of the central constriction were substituted with glycine using SCWRL4 (Krivov et al., 2009). Following steepest descents energy minimization and 1 ns equilibration, production runs were carried out until all of the secondary structure within this region had been abolished, according to DSSP analysis (Joosten et al., 2011); ~110 ns for the ATP-bound complex and ~80 ns for the ADP-bound complex (Figure 2—figure supplement 3). For each simulation, three time points (70, 80 and 110 ns for ATP and 50, 65 and 80 ns for ADP, with only the 50 and 65 ns simulations run for the tandem sequence (see below) simulations) were chosen which displayed no structure at all according to DSSP analysis (Figure 2—figure supplement 3). This represents a reasonably broad configurational sampling, and prevents the analyses being too biased by a specific starting configuration.
For each time point, an 18 residue region of pre-protein through the SecY channel was modelled in one of two different ways, with the first nine residues (ELERQHTFA) duplicated either in a mirrored (ELERQHTFAAFTHQRELE) or tandem (ELERQHTFAELERQHTFA) manner, using SCWRL4 (Krivov et al., 2009). These were then minimised using steepest descents, equilibrated for 1 ns and simulated for between 90 and 110 ns. The final 75 ns of each simulation run were then analysed for secondary structure. Simulations were run on EPCC’s Cirrus HPC Service.
Due to sampling restrictions, we were realistically only able to analyse a single nine residue sequence. We opted to use the sequence from the original crystal structure, as this has already been shown to be stable in the SecY channel, thereby reducing the potential for artefacts arising from changes of sequence. Additionally, the sequence has already been shown to fold within a reasonable time frame (Figure 2—figure supplement 2B).
Water dynamics simulations
Request a detailed protocolTo model the water dynamics in and around the translocon, a series of simulations were run using structural snapshots of the 1 µs SecA-SecYE-PP ATP simulation from Figure 2 as starting points (500, 502, 504, 506, 508, 510, 600, 602, 604, 606, 608, 610, 700, 702, 704, 706, 708, 710, 800, 802, 804, 806, 808, 810, 900, 902, 904, 906, 908, 910, 1000 ns). Water dynamics simulations were run in the NVE ensemble, meaning a constant number of particles, volume and energy were maintained throughout. The advantage of this is to avoid introduction of artificial perturbations through use of a thermostat or barostat. Simulations were run with a 1 fs time step, writing coordinates every 5 fs. The Verlet cutoff scheme was used with a buffer size of 0.001 kJ mol−1 ps−1 to achieve proper energy conservation, as monitored through following the energy of the simulations.
Mean squared displacement calculations
Request a detailed protocolTo model the translational dynamics of the waters in the system, mean squared displacement (MSD) calculations were employed. MSD is a common statistical mechanics measure of translational motion. It follows the progression of an atom in relation to a reference position as a product of time, revealing the extent of its exploration of 3D space. Analyses were carried out as described previously (Capponi et al., 2015). Briefly, a box of 5×5×12 nm was built around SecY, with the geometric centre at the pore ring. This box was subdivided into 24 slices 0.5 nm thick, and the waters from each slice were analysed separately. The MSD of the waters in each slice was plotted, and the last 25 ps was fitted to power law equations (Equation 1).
Where k is a fitting parameter, and the exponent A provides information on the molecular diffusivity. From this, the A value was plotted for each slice.
HOLE
Request a detailed protocolTo analyse the SecY cavity sizes, snapshots were taken from the 1 µs SecA-SecYE-PP ATP and ADP simulations at 500, 502, 504, 506, 508, 510, 600, 602, 604, 606, 608, 610, 700, 702, 704, 706, 708, 710, 800, 802, 804, 806, 808, 810, 900, 902, 904, 906, 908, 910 and 1000 ns. In addition, snapshots were taken at the same time points from simulations previously published (Allen et al., 2016) of the SecA-SecYEG complex without pre-protein, with ATP or ADP bound. Snapshots were also taken from simulations of T. thermophilus SecYEG (Tanaka et al., 2015) in a resting and open conformation (see Targeted MD section above).
For each snapshot, the centre of the SecY channel was initially defined as the geometric centre of the pore ring residues, and this was used to seed the HOLE calculations. The correct siting of the cavity was determined with visual inspection using VMD (Humphrey et al., 1996). For the successful calculations, data were extracted for 6.5 Å on either side of the pore ring. The area under the curves were computed and integrated using the trapezoidal rule with partitions of 0.5 Å.
EPR DEER
Request a detailed protocolProteins for spin-labelling were produced as described previously (Collinson et al., 2001). Complexes were labelled by incubation with a 20 fold excess of MTS and 40 μM cardiolipin at room temperature for 30 min, in buffer A (20 mM Tris pH 8, 130 mM NaCl, 10% glycerol with 0.02% by volume n-Dodecyl-β-D-Maltoside [DDM]). Following this, excess (unbound) cardiolipin and MTS were removed using size exclusion chromatography with a superose 6 10/300 column at 4°C. Purified fractions were concentrated and stored at −80°C. SecA was purified as described previously (Gold et al., 2007) and stored in buffer B (20 mM Tris pH 8, 50 mM KCl, 2 mM MgCl2).
The samples for EPR were prepared using these stock solutions. The proteins were mixed and diluted into the deuterated buffer C (20 mM Tris pH 8, 50 mM KCl, 2 mM MgCl2 in D2O) with the required amount of DDM to keep the concentration at 0.02% and MgCl2 was kept constant at 2 mM. SecYEG was added to 10 μM, with SecA in 5-fold excess and nucleotide (ADP or the non-hydrolysable AMPPNP) added to 2 mM final concentration. Glycerol-d8 30% by volume was added to give a final volume of 60 μl.
The mixtures were loaded into 3 mm OD suprasil EPR tubes and flash frozen. The samples were loaded into the high-power Q-band (34 GHz, ER 5106QT-2w cylindrical resonator) Bruker Elexsys E580 spectrometer at the University of St Andrews.
The double electron electron resonance (DEER) experiments were carried out at 50 K with a four pulse, deadtime free sequence (Milov et al., 1981; Martin et al., 1998; Jeschke, 2012). The pump pulse (14 ns) was set to the maximum of the echo detected field swept spectrum of the nitroxide. The observer sequence was set 80 MHz offset from the pump frequency and used 32 ns pulses. The shot repetition time was optimised to measure at least 70% of the echo height (3 ms) and for each scan, 30 shots were taken at each time point. Furthermore, a two-step phase cycle with nuclear modulation averaging (5 steps of 24 ns separation with an initial separation of 200 ns) was applied. Using DeerAnalysis2016, the data had 800 ns cut to remove artefacts at the end of the data, with a 129 ns zerotime. It was found necessary to use a fifth-order polynomial background correction (Jeschke et al., 2006). Distance distributions were calculated by Tikhonov Regularisation with the Regularisation parameter determined by the programme.
Note that for a fully spin-labelled SecYEG, the modulation depth would be expected to be about 0.3 rather than values in the approximate range 0.1 to 0.15 (Figure 5E) and thus the sample may not have been fully labelled. The differences between conditions were not experimentally reproducible with alternative samples and therefore may be due to small changes in the experimental conditions. The higher frequency of the measured dipole-dipole coupling in the case of spin-labelled SecYEG with SecA and AMPPNP indicates an alteration of the cavity dimensions. More details and experiments are shown in the Supporting Information (Figure 5—figure supplement 2).
The research data supporting the DEER data can be accessed at https://doi.org/10.17630/0fedaeec-7e27-4876-a6d1-cda2d3a6799c
Hydrogen deuterium exchange mass spectrometry (HDX-MS)
Request a detailed protocolPeptide identification and peptide coverage were optimized for SecYEG and SecA using un-deuterated controls. The sample workflow for HDX-MS involved mixing 10 μM of SecYEG and 15 μM of SecA and incubating for 10 min on ice. 5 μl of SecYEG-SecA complex was diluted into 95 μL of equilibration buffer C (20 mM Tris pH 8, 50 mM KCl, 2 mM MgCl2 and 0.02% DDM in H2O) or with deuterated buffer C (20 mM Tris pH 8, 50 mM KCl, 2 mM MgCl2 and 0.02% DDM in D2O) at 25°C. For experiments analysing the dynamics of the SecYEG-SecA in presence of either AMPPNP or ADP, 1 mM of nucleotide was added to the protein mixture and to buffer A.
The deuterated samples were then incubated for 0.25, 1, 5, and 30 min at 25°C. All samples were quenched with 100 μL of quench buffer (0.8% formic acid and 0.1% DDM). The mixed solution was at pH 2.4. Peptide digestion was then performed on-line using a Enzymate online digestion column (Waters) in 0.1% formic acid in water at 20°C and with a flow rate of 200 μL/min. The pepsin column was washed with pepsin wash recommended by the manufacturer (0.8% formic acid, 1.5 M Gu-HCl, 4% MeOH). To prevent significant peptide carry-over from the pepsin column, a blank run was performed between sample runs.
Peptide fragments were then trapped using an Acquity BEH C18 1.7 μM VANGUARD chilled pre-column for 3 min. Peptides were then eluted into a chilled Acquity UPLC BEH C18 1.7 μM 1.0 × 100 mm column using an 8–40% gradient of 0.1% formic acid in acetonitrile at 40 μL/min. Peptide fragments were ionized by positive electrospray into a Synapt G2-Si mass spectrometer (Waters). For high-energy acquisition of product ions, a 20 to 30 V trap collision energy ramp was used to acquire the MSE data. Mass accuracy correction was performed using Leucine Enkephalin (LeuEnk) as a lock mass and iodide was used for mass spectrometry calibration.
All deuterated time points and un-deuterated controls were carried out in triplicate. MSE data from un-deuterated controls samples of SecYEG and SecA were used for sequence identification using the Waters ProteinLynx Global Server 2.5.1 (PLGS) and filtered using DynamX (v. 3.0). Filtering parameters used were a minimum and maximum peptide sequence length of 4 and 25, respectively, minimum intensity of 1000, minimum MS/MS products of 2, minimum products per amino acid of 0.2, and a maximum MH +error threshold of 5 ppm. Furthermore, all the spectra were examined and checked visually and only peptides with a high signal to noise ratios were used for HDX-MS analysis.
Steered MD
Request a detailed protocol1 µs snapshots were taken from the SecA-SecYE-PP ADP- and ATP-bound simulations, in which four pre-protein residues in the exterior SecY cavity had formed an α-helical configuration (Figure 2A). From these, SecA was removed and the backbone of the pre-protein was broken at a position equivalent to proOmpA residue 25, separating the signal sequence and rest of the pre-protein. All residues after the equivalent of proOmpA residue 40 were removed (Figure 7A), and simulated for 5 ns. The α-helix was then either forced into an unstructured conformation using steered MD with a pulling force of 1000 kJ mol−1 nm−1 for 36 ps on the C-terminal residue, or stabilised using a constraint of 0.2–0.25 nm between the backbones of the i and i + 4 helix residues, followed by five ns unbiased simulation. Then, 10 independent steered MD simulations were run for each state, where the substrate was pulled from the C-terminal residue in a z-axis direction using a force constant of 600 kJ mol−1 nm−1. For each repeat, the time take for the helical/non-helical region to cross the pore was measured, based on the centre-of-mass distance between the region and the SecY pore residues, Ile-78, Ile-183, Ile-275 and Ile-404 (G. thermodentrificans numbering).
Data availability
All data generated during this study are included in the Figures of the mansucript. EPR data is available at https://doi.org/10.17630/0fedaeec-7e27-4876-a6d1-cda2d3a6799c.
-
University of St AndrewsEPR data from ATP-induced asymmetric pre-protein folding as a driver of protein translocation through the Sec machinery.https://doi.org/10.17630/0fedaeec-7e27-4876-a6d1-cda2d3a6799c
References
-
GROMACS: A message-passing parallel molecular dynamics implementationComputer Physics Communications 91:43–56.https://doi.org/10.1016/0010-4655(95)00042-E
-
The twin-arginine protein translocation pathwayAnnual Review of Biochemistry 84:843–864.https://doi.org/10.1146/annurev-biochem-060614-034251
-
A novel reversible thiol-specific spin label: papain active site labeling and inhibitionAnalytical Biochemistry 119:450–455.https://doi.org/10.1016/0003-2697(82)90612-1
-
Optimization of the additive CHARMM all-atom protein force field targeting improved sampling of the backbone φ, ψ and side-chain χ(1) and χ(2) dihedral anglesJournal of Chemical Theory and Computation 8:3257–3273.https://doi.org/10.1021/ct300400x
-
The sec and prl genes of Escherichia coliJournal of Bioenergetics and Biomembranes 22:291–310.https://doi.org/10.1007/BF00763169
-
Canonical sampling through velocity rescalingThe Journal of Chemical Physics 126:014101.https://doi.org/10.1063/1.2408420
-
Unlocking the Bacterial secy transloconStructure 24:518–527.https://doi.org/10.1016/j.str.2016.02.001
-
SoftwareThe PyMOL Molecular Graphics SystemThe PyMOL Molecular Graphics System.
-
Analysis of protein conformation and dynamics by hydrogen/deuterium exchange MSAnalytical Chemistry 81:7870–7875.https://doi.org/10.1021/ac901154s
-
Hydrogen exchange and structural dynamics of proteins and nucleic acidsQuarterly Reviews of Biophysics 16:521–655.https://doi.org/10.1017/S0033583500005217
-
Analysis of polypeptide movement in the SecY channel during SecA-mediated protein translocationJournal of Biological Chemistry 283:15709–15715.https://doi.org/10.1074/jbc.M710356200
-
Allosteric Regulation of SecAJournal of Biological Chemistry 282:17424–17432.https://doi.org/10.1074/jbc.M702066200
-
The dynamic action of SecA during the initiation of protein translocationBiochemical Journal 449:695–705.https://doi.org/10.1042/BJ20121314
-
The Sec61/SecY complex is inherently deficient in translocating intrinsically disordered proteinsJournal of Biological Chemistry 292:21383–21396.https://doi.org/10.1074/jbc.M117.788067
-
Free energy of nascent-chain folding in the transloconJournal of the American Chemical Society 133:7602–7607.https://doi.org/10.1021/ja2019299
-
Folding of a nascent peptide on the ribosomeProgress in Nucleic Acid Research and Molecular Biology 66:41–66.https://doi.org/10.1016/S0079-6603(00)66026-9
-
Tuning the properties of nitroxide spin labels for use in electron paramagnetic resonance spectroscopy through chemical modification of the nitroxide frameworkElectron Paramagnetic Resonance 25:1.
-
SecA protein needs both acidic phospholipids and SecY/E protein for functional high-affinity binding to the Escherichia coli plasma membraneThe Journal of Biological Chemistry 266:24596–24600.
-
Analysis of domain motions by approximate normal mode calculationsProteins: Structure, Function, and Genetics 33:417–429.https://doi.org/10.1002/(SICI)1097-0134(19981115)33:3<417::AID-PROT10>3.0.CO;2-8
-
VMD: visual molecular dynamicsJournal of Molecular Graphics 14:33–38.https://doi.org/10.1016/0263-7855(96)00018-5
-
DeerAnalysis2006—a comprehensive software package for analyzing pulsed ELDOR dataApplied Magnetic Resonance 30:473–498.https://doi.org/10.1007/BF03166213
-
DEER distance measurements on proteinsAnnual Review of Physical Chemistry 63:419–446.https://doi.org/10.1146/annurev-physchem-032511-143716
-
A series of PDB related databases for everyday needsNucleic Acids Research 39:D411–D419.https://doi.org/10.1093/nar/gkq1105
-
Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic LiquidsJournal of the American Chemical Society 118:11225–11236.https://doi.org/10.1021/ja9621760
-
Improved prediction of protein side-chain conformations with SCWRL4Proteins: Structure, Function, and Bioinformatics 77:778–795.https://doi.org/10.1002/prot.22488
-
Bacterial Sec protein transport is rate-limited by precursor length: a single turnover studyMolecular Biology of the Cell 20:4256–4266.https://doi.org/10.1091/mbc.e09-01-0075
-
Folding zones inside the ribosomal exit tunnelNature Structural & Molecular Biology 12:1123–1129.https://doi.org/10.1038/nsmb1021
-
The MARTINI force field: coarse grained model for biomolecular simulationsThe Journal of Physical Chemistry B 111:7812–7824.https://doi.org/10.1021/jp071097f
-
Application of ELDOR in electron-spin echo for paramagnetic center space distribution in solidsFizika Tverdogo Tela 23:975–982.
-
The MARTINI Coarse-Grained Force Field: Extension to ProteinsJournal of Chemical Theory and Computation 4:819–834.https://doi.org/10.1021/ct700324x
-
Protein traffic in bacteria: multiple routes from the ribosome to and across the membraneProgress in Nucleic Acid Research and Molecular Biology 66:107–157.https://doi.org/10.1016/S0079-6603(00)66028-2
-
Constant pressure molecular dynamics for molecular systemsMolecular Physics 50:1055–1076.https://doi.org/10.1080/00268978300102851
-
Polymorphic transitions in single crystals: A new molecular dynamics methodJournal of Applied Physics 52:7182–7190.https://doi.org/10.1063/1.328693
-
A large conformational change couples the ATP binding site of SecA to the SecY protein channelJournal of Molecular Biology 374:965–976.https://doi.org/10.1016/j.jmb.2007.09.086
-
Comparative protein modelling by satisfaction of spatial restraintsJournal of Molecular Biology 234:779–815.https://doi.org/10.1006/jmbi.1993.1626
-
HOLE: a program for the analysis of the pore dimensions of ion channel structural modelsJournal of Molecular Graphics 14:354–360.https://doi.org/10.1016/S0263-7855(97)00009-X
-
ACPYPE - antechamber python parser interfaceBMC Research Notes 5:367.https://doi.org/10.1186/1756-0500-5-367
-
From Coarse Grained to Atomistic: A Serial Multiscale Approach to Membrane Protein SimulationsJournal of Chemical Theory and Computation 7:1157–1166.https://doi.org/10.1021/ct100569y
-
Site-directed mutagenesis of colicin E1 provides specific attachment sites for spin labels whose spectra are sensitive to local conformationProteins: Structure, Function, and Genetics 6:294–305.https://doi.org/10.1002/prot.340060312
-
SecA supports a constant rate of preprotein translocationJournal of Biological Chemistry 281:15709–15713.https://doi.org/10.1074/jbc.M600205200
-
United atom lipid parameters for combination with the optimized potentials for liquid simulations all-atom force fieldJournal of Chemical Theory and Computation 5:1803–1813.https://doi.org/10.1021/ct900086b
-
Hydrogen exchange mass spectrometry for the analysis of protein dynamicsMass Spectrometry Reviews 25:158–170.https://doi.org/10.1002/mas.20064
-
Stabilization of proteins in confined spacesBiochemistry 40:11289–11293.https://doi.org/10.1021/bi0155504
Article and author information
Author details
Funding
Biotechnology and Biological Sciences Research Council (BB/M003604/1)
- Robin A Corey
- Ian Collinson
Wellcome (109854/Z/15/Z)
- Zainab Ahdash
- Argyris Politis
Wellcome (099149/Z/12/Z)
- Anokhi Shah
- Janet E Lovett
Engineering and Physical Sciences Research Council (ep/m508214/1)
- Anokhi Shah
Biotechnology and Biological Sciences Research Council (BB/I008675/1)
- William J Allen
- Ian Collinson
Wellcome (104632)
- William J Allen
- Ian Collinson
European Regional Development Fund (CZ.02.1.01/0.0/0.0/15_003/0000441)
- Tomas Fessl
Royal Society (University Research Fellowship)
- Janet E Lovett
Biotechnology and Biological Sciences Research Council (BB/N015126/1)
- Ian Collinson
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was funded by the BBSRC: BB/M003604/1, BB/I008675/1 and BB/N015126/1 to RAC, WJA and IC, and Wellcome: 104632 to IC and WJA and 109854/Z/15/Z to AP and ZA. JEL thanks the Royal Society for a University Research Fellowship and Wellcome for a Multi-User Equipment Grant (099149/Z/12/Z). TF is supported from European Regional Development Fund-Project ‘Mechanisms and dynamics of macromolecular complexes: from single molecules to cells’ (No. CZ.02.1.01/0.0/0.0/15_003/0000441). EP is the recipient of an Imperial College London Institute of Chemical Biology EPSRC CDT studentship. AS is funded by the EPSRC (ep/m508214/1). This work was carried out using the computational facilities of the Advanced Computing Research Centre, University of Bristol (http://www.bris.ac.uk/acrc/). Additional simulations were carried out using computer time on EPCC’s Cirrus HPC Service (https://www.epcc.ed.ac.uk/cirrus) and on the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk), provided by HECBioSim, the UK High End Computing Consortium for Biomolecular Simulation (http://www.hecbiosim.ac.uk/), supported by the EPSRC.
Copyright
© 2019, Corey 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
-
- 2,908
- views
-
- 471
- downloads
-
- 32
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
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)
Further reading
-
- Biochemistry and Chemical Biology
The biosynthesis of cyclic 3′,5′-adenosine monophosphate (cAMP) by mammalian membrane-bound adenylyl cyclases (mACs) is predominantly regulated by G-protein-coupled receptors (GPCRs). Up to now the two hexahelical transmembrane domains of mACs were considered to fix the enzyme to membranes. Here, we show that the transmembrane domains serve in addition as signal receptors and transmitters of lipid signals that control Gsα-stimulated mAC activities. We identify aliphatic fatty acids and anandamide as receptor ligands of mAC isoforms 1–7 and 9. The ligands enhance (mAC isoforms 2, 3, 7, and 9) or attenuate (isoforms 1, 4, 5, and 6) Gsα-stimulated mAC activities in vitro and in vivo. Substitution of the stimulatory membrane receptor of mAC3 by the inhibitory receptor of mAC5 results in a ligand inhibited mAC5–mAC3 chimera. Thus, we discovered a new class of membrane receptors in which two signaling modalities are at a crossing, direct tonic lipid and indirect phasic GPCR–Gsα signaling regulating the biosynthesis of cAMP.
-
- Biochemistry and Chemical Biology
The conformational ensemble and function of intrinsically disordered proteins (IDPs) are sensitive to their solution environment. The inherent malleability of disordered proteins, combined with the exposure of their residues, accounts for this sensitivity. One context in which IDPs play important roles that are concomitant with massive changes to the intracellular environment is during desiccation (extreme drying). The ability of organisms to survive desiccation has long been linked to the accumulation of high levels of cosolutes such as trehalose or sucrose as well as the enrichment of IDPs, such as late embryogenesis abundant (LEA) proteins or cytoplasmic abundant heat-soluble (CAHS) proteins. Despite knowing that IDPs play important roles and are co-enriched alongside endogenous, species-specific cosolutes during desiccation, little is known mechanistically about how IDP-cosolute interactions influence desiccation tolerance. Here, we test the notion that the protective function of desiccation-related IDPs is enhanced through conformational changes induced by endogenous cosolutes. We find that desiccation-related IDPs derived from four different organisms spanning two LEA protein families and the CAHS protein family synergize best with endogenous cosolutes during drying to promote desiccation protection. Yet the structural parameters of protective IDPs do not correlate with synergy for either CAHS or LEA proteins. We further demonstrate that for CAHS, but not LEA proteins, synergy is related to self-assembly and the formation of a gel. Our results suggest that functional synergy between IDPs and endogenous cosolutes is a convergent desiccation protection strategy seen among different IDP families and organisms, yet the mechanisms underlying this synergy differ between IDP families.