γ-Secretase Studied by Atomistic Molecular Dynamics Simulations: Global Dynamics, Enzyme Activation, Water Distribution and Lipid Binding

Emerging structures and dynamic mechanisms of γ-secretase for Alzheimer's disease

γ-Secretase, called "the proteasome of the membrane," is a membrane-embedded protease complex that cleaves 150+ peptide substrates with central roles in biology and medicine, including amyloid precursor protein and the Notch family of cell-surface receptors. Mutations in γ-secretase and amyloid precursor protein lead to early-onset familial Alzheimer's disease. γ-Secretase has thus served as a critical drug target for treating familial Alzheimer's disease and the more common late-onset Alzheimer's disease as well. However, critical gaps remain in understanding the mechanisms of processive proteolysis of substrates, the effects of familial Alzheimer's disease mutations, and allosteric modulation of substrate cleavage by γ-secretase. In this review, we focus on recent studies of structural dynamic mechanisms of γ-secretase. Different mechanisms, including the "Fit-Stay-Trim," "Sliding-Unwinding," and "Tilting-Unwinding," have been proposed for substrate proteolysis of amyloid precursor protein by γ-secretase based on all-atom molecular dynamics simulations. While an incorrect registry of the Notch1 substrate was identified in the cryo-electron microscopy structure of Notch1-bound γ-secretase, molecular dynamics simulations on a resolved model of Notch1-bound γ-secretase that was reconstructed using the amyloid precursor protein-bound γ-secretase as a template successfully captured γ-secretase activation for proper cleavages of both wildtype and mutant Notch, being consistent with biochemical experimental findings. The approach could be potentially applied to decipher the processing mechanisms of various substrates by γ-secretase. In addition, controversy over the effects of familial Alzheimer's disease mutations, particularly the issue of whether they stabilize or destabilize γ-secretase-substrate complexes, is discussed. Finally, an outlook is provided for future studies of γ-secretase, including pathways of substrate binding and product release, effects of modulators on familial Alzheimer's disease mutations of the γ-secretase-substrate complexes. Comprehensive understanding of the functional mechanisms of γ-secretase will greatly facilitate the rational design of effective drug molecules for treating familial Alzheimer's disease and perhaps Alzheimer's disease in general.

Selective Protonation of Catalytic Dyad for γ-Secretase-Mediated Hydrolysis Revealed by Multiscale Simulations

γ-Secretase plays a crucial role in producing disease-related amyloid-β proteins by cleaving the amyloid precursor protein (APP). The enzyme employs its catalytic dyad containing two aspartates (Asp257 and Asp385) to hydrolyze the substrate by a general acid-base catalytic mechanism, necessitating monoprotonation of the two aspartates for efficient hydrolysis. However, the precise protonation states of the aspartates remain uncertain. In this study, we employed a multiscale computational approach to investigate the dependence of the catalytic efficiency of γ-secretase on the protonation states of its catalytic dyad. Over 200 ms unbiased atomistic simulations of the substrate-enzyme complex reveal diverse orientations of the scissile bond of the bound substrate and accessible structural ensembles of the catalytic dyad with Asp257-Asp385 distances fluctuating between 4 and 10 Å. With a quantum mechanics/molecular mechanics (QM/MM) approach accelerated by enhanced sampling techniques, we find that the first step of the hydrolysis reaction, i.e., the formation of a gem-diol intermediate, experiences a higher reaction barrier by ∼2 kcal/mol when Asp385 is protonated. Furthermore, we find that Arg269 of the enzyme is most likely responsible for this preference of the protonation state: its basic side chain is spatially close to that of Asp257 and specifically stabilizes the transition state electrostatically when Asp257 is protonated. Collectively, our study suggests that Asp257 is likely the favored protonation site for APP cleavage by γ-secretase.

Enzyme-substrate hybrid β-sheet controls geometry and water access to the γ-secretase active site

γ-Secretase is an aspartyl intramembrane protease that cleaves the amyloid precursor protein (APP) involved in Alzheimer's disease pathology and other transmembrane proteins. Substrate-bound structures reveal a stable hybrid β-sheet immediately following the substrate scissile bond consisting of β1 and β2 from the enzyme and β3 from the substrate. Molecular dynamics simulations and enhanced sampling simulations demonstrate that the hybrid β-sheet stability is strongly correlated with the formation of a stable cleavage-compatible active geometry and it also controls water access to the active site. The hybrid β-sheet is only stable for substrates with 3 or more C-terminal residues beyond the scissile bond. The simulation model allowed us to predict the effect of Pro and Phe mutations that weaken the formation of the hybrid β-sheet which were confirmed by experimental testing. Our study provides a direct explanation why γ-secretase preferentially cleaves APP in steps of 3 residues and how the hybrid β-sheet facilitates γ-secretase proteolysis.

Effects of presenilin-1 familial Alzheimer's disease mutations on γ-secretase activation for cleavage of amyloid precursor protein

Presenilin-1 (PS1) is the catalytic subunit of γ-secretase which cleaves within the transmembrane domain of over 150 peptide substrates. Dominant missense mutations in PS1 cause early-onset familial Alzheimer's disease (FAD); however, the exact pathogenic mechanism remains unknown. Here we combined Gaussian accelerated molecular dynamics (GaMD) simulations and biochemical experiments to determine the effects of six representative PS1 FAD mutations (P117L, I143T, L166P, G384A, L435F, and L286V) on the enzyme-substrate interactions between γ-secretase and amyloid precursor protein (APP). Biochemical experiments showed that all six PS1 FAD mutations rendered γ-secretase less active for the endoproteolytic (ε) cleavage of APP. Distinct low-energy conformational states were identified from the free energy profiles of wildtype and PS1 FAD-mutant γ-secretase. The P117L and L286V FAD mutants could still sample the "Active" state for substrate cleavage, but with noticeably reduced conformational space compared with the wildtype. The other mutants hardly visited the "Active" state. The PS1 FAD mutants were found to reduce γ-secretase proteolytic activity by hindering APP residue L49 from proper orientation in the active site and/or disrupting the distance between the catalytic aspartates. Therefore, our findings provide mechanistic insights into how PS1 FAD mutations affect structural dynamics and enzyme-substrate interactions of γ-secretase and APP.

Elucidating the Protonation State of the γ-Secretase Catalytic Dyad

γ-Secretase (GS) is an intramembrane aspartyl protease that participates in the sequential cleavage of C99 to generate different isoforms of the amyloid-β (Aβ) peptides that are associated with the development of Alzheimer's disease. Due to its importance in the proteolytic processing of C99 by GS, we performed pH replica exchange molecular dynamics (pH-REMD) simulations of GS in its apo and substrate-bound forms to sample the protonation states of the catalytic dyad. We found that the catalytic dyad is deprotonated at physiological pH in our apo form, but the presence of the substrate at the active site displaces its monoprotonated state toward physiological pH. Our results show that Asp257 acts as the general base and Asp385 as the general acid during the cleavage mechanism. We identified different amino acids such as Lys265, Arg269, and the PAL motif interacting with the catalytic dyad and promoting changes in its acid-base behavior. Finally, we also found a significant pK_a shift of Glu280 related to the internalization of TM6-CT in the GS-apo form. Our study provides critical mechanistic insight into the GS mechanism and the basis for future research on the genesis of Aβ peptides and the development of Alzheimer's disease.

A thermodynamic investigation of amyloid precursor protein processing by human γ-secretase

Human γ-secretase cleaves the transmembrane domains (TMDs) of amyloid precursor protein (APP) into pathologically relevant amyloid-β peptides (Aβs). The detailed mechanisms of the unique endoproteolytic cleavage by the presenilin 1 domain (PS1) of γ-secretase are still poorly understood. Herein, we provide thermodynamic insights into how the α-helical APP TMD is processed by γ-secretase and elucidate the specificity of Aβ48/Aβ49 cleavage using unbiased molecular dynamics and bias-exchange metadynamics simulations. The thermodynamic data show that the unwinding of APP TMD is driven by water hydration in the intracellular pocket of PS1, and the scissile bond T32-L33 or L33-V34 of the APP TMD can slide down and up to interact with D257/D385 to achieve endoproteolysis. In the wild-type system, the L33-V34 scissile bond is more easily hijacked by D257/D385 than T32-L33, resulting in higher Aβ49 cleavage, while the T32N mutation on the APP TMD decreases the energy barrier of the sliding of the scissile bonds and increases the hydrogen bond occupancy for Aβ48 cleavage. In summary, the thermodynamic analysis elucidates possible mechanisms of APP TMD processing by PS1, which might facilitate rational drug design targeting γ-secretase.

Thermodynamic analysis from unbiased molecular dynamics and bias-exchange metadynamics simulations reveals possible mechanisms on how γ-secretase cleaves the transmembrane domains of amyloid precursor protein into amyloid-β peptides.

Structure and mechanism of the γ-secretase intramembrane protease complex

γ-Secretase is a membrane protein complex that proteolyzes within the transmembrane domain of >100 substrates, including those derived from the amyloid precursor protein and the Notch family of cell surface receptors. The nine-transmembrane presenilin is the catalytic component of this aspartyl protease complex that carries out hydrolysis in the lipid bilayer. Advances in cryoelectron microscopy have led to the elucidation of the structure of the γ-secretase complex at atomic resolution. Recently, structures of the enzyme have been determined with bound APP- or Notch-derived substrates, providing insight into the nature of substrate recognition and processing. Molecular dynamics simulations of substrate-bound enzymes suggest dynamic mechanisms of intramembrane proteolysis. Structures of the enzyme bound to small-molecule inhibitors and modulators have also been solved, setting the stage for rational structure-based drug discovery targeting γ-secretase.

Mechanism of Tripeptide Trimming of Amyloid β-Peptide 49 by γ-Secretase

The membrane-embedded γ-secretase complex processively cleaves within the transmembrane domain of amyloid precursor protein (APP) to produce 37-to-43-residue amyloid β-peptides (Aβ) of Alzheimer's disease (AD). Despite its importance in pathogenesis, the mechanism of processive proteolysis by γ-secretase remains poorly understood. Here, mass spectrometry and Western blotting were used to quantify the efficiency of tripeptide trimming of wild-type (WT) and familial AD (FAD) mutant Aβ49. In comparison to WT Aβ49, the efficiency of tripeptide trimming was similar for the I45F, A42T, and V46F Aβ49 FAD mutants but substantially diminished for the I45T and T48P mutants. In parallel with biochemical experiments, all-atom simulations using a novel peptide Gaussian accelerated molecular dynamics (Pep-GaMD) method were applied to investigate the tripeptide trimming of Aβ49 by γ-secretase. The starting structure was the active γ-secretase bound to Aβ49 and APP intracellular domain (AICD), as generated from our previous study that captured the activation of γ-secretase for the initial endoproteolytic cleavage of APP (Bhattarai, A., ACS Cent. Sci. 2020, 6, 969-983). Pep-GaMD simulations captured remarkable structural rearrangements of both the enzyme and substrate, in which hydrogen-bonded catalytic aspartates and water became poised for tripeptide trimming of Aβ49 to Aβ46. These structural changes required a positively charged N-terminus of endoproteolytic coproduct AICD, which could dissociate during conformational rearrangements of the protease and Aβ49. The simulation findings were highly consistent with biochemical experimental data. Taken together, our complementary biochemical experiments and Pep-GaMD simulations have enabled elucidation of the mechanism of tripeptide trimming of Aβ49 by γ-secretase.

Mechanistic regulation of γ-secretase by their substrates

γ-Secretase (GS) is a transmembrane (TM) enzyme that plays important roles in the processing of approximately 90 substrates.

Characterizing the Chemical Space of γ-Secretase Inhibitors and Modulators

γ-Secretase (GS) is one of the most attractive molecular targets for the treatment of Alzheimer's disease (AD). Its key role in the final step of amyloid-β peptides generation and its relationship in the cascade of events for disease development have caught the attention of many pharmaceutical groups. Over the past years, different inhibitors and modulators have been evaluated as promising therapeutics against AD. However, despite the great chemical diversity of the reported compounds, a global classification and visual representation of the chemical space for GS inhibitors and modulators remain unavailable. In the present work, we carried out a two-dimensional (2D) chemical space analysis from different classes and subclasses of GS inhibitors and modulators based on their structural similarity. Along with the novel structural information available for GS complexes, our analysis opens the possibility to identify compounds with high molecular similarity, critical to finding new chemical structures through the optimization of existing compounds and relating them with a potential binding site.

Role of cholesterol in substrate recognition by [Formula: see text]-secretase

[Formula: see text]-Secretase is an enzyme known to cleave multiple substrates within their transmembrane domains, with the amyloid precursor protein of Alzheimer's Disease among the most prominent examples. The activity of [Formula: see text]-secretase strictly depends on the membrane cholesterol content, yet the mechanistic role of cholesterol in the substrate binding and cleavage remains unclear. In this work, we used all-atom molecular dynamics simulations to examine the role of cholesterol in the initial binding of a direct precursor of [Formula: see text]-amyloid polypeptides by [Formula: see text]-secretase. We showed that in cholesterol-rich membranes, both the substrate and the enzyme region proximal to the active site induce a local membrane thinning. With the free energy methods we found that in the presence of cholesterol the substrate binds favorably to the identified exosite, while cholesterol depletion completely abolishes the binding. To explain these findings, we directly examined the role of hydrophobic mismatch in the substrate binding to [Formula: see text]-secretase, showing that increased membrane thickness results in higher propensity of the enzyme to bind substrates. Therefore, we propose that cholesterol promotes substrate binding to [Formula: see text]-secretase by increasing the membrane thickness, which leads to the negative hydrophobic mismatch between the membrane and binding partners.

The Role of Cholesterol in Amyloidogenic Substrate Binding to the γ-Secretase Complex

Alzheimer’s disease is the most common progressive neurodegenerative disorder and is characterized by the presence of amyloid β (Aβ) plaques in the brain. The γ-secretase complex, which produces Aβ, is an intramembrane-cleaving protease consisting of four membrane proteins. In this paper we investigated the amyloidogenic fragments of amyloid precursor protein (substrates Aβ₄₃ and Aβ₄₅, leading to less amyloidogenic Aβ₄₀ and more amyloidogenic Aβ₄₂, respectively) docked to the binding site of presenilin, the catalytic subunit of γ-secretase. In total, we performed 9 μs of all-atom molecular dynamics simulations of the whole γ-secretase complex with both substrates in low (10%) and high (50%) concentrations of cholesterol in the membrane. We found that, at the high cholesterol level, the Aβ₄₅ helix was statistically more flexible in the binding site of presenilin than Aβ₄₃. An increase in the cholesterol concentration was also correlated with a higher flexibility of the Aβ₄₅ helix, which suggests incompatibility between Aβ₄₅ and the binding site of presenilin potentiated by a high cholesterol level. However, at the C-terminal part of Aβ₄₅, the active site of presenilin was more compact in the case of a high cholesterol level, which could promote processing of this substrate. We also performed detailed mapping of the cholesterol binding sites at low and high cholesterol concentrations, which were independent of the typical cholesterol binding motifs.

Protein Predictive Modeling and Simulation of Mutations of Presenilin-1 Familial Alzheimer's Disease on the Orthosteric Site

Alzheimer's disease pathology is characterized by β-amyloid plaques and neurofibrillary tangles. Amyloid precursor protein is processed by β and γ secretase, resulting in the production of β-amyloid peptides with a length ranging from 38 to 43 amino acids. Presenilin 1 (PS1) is the catalytic unit of γ-secretase, and more than 200 PS1 pathogenic mutations have been identified as causative for Alzheimer's disease. A complete monocrystal structure of PS1 has not been determined so far due to the presence of two flexible domains. We have developed a complete structural model of PS1 using a computational approach with structure prediction software. Missing fragments Met1-Glut72 and Ser290-Glu375 were modeled and validated by their energetic and stereochemical characteristics. Then, with the complete structure of PS1, we defined that these fragments do not have a direct effect in the structure of the pore. Next, we used our hypothetical model for the analysis of the functional effects of PS1 mutations Ala246GLu, Leu248Pro, Leu248Arg, Leu250Val, Tyr256Ser, Ala260Val, and Val261Phe, localized in the catalytic pore. For this, we used a quantum mechanics/molecular mechanics (QM/MM) hybrid method, evaluating modifications in the topology, potential surface density, and electrostatic potential map of mutated PS1 proteins. We found that each mutation exerts changes resulting in structural modifications of the active site and in the shape of the pore. We suggest this as a valid approach for functional studies of PS1 in view of the possible impact in substrate processing and for the design of targeted therapeutic strategies.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

How Mutations Perturb γ-Secretase Active Site Studied by Free Energy Simulations

γ-Secretase is involved in processing of the amyloid precursor protein (APP) and generation of short Aβ peptides that may play a key role in neurodegenerative diseases such as Alzheimer's disease (AD). Several mutations in γ-secretase influence its activity, resulting in early AD onset (Familial AD or FAD mutations). The molecular details of how mutations, not located close to the active site, can affect enzyme activity is not understood. In molecular dynamics simulations of γ-secretase in the absence of substrate (apo), we identified two active site conformational states characterized by a direct contact between catalytic Asp residues (closed state) and an open water-bridged state. In the presence of substrate, only conformations compatible with the open active site geometry are accessible. Systematic free energy simulations on wild type and FAD mutations indicate a free energy difference between closed and open states that is significantly modulated by FAD mutations and correlates with the corresponding experimental activity. For mutations with reduced activity, an increased penalty for open-state transitions was found. Only for two mutations located at the active site a direct perturbation of the open-state geometry was observed that could directly explain the drop of enzyme activity. The simulations suggest that modulation of the closed/open equilibrium and perturbation of the open (active) catalytic geometry are possible mechanisms of how FAD mutations affect γ-secretase activity. The results also offer an explanation for the experimental finding that FAD mutations, although not located at the interface to the substrate, mainly destabilize the enzyme-substrate complex.

Membrane dynamics of γ-secretase with the anterior pharynx-defective 1B subunit

The four-subunit protease complex γ-secretase cleaves many single-pass transmembrane (TM) substrates, including Notch and β-amyloid precursor protein to generate amyloid-β (Aβ), central to Alzheimer's disease. Two of the subunits anterior pharynx-defective 1 (APH-1) and presenilin (PS) exist in two homologous forms APH1-A and APH1-B, and PS1 and PS2. The consequences of these variations are poorly understood and could affect Aβ production and γ-secretase medicine. Here, we developed the first complete structural model of the APH-1B subunit using the published cryo-electron microscopy (cryo-EM) structures of APH1-A (Protein Data Bank: 5FN2, 5A63, and 6IYC). We then performed all-atom molecular dynamics simulations at 303 K in a realistic bilayer system to understand both APH-1B alone and in γ-secretase without and with substrate C83-bound. We show that APH-1B adopts a 7TM topology with a water channel topology similar to APH-1A. We demonstrate direct transport of water through this channel, mainly via Glu84, Arg87, His170, and His196. The apo and holo states closely resemble the experimental cryo-EM structures with APH-1A, however with subtle differences: The substrate-bound APH-1B γ-secretase was quite stable, but some TM helices of PS1 and APH-1B rearranged in the membrane consistent with the disorder seen in the cryo-EM data. This produces different accessibility of water molecules for the catalytic aspartates of PS1, critical for Aβ production. In particular, we find that the typical distance between the catalytic aspartates of PS1 and the C83 cleavage sites are shorter in APH-1B, that is, it represents a more closed state, due to interactions with the C-terminal fragment of PS1. Our structural-dynamic model of APH-1B alone and in γ-secretase suggests generally similar topology but some notable differences in water accessibility which may be relevant to the protein's existence in two forms and their specific function and location.

Side-by-side comparison of Notch- and C83 binding to γ-secretase in a complete membrane model at physiological temperature

γ-Secretase cleaves the C99 fragment of the amyloid precursor protein, leading to formation of aggregated β-amyloid peptide central to Alzheimer's disease, and Notch, essential for cell regulation. Recent cryogenic electron microscopy (cryo-EM) structures indicate major changes upon substrate binding, a β-sheet recognition motif, and a possible helix unwinding to expose peptide bonds towards nucleophilic attack. Here we report side-by-side comparison of the 303 K dynamics of the two proteins in realistic membranes using molecular dynamics simulations. Our ensembles agree with the cryo-EM data (full-protein Cα-RMSD = 1.62–2.19 Å) but reveal distinct presenilin helix conformation states and thermal β-strand to coil transitions of C83 and Notch100. We identify distinct 303 K hydrogen bond dynamics and water accessibility of the catalytic sites. The RKRR motif (1758–1761) contributes significantly to Notch binding and serves as a “membrane anchor” that prevents Notch displacement. Water that transiently hydrogen bonds to G1753 and V1754 probably represents the catalytic nucleophile. At 303 K, Notch and C83 binding induce different conformation states, with Notch mostly present in a closed state with shorter Asp–Asp distance. This may explain the different outcome of Notch and C99 cleavage, as the latter is more imprecise with many products. Our identified conformation states may aid efforts to develop conformation-selective drugs that target C99 and Notch cleavage differently, e.g. Notch-sparing γ-secretase modulators.

Distinct membrane dynamics and conformations of C83- and Notch-bound γ-secretase may aid the development of Notch-sparing treatments of Alzheimer's disease.

Mechanisms of γ-Secretase Activation and Substrate Processing

Amyloid β-peptide, the principal component of characteristic cerebral plaques of Alzheimer's disease (AD), is produced through intramembrane proteolysis of the amyloid precursor protein (APP) by γ-secretase. Despite the importance in the pathogenesis of AD, the mechanisms of intramembrane proteolysis and substrate processing by γ-secretase remain poorly understood. Here, complementary all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method and biochemical experiments were combined to investigate substrate processing of wildtype and mutant APP by γ-secretase. The GaMD simulations captured spontaneous activation of γ-secretase, with hydrogen bonded catalytic aspartates and water poised for proteolysis of APP at the ε cleavage site. Furthermore, GaMD simulations revealed that familial AD mutations I45F and T48P enhanced the initial ε cleavage between residues Leu49-Val50, while M51F mutation shifted the ε cleavage site to the amide bond between Thr48-Leu49. Detailed analysis of the GaMD simulations allowed us to identify distinct low-energy conformational states of γ-secretase, different secondary structures of the wildtype and mutant APP substrate, and important active-site subpockets for catalytic function of the enzyme. The simulation findings were highly consistent with experimental analyses of APP proteolytic products using mass spectrometry and Western blotting. Taken together, the GaMD simulations and biochemical experiments have enabled us to elucidate the mechanisms of γ-secretase activation and substrate processing, which should facilitate rational computer-aided drug design targeting this functionally important enzyme.

Uncovering the Binding Mode of γ -Secretase Inhibitors

Knowledge of how transition state inhibitors bind to γ-secretase is of major importance for the design of new Alzheimer's disease therapies. On the basis of the known structure of γ-secretase in complex with a fragment of the amyloid precursor protein, we generated a structural model of γ-secretase in complex with the effective L-685,458 transition state inhibitor. The predicted binding mode is in excellent agreement with experimental data, mimicking all enzyme-substrate interactions at the active site and forming the relevant transition state geometry with the active site aspartate residues. The model also indicates the possible location and nature of the amino acid residues forming the proposed binding pockets S1', S2', and S3' near the active site that are occupied by chemical groups of the inhibitor. In addition, we found that the stability of the complex is very likely sensitive to the pH value. Comparative simulations on the binding of L-685,458 and the epimer L682,679 allowed us to explain the strongly reduced affinity of the epimer for γ-secretase. The structural model could form a valuable basis for the design of new or modified γ-secretase inhibitors.

Structure and dynamics of γ-secretase with presenilin 2 compared to presenilin 1

Severe early-onset familial Alzheimer's disease (FAD) is caused by more than 200 different mutations in the genes coding for presenilin, the catalytic subunit of the 4-subunit protease complex γ-secretase, which cleaves the C99 fragment of the amyloid precursor protein (APP) to produce Aβ peptides. γ-Secretase exists with either of two homologues, PS1 and PS2. All cryo-electron microscopic structures and computational work has so far focused on γ-secretase with PS1, yet PS2 mutations also cause FAD. A central question is thus whether there are structural and dynamic differences between PS1 and PS2. To address this question, we use the cryo-electron microscopic data for PS1 to develop the first structural and dynamic model of PS2-γ-secretase in the catalytically relevant mature membrane-bound state at ambient temperature, equilibrated by three independent 500 ns molecular dynamics simulations. We find that the characteristic nicastrin extra-cellular domain breathing mode and major movements in the cytosolic loop between TM6 and TM7 occur in both PS2- and PS1-γ-secretase. The overall structures and conformational states are similar, suggesting similar catalytic activities. However, at the sequence level, charge-controlled membrane-anchoring is extracellular for PS1 and intracellular for PS2, which suggests different subcellular locations. The tilt angles of the TM2, TM6, TM7 and TM9 helices differ in the two forms of γ-secretase, suggesting that the two proteins have somewhat different substrate processing and channel sizes. Our MD simulations consistently indicated that PS2 retains several water molecules near the catalytic site at the bilayer, as required for catalysis. The possible reasons for the differences of PS1 and PS2 are discussed in relation to their location and function.

APH-1A Component of γ-Secretase Forms an Internal Water and Ion-Containing Cavity

Anterior pharynx-defective 1A (APH-1A) is a seven transmembrane component of γ-secretase (GS), an aspartyl protease enzyme involved in the production of toxic amyloid-β peptides in Alzheimer's disease patients. Cryo-electron microscopy structures of the enzyme complex revealed a central cavity in its APH-1A component, similar to water-containing cavities in G-protein coupled receptors (GPCRs). In this work, we performed molecular dynamics and umbrella sampling simulations to understand the role of the APH-1A cavity in the GS complex. Our results suggest that APH-1A is able to store water molecules in its inner cavity and transport some of them between cell spaces. Additionally, APH-1A allows the influx of extracellular cations into a central hydrophilic cavity but cannot transport them into the intracellular space. Overall, this study seeks to describe an alternative APH-1A function in GS besides its complex stabilization role and provide novel approaches to understand the functioning of the GS enzyme.

Modulating Hinge Flexibility in the APP Transmembrane Domain Alters γ-Secretase Cleavage

Intramembrane cleavage of the β-amyloid precursor protein C99 substrate by γ-secretase is implicated in Alzheimer's disease pathogenesis. Biophysical data have suggested that the N-terminal part of the C99 transmembrane domain (TMD) is separated from the C-terminal cleavage domain by a di-glycine hinge. Because the flexibility of this hinge might be critical for γ-secretase cleavage, we mutated one of the glycine residues, G38, to a helix-stabilizing leucine and to a helix-distorting proline. Both mutants impaired γ-secretase cleavage and also altered its cleavage specificity. Circular dichroism, NMR, and backbone amide hydrogen/deuterium exchange measurements as well as molecular dynamics simulations showed that the mutations distinctly altered the intrinsic structural and dynamical properties of the substrate TMD. Although helix destabilization and/or unfolding was not observed at the initial ε-cleavage sites of C99, subtle changes in hinge flexibility were identified that substantially affected helix bending and twisting motions in the entire TMD. These resulted in altered orientation of the distal cleavage domain relative to the N-terminal TMD part. Our data suggest that both enhancing and reducing local helix flexibility of the di-glycine hinge may decrease the occurrence of enzyme-substrate complex conformations required for normal catalysis and that hinge mobility can thus be conducive for productive substrate-enzyme interactions.

Structural Modeling of γ-Secretase Aβ n Complex Formation and Substrate Processing

The intramembrane aspartyl protease γ-secretase (GSEC) cleaves single-span transmembrane helices including the C-terminal fragment of the amyloid precursor protein (APP). This substrate is initially cleaved at the ϵ-site followed by successive processing (trimming) events mostly in steps of three amino acids. GSEC is responsible for the formation of N-terminal APP amyloid-β (Aβ) peptides of different length (e.g., Aβ₄₂) that can form aggregates involved in Alzheimer's disease pathogenesis. The molecular mechanism of GSEC-APP substrate recognition is key for understanding how different peptide products are formed and could help in designing APP-selective modulators. Based on the known structure of apo GSEC and the APP-C99 fragment we have generated putative structural models of the initial binding in three different possible modes using extensive molecular dynamics (MD) simulations. The binding mode with the substrate helix located in a cleft between the transmembrane helices 2 and 3 of the presenilin subunit was identified as a most likely binding mode. Based on this arrangement, the processing steps were investigated using restraint MD simulations to pull the scissile bond (for each processing step) into a transition like (cleavable) state. This allowed us to analyze in detail the motions and energetic contributions of participating residues. The structural model agrees qualitatively well with the influence of many mutations in GSEC and C99. It also explains the effects of inhibitors, cross-linking, as well as spectroscopic data on GSEC substrate binding and can serve as working model for the future planning of structural and biochemical studies.