Active site geometry stabilization of a presenilin homolog by the lipid bilayer promotes intramembrane proteolysis

  1. Lukas P Feilen
  2. Shu-Yu Chen
  3. Akio Fukumori
  4. Regina Feederle
  5. Martin Zacharias
  6. Harald Steiner  Is a corresponding author
  1. German Center for Neurodegenerative Diseases, Germany
  2. Center of Functional Protein Assemblies and Physics Department T38, Technical University of Munich, Germany
  3. Department of Pharmacotherapeutics II, Osaka Medical and Pharmaceutical University, Japan
  4. Institute for Diabetes and Obesity, Monoclonal Antibody Core Facility, Helmholtz Munich, German Research Center for Environmental Health, Germany
  5. Biomedical Center (BMC), Division of Metabolic Biochemistry, Faculty of Medicine, Germany

Abstract

Cleavage of membrane proteins in the lipid bilayer by intramembrane proteases is crucial for health and disease. Although different lipid environments can potently modulate their activity, how this is linked to their structural dynamics is unclear. Here, we show that the carboxy-peptidase-like activity of the archaeal intramembrane protease PSH, a homolog of the Alzheimer’s disease-associated presenilin/γ-secretase is impaired in micelles and promoted in a lipid bilayer. Comparative molecular dynamics simulations revealed that important elements for substrate binding such as transmembrane domain 6a of PSH are more labile in micelles and stabilized in the lipid bilayer. Moreover, consistent with an enhanced interaction of PSH with a transition-state analog inhibitor, the bilayer promoted the formation of the enzyme’s catalytic active site geometry. Our data indicate that the lipid environment of an intramembrane protease plays a critical role in structural stabilization and active site arrangement of the enzyme-substrate complex thereby promoting intramembrane proteolysis.

Editor's evaluation

This work provides a strong contribution to our understanding of intramembrane proteolysis and in particular the subtle structural but significant influence of the lipid bilayer on proteolytic activity and coordination of the active site geometry.

https://doi.org/10.7554/eLife.76090.sa0

eLife digest

Cutting proteins into pieces is a crucial process in the cell, allowing several important processes to take place, including cell differentiation (which allows cells to develop into specific types), cell death, protein quality control, or even where in the cell a protein will end up. However, the specialized proteins that carry out this task, known as proteases, can also be involved in the development of disease. For example, in the brain, a protease called γ-secretase cuts up the amyloid-β protein precursor, producing toxic forms of amyloid-β peptides that are widely believed to cause Alzheimer’s disease.

Proteases like γ-secretase carry out their role in the membrane, the layer of fats (also known as lipids) that forms the outer boundary of the cell. The environment in this area of the cell can influence the activity of proteases, but it is poorly understood how this happens.

One way to address this question would be to compare the activity of γ-secretase in the lipid environment of the membrane to its activity when it is entirely surrounded by different molecules, such as detergent molecules. Unfortunately, γ-secretase is not active when it is removed from its lipid environment by a detergent, making it difficult to perform this comparison. To overcome this issue, Feilen et al. chose to study PSH, a protease similar to γ-secretase that produces the same amyloid-β peptides but remains active in detergent.

When Feilen et al. mixed PSH with lipid molecules like those found in the membrane and amyloid-β precursor protein, PSH produced amyloid-β peptides including those that are thought to cause Alzheimer’s. However, when a detergent was substituted for the lipid molecules this led to longer amyloid-β peptides than usual, indicating that PSH was not able to cut proteins as effectively. The change in environment appeared to reduce PSH’s ability to progressively trim small segments from the peptides.

Computer modelling of the protease’s structure in lipids versus detergent supported the experimental findings: the model predicted that the areas of PSH important for recognizing and cutting other proteins would be more stable in the membrane compared to the detergent.

These results indicate that the cell membrane plays a vital role in the stability of the active regions of proteases that are cleaving in this environment. In the future, this could help to better understand how changes to the lipid molecules in the membrane may contribute to the activity of γ-secretase and its role in Alzheimer’s disease.

Introduction

Intramembrane proteolysis is a crucial cellular mechanism underlying many fundamental physiological processes (Erez et al., 2009; Beard et al., 2019). It is also involved in pathological conditions, most prominently in Alzheimer´s disease (AD). Here, intramembrane cleavage within the transmembrane domain (TMD) of the amyloid precursor protein (APP) derived C99 substrate by γ-secretase results in the release of a variety of amyloid β-peptide (Aβ) species (Steiner et al., 2018). The longer forms, Aβ42 and Aβ43, are toxic to neurons and believed to trigger the onset of AD (Selkoe and Hardy, 2016). γ-Secretase is a membrane-embedded protein complex consisting of four components (Yang et al., 2017). The catalytic subunit presenilin is an aspartyl intramembrane protease (Wolfe et al., 1999; Li et al., 2000; Steiner et al., 2000; Steiner et al., 1999; Kimberly et al., 2000) present in the mammalian γ-secretase complexes as either presenilin-1 (PS1) or presenilin-2 variant (Yu et al., 1998; Saura et al., 1999). Mutations in PS1 are the major cause of familial AD (FAD) and cause an imbalance in the production of Aβ species that leads to relative increases of the longer forms over the normally major form Aβ40 (Steiner et al., 2018). Presenilins are evolutionary highly conserved proteins and related to the signal peptide peptidase (SPP) family of intramembrane proteases (Ponting et al., 2002; Weihofen et al., 2002). Ancestral precursors of presenilin and SPP exist in several archaea (Torres-Arancivia et al., 2010) and share key signature motifs including the protease family-defining GxGD active site motif (Steiner et al., 2000) with presenilin and SPP. The archaeal homolog from Methanoculleus marisnigri JR1 termed presenilin/SPP homolog (PSH) is capable of cleaving C99 and several other substrates (Torres-Arancivia et al., 2010; Dang et al., 2015; Naing et al., 2015; Naing et al., 2018). Similar to presenilin in the γ-secretase complex, PSH appears to cleave C99 in a sequential manner starting by initial ε-site cleavages between L49 and V50 (ε49) or T48 and L49 (ε48) followed by the release of various Aβ species from stepwise carboxy-terminal trimming cleavages (Dang et al., 2015; Takami et al., 2009). However, in contrast to presenilin, which requires complex formation with the other γ-secretase complex components for activity (Takasugi et al., 2003; Edbauer et al., 2003; Kimberly et al., 2003), PSH is active without accessory components. The crystal structure of PSH revealed first important insights into aspartyl intramembrane proteases showing that the two catalytic aspartate residues of the active site in TMD6 and TMD7 directly face each other and locate in a water-accessible cavity (Li et al., 2012). Subsequent cryo-electron microscopy (cryo-EM) structural analysis of γ-secretase showed that presenilin adopts a structure in the complex very similar to that of PSH (Sun et al., 2015) with nearly identical positions of the catalytic residues (Bai et al., 2015b). Further cryo-EM studies showed that binding of APP and Notch substrates causes major conformational changes in both enzyme and substrate (Zhou et al., 2019; Yang et al., 2019). These led to an enzyme-substrate complex (E-S) with an extended TMD6 by formation of a new and stable TMD6a helix as well as a hybrid β-sheet between enzyme and substrate that causes unfolding of the ε-cleavage site region in the substrate (Zhou et al., 2019; Yang et al., 2019). Interestingly, formation of the TMD6a helix was also observed by cryo-EM upon inhibitor binding thus partially mimicking the substrate-bound state (Bai et al., 2015a; Yang et al., 2021).

The very similar structural folds of presenilin and PSH and the ability to cleave C99 in the TMD at the same sites as γ-secretase (Torres-Arancivia et al., 2010; Dang et al., 2015) make PSH an attractive model for the intrinsic protease activity of presenilin. To gain basic insights into the enzymatic workings of presenilin proteases, we thus set out to characterize the influence of two fundamentally different hydrophobic environments on the activity of PSH and asked if cleavage of C99 by the solubilized enzyme in detergent micelles would differ from a lipid-reconstituted state and if so, whether such differences could be correlated with the structural dynamics of this prototype presenilin protease or its E-S. Although the influence of lipids on the activity of presenilin and other intramembrane proteases is well documented (Paschkowsky et al., 2018), there are so far no studies in which biochemically determined activities of these proteases were linked with structural information that could explain how lipids, in particular a membrane bilayer environment, affect intramembrane protease structural dynamics and enzyme function. Since presenilins are not active in detergent micelles without lipids (Zhou et al., 2010), this critical question can however not be addressed for γ-secretase directly and requires a suitable model protease such as PSH. We found that detergent-solubilized PSH has a reduced carboxy-terminal trimming activity, that is processivity, compared to γ-secretase giving rise to an increased production of very long Aβ species such as Aβ46. Strikingly, the reconstitution of PSH into a lipid bilayer strongly promoted the protease processivity to shorter Aβ species such as Aβ38 highlighting the important role of the lipid membrane environment for intramembrane proteolysis. Furthermore, it enhanced the binding of a transition-state analog (TSA) γ-secretase inhibitor (GSI) affinity probe suggesting a more stable active site conformation in the lipid bilayer. These biochemical studies were accompanied by comparative modeling and molecular dynamics (MD) simulations to study the effect of detergent micelle and membrane lipid environment on substrate-bound PSH. In good agreement with the experimental data, the computational data suggest that the stabilization of TMD6a and the active site can explain the increased processivity and inhibitor binding in the membrane bilayer. Mutational analysis confirmed the assumed critical functional role of β-sheet and TMD6a corroborating the computational analysis of substrate-bound PSH. Collectively, these data provide insights into how structural adaptations occurring in response to changes in the hydrophobic environment from a micellar membrane mimetic to a real lipid bilayer translate into activity changes of an intramembrane protease. Moreover, with general implications for intramembrane proteolysis, they show how a lipid bilayer allows the formation of a stabilized active site geometry poised for substrate cleavage.

Results

PSH cleaves APP C99 to longer Aβ species

To get insights into the intrinsic protease activity of presenilin, we set out to further characterize the intramembrane cleavage of C99 by PSH (Figure 1A). Consistent with previous findings (Dang et al., 2015), n-dodecyl β-D-maltoside (DDM)-solubilized, His-affinity-purified PSH could cleave the C99-based APP C100-His6 substrate (Edbauer et al., 2003) as demonstrated by the generation of the APP intracellular domain (AICD) and Aβ cleavage products (Figure 1B). Cleavage was inhibited by the TSA GSI L-685,458 (Shearman et al., 2000; Figure 1B) although much higher, micromolar concentrations were needed for efficient inhibition compared to those known for γ-secretase (Li et al., 2000; Shearman et al., 2000). Analysis of the Aβ profile using Tris-Bicine urea SDS-PAGE (Figure 1C) and MALDI-TOF mass spectrometry (Figure 1D, Figure 1—figure supplement 1) showed that Aβ40 and Aβ42 were the major Aβ species produced with a preference of Aβ42 over Aβ40. Interestingly, besides the increased generation of Aβ42 even longer Aβ species such as Aβ46 were relatively abundant. This suggests that PSH cleaves C99 at the same sites as γ-secretase but with reduced processivity.

Figure 1 with 1 supplement see all
Cleavage of APP C99 by PSH.

(A) Schematic illustration of APP C99 cleavage by PSH. PSH cleaves C99 and releases an AICD fragment and Aβ peptides. The epitope of the PSH specific antibody 6F4 in the loop between TMD6 and TMD7 is indicated. (B) Analysis of PSH activity in DDM micelles after incubation with C100-His6 substrate overnight at 37 °C by immunoblotting for AICD (Y188) and Aβ (2D8). Specificity of substrate cleavage by PSH in the assay was controlled by sample incubation at 4 °C or 37 °C in presence of the GSI L-685,458 (20 µM). Immunoblotting of PSH (6F4) was performed to control for PSH levels. (C) Aliquot of samples from (B) separated by Tris-Bicine urea SDS-PAGE for identifying Aβ species produced by PSH in DDM micelles and analysis by immunoblotting (2D8). In (B) and (C), representative immunoblots from three to six independent biological replicates (i.e. independent protease preparations) are shown. (D) Representative MALDI-TOF MS spectrum of Aβ profile generated by PSH in DDM micelles from four independent biological replicates. The intensity of the highest peak was set to 100%. A GSI control is shown in Figure 1—figure supplement 1 and observed masses for identified Aβ species are shown in Figure 1—source data 1.

Figure 1—source data 1

Immunoblot images (raw and annotated) of cleavage assay (Source data for Figure 1B, C).

https://cdn.elifesciences.org/articles/76090/elife-76090-fig1-data1-v2.zip
Figure 1—source data 2

Calculated and observed masses for Aβ species in MALDI-TOF mass spectrometry (Source data for Figure 1D).

https://cdn.elifesciences.org/articles/76090/elife-76090-fig1-data2-v2.xlsx

Lipid membrane enhances the processivity of PSH

We next investigated how a membrane environment of PSH influences the cleavage and processivity of the protease. Previous studies have shown that γ-secretase activity is dependent on the membrane environment as modulations of the lipid composition and/or bulk membrane properties in cell-free assays affected total activity as well as the ratios of the Aβ species generated (Osenkowski et al., 2008; Osawa et al., 2008; Holmes et al., 2012; Winkler et al., 2012). Furthermore, it was also shown that varying the pH in cell-free assays can modulate the total activity as well as the processivity of γ-secretase (Quintero-Monzon et al., 2011). To investigate whether and how PSH cleavage of C99 would respond to a change from the micellar environment in DDM to a lipid bilayer environment, PSH was reconstituted in defined small unilamellar vesicles (SUVs) composed of palmitoyl-oleoyl PC (POPC), the most abundant phospholipid of biological membranes (Figure 2—figure supplement 1A, B). We then performed PSH in vitro assays with DDM-solubilized or POPC-reconstituted PSH in a pH range from 5.5 to 9.0. As shown in Figure 2A, the total activity was highest in the mild acidic to mild alkaline pH range and sharply dropped at pH values above 8.0. The pH optima for both conditions were very similar and lying around pH 7.0. However, compared to the DDM-solubilized enzyme, the processivity of PSH was strongly promoted in the lipid bilayer environment of the POPC SUVs as seen by a strongly increased production of Aβ38 and Aβ40 and the strong reduction of longer Aβ species (Figure 2B). Interestingly, the processivity of the reconstituted PSH appeared to be more reduced at alkaline pH values as judged from the appearance of Aβ species longer than Aβ42 (Figure 2B) at pH 7.5 and higher. A direct comparison of the Aβ profiles at pH 7.0 confirmed the increased processivity of PSH in the POPC lipid bilayer (Figure 2C, Figure 2—figure supplement 2A). Finally, we investigated the initial ε-site cleavages of C99 by PSH in DDM micelles or POPC vesicles at this pH. Mass spectrometry analysis showed that C99 was cleaved in both conditions at the ε49 and ε48 cleavage sites resulting in the release of the two N-terminally distinct AICD50 (ε49) and AICD51 (ε48) species (Figure 2D, Figure 2—figure supplement 2B). Collectively, these data show that the lipid environment increases the processivity of PSH in cleaving C99.

Figure 2 with 2 supplements see all
Comparison of PSH cleavage activity and processivity in DDM micelles and POPC bilayer.

(A) Analysis of PSH activity in DDM micelles and POPC vesicles after incubation with C100-His6 substrate at 37 °C overnight by immunoblotting for AICD (Y188) and Aβ (2D8). Immunoblotting of PSH (6F4) was performed to control for PSH levels. (B) Separation of Aβ species produced by PSH in DDM micelles and POPC vesicles by Tris-Bicine urea SDS-PAGE and analysis by immunoblotting for Aβ (2D8). In (A) and (B), representative immunoblots from six independent biological replicates are shown. Confirmation of PSH reconstitution in POPC SUVs is shown in Figure 2—figure supplement 1. (C, D) MALDI-TOF MS analysis of Aβ (C) and AICD (D) species generated by PSH in DDM micelles and POPC vesicles at pH 7.0. Representative mass spectra from four independent biological replicates are shown. The intensity of the highest peak was set to 100%. GSI controls are shown in Figure 2—figure supplement 2 and observed masses for identified Aβ and AICD species are shown in Figure 2—source data 1.

Figure 2—source data 1

Immunoblot images (raw and annotated) of cleavage assays (Source data for Figure 2A, B).

https://cdn.elifesciences.org/articles/76090/elife-76090-fig2-data1-v2.zip
Figure 2—source data 2

Calculated and observed masses for Aβ, and AICD species in MALDI-TOF mass spectrometry (Source data for Figure 2C and D).

The differences between calculated and observed masses for AICD species arise from different adducts as specified in the respective tables. The GSI controls in Figure 2—figure supplement 2B show that the observed peaks are specific for AICD species derived from PSH cleavage.

https://cdn.elifesciences.org/articles/76090/elife-76090-fig2-data2-v2.xlsx

Enhanced processivity of PSH is independent of the APP substrate N-terminus

Since the activities of DDM-solubilized and reconstituted PSH clearly differed, particularly in the processivity, we next sought to understand the underlying basis for this behavior at the level of its structural dynamics. Since there is no structure of PSH in complex with C99 available and structural investigations on γ-secretase in complex with an APP substrate were so far only performed with C83, an N-terminally shorter alternative C-terminal APP fragment generated by α-secretase (Lichtenthaler et al., 2011), we first tested whether C83 is processed similarly to C99. We thus analysed its cleavage by PSH in DDM micelles and the POPC bilayer at pH 7.0. C83 was cleaved by PSH in both conditions resulting in the generation of an AICD and the Aβ-equivalent cleavage product p3 (Lichtenthaler et al., 2011; Figure 3A). As judged from the processivity-reflecting ratios of p3 species ending at position 40 and 42, in contrast to the DDM micelle environment, the POPC bilayer enhanced the processivity and caused an increase in the relative production of shorter p3 species (Figure 3B). The increase in processivity was similar to that observed under these conditions for the corresponding Aβ species produced from C99 (Figure 3C). These data show that the increased processivity in the POPC environment is independent of the N-terminus of the substrate and that C83 and C99 behave comparable in both environments, so that C83 should be suitable as C99 surrogate for structural modeling and molecular dynamics simulations of PSH in complex with a substrate using information from the C83-bound γ-secretase.

Cleavage of APP C83 by PSH.

(A) Analysis of PSH activity in DDM and POPC environment after incubation with C83-His6 and C100-His6 substrates at 37 °C and pH 7.0 by immunoblotting for AICD (penta-His) and p3 (Aβ (22-35)). Immunoblotting of PSH (6F4) was performed to control for PSH levels. (B, C) p3-40/p3-42 ratio (B) and Aβ40/Aβ42 ratio (C) from PSH activity assays in DDM (red) and POPC (blue) environment analyzed by ECL-IA. Quantitative data are represented as mean ± standard deviation (SD) (n=3 biological replicates). Source data are shown in Figure 3—source data 1.

Figure 3—source data 1

Immunoblot images (raw and annotated) of cleavage assays (Source data for Figure 3A).

https://cdn.elifesciences.org/articles/76090/elife-76090-fig3-data1-v2.zip
Figure 3—source data 2

Raw values of p3 and Aβ concentrations measured in the ECL-IA and calculated p3-40/p3-42 and Aβ40/Aβ42 ratios (Source data for Figure 3B, C).

https://cdn.elifesciences.org/articles/76090/elife-76090-fig3-data2-v2.xlsx

Structural modeling shows key features of substrate-bound γ-secretase in PSH

Since no experimental structure of the substrate-bound (holo) PSH is available, we generated 3 different starting structures, models 1, 2 and 3, for the PSH holo form using template-based modeling by assuming that PSH binds its substrate in a similar way as PS1. Since the PSH crystal structure (Li et al., 2012) (PDB 4HYG) misses several residues and loop segments in its substrate-free (apo) form, it is necessary to include the cryo-EM structure of holo PS1 bound to the C83 γ-secretase substrate (Zhou et al., 2019) (PDB 6IYC) as an additional template. In model 1, we used the entire holo PS1 including the C83 substrate as template whereas in models 2 and 3 most of the apo PSH structure was included and only different parts of missing loops segments were modeled based on the 6IYC template (including the C83 substrate, see Materials and methods for details). The overall structures of our holo PSH models are close to the experimentally determined apo PSH crystal structure (Li et al., 2012) but include several residues and loop segments that are missing in the apo crystal structure (Table 1, Figure 4—figure supplement 1A, B, Figure 4—figure supplement 2). In particular, as exemplified for model 2 (Figure 4A), it includes the TMD6a helix (residues H171-E177) C-terminal of TMD6 and the hybrid β-sheet formed between the β2-strand (A213-G217) and the β3-strand (V50-K54) of the substrate, which is stabilized by backbone interactions with residue Q272 preceding TMD9 (Figure 4B). These two structural elements were also found in the γ-secretase complexes with bound C83 (Zhou et al., 2019) or Notch1 (Yang et al., 2019) (PDB 6IDF) and were not present in the substrate-free γ-secretase complex. In addition, a salt bridge between R70 and E181 in the homology-modeled PSH (Figure 4C) replaced the hydrogen bond partners Y159 and R278 of the substrate-bound γ-secretase structure (Figure 4D). Thus, key features for substrate interaction (TMD6a and hybrid β-sheet) known from the γ-secretase–substrate complexes are analogously found in our models of the PSH–C83 complex.

Figure 4 with 2 supplements see all
PSH homology model.

(A) Alignment of the modeled holo form (model 2) of PSH (blue) with APP C83 substrate (orange) and the crystal structure of PSH (PDB 4HYG) in the apo form (green) in side view (left panel) and top view (right panel). An overlay of all three models and their RMSD values are provided in Figure 4—figure supplement 1A and B. (B) Schematic representation of the hydrogen bonds formed between β3 of the substrate (orange), β2, and Q272, respectively, of PSH (blue). (C) Interaction of TMD3 and TMD6a in the C83-bound PSH model through residues R70 and E181. (D) Interaction of TMD3 and TMD6a through residues Y159 and R278 in the C83-bound γ-secretase cryo-EM structure (PDB 6IYC). (E) Analysis of WT and mutant PSH activity in DDM and POPC environment after incubation with C100-His6 substrate at 37 °C overnight by immunoblotting for AICD (Y188) and Aβ (2D8). Immunoblotting of PSH (6F4) was performed to control for PSH levels. The asterisks mark two substrate degradation bands, which are independent of PSH cleavage.

Figure 4—source data 1

Immunoblot images (raw and annotated) of cleavage assays (Source data for Figure 4E).

https://cdn.elifesciences.org/articles/76090/elife-76090-fig4-data1-v2.zip
Table 1
Templates used for model building of PSH in complex with C83.

Residues of PSH used for model building are indicated.

Template
6IYC(PS1, C83)4HYG(Chain B)Model 1
Model 1*
Model 2L7-D162, D220-L292complete
Model 3L7-A176, E210-A293complete
  1. *

    PSH residues L7-R193 and E210-A293 were modeled based on the template.

Verification of the β2-strand in the substrate-bound PSH models

Comparative modeling is sensitive to the choice of the templates and how the sequences are aligned together. A template with low quality or a sequence alignment with high uncertainties can lead to an unrealistic protein structure. Therefore, it is necessary to verify the presence of specific structural features in our substrate-bound PSH models that were observed in the substrate-bound PS1 structure. In case of γ-secretase, the deletion of the β2-strand (R377-L381) impaired the activity of the enzyme towards C83 and Notch1 (Zhou et al., 2019; Yang et al., 2019). To investigate whether this structural element is of similar functional importance in PSH, we mutated amino acid residues A213, F214, and V215 within the putative PSH β2-strand to prolines. Similar as done for PS1 (Zhou et al., 2019; Yang et al., 2019), we also deleted amino acid residues A213 to M216 and A213 to G217. When assessed for their enzymatic activity, the three mutants as well as the two deletion constructs showed clearly, and mainly strongly, decreased activities compared to wild type (WT) PSH in both DDM micelle or POPC bilayer conditions (Figure 4E). These results suggest that as the residues 213 to 216/217 are important for the activity of the protease they might indeed form the β2-strand observed in our structural models.

Comparative molecular dynamics simulations of PSH in micelle and membrane environment reveals reduced PSH flexibility in the lipid bilayer

To get insight into the molecular details on how micelle and membrane environments might influence PSH conformational dynamics, the constructed C83-bound PSH models were embedded in DDM micelle (150 DDM molecules) or POPC bilayer (302 POPC molecules) environments and in each case three simulations with different distribution of starting velocities were performed (each simulation length: 0.6 µs). In total, six systems were constructed and 18 trajectories were generated in silico.

Snapshots of the PSH holo form in both environments are shown for model 2 in Figure 5A. All simulated systems stayed overall close to the starting structures with root-mean-square deviations (RMSD) relative to the start structure of about or less than 4.0 Å and an overall lower RMSD for all model 2 simulations (Figure 5—figure supplement 1A). Larger deviations were observed for all model 1 simulations and for one trajectory for model 3. It is likely that model 1 is less realistic than models 2 and 3 because it is entirely based on the PS1 template structure and is missing structural information from PSH.

Figure 5 with 3 supplements see all
MD simulations of holo PSH forms in DDM micelle or POPC bilayer.

(A) PSH with bound C83 substrate (model 2) embedded in a DDM micelle environment (upper panel) and a POPC bilayer (lower panel). (B) The average number of H-bonds formed between the β3-strand of C83 and the β2-strand of PSH. Each data point stands for the average value throughout one trajectory and the error bars represent the SD of the mean of three data points. (C) The backbone RMSF of C83 of different models in DDM (red) and POPC (blue) environments averaged over three trajectories. The shaded areas represent the SD of the mean. (D) Water accessibility along the substrate TMD residues extracted from the simulations of the holo PSH in DDM (red) and POPC (blue) environments (water accessibility for a residue is obtained as the mean number of water molecules within 5 Å of any atom of the residue). The red arrows indicate the position of the two ε-cleavage sites. The error bars represent the SD of the mean (n=3 trajectories). (E) The backbone RMSF of PSH of different models in DDM (red) and POPC (blue), environments averaged over three trajectories (note that residues 194–209 are not included in our PSH models). The gray boxes highlight TMD6a and the shaded areas represent the SD of the mean. Enlarged views on backbone RMSF of residues K170 to P185 (including TMD6a) are shown in Figure 5—figure supplement 1B.

Figure 5—source data 1

Raw values of simulation data analysis (Source data for Figure 5B–E).

https://cdn.elifesciences.org/articles/76090/elife-76090-fig5-data1-v2.zip

We first analyzed the substrate mobility and interaction with PSH during the simulations in both environments. Stable substrate binding near the active site region involves the interaction of the β2-strand with the β3-strand at the C-terminus of the C83 substrate, which is required for substrate cleavage. The overall β-sheet interaction in terms of hydrogen bonds (H-bonds) was found to be similar in both the micelles and membrane environments except for model 3 where the H-bonds are more frequently formed in the bilayer (Figure 5B). We next investigated the mobility of the individual residues of C83 and PSH by calculating the root-mean-square fluctuation (RMSF) as well as the water accessibility of residues in the C83 substrate TMD by counting the average number of water molecules within 5 Å of the residue of interest. Both RMSF and water accessibility for each amino acid were not much different in both environments (Figure 5C and D). Notably, the substrate remained in a dry region from G37 to V46 and abruptly gained an increase in water accessibility around T48 and L49, which correspond to the initial ε-cleavage sites of C83 (and C99).

To quantitatively evaluate the flexible regions of PSH, RMSF profiles of each residue along the PSH sequence were calculated (Figure 5E). In all three models, larger fluctuations of hydrophilic loop 1 (HL1, between TMD1 and 2) were observed in DDM compared to POPC (Figure 5E). Similarly, the atomic fluctuations of TMD6a in models 1 and 2 were reduced in POPC versus DDM but no differences in the two environments were observed for TMD6a in model 3. In both environments, the atomic fluctuations of residues C-terminal of TMD6a stayed reduced in model 1 while they increased in model 3. In contrast to models 1 and 3, the atomic fluctuations of these residues stayed reduced in POPC and increased in DDM in model 2 (Figure 5E, Figure 5—figure supplement 1B). Furthermore, secondary structure analysis showed that TMD6a is mostly unfolded in model 1 and stable in model 3 (Figure 5—figure supplement 2A, C). Strikingly, when model 2 is placed in the micelle environment, TMD6a of substrate-bound PSH underwent a conformational transition between an α-helix and a loop structure (Figure 5—figure supplement 2B), while this transition was not observed for TMD6a in the POPC bilayer showing that the membrane environment stabilizes TMD6a (Figure 5—figure supplement 2B).

Despite the high structural similarity between the starting structures of model 2 and model 3 (RMSD = 0.163, Figure 4—figure supplement 1B), a difference in TMD6a positioning was observed during the simulations of these two models in the POPC environment. While TMD6a is located closer to C83 in model 2, it is located further away from C83 in model 3 (Figure 5—figure supplement 3A). In addition, when we directly compared the RMSF profiles of models 2 and 3 (Figure 5E) with each other, we observed that the residues immediately C-terminal of TMD6a are slightly more mobile in model 3 (Figure 5—figure supplement 3B, C). Because these residues are spatially close to HL4 between TMD4 and TMD5, HL4 also becomes more flexible in model 3 (Figure 5—figure supplement 3B). In fact, the conformational discrepancy in these regions arise from the model building of the two models which were built differently for residues A163 to R193 and E210 to G219 (Table 1, see Materials and methods for details). With a closer contact between TMD6a and C83 as well as an overall lower RMSD for model 2, we justified that model 2 may describe the dynamics of the PSH-C83 complex best so that model 2 is therefore used for our following analysis.

DDM insertion leads to an unwinding of PSH TMD6a

The biochemical cleavage data as well as the fluctuations in the MD simulations indicate a potential weakening of the E-S interaction in the micelle environment. This might explain the remarkable shift in processivity of PSH in the presence of a membrane lipid bilayer. In our PSH models, TMD6a creates a hydrophobic patch (formed by residues M172, I173, L175, and A176) that contacts the C83 substrate in the ε-cleavage site region (e.g. V50 and L52, Figure 6A). Similar interactions are found in the experimentally resolved C83-bound and Notch1-bound γ-secretase structures (Zhou et al., 2019; Yang et al., 2019), as well in the GSI-bound γ-secretase structures (Yang et al., 2021; Figure 6—figure supplement 1A-D). To gain a more detailed mechanistic view on how DDM and POPC molecules modulate the E-S stability, we calculated the structural properties of these molecules. Furthermore, we investigated how these molecules interact with PSH and C83. In the bilayer environment, the simulations indicate that POPC molecules are well ordered as indicated by a high lipid order parameter SCH (Figure 6—figure supplement 2A, B) and the computed area per lipid of ~68 Å2 is close to the experimentally determined value of 64.3 Å2 (Kučerka et al., 2011; Figure 6—figure supplement 2C). In contrast, in the micelle environment, DDM molecules are more mobile and can change their orientation more freely, as indicated by the lower lipid order parameter SCH for DDM compared to POPC (Figure 6—figure supplement 2A, B). In a larger micelle environment, with 50% more DDM molecules (225 DDM molecules, Figure 6—figure supplement 3A) the lipid order parameter of DDM improved (Figure 6—figure supplement 3B). Nevertheless, the number of H-bonds in the β-sheet (Figure 6—figure supplement 3C), the RSMF profile (Figure 6—figure supplement 3D) and the water accessibility in the C83 substrate (Figure 6—figure supplement 3E) did not differ from the respective values in the smaller DDM micelle (150 DDM molecules). The atomic fluctuations also did not differ largely between the different micelle sizes except for strongly reduced fluctuations of HL1 in the larger micelle (Figure 6—figure supplement 3F) originating from the more ordered DDM molecules around HL1 (Figure 6—figure supplement 3H). Furthermore, TMD6a and residues C-terminal of TMD6a were found to be also flexible in the larger DDM micelle as observed for the smaller micelle (Figure 6—figure supplement 3F, G), indicating that the size of the DDM micelle does not influence the observed differences in flexibility of TMD6a between DDM and POPC environment.

Figure 6 with 3 supplements see all
Destabilization of PSH TMD6a in a DDM micelle.

(A) Hydrophobic interactions between PSH regions and the C83 substrate. The right panel shows an enlarged view of the interaction of a hydrophobic patch (gray surface) of TMD6a with V50 and L52 of the substrate in the POPC bilayer. (B, C) Snapshots of the DDM insertion between TMD2 and TMD6 in the first (B) and second (C) run of the simulations in DDM environment at 400 ns. The right panels show an enlarged view of the unspecific hydrogen bonding interactions between the DDM molecule and the TMD6a amino acid backbones. (D) Snapshot of the DDM insertion between TMD3 and TMD4 of the simulation in DDM environment at 430 ns. The right panel shows an enlarged view of the unspecific hydrogen bonding interactions between the DDM molecule and the amino acid backbones of the residues immediately C-terminal of TMD6a.

Although most of the DDM molecules are well aligned to the membrane normal at the protein periphery, some can transiently flip to a direction perpendicular to the membrane normal near to the gaps between TMDs (such as TMD2-TMD6 and TMD3-TMD4) (Figure 6B–D, Videos 1 and 2). When inserted into the intramolecular gaps between TMD2 and TMD6, the DDM molecule perturbs intramolecular interactions by forming unspecific hydrogen bonds with the adjacent amino acid backbones, and thus destabilizes TMD6a (Figure 6B–C). In addition, DDM inserts between TMD3 and TMD4 of PSH in the micelle environment and interacts with the loop C-terminal of TMD6a, corresponding to the higher RMSF observed for these residues (Figures 5E and 6D). In contrast, the well-ordered POPC molecules do not enter into the gaps between TMDs (Video 3) and do, therefore, not disturb intramolecular interactions. Collectively, these data suggest that a membrane lipid environment promotes the formation of a stabilized E-S of PSH with the APP C83 substrate by the stabilization of TMD6a, an important structural element involved in substrate stabilization.

Video 1
Detergent-enzyme interaction in the DDM environment.

600 ns trajectories of PSH (red) in complex with C83 (orange) and nearby DDM molecules. A DDM molecule enters into the gap between TMD2 and TMD6. TMD6a switches between a helical and a loop conformation with the interference of the disordered DDM molecule.

Video 2
Detergent-enzyme interaction in the DDM environment.

600 ns trajectories of PSH (red) in complex with C83 (orange) and nearby DDM molecules. A DDM molecule enters into the gap between TMD3 and TMD4. TMD6a switches between a helical and a loop conformation with the interference of the disordered DDM molecule.

Video 3
Lipid-enzyme interaction in the POPC environment.

600 ns trajectory of PSH (blue) in complex with C83 (orange) and nearby POPC molecules. POPC molecules do not enter between TMD gaps and TMD6a remains a stable helix throughout the whole trajectory.

Lysine mutations in TMD6a lead to helix unwinding and reduced activity

Our computational results suggest that the TMD6a helix plays an important role for substrate binding of PSH. Furthermore, the residues of the hydrophobic patch in TMD6a (M172, I173, L175, and A176) of PSH correspond to a homologous hydrophobic patch in the TMD6a of PS1 (L271, V272, T274, and A275), which is also affected by FAD mutations (Steiner et al., 2018). Some of them display a strong loss of function such as L271V and T274R (Sun et al., 2017) supporting the idea that TMD6a has an important function in substrate cleavage. To investigate the functional role that TMD6a plays in substrate stabilization, we performed additional MD simulations of in silico generated lysine mutations of M172, I173, L175, and A176 in the TMD6a hydrophobic patch in order to disrupt its nonpolar character. Because TMD6a is already unstable in DDM, the simulations were performed in the POPC environment only. The RMSD plots indicate that in most of the simulations, mutated PSH remained in an overall stable structure with RMSDs ~3.0 Å (similar to WT) (Figure 7A). In addition, no significant difference was found in C83 RMSF (Figure 7B), residue-wise water distribution (Figure 7C) and the hydrogen-bonding pattern of the β-sheet C-terminal of the ε-cleavage site (Figure 7D).

Figure 7 with 1 supplement see all
Impact of PSH TMD6a mutations on PSH structural dynamics and activity.

(A) RMSD of the WT (blue) and the mutated systems M172K (violet), I173K (green), L175K (gray), and A176K (orange) in the POPC bilayer environment. The solid, dashed, and dotted lines represent three different simulations with random initial velocities. The black dashed line indicates an RMSD of 4 Å. (B) The backbone RMSF of C83 of WT PSH and different TMD6a lysine mutants in POPC environment averaged over three trajectories. The shaded areas represent the SD of the mean. (C) Water accessibility along the substrate TMD residues extracted from the simulations of the four lysine-mutant holo-form PSH systems in the POPC bilayer environments (water accessibility for a residue is obtained as the mean number of water molecules within 5 Å of any atom of the residue). The red arrows indicate the position of the two ε-cleavage sites. The error bars represent the SD of the mean (n=3 trajectories). (D) The average H-bond formed between the β3-strand at the C-terminus of C83 and the β2-strand of WT and lysine-mutated PSH. Each data point stands for the average value throughout one trajectory and the error bars represent the SD of the mean of three data points. (E) The backbone RMSF of WT (blue) and M172K (violet), I173K (green), L175K (gray), and A176K (orange) mutated PSH in POPC. The box highlights TMD6a and the shaded areas represent the SD of the mean (n=3 trajectories). Larger RMSF of the A176K mutant observed in residues 235–243 correspond to a folding-unfolding event in the mobile loop between TMD7 and TMD8 in the third trajectory (see Figure 7—source data 1 and Figure 7—source data 2). (F) Analysis of WT and lysine-mutant PSH activity in DDM and POPC environment after incubation with C100-His6 substrate at 37 °C overnight by immunoblotting for AICD (Y188) and Aβ (2D8). Immunoblotting of PSH (6F4) was performed to control for PSH levels. The asterisks marks substrate degradation bands, which are independent of PSH cleavage.

Figure 7—source data 1

Raw values of simulation data analysis (Source data for Figure 7A–E).

https://cdn.elifesciences.org/articles/76090/elife-76090-fig7-data1-v2.zip
Figure 7—source data 2

Immunoblot images (raw and annotated) of cleavage assays (Source data for Figure 7F).

https://cdn.elifesciences.org/articles/76090/elife-76090-fig7-data2-v2.zip

The PSH RMSF plots show that all mutations destabilize TMD6a to different degrees while having only smaller effects in other regions in comparison to WT (Figure 7E). It is worth noting that L175K introduces the largest TMD6a fluctuation, compared to the WT and the other three lysine mutants. Nevertheless, all four mutations distorted the helical structure of TMD6a in at least one of the simulations (Figure 7—figure supplement 1). Finally, to experimentally validate these structural predictions, we analyzed the cleavage of C99 by these mutant forms of PSH. All four mutants showed a strongly decreased, nearly abolished cleavage of C99 compared to WT PSH in both DDM micelle or POPC bilayer conditions (Figure 7F) suggesting that TMD6a and its hydrophobic patch is an important structural element of PSH.

Lipid membrane environment stabilizes the active site geometry of PSH

The proteolysis reaction requires a specific geometry of all elements that form the active site. A critical issue is the distance between the two catalytic aspartic acids D162 and D220. Free energy calculation along the distance between D257 and D385 of PS1 has suggested that mutations disturbing the active site geometry and alter the distance between D257 and D385 correlate with changes in enzyme activity (Chen and Zacharias, 2020). In PSH, the geometry is characterized by the Cγ-Cγ distances between the D162 and D220. The distances appeared more frequently around 6.8 Å in our model when placed in a bilayer environment (Figure 8A, Figure 8—figure supplement 1). These distances correspond to a potentially catalytically active site geometry that can also accommodate a water molecule between the catalytic aspartates and L49 essential for the proteolytic cleavage (Figure 8B). In contrast, larger distances are more frequently sampled in the micelle environment. When the distance is enlarged here, an increased number of water molecules can access the catalytic center and disturb the catalytic geometry (Figure 8B). Proteolysis-compatible Cγ-Cγ distances below 7.0 Å were more frequently sampled in the membrane environment (~76%) compared to simulations in the DDM environment (~63%) (Figure 8A, Figure 8—figure supplement 1). Detailed geometries at the catalytic site of the E-S for a smaller and a larger Cγ-Cγ distance are depicted in Figure 8—figure supplement 2A, B.

Figure 8 with 3 supplements see all
Stabilization of the PSH active site geometry in a POPC bilayer.

(A) Histograms of the Cγ-Cγ distances between the D162 and D220 of PSH measured in DDM micelle (red) and POPC bilayer (blue) environments. The dashed line indicates the distance of 7 Å. The measured distances over time are shown in Figure 8—figure supplement 1C. (B) Snapshot of the catalytic cavity in DDM (left panel) and POPC (right panel) environment. The Cγ-Cγ distance between the two catalytic aspartates D162 and D220 is larger in DDM micelles and more water molecules enter the catalytic cavity between D162 and the substrate. Detailed geometries of these two active site conformations are depicted in Figure 8—figure supplement 2. (C) Immunoblot analysis of TSA-inhibitor binding to PSH in DDM micelles or POPC vesicles. PSH was affinity-precipitated by Merck C (a biotinylated derivative of L-685,458; 20 µM). To control for background binding and binding specificity, the affinity precipitation was assessed in the absence of Merck C as well as in the presence of excess amounts of the parental compound L-685,458 (2 mM) as competitor. The input represents 2.5% of the total sample used for the affinity precipitation. A representative immunoblot from four independent biological replicates is shown. (D) Quantitation of PSH binding by Merck C. Specific binding was defined as difference of PSH signals in the absence or presence of L-685,458 after additional subtraction of unspecific background binding signals. Quantitative data are represented as mean ± SD (n=4 biological replicates). The source data are shown in Figure 8—source data 1. (E, F) Inhibition assay of PSH in DDM micelles and POPC vesicles with increasing concentrations of L-685,458 (E) or Merck C (F), respectively. PSH activity was analyzed by immunoblotting for AICD (Y188) and Aβ (2D8) following incubation with C100-His6 substrate at 37 °C overnight. Representative immunoblots from three independent biological replicates are shown. The asterisks mark two substrate degradation bands, which are independent of PSH cleavage. (G) Inhibition assay of PSH reconstituted in POPC vesicles in the presence of 20 µM TSA and non-TSA γ-secretase inhibitors. PSH activity was analyzed by immunoblotting for AICD (Y188) and Aβ (2D8) following incubation with C100-His6 substrate at 37 °C overnight. Representative immunoblots from three independent biological replicates are shown.

Figure 8—source data 1

Raw values of simulation data analysis (Source data for Figure 8A and Figure 8—figure supplement 1A-C).

https://cdn.elifesciences.org/articles/76090/elife-76090-fig8-data1-v2.xlsx
Figure 8—source data 2

Immunoblot images (raw and annotated) of inhibitor precipitation assay and cleavage assays (Source data for Figure 8C, E-G).

https://cdn.elifesciences.org/articles/76090/elife-76090-fig8-data2-v2.zip
Figure 8—source data 3

Raw values of immunoblot quantitation (Source data for Figure 8D).

https://cdn.elifesciences.org/articles/76090/elife-76090-fig8-data3-v2.xlsx

To experimentally test whether the lipid environment influences the active site geometry of PSH, we used the L-685,458-based biotinylated affinity ligand Merck C (Beher et al., 2003) to capture PSH in DDM micelles and in POPC vesicles. As shown in Figure 8C, Merck C was able to capture specifically PSH as judged from binding competition in the presence of excess amounts of the parental compound L-685,458. In all these experiments, binding competition was stronger in the POPC vesicles than in DDM micelles indicating that the more labile DDM environment also weakens the competition of binding with the parental compound. In addition, also the background level of unspecific binding was higher in the latter environment contributing to the higher levels of unspecific PSH capture. In agreement with these observations, quantitation of specifically Merck C-bound PSH showed that the capture was markedly enhanced for the protease in the POPC bilayer (Figure 8D). Additional enzyme inhibition experiments further showed that both inhibitors, the parental L-685,458 as well as Merck C, inhibited PSH less well in DDM than in POPC (Figure 8E and F). All in all, these findings support the interpretation that the POPC bilayer stabilizes the active site, whereas it is destabilized in DDM micelles. Thus, these data suggest that the lipid environment stabilizes the geometry of the active site, which translates into the increased processivity of PSH in POPC vesicles.

Despite the improved binding and capture of PSH by the L-685,458-derived affinity probe in the POPC membrane environment, a potent inhibition of PSH by L-685,458 required rather high micromolar concentrations of this GSI. We thus finally tested whether other known GSIs would be more effective in inhibiting PSH cleavage of C99. Besides L-685,458, another TSA inhibitor (III-31C Esler et al., 2002, Figure 8—figure supplement 3A, B) and four non-TSA inhibitors with comparable potency (DAPT Dovey et al., 2001, LY411575 Lanz et al., 2004, Begacestat Mayer et al., 2008 and MRK-560 Best et al., 2006, Figure 8—figure supplement 3C-F) were tested for their potential to inhibit reconstituted PSH at pH 7.0. Remarkably, only the TSA inhibitor III-31C was able to inhibit C99 cleavage by PSH, whereas the non-TSA inhibitors were largely ineffective, even when used at the same high concentrations as for L-685,458 (20 µM) (Figure 8G). Overall, these inhibition data suggest that the stabilizing interactions of the TSA inhibitors in the PSH active site region may be different from that of γ-secretase (Yang et al., 2021; Hitzenberger and Zacharias, 2019). Such differences in stabilizing interactions might also affect the other GSIs that may bind too weakly to inhibit the enzyme.

Taken together, the MD simulations indicate a more stable, that is less fluctuating, geometry of the enzyme-substrate binding state around the enzyme active site in a POPC membrane environment compared to a DDM micelle environment. The DDM environment appears to destabilize important structural elements such as the TMD6a that is required for a stable enzyme-substrate interaction. These results can qualitatively explain the experimentally observed reduced processivity of PSH in DDM micelles and its boost in the POPC bilayer. Further, they emphasize the critical importance of the membrane environment for the formation of a conformationally stable active site geometry in the E-S complex, which is key for the efficient operation of intramembrane proteases in general.

Discussion

It has previously been demonstrated that the archaeal intramembrane protease PSH cleaves the APP substrate C99 into Aβ40 and Aβ42 in a manner very similar to γ-secretase (Dang et al., 2015). PSH can thus be used as a surrogate for γ-secretase allowing to study the proteolytic activity of its catalytic presenilin subunit in the absence of its complex partners. Here, we confirm and extend these prior findings by a further, more in depth characterization of C99 processing by PSH. We first found that detergent-solubilized PSH cleaves C99 in DDM micelles with a reduced processivity as evident from higher amounts of Aβ42 than Aβ40. The reduced processivity of PSH under these conditions was supported by the identification of longer Aβ species such as Aβ46, cleavage products that were not identified in the previous study. Strikingly, we found in our assay system that the processivity was strongly enhanced in a membrane bilayer when the enzyme was reconstituted into POPC SUVs. Under these conditions, PSH processivity was strongly promoted as seen by the increased production of Aβ38. The protease was pH-dependent and showed the highest activity in the mild acidic to mild alkaline pH range in both micelle and bilayer conditions. In the POPC bilayer, the processivity of PSH was increased up to neutral pH before it rapidly dropped in the alkaline pH range of 7.5–8.5, where longer Aβ species started to accumulate eventually remaining unprocessed. Although the pH/activity profile of PSH showed overall similarities to that of γ-secretase, there were some notable differences. Compared to γ-secretase, which has a pH optimum of 6.5 (Quintero-Monzon et al., 2011), that of PSH was shifted to neutral pH. Moreover, although the processivity of γ-secretase was increasingly impaired toward more alkaline conditions like for PSH shown here, paradoxically, relative increases of Aβ38 were observed for γ-secretase in parallel (Quintero-Monzon et al., 2011). Clearly, the most remarkable observation, however, was the rise in processivity when PSH was reconstituted into POPC membranes.

In searching for the underlying basis of these dramatic activity changes when changing from a micelle to a bilayer environment, we asked if these could be due to potential structural rearrangements that PSH undergoes in these two different environments. To investigate this possibility, PS1-based homology models were generated for PSH in the APP C83 substrate-bound holo form and the substrate-free apo form. The models revealed both β2-strand and TMD6a as structural elements, which we found by mutational analysis to be functionally highly critical for substrate cleavage by PSH. As observed for substrate-bound PS1 (Zhou et al., 2019; Yang et al., 2019), this suggests that they constitute important structural elements for substrate binding also for PSH.

The conformational dynamics of these structural models was evaluated in MD simulations to test whether structural changes might be observable that could explain the activity changes. We note that in previous simulations the atomistic dynamics of γ-secretase and the interaction with C83 have been studied (Mehra et al., 2020; Bhattarai et al., 2020), including also its activated state poised for ε-site cleavage (Bhattarai et al., 2020). However, comparative simulations in detergent micelles and lipid bilayer have so far not been performed. Among three models, model 2 is considered as the most realistic model and was chosen as working model for the E-S in our study. Our simulation results on PSH clearly showed more structural fluctuations of the protease in the micelle environment than in the membrane environment. These translated into less stable interactions with the substrate in the micelle compared to that in the bilayer, particularly of TMD6a with C83 in the active site region. In line with our mutational analysis, changing residues within the hydrophobic patch of TMD6a, which interacts with residues V50 and L52 at or near, respectively, the ε49-site of the substrate, disrupted interactions with the substrate in the MD simulations and strongly interfered with substrate cleavage in the PSH cleavage assays thus linking functional biochemical data with structural dynamics of PSH. TMD6a thus further emerges as an important structural element for substrate interaction that appears to be able to sense changes in the hydrophobic environment of the protease. The increased stability of TMD6a in POPC reflects a stabilized enzyme-substrate interaction that could likely translate into the enhanced processivity observed for the membrane environment. A more stable interaction increases the substrate residence time at the enzyme that has been shown for γ-secretase to be key for its processivity (Okochi et al., 2013) and that is also supported by MD simulations of a γ-secretase–C99 complex (Dehury et al., 2019). Moreover, a closer distance of the catalytic aspartate residues was much more frequently observed for the substrate-bound PSH holo form in the POPC bilayer suggesting that the formation of an active site geometry capable of peptide bond hydrolysis is promoted in the membrane environment. Consistent with these data, we found that binding of PSH to Merck C, a biotinylated derivative of L-685,458, was increased in the POPC bilayer. As shown previously for γ-secretase, L-685,458 interacts with the same subsite pockets as C83 and occupies a position of the substrate in the active site region close to where also the ε-cleavage sites of C83 become exposed and are unfolded (Yang et al., 2021; Hitzenberger and Zacharias, 2019).

Despite the demonstration of direct binding of the L-685,458 lead structure to PSH using the Merck C affinity ligand, L-685,458 inhibited PSH much less efficiently than γ-secretase, that is micromolar concentrations were needed to inhibit PSH compared to nanomolar concentrations known to inhibit γ-secretase. Likewise, and consistent with previous results (Dang et al., 2015), the related TSA inhibitor III-31C could inhibit C99 cleavage of PSH but again at micromolar concentrations. Since other non-TSA GSIs failed to inhibit PSH, only TSA inhibitors can interact with PSH and effectively inhibit the enzyme. This suggests that the binding sites for non-TSA GSIs are different or, more likely, that their interactions with PSH are too weak to inhibit the enzyme. As shown previously, the binding sites of the non-TSA GSIs Avagacestat and Semagacestat are similar to the binding site of the TSA GSI L-685,458 (Yang et al., 2021). The non-TSA GSIs occupy the position of the β-strand of the substrate but do not protrude to the catalytic site resulting in decreased interactions with γ-secretase compared to the TSA GSI (Yang et al., 2021).

Both our experimental studies and the corresponding comparative MD simulations therefore suggest that the higher conformational flexibility of PSH in micelles causes destabilized interaction with C83 and C99 and consequently a reduced processivity. In contrast, a lipid bilayer induces a less flexible conformation of PSH that allows a more stable interaction with substrate, thereby promoting the processivity of PSH. Our data support recent findings for γ-secretase that showed differences in the processivity in a phospholipid/detergent-based versus a lipid raft-like membrane environment (Szaruga et al., 2017) and now provide an underlying molecular basis for this behavior. Similar to an artificial destabilization of the PSH/presenilin fold in detergent micelles, computational analyses suggest that FAD mutations in presenilin cause structural destabilizations (Chen and Zacharias, 2020; Somavarapu and Kepp, 2016; Tang et al., 2019) which are consistent with the experimentally observed impact on substrate interactions of these mutants (Fukumori and Steiner, 2016; Trambauer et al., 2020) and their alteration of APP/Aβ E-S stabilities resulting in processivity impairments (Szaruga et al., 2017). We also note that a less stable E-S in detergent micelles might account for differences in cleavage site usage and in inhibition profiles for diverse C99-based substrates that were observed in previous PSH assays (Torres-Arancivia et al., 2010; Dang et al., 2015; Naing et al., 2018).

Presumably, due to their non-native environment, the available structures of GxGD-type proteases show catalytically inactive conformations with too distant catalytic residues. The large distance of the catalytic aspartates of PSH in the substrate-free apo form is also seen for γ-secretase (10.6 Å, Bai et al., 2015b) as well as in different GxGD-type aspartyl proteases like FlaK (12 Å, Hu et al., 2011) and seems to represent their inactive form. Upon substrate interaction, this distance is decreased bringing the two catalytic aspartates closer to the initial cleavage sites (Zhou et al., 2019; Yang et al., 2019). As now shown in our study, a lipid bilayer environment promotes the formation of a stable active-site geometry by bringing the catalytic residues, water and the substrate scissile bonds into a conformation that allows proteolysis to proceed more efficiently. As a general implication for intramembrane proteolysis, our data suggest that a lipid bilayer-mediated stabilization of the active-site geometry might also be observable for other intramembrane proteases of different catalytic types.

Taken together, in good correlation between experimental and simulation data, our results with PSH as a model intramembrane protease highlight an important role of the membrane lipid environment in providing a stabilized E-S conformation that is crucial for substrate processing in intramembrane proteolysis. Our data further underscore a key role of the conformational flexibility of presenilin/PSH TMD6a for substrate interactions and proteolytic cleavage of presenilin-type proteases. Most importantly, they provide evidence that the lipid bilayer promotes the formation of a conformationally stable active site geometry, which is of general importance for an efficient catalytic operation of intramembrane proteases.

Materials and methods

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Strain, strain background (Escherichia coli)BL21(DE3)RILAgilent TechnologiesCat# 230245
Recombinant DNA reagentpQE60-C100-His6Edbauer et al., 2003N/A
Recombinant DNA reagentpQE60-C83-His6This studyN/A
Recombinant DNA reagentpET21b-PSHLi et al., 2012N/AGift from Yigong Shi
AntibodyAnti-APP (C-terminus) Y188
(rabbit monoclonal)
AbcamCat# ab32136IB (immunoblot) (1:5000)
AntibodyAnti-APP (C-terminus) 6687
(rabbit polyclonal)
Steiner et al., 2000N/AIP (immunoprecipitation)
(1:100–1:200)
AntibodyAnti-APP (Aβ22–35) Aβ (22-35)
(rabbit polyclonal)
Sigma-AldrichCat# A3356IB (1:1000)
AntibodyAnti-APP (Aβ1–16) 2D8
(mouse monoclonal)
Shirotani et al., 2007N/AIB (3 µg/ml)
AntibodyAnti-APP (Aβ17–24) 4G8
(mouse monoclonal)
BioLegendCat# 800702IB (1:500-1:2500)
AntibodyAnti-PSH (residues 192–204)
6F4
(rat monoclonal)
This studyN/AIB (3 µg/ml), generation of antibody described further below
AntibodyAnti-Penta-His
(mouse monoclonal)
QiagenCat# 34660IB (1:1000)
Chemical compound, drugNi-NTA AgaroseQiagenCat# 30210
Chemical compound, drugCalbiosorb Adsorbent beadsCalbiochemCat# 206550Discontinued
Chemical compound, drugPOPCAvanti Polar LipidsCat# 850457PPowder
Chemical compound, drugRhodamine-DHPEInvitrogenCat# L1392
Chemical compound, drugSephacryl S-200 HRGE HealthcareCat# 17058410
Chemical compound, drugStreptavidin SepharoseGE HealthcareCat# 17511301
Chemical compound, drugL-685,458Sigma-AldrichCat# 565771InSolution γ-Secretase Inhibitor X, used in cleavage assays
Chemical compound, drugL-685,458Sigma-AldrichCat# L1790Powder, dissolved in DMSO and used in inhibitor affinity precipitation experiments
Chemical compound, drugMerck CTaros ChemicalsN/ABiotinylated L-685,458
Chemical compound, drugn-Dodecyl β-D-maltoside (DDM)MilliporeCat# 324355
Chemical compound, drugProtein G SepharoseCytivaCat# 17061801
Chemical compound, drugProtein A SepharoseCytivaCat# 17528001
Chemical compound, drugTropix
I-BLOCK
InvitrogenCat# T2015
Chemical compound, drugIII-31CSigma-AldrichCat# C0619
Chemical compound, drugDAPTBoehringer Ingelheim Pharma KGN/A
Chemical compound, drugLY411575Karlheinz BaumannN/A
Chemical compound, drugBegacestatKarlheinz BaumannN/A
Chemical compound, drugMRK-560Karlheinz BaumannN/A
Commercial assay or kitV-PLEX Plus Aβ Peptide Panel 1 (4G8) KitMeso Scale DiscoveryCat# K15199G
Commercial assay or kitNativePAGE 4 to 16%, Bis-Tris, 1.0 mm, Mini Protein Gels, 10 wellsInvitrogenCat#
BN1002BOX
Software, algorithmGelAnalyzer 19.1Istvan Lazar Jr., PhD Istvan Lazar Sr., PhD, CScN/Ahttp://www.gelanalyzer.com
Software, algorithmAMBER18Case et al., 2005N/A
Software, algorithmCHARMM-GUIJo et al., 2008N/A
Software, algorithmSWISS-MODELWaterhouse et al., 2018N/A
Software, algorithmPROPKA3.1Olsson et al., 2011; Sondergaard et al., 2011N/A
Software, algorithmDSSPKabsch and Sander, 1983; Touw et al., 2015N/A

Monoclonal antibody generation

Request a detailed protocol

Monoclonal antibody 6F4 (IgG2b/k) to PSH was raised in Wistar rat against amino acid residues 192–204 (KRADYSFRKEGLN) of PSH from Methanoculleus marisnigri.

PSH constructs

Request a detailed protocol

All constructs are based on the PSH expression construct in pET-21b used for structure determination (Li et al., 2012). PSH point mutations and deletion were generated using site-directed mutagenesis. DNA sequencing of the newly generated plasmid confirmed successful mutagenesis.

PSH expression and purification

Request a detailed protocol

Expression and purification of WT and mutant PSH was adopted from the published protocol (Li et al., 2012). In brief, E. coli BL21(DE3)RIL cells transformed with the pET-21b vector harboring an N-terminal 8 x His-tagged PSH were grown in LB medium to an optical density of 1.5 and expression was induced with 0.2 mM isopropyl β-D-1-thiogalactopyranoside (IPTG). PSH was expressed at 22 °C overnight and harvested cells were resuspended in resuspension buffer (25 mM Tris-HCl, pH 8.0, 150 mM NaCl). Cells were lysed by sonication; cell debris was removed by centrifugation and membranes were collected by ultracentrifugation at 150,000 x g for 1 hr. Membranes were solubilized in resuspension buffer containing 2% DDM by rocking at 4 °C for 2 hr. After ultracentrifugation at 150,000 x g for 30 min, the supernatant was incubated with Ni-NTA agarose beads (Qiagen) for 2 hr at room temperature. Beads were then washed with resuspension buffer containing 20 mM imidazole and 0.6% (w/v) DDM. PSH was eluted with resuspension buffer containing 250 mM imidazole and 0.6% (w/v) DDM. Correct folding of WT and mutant PSH was confirmed by dynamic light scattering (DLS, Malvern Instruments High Performance Particle Sizer) (Appendix 1—figure 1A), Blue Native (BN)-PAGE (Appendix 1—figure 1B) and nano differential scanning fluorimetry (nanoDSF, NanoTemper Tycho) (Appendix 1—figure 1C). For DLS, protein samples (25 µM) were analyzed in a Hellma Analytics High Precision Cell. For BN-PAGE, samples were prepared as described (Schägger and von Jagow, 1991) and separated using a Novex NativePAGE 4–16% Bis-Tris gel. Following electrophoresis, the gel was prepared for blotting as described (Winkler et al., 2009) and WT and mutant PSH were subjected to immunoblot analysis with antibody 6F4. For nanoDSF, protein samples (25 µM) were loaded into NanoTemper Tycho NT.6 capillaries, unfolding profiles of WT and mutant PSH were recorded, and the inflection temperatures (Ti) were obtained by automated data analysis.

PSH reconstitution in POPC vesicles

Request a detailed protocol

PSH reconstitution into POPC SUVs was based on the reconstitution of γ-secretase into SUVs (Winkler et al., 2012). SUVs were prepared in a low citrate buffer (5 mM sodium citrate, 3.5% glycerol, pH 6.4) and diluted 2.5 times with buffer (5 mM sodium citrate, 3.5% glycerol, 30 mM DTT, pH 6.4). One volume of purified PSH and four volumes of the vesicle preparation were mixed in the presence of an excess of Calbiosorb adsorbent beads (Calbiochem) and incubated at 4 °C overnight to allow the formation of proteoliposomes.

Validation of PSH reconstitution into POPC SUVs

Request a detailed protocol

To validate the incorporation of PSH into POPC vesicles, PSH was reconstituted into POPC vesicles containing the fluorescent marker lipid rhodamine-DHPE. These vesicles were then subjected to a small gel filtration column packed with Sephacryl S-200 HR to separate vesicles and free PSH. The vesicle content of each fraction was analyzed by measuring fluorescence (λex 530 nm, λem 590 nm) with Fluoroscan Asket Fl (Labsystems) and the presence of PSH in the fractions was analyzed by immunoblotting with antibody 6F4.

APP substrate constructs

Request a detailed protocol

Recombinant APP substrate C100-His6 was described before (Edbauer et al., 2003). The corresponding C83-His6 (containing an N-terminal methionine) was generated by PCR and cloned into pQE60.

Expression and purification of APP-based substrates

Request a detailed protocol

C100-His6 and C83-His6 were expressed in E. coli BL21(DE3)RIL cells after induction with IPTG at 37 °C for 4 hr. Cell pellets were resuspended in TE buffer (20 mM Tris (pH 7.5), 1 mM EDTA), sonified and inclusion bodies were collected by centrifugation. Inclusion bodies were lysed overnight at 4 °C in 20 mM Tris (pH 8.5), 6 M urea,1 mM CaCl2, 100 mM NaCl, 1% (w/v) SDS and 1% (v/v) Triton X-100 by rotation. The lysate was diluted 1:5 with 20 mM Tris (pH 7.5) and 150 mM NaCl and then incubated with Ni-NTA agarose for 2 hr at room temperature. Ni-NTA beads were washed extensively with TX-wash buffer (50 mM Tris (pH 8.5), 300 mM NaCl, 1% (v/v) Triton X-100), SDS-wash buffer (50 mM Tris (pH 8.5), 300 mM NaCl, 0.2% (w/v) SDS) and imidazole wash buffer (50 mM Tris (pH 8.5), 300 mM NaCl, 0.2% (w/v) SDS, 20 mM imidazole) before the elution of bound protein with elution buffer (50 mM Tris (pH 8.5), 300 mM NaCl, 0.2% (w/v) SDS, 150 mM imidazole).

PSH in vitro assay

Request a detailed protocol

The in vitro assays using recombinant APP substrates C100-His6 and C83-His6 were performed with either DDM-solubilized PSH or reconstituted PSH. To test PSH activity over a range of different pH values, the in vitro assay was performed in the presence of a master buffer (20 mM Bicine, 20 mM HEPES, 20 mM MES) adjusted to different pH values (pH 5.5–9.0). 1–2 µM PSH was incubated with 0.5 µM substrate overnight at 37 °C. The final DDM concentration was 0.02% for the assays in DDM micelles. For the assays with PSH reconstituted in POPC vesicles a small amount of DDM was added just below the critical micelle concentration (CMC) of 0.0087%. As for γ-secretase, detergent addition below the CMC is necessary to achieve enzyme activity after reconstitution (Winkler et al., 2012). Generated AICD, Aβ and p3 were analyzed by immunoblotting and in addition, Aβ and AICD species were determined by MALDI-TOF mass spectrometry (MS) analysis as described previously (Winkler et al., 2009; Page et al., 2008; Ebke et al., 2011). In brief, samples were diluted with IP-MS buffer (10 mM Tris (pH 8.0), 140 mM NaCl, 0.5 mM EDTA, 0.1% n-octyl-glucopyranoside) and immunoprecipitated for 16 hr at 4 °C with antibody 4G8 and protein G Sepharose for Aβ species or with antibody 6687 and protein A Sepharose for AICD species.

Electrochemiluminescence immunoassay (ECL-IA)

Request a detailed protocol

Aβ and p3 species from PSH in vitro assays were analyzed with the V-PLEX Plus Aβ Peptide Panel 1 (4G8) Kit (Meso Scale Discovery, K15199G) using SULFO-tagged anti-Aβ antibody (4G8) in a 1:50 dilution. Samples were centrifuged for 30 min at 21,000 x g and then diluted 1:25 to reach a sample concentration in the linear detection range. The ECL-IA was performed following the manufacturer’s protocol.

Inhibitor affinity precipitations of PSH

Request a detailed protocol

Streptavidin Sepharose beads (GE Healthcare) were washed three times with PBS and then blocked overnight with 1% I-Block in PBS at 4 °C and additionally for 30 min at room temperature. PSH in DDM micelles or in POPC vesicles was diluted with MES-buffer (50 mM MES pH 6.0, 500 mM NaCl, 5 mM MgCl2, 1 x PI mix complete (Roche)) to 1–2 µM PSH. To mimic the cleavage assay condition, a small amount of DDM was added just below the CMC of DDM (0.0087%) to the POPC reconstituted PSH. The diluted PSH solution was precleared with streptavidin beads for 30 min at 4 °C. To capture PSH, the precleared solution was incubated with 20 µM L-685,458-based biotinylated TSA inhibitor Merck C (Beher et al., 2003) (Taros Chemicals) in the presence of streptavidin beads for 2 hr at room temperature. To analyze non-specific binding, Merck C was omitted or a 100-fold molar excess of the parental compound L-685,458 was added. To quantify Merck C GSI binding to PSH, the chemiluminescence signal of the respective immunoblots were quantified using the LAS−4000 image reader (Fujifilm Life Science) and GelAnalyzer 19.1 software (http://www.gelanalyzer.com). Specific binding was calculated as the difference between the binding of PSH and the binding of PSH in presence of the competitor L-685,458 after subtraction of unspecific PSH binding to the beads.

Molecular dynamics simulations

Request a detailed protocol

The available crystal structures of PSH, PDB 4HYG (Li et al., 2012) and PDB 4Y6K (Dang et al., 2015) include several amino acid substitutions compared to the WT PSH sequence and in addition, several important loop segments are missing. In both crystals, the enzyme forms a tetramer that may also stabilize a structure different from the solution and substrate-bound conformation. Under the assumption that a substrate is bound to PSH in an analogous fashion as in the PS1 homolog, we used the option to generate a comparative model structure using the C83 substrate-bound γ-secretase structure (PDB 6IYC) (Zhou et al., 2019) as a template. Three models were generated using the SWISS online server (Waterhouse et al., 2018) and the program MODELLER (Eswar et al., 2006) with a same sequence alignment strategy (Figure 4—figure supplement 2) but different homology protocols. Model 1 was built by taking only PDB 6IYC as the template as generated by the SWISS online server (Table 1). Model 2 was generated by using residues ranging from L7 to D162 and D220 to L292 from chain B of PDB 4HYG and model 1 as templates with the MODELLER multi-template method (Table 1). Model 3 was built by taking all residues resolved in chain B of PDB 4HYG and model 1 with the MODELLER multi-template method (Table 1). Similar to the PS1 template structure, the final holo PSH structure is composed of two fragments with an N-terminal fragment from L7-R193 and a C-terminal fragment from E210 to A293. The generated holo-state PSH structures were then embedded in two different environments: micelle capsules consisting of 150 DDM molecules and a membrane bilayer with 302 POPC molecules. An additional, larger micelle capsule system of model 2 was constructed with 50% more, namely 225, DDM, molecules. Lysine mutations M172K, I173K, L175K and A176K were constructed based on model 2 with RMSD (WT vs mutant) < 0.1 Å, and embedded in a membrane bilayer system with 302 POPC molecules using the CHARMM-GUI online server (Cheng et al., 2013; Wu et al., 2014). All 11 systems were prepared and solvated in explicit TIP3P water (Jorgensen et al., 1983) at a salt concentration of 0.15 M KCl using the CHARMM-GUI online server .

The interaction of proteins, lipid, and micelles is described by the charmm36m force field (Huang et al., 2017). Each system was simulated using the AMBER18 pmemd GPU accelerated version (Case et al., 2005) in combination with a Berendsen barostat (1 bar) and a Langevin thermostat (303.15 K). The hydrogen mass repartitioning method was used allowing a time step of 4 fs. Three simulations with 600 ns each were performed for each system, in total 33 NPT trajectories were generated for further analysis. Non-bonded cutoff was set to 12 Å with a force-based switching distance of 10 Å. D220 was selected to be protonated while D162 was unprotonated according to the pKa prediction on the existing PSH and PS1 structures by PROPKA3.1 (Olsson et al., 2011; Sondergaard et al., 2011; Table 2).

Table 2
pKa predictions as calculated by PROPKA3.1 for published PSH and γ-secretase structures.

pKa values of the catalytic aspartate residue which is most likely protonated are indicated in red.

PDB IDEnyzmeLigandpKa (D162, D220)
4HYGPSHNone5.04, 6.62
4Y6KPSHIII-31-C5.63, 7.52
PDB IDEnyzmeLigandpKa (D257, D385)
4UISPS1None3.18, 6.16
5A63PS1None4.42, 6.16
5FN5PS1None4.98, 3.63
5FN4PS1Unknown helix4.70, 4.71
5FN3PS1Unknown helix4.90, 7.13
5FN2PS1DAPT5.13, 9.93
6IYCPS1C836.39, X*
6IDFPS1Notch16.21, X*
6LR4PS1Semagacestat6.12, 7.94
6LQGPS1Avagacestat6.08, 7.22
7V9IPS1L-685,4587.11, 8.90
7D8XPS1L-685,458 and E20127.01, 8.69
  1. *

    For structure determination D385 was mutated to alanine and therefore no pKa value is given.

The lipid tail order parameter SCH was computed in model 2 in DDM and POPC environments to show the orientation and the ordering of the concerning CH vector (Tieleman et al., 1997; Vermeer et al., 2007) with respect to the protein principle axis, which was aligned to the lipid normal in the POPC environment. In addition, the area per lipid was computed in model 2 on both leaflets to verify the reliability of our POPC lipid model.

RMSF of PSH and C83 curves were calculated by taking the last 200 ns with the time-average PSH structure of each simulation as the reference and only taking the backbone atoms for the calculation. Secondary structure of PSH TMD6a was calculated using the DSSP method (Kabsch and Sander, 1983; Touw et al., 2015).

Appendix 1

Protein quality control

To check for the quality of the protein preparations, WT and mutant PSH were analyzed by DLS, BN-PAGE and nanoDSF. In DLS experiments, the Z-Average (Z-Ave) value for WT and mutant PSH was below 100 nm, which indicates that the different PSH preparations were not aggregated (Appendix 1—figure 1A). BN-PAGE immunoblot analysis showed a band for monomeric and dimeric PSH for all constructs but no higher molecular weight aggregate formation was observed (Appendix 1—figure 1B). nanoDSF experiments showed an inflection temperature (Ti) of round about 72.5 °C for all constructs (Appendix 1—figure 1C) indicating that the introduction of single point mutations into PSH does not influence the thermal stability and folding of the protein.

Appendix 1—figure 1
Quality control of WT and mutant PSH.

(A) Analysis of protein aggregation of WT and mutant PSH by DLS. (B) Analysis of protein aggregation of WT and mutant PSH by BN-PAGE followed by immunoblotting for PSH (6F4). (C) Analysis of protein misfolding of WT and mutant PSH by nanoDSF.

Data availability

For all figures the source data are provided in the respective source data files. The coordinate and trajectory files of all simulations can be accessed at Zenodo: https://doi.org/10.5281/zenodo.6487373.

The following previously published data sets were used
    1. Zhou R
    2. Yang G
    3. Guo X
    4. Zhou Q
    5. Lei J
    6. Shi Y
    (2019) RCSB Protein Data Bank
    ID 6IYC. Recognition of the amyloid precursor protein by human γ-secretase.
    1. Li X
    2. Dang S
    3. Yan C
    4. Gong X
    5. Wang J
    6. Shi Y
    (2013) RCSB Protein Data Bank
    ID 4HYG. Structure of a presenilin family intramembrane aspartate protease.

References

Decision letter

  1. M Joanne Lemieux
    Reviewing Editor; University of Alberta, Canada
  2. José D Faraldo-Gómez
    Senior Editor; National Heart, Lung and Blood Institute, National Institutes of Health, United States
  3. M Joanne Lemieux
    Reviewer; University of Alberta, Canada
  4. Sean Workman
    Reviewer; University of Regina, Canada

In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.

Decision letter after peer review:

[Editors’ note: the authors submitted for reconsideration following the decision after peer review. What follows is the decision letter after the first round of review.]

Thank you for submitting your work entitled "Active site geometry stabilization of a presenilin homolog by the lipid bilayer promotes intramembrane proteolysis" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, including Joanne Lemieux as the Reviewing Editor and Reviewer #1, and the evaluation has been overseen by a Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Michael S. Wolfe (Reviewer #3).

Our decision has been reached after a consultation among editors and reviewers. Based on these discussions and the individual reviews below, we regret to inform you that eLife will not be considering this submission for publication. We recognize the importance of this work and the progress made towards a better understanding of the factors that influence intramembrane proteolysis. Indeed, the reviewers' evaluation of the experimental work was positive overall. However, the reviewers also identified critical problems in the computational component. Although it seems possible to address these problems, it also seems that the revision would exceed the two-month period that is typically allowed. Hence, we believe the appropriate action is to reject the current version of your manuscript and allow you to consider your options. Should you decide to address all the concerns raised by the reviewers, we would be willing to examine a new version of your manuscript – but please note it would be considered a new submission, and therefore, that it might not be sent for review or be evaluated by the same reviewers.

Reviewer #1:

The activity of the archeal presenilin homolog (PSH) is assessed using a portion of the amyloid precursor substrate C99 substrate in DDM detergent and also in a POPC liposome environment. Mass spectrometry was used to precise indicate cleavage species as well as SDS. The results clearly show lipid influence both the rate of cleavage and the species generates. Modelling was used to determine structures of the PSH with the substrate in both app and holo forms. Similar to the cryoEM structures of PS with Notch and APP, a hybrid enzyme-substrate beta-sheet is observed between substrate and active site, suggesting a valid model for MD situations. These were used in MD simulations in both detergent and lipid environments. Importantly protonation states were analyzed for this Aspartyl protease, in which the two Asp residues were evaluated independently. D220 was determined to be the protonated Asp with D162 left charged in the simulations in different environments. The MD simulations show there were less structural changes in the DOPC embedded protein compared to the non-lipid sample, with a notable shift in TMD6a only in the DDM sample . MD simulations support the above experimentally observations that the lipid stabilize the E-S complex.

Overall this is an interesting paper that combines both experimental and structural studies to rationale differences in membrane environment and protease activity.

Reviewer #2:

I have several comments and suggestions for the authors to consider improving the quality of the paper:

1. Utilization of C100 for enzymatic analysis.

Authors utilized the recombinant C100-His6 as a substrate for the enzymatic assay. In contrast, they modeled the complex structure of PSH with C83, as the template cryo-EM structure was PS1-C83 complex. Essentially γ-secretase cleaves several substrates regardless the sequence. However, N-terminal length of the substrate might affect the proteolytic efficacy of the γ-secretase. In fact, Funamoto et al. (Nat Commun 2013) reported that AICD production from C83-FLAG was much greater than that from C99-FLAG (i.e., distinct Km and Vmax values). These data suggest the possibility that the formation of E-S complex is regulated by the N-terminal length of the substrate in the proteolytic mechanism of the γ-secretase. Thus, I would recommend the authors to compare the cleavage and trimming patterns of C83 by PSH in either DDM or POPC, as shown in the modeling analysis and the MD simulation. Such comparison strengthens the author's conclusion that the stabilization of E-S complex is critical to the intramembrane proteolysis.

2. Importance of TMD6a and hybrid β-sheet in the proteolysis by PSH.

Modeled structure of PSH-C83 strongly implicates the importance of these two structural elements in the proteolysis of PSH. However, it remains unclear whether these elements are truly required for the proteolytic activity of PSH. The possibility that these structures artificially modeled because the authors used the γ-secretase-C83 structure as a template. Thus, authors should test the proteolytic activity of mutant PSH that carries amino acid substitutions in TMD6a or beta2-strand to abolish the interaction with the substrate.

3. Figure presentation

In figure 5A, it is difficult to understand the difference of TMD6a conformation in DDM/POPC because of superimposed structure. They should be presented separately.

Reviewer #3:

In this manuscript, Feilen et al. combined biochemical experiments and molecular dynamics (MD) simulations to investigate the effects of detergent solubilization versus lipid vesicle reconstitution on intramembrane protease activity of an archaeal presenilin homolog (PSH). This PSH has been previously reported to process amyloid precursor (APP)-based substrate to amyloid β-peptides (A-β) in a manner closely similar to that accomplished by the presenilin-containing γ-secretase complex. The archaeal PSH is employed here as a surrogate to gain mechanistic insight into γ-secretase.

The authors showed that the carboxypeptidase-like activity of PSH is impaired in DDM micelles compared to a POPC lipid bilayer. Evidence is also provided suggesting that DDM-solubilized PSH binds more weakly to a transition-state analog inhibitor compared with POPC-bilayer PSH. Comparative MD simulations suggested that the lipid bilayer stabilized relevant structural elements of PSH for substrate binding (notably TMD6a) and formation of the enzyme active site geometry for proteolysis. A number of suggestions that could help improving the manuscript include the following:

1. While it is understood that archaeal PSH is taken as a surrogate for presenilin in the γ-secretase complex, the authors should explain why this is necessary. All the described experiments, including the MD simulations, could have been conducted with γ-secretase itself. In the discussion, the specific implications for γ-secretase (especially FAD mutations) should be de-emphasized and more emphasis put on the implications for intramembrane proteolysis in general.

2. For the MS analysis of A-β peptide products in Figure 1D, a table of observed vs. calculated m/z should be provided. Some of these peaks (A-β-43, -45, and -46) are quite weak. The same is true of AICD MS analysis in Figure 2D.

3. Given crystal structures of the PSH, is it more reasonable to preserve protein coordinates in these crystal structures, but only add the missing residues (e.g., TMD6a) using the PS1/γ-secretase cryo-EM structures as template? The Results section is vague about how the homology modeling of enzyme-substrate complex was generated; an additional sentence or two should be provided so the reader does not have to refer to the experimental section to answer this basic question.

4. Assuming pKa calculations depend on local geometry of the protonation site, the protein structure(s) used for the calculations need to be described clearly. Moreover, how do the residue pKa's depend on the protein structures (e..g, 4Y6K and 4HYG crystal structures of PSH, the apo and holo cryo-EM structures of PS1/γ-secretase, and simulation equilibrated protein structures in notably two different conformations with the D162-D220 distance centered around ~6.5 and ~8.2 Angstroms in Figure 6A).

5. It is unclear why the Amber force field was used in simulations of holo PSH with two different protonation states, but CHARMM36m in simulations of the apo and holo PSH in different membrane environments. It would help to evaluate the force field differences and potential effects by adding simulations using CHARMM36m to the holo PSH with different protonation states or simulations using Amber to the apo and holo PSH in different membrane environments.

6. In Figure 3A, what is the RMSD between the PSH homology model and its crystal structures?

7. In Figures 4C-4E, it would help to add error bars (standard deviations) to examine what differences are significant. The authors calculated RMSFs of PSH in the apo and holo forms to describe its stability in different lipid environments. It is better to show error bars (with total simulation times across different replicates) to describe RMSF differences in notably, the TMD6a, TMD4 and TM2-TM3 loop. There is a notable RMSF difference just beyond TMD6a, which should be mentioned and explained. In addition to RMSF, further simulation analysis such as comparison of the distances between the catalytic aspartate and scissile peptide bond in C83 could provide more insights.

8. In Figure 5A, it would help to quantitatively calculate and plot the helicity of TMD6a and/or secondary structures of residues in TMD6a as a function of simulation time and compare these quantities between the different simulated systems.

9. In Figure 5C-5D, the authors described that the DDM molecule can insert itself between TMD2 and TMD6 and intervenes intra- and intermolecular interactions and thus destabilize TMD6a. It would be better to have more explanation if these insertions are with just one lipid molecule or there are multiple molecules involved during different time frames of the simulations. In addition, it would be more convincing to see atomic detailed interaction between the DDM molecule and the protein, especially because DDM is a nonionic detergent. Moreover, it could help to calculate -SCD order parameters that are usually obtained from NMR experiments to measure orientational anisotropy of the C-H bonds lipid chains and quantify differences of the lipid orientations.

10. For Figure 6A, it would help to plot the D162-D220 distance vs. time in simulations of the different systems also. What are the corresponding distance values in the PSH crystal and PS1 cryo-EM structures? In addition to a main peak at ~6.5 Å distance, there seems another peak at ~8.2 Å distance between D162-D220; what is this conformational state?

11. In Figure 6B, there is substantial non-specific binding of the biotinylated Merck C to PSH in DDM; the parent inhibitor is essentially not competing, even at 20 microM. This makes the interpretation that specific binding is stronger in POPC vesicles less convincing. Some comment to this effect should be added to the Results section.

Reviewer #4:

Feilen et al. address an interesting mechanistic and biophysical question, namely the significance of the lipid environment for intramembrane enzymatic activity – through biochemical experiments and MD simulations. While the experimental component seems compelling, my opinion is that the structural/computational element is unsuitable for publication in eLife. It is well known that the construction of homology models relies on multiple arbitrary decisions, from the choice of template structure to the sequence alignments to the scoring function – which together imply a degree of inaccuracy that is simply unknown a priori. When a complex is modeled, the uncertainties accumulate. At the same time, MD simulations are, by design, highly sensitive to the details of the input structure. It follows, that systematic inaccuracies in the input model for an MD simulation might result in observations that have no mechanistic significance. At the level of an eLife publication, therefore, it is essential that the authors demonstrate that their conclusions are robust and independent from those built-in uncertainties. A possible way forward would be to carry out simulations of existing experimental structures that might be relevant (with or without minor modifications). Alternatively, or in addition, the authors could consider equally plausible but meaningfully different homology models, constructed on the basis of different assumptions. Either way, I do not believe the manuscript can move forward to publication without an extensive overhaul of the computational section (or its removal), so as to clearly ascertain the conclusions are indeed significant and robust.

A related issue pertains to the comparison of PC bilayers vs DDM micelles. The authors construct and examine one micelle system with a specific number of DDM molecules solubilizing the protein in a specific volume. This choice seems again arbitrary – or at least it is not explained. While the characteristics of lipid bilayer models in odellingn have been extensively studied and optimized (area per lipid, bending modulus, etc), I am unclear the same applies to DDM micelles. What are the observables that give the authors confidence that the structural and elastic properties of their micelle model are realistic? Given that the differences between the POPC and DDM simulations are ultimately modest, this comparative analysis needs to be much more systematic than it currently is to merit publication in eLife. As with the homology odelling, the question is whether different micelles or different DDM models would lead to alternative conclusions.

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

Thank you for resubmitting your work entitled “Active site geometry stabilization of a presenilin homolog by the lipid bilayer promotes intramembrane proteolysis” for further consideration by eLife. Your revised article has been evaluated by 3 reviewers, including Joanne Lemieux as the Reviewing Editor and reviewer #1 and overseen by José Faraldo-Gómez as Senior Editor.

Reviewers and editors agree that the manuscript has been improved but there are some remaining issues that need to be addressed, outlined below. While recognizing the merits of the work, all agree the authors at times overstate the novelty of their findings and that the narrative must be toned down on account of the published work in this area. Furthermore, the reviewers add some insight that could be brought into the discussion.

Reviewer #1:

The authors provide an exciting body of work to support the role of the lipid bilayer in stabilizing the coordination of active site residues. The PSH protein is a suitable model since it does not contain any co-factors and is an active protease on its own. The authors demonstrate that the lipid bilayer enhances processivity of PSH. The data is supported by in vitro cleavage assays in detergent and lipid vesicles. Furthermore, in this revised version there are MD analyses conducted of several apo- and substrate: protease complexes, which provide insight into the protease dynamics and interactions with substrates. The revisions add depth to the paper and is suitable for publication.

Reviewer #2:

Comments and suggestions that could improve the manuscripts are below:

1. In the case of the deletions/proline mutations in the hybrid β-sheet and the lysine mutations in TMD6a, I would have liked to see some type of experiment (DLS, SEC, etc.) to show that the mutations did not simply result in aggregated or misfolded protein. I find it troublesome that it is rare to see confirmation that point mutations to key residues did not completely disrupt protein folding, particularly when detergents are being used for solubilization.

2. It would be worthwhile to consider the possibility that the inhibitors partition into the DDM micelle in the detergent environment. This is a known phenomenon and could explain the poor inhibition of PSH by both L-685,458 and Merck C in the DDM environment, particularly if the DDM that was added to the POPC experiments was to aid in the solubilization of what I gather is a highly insoluble molecule. In general, more detail in the experimental section regarding the in vitro cleavage assay and affinity precipitations would be helpful.

3. The manuscript suggests that it is remarkable to see a rise in processivity when PSH is reconstituted in POPC membranes. In fact, highly inconsistent enzymatic activities (and structures!) have been observed for numerous enzymes in detergent vs. lipid bilayers (e.g. MsbA, MalFGK). It could be beneficial to consider these examples in the discussion and tone down the language, as it is not particularly surprising to see this sort of inconsistency in these different environments.

4. Both the in vitro and in silico experiments in the lipid environment were carried out in POPC, but PSH is an archaeal homologue. The lipids found in archaea are very different to other membranes, and while it may be outside the scope of this study to carry out the experiments in the natural PSH lipid environment (or not possible due to availability), it may be worthwhile running MD simulations a more representative environment, particularly because it well established that the identity of the annular lipids around an enzyme can significantly affect its activity.

5. Related to the above point, I realize that working with a single lipid simplifies things, but could or should the experiments in lipid bilayer not have been done in a brain lipid extract to more accurately recapitulate the environment of the enzyme that PSH is meant to be a surrogate for? Or as above, at the very least the MD simulations could have been done in a more representative environment.

Reviewer #3:

The manuscript by Feilen et al. underwent significant revisions and addressed the reviewer's concerns adequately. The manuscript focuses on the importance of the lipid environment for intramembrane proteolysis. Before I list a number of only smaller suggestions for this manuscript, I would like to mention that while I fully understand that this is beyond the scope of the present manuscript, it would have been interesting if the authors could have enforced the aspect of direct lipid-enzyme interactions, about which very little is known. For example, could the authors deduce amino acids in PSH that are important for the interaction with POPC molecules, and mutate those? Or is the stabilizing effect of POPC solely conveyed through the self-ordering of membrane molecules?

https://doi.org/10.7554/eLife.76090.sa1

Author response

[Editors’ note: the authors resubmitted a revised version of the paper for consideration. What follows is the authors’ response to the first round of review.]

Reviewer #2:

I have several comments and suggestions for the authors to consider improving the quality of the paper:

1. Utilization of C100 for enzymatic analysis.

Authors utilized the recombinant C100-His6 as a substrate for the enzymatic assay. In contrast, they modeled the complex structure of PSH with C83, as the template cryo-EM structure was PS1-C83 complex. Essentially γ-secretase cleaves several substrates regardless the sequence. However, N-terminal length of the substrate might affect the proteolytic efficacy of the γ-secretase. In fact, Funamoto et al. (Nat Commun 2013) reported that AICD production from C83-FLAG was much greater than that from C99-FLAG (i.e., distinct Km and Vmax values). These data suggest the possibility that the formation of E-S complex is regulated by the N-terminal length of the substrate in the proteolytic mechanism of the γ-secretase. Thus, I would recommend the authors to compare the cleavage and trimming patterns of C83 by PSH in either DDM or POPC, as shown in the modeling analysis and the MD simulation. Such comparison strengthens the author's conclusion that the stabilization of E-S complex is critical to the intramembrane proteolysis.

We thank the reviewer for the recommendations on performing cleavage assays with APP C83. In the new version of the manuscript, we now also analyzed the cleavage of C83 by PSH in the DDM and POPC environments to test whether substrate N-terminus length might affect our results. As shown in Figure 3A, we observed robust cleavage of C83 also by PSH in both environments. In line with our observations for C99 processing by PSH in the different environments, POPC also increased the processive cleavage of C83 resulting in the production of shorter p3 species (Figure 3B and Figure 3 – source data 1). The similarity of the C83 (Figure 3B and Figure 3 – source data 1) and C99 processivity data (Figure 3C and Figure 3 – source data 1) further strengthen our conclusion that the stabilization of the E-S complex is indeed critical for intramembrane proteolysis. Furthermore, these data inform that the differences in the processivity of PSH observed for C99 in detergent micelle vs lipid bilayer environments should be explorable at the structural level with models of a PSH–C83 E-S complex.

2. Importance of TMD6a and hybrid β-sheet in the proteolysis by PSH.

Modeled structure of PSH-C83 strongly implicates the importance of these two structural elements in the proteolysis of PSH. However, it remains unclear whether these elements are truly required for the proteolytic activity of PSH. The possibility that these structures artificially modeled because the authors used the γ-secretase-C83 structure as a template. Thus, authors should test the proteolytic activity of mutant PSH that carries amino acid substitutions in TMD6a or beta2-strand to abolish the interaction with the substrate.

We agree with the reviewer that the two important structural elements TMD6a and 2-strand might have arisen from the PS1 template of our homology modeling and therefore might be only structural artefacts.

Therefore, we performed mutational analysis of these structural elements and assed their proteolytic activity. As shown in Figure 4E and Figure 7F, mutations within TMD6a and 2-strand abolished the activity of PSH towards C99. Therefore we believe that both structural elements are functionally present in PSH and do not represent an artefact from model building.

3. Figure presentation

In figure 5A, it is difficult to understand the difference of TMD6a conformation in DDM/POPC because of superimposed structure. They should be presented separately.

We agree with the reviewer. In our new version of the manuscript, we now provide in Figure 6 separate figures for the PSH structures in POPC (Figure 6A) and DDM (Figure 6B-D). This makes it easier to differ between the stable TMD6a conformation in POPC and the unfolded conformation in DDM.

Reviewer #3:

In this manuscript, Feilen et al. combined biochemical experiments and molecular dynamics (MD) simulations to investigate the effects of detergent solubilization versus lipid vesicle reconstitution on intramembrane protease activity of an archaeal presenilin homolog (PSH). This PSH has been previously reported to process amyloid precursor (APP)-based substrate to amyloid β-peptides (A-β) in a manner closely similar to that accomplished by the presenilin-containing γ-secretase complex. The archaeal PSH is employed here as a surrogate to gain mechanistic insight into γ-secretase.

The authors showed that the carboxypeptidase-like activity of PSH is impaired in DDM micelles compared to a POPC lipid bilayer. Evidence is also provided suggesting that DDM-solubilized PSH binds more weakly to a transition-state analog inhibitor compared with POPC-bilayer PSH. Comparative MD simulations suggested that the lipid bilayer stabilized relevant structural elements of PSH for substrate binding (notably TMD6a) and formation of the enzyme active site geometry for proteolysis. A number of suggestions that could help improving the manuscript include the following:

1. While it is understood that archaeal PSH is taken as a surrogate for presenilin in the γ-secretase complex, the authors should explain why this is necessary. All the described experiments, including the MD simulations, could have been conducted with γ-secretase itself. In the discussion, the specific implications for γ-secretase (especially FAD mutations) should be de-emphasized and more emphasis put on the implications for intramembrane proteolysis in general.

As recommended by the reviewer, we have now more clearly described in the introduction why PSH is needed as presenilin/-secretase surrogate for our study. The principal problem is that -secretase, unlike PSH, is not active without lipids (e.g., Zhou et al. 2010), so that a comparison of processivity in lipid bilayer versus detergent-micelle environment is not possible. This type of analysis can therefore not be done for -secretase itself. As further recommended, we also followed the reviewer´s suggestion to de-emphasize the implications for FAD mutations in -secretase and rewrote the respective part in the Discussion, highlighting more the implications of our findings for intramembrane proteolysis in general.

References:

Zhou et al. (2010) Dependency of -secretase complex activity on the structural integrity of the bilayer, Biochem. Biophys. Res. Commun., 12, 402(2), 291-296.

2. For the MS analysis of A-β peptide products in Figure 1D, a table of observed vs. calculated m/z should be provided. Some of these peaks (A-β-43, -45, and -46) are quite weak. The same is true of AICD MS analysis in Figure 2D.

As requested by the reviewer, we now provide the calculated and observed masses for the peaks of our MS analysis. These are found in the corresponding source data files (Figure 1 – source data 1 and Figure 2 – source data 1). We have worked hard to get better spectra for the weaker peaks but were so far not successful and could thus not replace this figure with better quality data. Besides matching masses, the peaks are specific as demonstrated by their absence in inhibitor controls (Figure 1 —figure supplement 1 and Figure 2 —figure supplement 2).

3. Given crystal structures of the PSH, is it more reasonable to preserve protein coordinates in these crystal structures, but only add the missing residues (e.g., TMD6a) using the PS1/γ-secretase cryo-EM structures as template? The Results section is vague about how the homology modeling of enzyme-substrate complex was generated; an additional sentence or two should be provided so the reader does not have to refer to the experimental section to answer this basic question.

We thank the reviewer for his recommendations on model building and also agree that they should be better described in the Results section. In the new version of the manuscript, for which two additional C83-bound PSH models were generated, this is now more clearly described in the main text. The models are generated by three different approaches, including one considering only the PS1 cryo-EM structure (model 1, Table 1) and two including also the coordinates of the PSH crystal structure (model 2 and model 3, Table 1). It is important to note that although the X-ray crystal structure of apo PSH is available (PDB 4HYG), residues 182 to 209 had been proteolytically removed and five stabilizing mutations were introduced. Furthermore, due to its tetrameric nature, the possible substrate entries between TMDs 2 and 6 and TMDs 2 and 3 are blocked implicating an artificial TMD arrangement in the crystal structure. Therefore, it is questionable if the 4HYG structure alone is a good starting model for a PSH-substrate complex. However, we like to emphasize that our model 2 (corresponding to the model used in the first version of the manuscript) is indeed based mostly on the 4HYG structure, but the missing relevant parts for substrate binding are replaced by the PS1 template (Table 1). Our third model considers also parts of the structural information of PS1/-secretase cryo-EM structures ( Table 1), aiming to generate a more accurate model that contains the reliable part of the apo-PSH structure and the critical features revealed in the PS1-C83 complex. As stated above, we also include a more detailed description of the homology modeling in the new version of the manuscript in the Results and Materials and methods sections.

4. Assuming pKa calculations depend on local geometry of the protonation site, the protein structure(s) used for the calculations need to be described clearly. Moreover, how do the residue pKa's depend on the protein structures (e..g, 4Y6K and 4HYG crystal structures of PSH, the apo and holo cryo-EM structures of PS1/γ-secretase, and simulation equilibrated protein structures in notably two different conformations with the D162-D220 distance centered around ~6.5 and ~8.2 Angstroms in Figure 6A).

We agree and have accordingly revised the relevant sections in our new manuscript providing additional information regarding the pKa values of the catalytic residues. The pKa prediction on the two aspartic acids of PS1/PSH are now listed in Table 2. For basically all available structures (except for PDB 5FN5, a structure of one of the PS1/-secretase apo-states (apo-state 3)), D220 of PSH (D385 in PS1) in TMD7 is predicted to have a higher pKa than the second aspartic acid residue in TMD6, D162 (D257 in PS1).

5. It is unclear why the Amber force field was used in simulations of holo PSH with two different protonation states, but CHARMM36m in simulations of the apo and holo PSH in different membrane environments. It would help to evaluate the force field differences and potential effects by adding simulations using CHARMM36m to the holo PSH with different protonation states or simulations using Amber to the apo and holo PSH in different membrane environments.

In our previous computational work, the AMBER force field is usually implemented since it is more compatible with our computational engine. However, the AMBER force field does not have the parameters describing the DDM molecules, and therefore we needed to switch the force field from AMBER to CHARMM36m. Since the parameter files with the CHARMM36m force field generated by CHARMM-GUI server cannot be modified easily, changing the protonation state and calculating the substrate-binding energy using the MMPBSA method might include more uncertainties such as the need to define the implicit radii of each atom etc. The verification of the protonation states of the catalytic residues using the AMBER force field is now removed in the new manuscript and is replaced by showing the list of pKa predictions in Table 2.

6. In Figure 3A, what is the RMSD between the PSH homology model and its crystal structures?

We give now more details of the RMSD values in Figure 4 —figure supplement 1. We show the RMSD values between our three models and the two published PSH crystal structures (PDB 4HYG, 4Y6K, chain B) as well as the RMSD between the homology models that we generated.

7. In Figures 4C-4E, it would help to add error bars (standard deviations) to examine what differences are significant. The authors calculated RMSFs of PSH in the apo and holo forms to describe its stability in different lipid environments. It is better to show error bars (with total simulation times across different replicates) to describe RMSF differences in notably, the TMD6a, TMD4 and TM2-TM3 loop. There is a notable RMSF difference just beyond TMD6a, which should be mentioned and explained. In addition to RMSF, further simulation analysis such as comparison of the distances between the catalytic aspartate and scissile peptide bond in C83 could provide more insights.

As recommended by the reviewer, standard deviations of the RMSFs across three different replicates are now included in the plots as shaded areas (Figure 5C and E, Figure 5 —figure supplement 1B, Figure 5 —figure supplement 3B and C, Figure 6 —figure supplement 3D, F and G, Figure 7B and E). The RMSF difference between bilayer and micelle environments at the region immediately C-terminal of TMD6a (T182A184) are now mentioned and discussed in new version. Similar as what was observed for TMD6a, these residues are also interacting with the DDM molecules inserted between TMD3 and TMD4.

In Figure 8 —figure supplement 2, we provide a more detailed geometric picture of the active site of the two sampled conformations shown in Figure 8B including also distances of the catalytic residues to the substrate scissile bond. Although both conformations are capable of forming hydrogen bonds with the protonated aspartic acid (D220), the geometry with a larger C-C distance between the catalytic aspartic acids is less likely for the enzyme to execute the hydrolysis reaction because of the larger distance between the scissile bond and the unprotonated aspartic acid (D162). Overall, our additional analysis of various relevant distances shown in Figure 8 —figure supplement 2 indicates that the CC distance should be the preferred criterion to characterize a functional active site geometry as compared to merely using H-bond distances.

8. In Figure 5A, it would help to quantitatively calculate and plot the helicity of TMD6a and/or secondary structures of residues in TMD6a as a function of simulation time and compare these quantities between the different simulated systems.

As recommended by the reviewer, evaluation of the secondary structure of TMD6a and surrounding residues is now shown in Figure 5 —figure supplement 2 and Figure 7 —figure supplement 1.

9. In Figure 5C-5D, the authors described that the DDM molecule can insert itself between TMD2 and TMD6 and intervenes intra- and intermolecular interactions and thus destabilize TMD6a. It would be better to have more explanation if these insertions are with just one lipid molecule or there are multiple molecules involved during different time frames of the simulations. In addition, it would be more convincing to see atomic detailed interaction between the DDM molecule and the protein, especially because DDM is a nonionic detergent. Moreover, it could help to calculate -SCD order parameters that are usually obtained from NMR experiments to measure orientational anisotropy of the C-H bonds lipid chains and quantify differences of the lipid orientations.

As recommended by the reviewer, the current manuscript now describes in more detail how DDM perturbs intramolecular PSH TMD interactions. Figure 6B and C shows that the inserted DDM molecule disturbs the TMD6a helix mainly by unspecific hydrogen bond interactions with the PSH backbone. This also happens on the residues immediately Cterminal to TMD6a (Figure 6D). During our simulation time, mostly only one DDM molecule was found to insert into the gap and interfere with the PSH backbone. As further recommended, we calculated lipid chain C-H bond order parameters to quantify differences of the lipid orientation. The SCH order parameters calculated are now included in Figure 6 —figure supplement 2B and demonstrate the average orientation of DDM and POPC molecules.

10. For Figure 6A, it would help to plot the D162-D220 distance vs. time in simulations of the different systems also. What are the corresponding distance values in the PSH crystal and PS1 cryo-EM structures? In addition to a main peak at ~6.5 Å distance, there seems another peak at ~8.2 Å distance between D162-D220; what is this conformational state?

As recommended by the reviewer, the D162-D220 C-C distance over simulation time is now shown in Figure 8 —figure supplement 1C. The C-C distances of the substrate-bound structures are not measurable since D385 is mutated to alanine. The C-C distances in other PDB structures of PSH and PS1/-secretase are 3.99 Å (4HYG), 3.99 Å (4UIS), 6.58 Å (5A63), 5.35 Å (5FN5), 11.48 Å (5FN4), 5.06 Å (5FN3), and 3.89 Å (5FN2). The D162-D220 distance at ~8.2 Å corresponds to the geometry where more than two water molecules are accommodated in the catalytic center.

11. In Figure 6B, there is substantial non-specific binding of the biotinylated Merck C to PSH in DDM; the parent inhibitor is essentially not competing, even at 20 microM. This makes the interpretation that specific binding is stronger in POPC vesicles less convincing. Some comment to this effect should be added to the Results section.

This effect puzzled us as well. We agree with the suggestion of the reviewer and have thus commented on this effect regarding its interpretation in the results. As seen in the calculations of binding presented in Figure 8 – source data 1, binding of PSH is, however, still inhibited by the parental compound L-685,458 in DDM, although only to a very minor extent and much less compared to POPC. There was generally also much stronger background binding to the streptavidin beads in the absence of Merck C in DDM, adding to the low amount of specifically captured PSH in this condition. We thus speculate that the more labile active site in DDM also weakens the competition of binding. The reduced levels of specific binding in DDM support the interpretation that the active site geometry is more stabilized in the POPC bilayer. In support of this view, additional enzyme inhibition experiments presented in Figure 8E and F showed that L-685,458 and Merck C both inhibit PSH less well in DDM than in POPC.

Reviewer #4:

Feilen et al. address an interesting mechanistic and biophysical question, namely the significance of the lipid environment for intramembrane enzymatic activity – through biochemical experiments and MD simulations. While the experimental component seems compelling, my opinion is that the structural/computational element is unsuitable for publication in eLife. It is well known that the construction of homology models relies on multiple arbitrary decisions, from the choice of template structure to the sequence alignments to the scoring function – which together imply a degree of inaccuracy that is simply unknown a priori. When a complex is modeled, the uncertainties accumulate. At the same time, MD simulations are, by design, highly sensitive to the details of the input structure. It follows, that systematic inaccuracies in the input model for an MD simulation might result in observations that have no mechanistic significance. At the level of an eLife publication, therefore, it is essential that the authors demonstrate that their conclusions are robust and independent from those built-in uncertainties. A possible way forward would be to carry out simulations of existing experimental structures that might be relevant (with or without minor modifications). Alternatively, or in addition, the authors could consider equally plausible but meaningfully different homology models, constructed on the basis of different assumptions. Either way, I do not believe the manuscript can move forward to publication without an extensive overhaul of the computational section (or its removal), so as to clearly ascertain the conclusions are indeed significant and robust.

The reviewer´s concerns on using comparative modeling are well taken – we are aware of the mentioned limitations and the uncertainties that can arise. As recommended by the reviewer, we have thus generated two additional model structures (Figure 4 —figure supplement 1) in the new version based on the same sequence alignment method but different in the ways of including different structural information from the existing structures (Figure 4 —figure supplement 2 and Table 1). It is important to note that although the X-ray crystal structure of apo PSH is available (PDB 4HYG), residues 182 to 209 are missing and five mutations are introduced. Furthermore, the tetrameric complex of four PSH molecules leads to a crystallographic artefact in the unit cell, which might block possible substrate entry sites between TMDs 2 and 6 and TMDs 2 and 3. Therefore, the available crystal structure is presumably far away from the substrate-bound structure. Our homology models consider also the structural information of PS1/-secretase cryo-EM structures in complex with the C83 substrate (Table 1), aiming to generate a more accurate model that contains the reliable part of apo PSH structure and the critical features revealed in the PS1-C83 complex.

We then performed comparative simulations on all model structures and based on the results argue that model 2 (which contains the reliable part of apo PSH and the substrate-binding region modelled based on the PS1 substrate-bound structure) is most realistic. Since our MD simulation studies on the model structures gave consistent results, we are confident that give likely and reliable explanations for the experimental results.

In addition, simulations can be used to create a hypothesis or to make predictions on the role of residues or protein segments. In this regard, our MD simulations are quite successful: based on the computational data, biochemical experiments were performed to evaluate the importance of the TMD6a and 2 structural features that were identified for holo PSH. By mutational analysis, we could indeed show that TMD6a and 2-strand are critical for substrate cleavage (Figure 4E and Figure 7F), strongly suggesting the reliability of homology model 2, which was selected as our working model. Therefore, we believe that our homology-modeling approach is suitable to study the dynamics of PSH embedded in two different environments at an atomic scale.

Taken together, we thus believe that the combination of our biochemical and computational approaches provides us with a picture of not only the unprecedented substrate-bound PSH structure but also with a molecular explanation of how a micelle environment destabilizes the active site geometry and interrupts the stability of the E-S complex.

A related issue pertains to the comparison of PC bilayers vs DDM micelles. The authors construct and examine one micelle system with a specific number of DDM molecules solubilizing the protein in a specific volume. This choice seems again arbitrary – or at least it is not explained. While the characteristics of lipid bilayer models in odellingn have been extensively studied and optimized (area per lipid, bending modulus, etc), I am unclear the same applies to DDM micelles. What are the observables that give the authors confidence that the structural and elastic properties of their micelle model are realistic? Given that the differences between the POPC and DDM simulations are ultimately modest, this comparative analysis needs to be much more systematic than it currently is to merit publication in eLife. As with the homology odelling, the question is whether different micelles or different DDM models would lead to alternative conclusions.

We agree with the reviewer that the chosen DDM micelle configuration could potentially be another concern and have thus followed the advice to control our data for such micelle-mediated effects. As recommended by the reviewer, our new manuscript now also includes an additional simulation system with 50% more DDM molecules (Figure 6 —figure supplement 3) The number of DDM molecules and micelle size were chosen based on Cheng et al. 2013, who parameterized the micelle parameters. By comparing the size of PSH to the size of proteins OmpA (PDB 1BXW) and OmpF (PDB 3POX) simulated by Cheng et al., our substrate-bound PSH is larger than OmpA but smaller than the OmpF trimer. With 80 DPC (n-dodecylphosphocholine) molecules in the former system and 160 in the latter, we consider 150 DDM molecules being a suitable number to investigate its dynamics in our system. In addition, the dynamics of the newly constructed holo PSH in the larger DDM micelle is consistent with the initial micelle size, which further consolidates the impact of DDM shown in our study. In the new manuscript version, we also include the area per lipid of POPC molecules (Figure 6 —figure supplement 2C) as well as SCH order parameters of POPC and DDM (Figure 6 —figure supplement 2B). Overall, the values of area per lipid and the SCH value for POPC agree well with the experimental data. Although the SCH order parameter for DDM is to our knowledge not available to date and was therefore computed from our simulations, our analysis shows how disordered the DDM molecules are and from our structural depiction, we could illustrate how the disordered character of DDM perturbs the intramolecular interactions of PSH, in particular with those of TMD6a.

In conclusion, while we acknowledge the possibility that different sizes of micelle or lipid bilayer models may give divergent results from our initially chosen simulation systems, our new data with a larger DDM micelle size now suggest that this is rather unlikely as our principal observations were recapitulated with the additional analyses of alternative detergent micelle configurations. These additional data are thus overall in line with what one would expect for these two different conditions as the basic geometrical properties of a micelle, a spherical amphiphilic environment with disordered detergent molecules, and those of a bilayer, a relatively flat environment with ordered phospholipids molecules, are principally also shared by other detergent or phospholipid molecules.

References:

Cheng et al. (2013) CHARMM-GUI micelle builder for pure/mixed micelle and protein/micelle complex systems, J. Chem. Inf. Model, 53(8), 2171-

2180.

[Editors’ note: what follows is the authors’ response to the second round of review.]

Reviewers and editors agree that the manuscript has been improved but there are some remaining issues that need to be addressed, outlined below. While recognizing the merits of the work, all agree the authors at times overstate the novelty of their findings and that the narrative must be toned down on account of the published work in this area. Furthermore, the reviewers add some insight that could be brought into the discussion.

Reviewer #2:

Comments and suggestions that could improve the manuscripts are below:

1. In the case of the deletions/proline mutations in the hybrid β-sheet and the lysine mutations in TMD6a, I would have liked to see some type of experiment (DLS, SEC, etc.) to show that the mutations did not simply result in aggregated or misfolded protein. I find it troublesome that it is rare to see confirmation that point mutations to key residues did not completely disrupt protein folding, particularly when detergents are being used for solubilization.

We agree with the reviewer that mutations might potentially induce aggregation and/or misfolding of the protein. To rule this out, we performed dynamic light scattering (DLS), Blue Native PAGE (BN-PAGE) and nano differential scanning fluorometry (nanoDSF) experiments with WT and mutant PSH. As shown in Appendix 1 – figure 1 A (DLS), B (BN-PAGE) and C (nanoDSF), all studied mutants behaved similar as WT PSH in these analyses suggesting that the mutants do not disrupt protein folding.

2. It would be worthwhile to consider the possibility that the inhibitors partition into the DDM micelle in the detergent environment. This is a known phenomenon and could explain the poor inhibition of PSH by both L-685,458 and Merck C in the DDM environment, particularly if the DDM that was added to the POPC experiments was to aid in the solubilization of what I gather is a highly insoluble molecule. In general, more detail in the experimental section regarding the in vitro cleavage assay and affinity precipitations would be helpful.

We thank the reviewer for pointing out this possibility, which we took into consideration.

However, in the POPC environment, DDM was not added to aid the solubilization of the L-685,458 γ-secretase inhibitor (which is added to the assays from a stock solution in DMSO), but to ensure enzyme activity of PSH. As we show in Author response image 1, PSH is not active in the POPC bilayer without the addition of small amounts of DDM (below the CMC). Interestingly, this seems to be a general phenomenon for presenilin-type intramembrane proteases and was already reported for γ-secretase in our Winkler et al. 2012 study. As shown in Supplementary figure 2B of this paper, also γ-secretase requires the addition of small amounts of CHAPSO detergent for activity when it is reconstituted in POPC vesicles. The reason for this requirement is however not understood yet. We included this information in the Materials and methods in our revised version and as suggested by the reviewer have also slightly revised the relevant text passages of our cleavage and inhibitor binding assays to provide the reader with more detail.

Author response image 1
Comparison of PSH activity in the POPC bilayer with and without the addition of DDM.

Analysis of PSH activity in POPC vesicles with and without the addition of 0.008% DDM after incubation with C100-His6 at 37 °C overnight by immunoblotting for AICD (Y188) and Aβ (2D8). The asterisks mark substrate degradation bands, which are independent of PSH cleavage.

Furthermore, while we cannot completely rule out the possibility of some inhibitor partition into DDM micelles, we do think that the inhibitors reach the active site in sufficient amounts due to their principal chemical properties. The mimicking of the transition state of a substrate peptide bond attacked by a presenilin-type aspartyl protease should direct these inhibitors to the active site. In fact, structural investigations of PSH with the γ-secretase inhibitor L-682,679, which is from the same chemical class as L-685,458 and Merck C, showed that the inhibitor enters the catalytical cleft (Dang et al. 2015). Also, for γ-secretase it was observed that L-685,458 binds directly at the active site (Yang et al. 2021). In addition, in MD simulations of γ-secretase with L-685,458 the inhibitor stays tightly bound to the active site (Hitzenberger et al. 2019).

References:

Winkler et al. (2012) Generation of Alzheimer disease-associated amyloid β42/43 peptide by γ-secretase can be inhibited directly by modulation of membrane thickness, J. Biol. Chem., 287(25), 21326-21334.

Dang et al. (2015) Cleavage of amyloid precursor protein by an archaeal presenilin homologue PSH, Proc. Natl. Acad. Sci. USA, 112(11), 3344-3349.

Yang et al. (2021) Structural basis of γ-secretase inhibition and modulation by small molecule drugs, Cell, 184(2), 521-533.e14.

Hitzenberger et al. (2019) Uncovering the binding mode of γ-secretase inhibitors, ACS Chem. Neurosci., 10(8), 3398-3403.

3. The manuscript suggests that it is remarkable to see a rise in processivity when PSH is reconstituted in POPC membranes. In fact, highly inconsistent enzymatic activities (and structures!) have been observed for numerous enzymes in detergent vs. lipid bilayers (e.g. MsbA, MalFGK). It could be beneficial to consider these examples in the discussion and tone down the language, as it is not particularly surprising to see this sort of inconsistency in these different environments.

We thank the reviewer for drawing our attention to studies on other membrane proteins that, for various reasons, show distinct changes in their activities in detergent micelles vs. membrane bilayers. However, in our study, we could consistently link activity differences in these two environments to structural changes of the enzyme-substrate complex and the active site geometry of the protease as determined by MD simulations that were in agreement with the biochemical analyses which included mutational analysis as well as the use of inhibitors probing the conformation of the active site. Nevertheless, since our findings may be misconceived in the light of studies with other membrane proteins in different environments, we followed the advice of the reviewer and toned down our wording by changing two sentences of the manuscript in the summarizing paragraph of the Discussion. This section is now modified as follows:

“Taken together, in good correlation between experimental and simulation data, our results with PSH as a model intramembrane protease highlight an important role of the membrane lipid environment in providing a stabilized E-S conformation that is crucial for substrate processing in intramembrane proteolysis.”

This rewording should capture the novelty of the insights provided by our study without overstatements.

4. Both the in vitro and in silico experiments in the lipid environment were carried out in POPC, but PSH is an archaeal homologue. The lipids found in archaea are very different to other membranes, and while it may be outside the scope of this study to carry out the experiments in the natural PSH lipid environment (or not possible due to availability), it may be worthwhile running MD simulations a more representative environment, particularly because it well established that the identity of the annular lipids around an enzyme can significantly affect its activity.

As already noted by the reviewer, archaeal lipids are not available for biochemical experiments, and we can therefore not perform experiments with the natural lipids for this protease. Additionally, to our knowledge nothing is known about the lipid composition of Methanoculleus marisnigri, the archaeal origin of PSH. Therefore, it will be difficult to perform biochemical experiments in a natural PSH environment and given this limitation, accompanying MD simulations will have only little if any informative value. Additionally, proper force fields for archaeal lipids are only limitedly available. Furthermore, the study represented here uses a non-natural substrate of PSH, since so far, no natural substrates are described. Even though we agree with the reviewer that the identity of the annular lipids around an enzyme can significantly affect its activity, we think that investigating PSH in a natural lipid environment would in this case also best be done with a natural substrate in order to study a homologous system. But we agree that while such a study is currently not doable and also outside the scope of this present study, it could be an interesting subject of future studies.

5. Related to the above point, I realize that working with a single lipid simplifies things, but could or should the experiments in lipid bilayer not have been done in a brain lipid extract to more accurately recapitulate the environment of the enzyme that PSH is meant to be a surrogate for? Or as above, at the very least the MD simulations could have been done in a more representative environment.

Simulating PSH in a bilayer composed of brain lipids is challenging due to the complex composition of brain lipids. From the available brain lipid extracts – also used in our experiments – only 40% of the lipids are known. This makes an accurate modeling of a brain lipid environment impossible. Furthermore, in our timeframe of 600 ns it will not be possible to sample all possibilities in a complex lipid system and influences of low abundance lipids might not be recognized at all. Nevertheless, we performed the cleavage assay with PSH reconstituted in vesicles prepared from a brain lipid extract and observed that this complex lipid environment also facilitates the processive cleavage of PSH (Author response image 2A–C) indicating that this observation is not limited to the lipid type used in our studies.

Author response image 2

Cleavage of C99 by PSH reconstituted in brain lipid vesicles. (A) Analysis of PSH activity in DDM micelles, POPC vesicles and brain lipid vesicles after incubation with C100-His6 substrate at 37 °C overnight by immunoblotting for AICD (Y188) and Aβ (2D8). The asterisks mark substrate degradation bands, which are independent of PSH cleavage. (B) Separation of Aβ species produced by PSH in DDM micelles, POPC vesicles and brain lipid vesicles by Tris-Bicine urea SDS-PAGE and analysis by immunoblotting for Aβ (2D8). (C) MALDI-TOF MS analysis of Aβ species generated by PSH in brain lipid vesicles at pH 7.0. The intensity of the highest peak was set to 100%.

Reviewer #3:

The manuscript by Feilen et al. underwent significant revisions and addressed the reviewer's concerns adequately. The manuscript focuses on the importance of the lipid environment for intramembrane proteolysis. Before I list a number of only smaller suggestions for this manuscript, I would like to mention that while I fully understand that this is beyond the scope of the present manuscript, it would have been interesting if the authors could have enforced the aspect of direct lipid-enzyme interactions, about which very little is known. For example, could the authors deduce amino acids in PSH that are important for the interaction with POPC molecules, and mutate those? Or is the stabilizing effect of POPC solely conveyed through the self-ordering of membrane molecules?

We agree with the reviewer that it would be interesting to understand whether the effects we observed in the bilayer system are caused by a direct lipid-protein interaction or solely conveyed through the self-ordering of membrane molecules. Therefore, we analyzed the residue-wise solvent residence time of POPC and DDM using PyLipID (Song et al. 2022) with a dual cutoff distance of 1.2 Å and 0.8 Å (Author response image 3A). Interestingly, this analysis shows that both POPC and DDM molecules reside longest in the gap between the Nterminal part of the TMD of C83 and the N-terminal part of TMD3 of PSH (Author response image 3B, upper panel). Furthermore, POPC molecules show additional binding to several other regions (e.g. the loop between TMD7 and TMD8) (Author response image 3B, lower panel). The long residence time between C83 and TMD3 of PSH is due to the large gap into which POPC and DDM can insert, whereas this entry is blocked by the hydrophobic interaction between I31 of C83 and I242 of nicastrin in the C83-bound γ-secretase complex (Author response image 3C). Strikingly, the POPC molecule that we identified residing at the surface between TMD1, TMD5, TMD7 and TMD8 binds in a conformation similar to the ones resolved in the C83-bound and Notch1-bound γ-secretase complex structures (Author response image 3D).

Author response image 3
Binding of POPC and DDM molecules on the PSH surface.

(A) The residue-wise averaged solvent residence time of POPC (blue) or DDM (red) on PSH surface. (B) Visualization of the DDM (middle) and POPC (right) residence time based on (A) from two different angles. Surfaces with long, medium and short lipid residence time are colored in blue, white, and red, respectively. Residence time is visualized in surface representation. (C) Binding of DDM (left) and POPC (middle) molecules in the gap between C83 and PSH TMD3 in the PSH simulations. In γ-secretase, this gap is blocked by the hydrophobic contact between nicastrin I242 (green) and C83 I31 (orange, right, 6IYC). (D)Binding of a POPC molecule to the PSH–C83 complex from MD simulations (left), Notchbound γ-secretase (middle, 6IDF) and C83-bound γ-secretase (right, 6IYC).

The similar lipid interaction sites in PSH and PS1 are very interesting and might strengthen the hypothesis that direct lipid-enzyme interactions might influence the enzyme activity. This is an interesting starting point for future studies on the lipid-enzyme interaction but is in our opinion out of scope for our present study. Of note, previous structural investigations in different intramembrane proteases have identified lipid molecules that were bound to the respective intramembrane protease. But none of these studies revealed a functional role of the enzymelipid interaction (Ben-Shem et al. 2007, Lemieux et al. 2007, Bai et al. 2015, Yang et al. 2019, Zhou et al. 2019).

References:

Song et al. (2022) PyLipID: A python package for analysis of protein-lipid interactions from molecular dynamics simulations. J. Chem. Theory Comput., 18(2), 1188-1201.

Ben-Shem et al. (2007) Structural basis for intramembrane proteolysis by rhomboid serine proteases, Proc. Natl. Acad. Sci. USA. 104(2), 462-466.

Lemieux et al. (2007) The crystal structure of the rhomboid peptidase from Haemophilus influenzae provides insight into intramembrane proteolysis, Proc. Natl. Acad. Sci. USA, 104(3), 750-754.

Bai et al. (2015) An atomic structure of human γ-secretase, Nature, 525(7568), 212-217.

Yang et al. (2019) Structural basis of Notch recognition by human γ-secretase, Nature, 565(7738), 192-197. Zhou et al., (2019) Recognition of the amyloid precursor protein by human γ-secretase, Science, 363(6428), eaaw0930.

https://doi.org/10.7554/eLife.76090.sa2

Article and author information

Author details

  1. Lukas P Feilen

    German Center for Neurodegenerative Diseases, Munich, Germany
    Contribution
    Conceptualization, Formal analysis, Investigation, Project administration, Visualization, Writing - review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8221-6742
  2. Shu-Yu Chen

    Center of Functional Protein Assemblies and Physics Department T38, Technical University of Munich, Garching, Germany
    Contribution
    Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft
    Competing interests
    No competing interests declared
  3. Akio Fukumori

    Department of Pharmacotherapeutics II, Osaka Medical and Pharmaceutical University, Takatsuki, Japan
    Contribution
    Construct generation, Resources, Writing - review and editing
    Competing interests
    No competing interests declared
  4. Regina Feederle

    1. German Center for Neurodegenerative Diseases, Munich, Germany
    2. Institute for Diabetes and Obesity, Monoclonal Antibody Core Facility, Helmholtz Munich, German Research Center for Environmental Health, Neuherberg, Germany
    Contribution
    Antibody generation, Resources
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3981-367X
  5. Martin Zacharias

    Center of Functional Protein Assemblies and Physics Department T38, Technical University of Munich, Garching, Germany
    Contribution
    Conceptualization, Funding acquisition, Resources, Supervision, Writing - review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5163-2663
  6. Harald Steiner

    1. German Center for Neurodegenerative Diseases, Munich, Germany
    2. Biomedical Center (BMC), Division of Metabolic Biochemistry, Faculty of Medicine, Munich, Germany
    Contribution
    Conceptualization, Funding acquisition, Project administration, Supervision, Writing – original draft
    For correspondence
    harald.steiner@med.uni-muenchen.de
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3935-0318

Funding

Deutsche Forschungsgemeinschaft (263531414 / FOR2290)

  • Martin Zacharias
  • Harald Steiner

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

Acknowledgements

We thank Yigong Shi for the PSH expression plasmid, Karlheinz Baumann for GSIs and Gabriele Basset, Frits Kamp and Alice Sülzen for technical assistance as well as Shibojyoti Lahiri, Ignasi Forné and Axel Imhof from the Protein Analysis Unit of the Biomedical Center Munich (BMC) and Michaela Smolle from the Biophysics Core Facility of the BMC for access to their instruments, helpful discussions and advice.

Senior Editor

  1. José D Faraldo-Gómez, National Heart, Lung and Blood Institute, National Institutes of Health, United States

Reviewing Editor

  1. M Joanne Lemieux, University of Alberta, Canada

Reviewers

  1. M Joanne Lemieux, University of Alberta, Canada
  2. Sean Workman, University of Regina, Canada

Version history

  1. Received: December 3, 2021
  2. Preprint posted: January 10, 2022 (view preprint)
  3. Accepted: May 16, 2022
  4. Accepted Manuscript published: May 17, 2022 (version 1)
  5. Version of Record published: July 14, 2022 (version 2)

Copyright

© 2022, Feilen 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

  • 849
    Page views
  • 198
    Downloads
  • 1
    Citations

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

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. Lukas P Feilen
  2. Shu-Yu Chen
  3. Akio Fukumori
  4. Regina Feederle
  5. Martin Zacharias
  6. Harald Steiner
(2022)
Active site geometry stabilization of a presenilin homolog by the lipid bilayer promotes intramembrane proteolysis
eLife 11:e76090.
https://doi.org/10.7554/eLife.76090

Further reading

    1. Biochemistry and Chemical Biology
    2. Cell Biology
    Maria Körner, Susanne R Meyer ... Alexander Buchberger
    Research Article Updated

    The ATPase p97 (also known as VCP, Cdc48) has crucial functions in a variety of important cellular processes such as protein quality control, organellar homeostasis, and DNA damage repair, and its de-regulation is linked to neuromuscular diseases and cancer. p97 is tightly controlled by numerous regulatory cofactors, but the full range and function of the p97–cofactor network is unknown. Here, we identify the hitherto uncharacterized FAM104 proteins as a conserved family of p97 interactors. The two human family members VCP nuclear cofactor family member 1 and 2 (VCF1/2) bind p97 directly via a novel, alpha-helical motif and associate with p97-UFD1-NPL4 and p97-UBXN2B complexes in cells. VCF1/2 localize to the nucleus and promote the nuclear import of p97. Loss of VCF1/2 results in reduced nuclear p97 levels, slow growth, and hypersensitivity to chemical inhibition of p97 in the absence and presence of DNA damage, suggesting that FAM104 proteins are critical regulators of nuclear p97 functions.

    1. Biochemistry and Chemical Biology
    2. Epidemiology and Global Health
    Takashi Sasaki, Yoshinori Nishimoto ... Yasumichi Arai
    Research Article

    Background: High levels of circulating adiponectin are associated with increased insulin sensitivity, low prevalence of diabetes, and low body mass index (BMI); however, high levels of circulating adiponectin are also associated with increased mortality in the 60-70 age group. In this study, we aimed to clarify factors associated with circulating high-molecular-weight (cHMW) adiponectin levels and their association with mortality in the very old (85-89 years old) and centenarians.

    Methods: The study included 812 (women: 84.4%) for centenarians and 1,498 (women: 51.7%) for the very old. The genomic DNA sequence data were obtained by whole genome sequencing or DNA microarray-imputation methods. LASSO and multivariate regression analyses were used to evaluate cHMW adiponectin characteristics and associated factors. All-cause mortality was analyzed in three quantile groups of cHMW adiponectin levels using Cox regression.

    Results: The cHMW adiponectin levels were increased significantly beyond 100 years of age, were negatively associated with diabetes prevalence, and were associated with SNVs in CDH13 (p = 2.21 × 10-22) and ADIPOQ (p = 5.72 × 10-7). Multivariate regression analysis revealed that genetic variants, BMI, and high-density lipoprotein cholesterol (HDLC) were the main factors associated with cHMW adiponectin levels in the very old, whereas the BMI showed no association in centenarians. The hazard ratios for all-cause mortality in the intermediate and high cHMW adiponectin groups in very old men were significantly higher rather than those for all-cause mortality in the low level cHMW adiponectin group, even after adjustment with BMI. In contrast, the hazard ratios for all-cause mortality were significantly higher for high cHMW adiponectin groups in very old women, but were not significant after adjustment with BMI.

    Conclusions: cHMW adiponectin levels increased with age until centenarians, and the contribution of known major factors associated with cHMW adiponectin levels, including BMI and HDLC, varies with age, suggesting that its physiological significance also varies with age in the oldest old.

    Funding: This study was supported by grants from the Ministry of Health, Welfare, and Labour for the Scientific Research Projects for Longevity; a Grant-in-Aid for Scientific Research (No 21590775, 24590898, 15KT0009, 18H03055, 20K20409, 20K07792, 23H03337) from the Japan Society for the Promotion of Science; Keio University Global Research Institute (KGRI), Kanagawa Institute of Industrial Science and Technology (KISTEC), Japan Science and Technology Agency (JST) Research Complex Program 'Tonomachi Research Complex' Wellbeing Research Campus: Creating new values through technological and social innovation (JP15667051), the Program for an Integrated Database of Clinical and Genomic Information from the Japan Agency for Medical Research and Development (No. 16kk0205009h001, 17jm0210051h0001, 19dk0207045h0001); the medical-welfare-food-agriculture collaborative consortium project from the Japan Ministry of Agriculture, Forestry, and Fisheries; and the Biobank Japan Program from the Ministry of Education, Culture, Sports, and Technology.