Characterizing the structural ensemble of γ-secretase using a multiscale molecular dynamics approach

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.

Molecular Dynamics Activation of γ-Secretase for Cleavage of the Notch1 Substrate

γ-Secretase is an intramembrane aspartyl protease complex that cleaves the transmembrane domain of over 150 peptide substrates, including amyloid precursor protein (APP) and the Notch family of receptors, via two conserved aspartates D257 and D385 in the presenilin-1 (PS1) catalytic subunit. However, while the activation of γ-secretase for cleavage of APP has been widely studied, the cleavage of Notch by γ-secretase remains poorly explored. Here, we combined Gaussian accelerated molecular dynamics (GaMD) simulations and mass spectrometry (MS) analysis of proteolytic products to present the first dynamic models for cleavage of Notch by γ-secretase. MS showed that γ-secretase cleaved the WT Notch at Notch residue G34, while cleavage of the L36F mutant Notch occurred at Notch residue C33. Initially, we prepared our simulation systems starting from the cryoEM structure of Notch-bound γ-secretase (PDB: 6IDF) and failed to capture the proper cleavages of WT and L36F Notch by γ-secretase. We then discovered an incorrect registry of the Notch substrate in the PS1 active site through alignment of the experimental structure of Notch-bound (PDB: 6IDF) and APP-bound γ-secretase (PDB: 6IYC). Every residue of the APP substrate was systematically mutated to the corresponding Notch residue to prepare a resolved model of Notch-bound γ-secretase complexes. GaMD simulations of the resolved model successfully captured γ-secretase activation for proper cleavages of both WT and L36F mutant Notch. Our findings presented here provided mechanistic insights into the structural dynamics and enzyme-substrate interactions required for γ-secretase activation for cleavage of Notch and other substrates.

Interaction of Substrates with γ-Secretase at the Level of Individual Transmembrane Helices-A Methodological Approach

Intramembrane proteases, such as γ secretase, typically recruit multiple substrates from an excess of single-span membrane proteins. It is currently unclear to which extent substrate recognition depends on specific interactions of their transmembrane domains (TMDs) with TMDs of a protease. Here, we investigated a large number of potential pairwise interactions between TMDs of γ secretase and a diverse set of its substrates using two different configurations of BLaTM, a genetic reporter system. Our results reveal significant interactions between TMD2 of presenilin, the enzymatic subunit of γ secretase, and the TMD of the amyloid precursor protein, as well as of several other substrates. Presenilin TMD2 is a prime candidate for substrate recruitment, as has been shown from previous studies. In addition, the amyloid precursor protein TMD enters interactions with presenilin TMD 4 as well as with the TMD of nicastrin. Interestingly, the Gly-rich interfaces between the amyloid precursor protein TMD and presenilin TMDs 2 and 4 are highly similar to its homodimerization interface. In terms of methodology, the economics of the newly developed library-based method could prove to be a useful feature in related future work for identifying heterotypic TMD-TMD interactions within other biological contexts.

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.

The Binding of Different Substrate Molecules at the Docking Site and the Active Site of γ-Secretase Can Trigger Toxic Events in Sporadic and Familial Alzheimer's Disease

Pathogenic changes in γ-secretase activity, along with its response to different drugs, can be affected by changes in the saturation of γ-secretase with its substrate. We analyze the saturation of γ-secretase with its substrate using multiscale molecular dynamics studies. We found that an increase in the saturation of γ-secretase with its substrate could result in the parallel binding of different substrate molecules at the docking site and the active site. The C-terminal domain of the substrate bound at the docking site can interact with the most dynamic presenilin sites at the cytosolic end of the active site tunnel. Such interactions can inhibit the ongoing catalytic activity and increase the production of the longer, more hydrophobic, and more toxic Aβ proteins. Similar disruptions in dynamic presenilin structures can be observed with different drugs and disease-causing mutations. Both, C99-βCTF-APP substrate and its different Aβ products, can support the toxic aggregation. The aggregation depends on the substrate N-terminal domain. Thus, the C99-βCTF-APP substrate and β-secretase path can be more toxic than the C83-αCTF-APP substrate and α-secretase path. Nicastrin can control the toxic aggregation in the closed conformation. The binding of the C99-βCTF-APP substrate to γ-secretase can be controlled by substrate channeling between the nicastrin and β-secretase. We conclude that the presented two-substrate mechanism could explain the pathogenic changes in γ-secretase activity and Aβ metabolism in different sporadic and familial cases of Alzheimer's disease. Future drug-development efforts should target different cellular mechanisms that regulate the optimal balance between γ-secretase activity and amyloid metabolism.

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.

Emergence of a recurrent insertion in the N-terminal domain of the SARS-CoV-2 spike glycoprotein

Tracking the evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) through genomic surveillance programs is undoubtedly one of the key priorities in the current pandemic situation. Although the genome of SARS-CoV-2 acquires mutations at a slower rate compared with other RNA viruses, evolutionary pressures derived from the widespread circulation of SARS-CoV-2 in the human population have progressively favored the global emergence, though natural selection, of several variants of concern that carry multiple non-synonymous mutations in the spike glycoprotein. These are often placed in key sites within major antibody epitopes and may therefore confer resistance to neutralizing antibodies, leading to partial immune escape, or otherwise compensate infectivity deficits associated with other non-synonymous substitutions. As previously shown by other authors, several emerging variants carry recurrent deletion regions (RDRs) that display a partial overlap with antibody epitopes located in the spike N-terminal domain (NTD). Comparatively, very little attention had been directed towards spike insertion mutations prior to the emergence of the B.1.1.529 (omicron) lineage. This manuscript describes a single recurrent insertion region (RIR1) in the N-terminal domain of SARS-CoV-2 spike protein, characterized by at least 49 independent acquisitions of 1-8 additional codons between Val213 and Leu216 in different viral lineages. Even though RIR1 is unlikely to confer antibody escape, its association with two distinct formerly widespread lineages (A.2.5 and B.1.214.2), with the quickly spreading omicron and with other VOCs and VOIs warrants further investigation concerning its effects on spike structure and viral infectivity.

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.

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.

Deletion of Dcf1 Reduces Amyloid-β Aggregation and Mitigates Memory Deficits

BACKGROUND

Alzheimer's disease (AD) is a progressive neurodegenerative disease. One of the pathologies of AD is the accumulation of amyloid-β (Aβ) to form senile plaques, leading to a decline in cognitive ability and a lack of learning and memory. However, the cause leading to Aβ aggregation is not well understood. Dendritic cell factor 1 (Dcf1) shows a high expression in the entorhinal cortex neurons and neurofibrillary tangles in AD patients.

OBJECTIVE

Our goal is to investigate the effect of Dcf1 on Aβ aggregation and memory deficits in AD development.

METHODS

The mouse and Drosophila AD model were used to test the expression and aggregation of Aβ, senile plaque formation, and pathological changes in cognitive behavior during dcf1 knockout and expression. We finally explored possible drug target effects through intracerebroventricular delivery of Dcf1 antibodies.

RESULTS

Deletion of Dcf1 resulted in decreased Aβ42 level and deposition, and rescued AMPA Receptor (GluA2) levels in the hippocampus of APP-PS1-AD mice. In Aβ42 AD Drosophila, the expression of Dcf1 in Aβ42 AD flies aggravated the formation and accumulation of senile plaques, significantly reduced its climbing ability and learning-memory. Data analysis from all 20 donors with and without AD patients aged between 80 and 90 indicated a high-level expression of Dcf1 in the temporal neocortex. Dcf1 contributed to Aβ aggregation by UV spectroscopy assay. Intracerebroventricular delivery of Dcf1 antibodies in the hippocampus reduced the area of senile plaques and reversed learning and memory deficits in APP-PS1-AD mice.

CONCLUSION

Dcf1 causes Aβ-plaque accumulation, inhibiting dcf1 expression could potentially offer therapeutic avenues.

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.

Cryo-temperature effects on membrane protein structure and dynamics

Cryo-electron structures revolutionize biology, yet cooling effects are unclear. Using a simulation protocol of hot, cold, and rapidly cooled γ-secretase we identify cryo-contraction and modes relevant to Aβ production and cryo-analysis in general.

Insights into membrane-bound presenilin 2 from all-atom molecular dynamics simulations

Presenilins 1 and 2 (PS1 or PS2) are main genetic risk factors of familial Alzheimer's disease (AD) that produce the β-amyloid (Aβ) peptides and also have important stand-alone functions related to, e.g. calcium signaling. Most work so far has focused on PS1, but humans carry both PS1 and PS2, and mutations in both cause AD. Here, we develop a computational model of PS2 in the membrane to address the question how pathogenic PS2 mutations affect the membrane-embedded protein. The models are based on cryo-electron microscopy structures of PS1 translated to PS2, augmented with missing residues and a complete all-atom membrane-water system, and equilibrated using three independent 500-ns simulations of molecular dynamics with a structure-balanced force field. We show that the nine-transmembrane channel structure is substantially controlled by major dynamics in the hydrophilic loop bridging TM6 and TM7, which functions as a 'plug' in the PS2 membrane channel. TM2, TM6, TM7 and TM9 flexibility controls the size of this channel. We find that most pathogenic PS2 mutations significantly reduce stability relative to random mutations, using a statistical ANOVA test with all possible mutations in the affected sites as a control. The associated loss of compactness may also impair calcium affinity. Remarkably, similar properties of the open state are known to impair the binding of substrates to γ-secretase, and we thus argue that the two mechanisms could be functionally related.Communicated by Ramaswamy H. Sarma.

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.

Molecular dynamics of C99-bound γ-secretase reveal two binding modes with distinct compactness, stability, and active-site retention: implications for Aβ production

The membrane protease γ-secretase cleaves the C99 fragment of the amyloid precursor protein, thus producing the Aβ peptides central to Alzheimer's disease. Cryo-electron microscopy has provided the topology but misses the membrane and loop parts that contribute to substrate binding. We report here an essentially complete atomic model of C99 within wild-type γ-secretase that respects all the experimental constraints and additionally describes loop, helix, and C99 substrate dynamics in a realistic all-atom membrane. Our model represents the matured auto-cleaved state required for catalysis. From two independent 500-ns molecular dynamic simulations, we identify two conformation states of C99 in equilibrium, a compact and a loose state. Our simulations provide a basis for C99 processing and Aβ formation and explain the production of longer and shorter Aβ, as the compact state retains C99 for longer and thus probably trims to shorter Aβ peptides. We expect pathogenic presenilin mutations to stabilize the loose over the compact state. The simulations detail the role of the Lys53-Lys54-Lys55 anchor for C99 binding, a loss of helicity of bound C99, and positioning of Thr48 and Leu49 leading to alternative trimming pathways on opposite sides of the C99 helix in three amino acid steps. The C99 binding topology resembles that of C83-bound γ-secretase without membrane but lacks a presenilin 1-C99 β-sheet, which could be induced by C83's stronger binding. The loose state should be selectively disfavored by γ-secretase modulators to increase C99 trimming and reduce the formation of longer Aβ, a strategy that is currently much explored but has lacked a structural basis.

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.

Quantifying correlations between mutational sites in the catalytic subunit of γ-secretase

Presenilin 1 (PS1) is the catalytic subunit of the γ-secretase complex which is involved in the generation of amyloid-β peptides (Aβ). Single point mutations in PS1 alter the cleavage pattern of the amyloid precursor protein (APP) and lead to the formation of aberrant Aβ peptides. To date, more than two hundred mutations distributed among almost a third of PS1's amino acids have been associated to the development of Alzheimer's disease (AD). Nevertheless, the mechanism by which mutations far from the catalytic site alter the γ-secretase's cleavage pattern remains unclear. In this work we analyzed correlated motions between amino acids in the wild type (WT) enzyme and 13 γ-secretase mutant models employing a multi-scale molecular dynamics approach. The effect of the protonation state of key catalytic residue Asp385 on the correlation networks was also evaluated. We observed that the strength and number of correlations is highly influenced in all mutant models in both protonation state models. The biggest changes were observed in mutants I83T, W165G, H214Y and L435F; the latest has been proved to drastically reduce γ-secretase activity. Finally, we made a classification of the studied mutations according to their correlation networks with amino acids at: (1) the interfaces with the other γ-secretase components, (2) the catalytic site, (3) the substrate entry site and (4) the substrate recognition site. Overall, this work provides insight into the allosteric communication networks of PS1.

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

γ-secretase, an intramembrane-cleaving aspartyl protease is involved in the cleavage of a large number of intramembrane proteins. The most prominent substrate is the amyloid precursor protein, whose proteolytic processing leads to the production of different amyloid Aβ peptides. These peptides are known to form toxic aggregates and may play a key role in Alzheimer's disease (AD). Recently, the three-dimensional structure of γ-secretase has been determined via Cryo-EM, elucidating the spatial geometry of this enzyme complex in different functional states. We have used molecular dynamics (MD) simulations to study the global dynamics and conformational transitions of γ-secretase, as well as the water and lipid distributions in and around the transmembrane domains in atomic detail. Simulations were performed on the full enzyme complex and on the membrane embedded parts alone. The simulations revealed global motions compatible with the experimental enzyme structures and indicated little dependence of the dynamics of the transmembrane domains on the soluble extracellular subunits. During the simulation on the membrane spanning part a transition between an inactive conformation (with catalytic residues far apart) toward a putatively active form (with catalytic residues in close proximity) has been observed. This conformational change is associated with a distinct rearrangement of transmembrane helices, a global compaction of the catalytically active presenilin subunit a change in the water structure near the active site and a rigidification of the protein fold. The observed conformational rearrangement allows the interpretation of the effect of several mutations on the activity of γ-secretase. A number of long-lived lipid binding sites could be identified on the membrane spanning surface of γ-secretase which may coincide with association regions of hydrophobic membrane helices to form putative substrate binding exosites.

Influence of membrane lipid composition on the structure and activity of γ-secretase

γ-Secretase (GS) is a multi-subunit membrane-embedded aspartyl protease that cleaves more than 80 integral membrane proteins, including the amyloid precursor protein (APP) to produce the amyloid-β (Aβ) peptide. Oligomerization and aggregation of the 42-amino acid length Aβ isoform in the brain has been associated with the development and progression of Alzheimer's disease (AD). Based on recent experimental structural studies and using multiscale computational modeling approaches, the conformational states and protein-membrane interactions of the GS complex embedded in six homogeneous and six heterogeneous lipid bilayers were characterized. In order to identify potential lipid and cholesterol binding sites, GS regions with high lipid/cholesterol occupancy values were analyzed using atomistic and coarse-grained simulations. Long lipid residence times were observed to be correlated with a large number of hydrogen bonds between the charged headgroups and key GS amino acids. This observation provides a plausible explanation for the inhibition of GS by charged lipids observed in previous experimental studies. Computed lateral pressure profiles suggest that higher transmembrane pressures favor active state conformations of the catalytic subunit. A probable mechanism for the regulation of the local stress response in cholesterol-rich multicomponent lipid bilayers is identified. Finally, it is demonstrated that interactions between the nicastrin extracellular domain and lipid headgroups leads to a compact structural conformation of the GS complex. Overall, this study provides valuable insight into the effect of bilayer lipid composition on the GS structural ensemble and its function.

Making the final cut: pathogenic amyloid-β peptide generation by γ-secretase

Alzheimer´s disease (AD) is a devastating neurodegenerative disease of the elderly population. Genetic evidence strongly suggests that aberrant generation and/or clearance of the neurotoxic amyloid-β peptide (Aβ) is triggering the disease. Aβ is generated from the amyloid precursor protein (APP) by the sequential cleavages of β- and γ-secretase. The latter cleavage by γ-secretase, a unique and fascinating four-component protease complex, occurs in the APP transmembrane domain thereby releasing Aβ species of 37-43 amino acids in length including the longer, highly pathogenic peptides Aβ42 and Aβ43. The lack of a precise understanding of Aβ generation as well as of the functions of other γ-secretase substrates has been one factor underlying the disappointing failure of γ-secretase inhibitors in clinical trials, but on the other side also been a major driving force for structural and in depth mechanistic studies on this key AD drug target in the past few years. Here we review recent breakthroughs in our understanding of how the γ-secretase complex recognizes substrates, of how it binds and processes β-secretase cleaved APP into different Aβ species, as well as the progress made on a question of outstanding interest, namely how clinical AD mutations in the catalytic subunit presenilin and the γ-secretase cleavage region of APP lead to relative increases of Aβ42/43. Finally, we discuss how the knowledge emerging from these studies could be used to therapeutically target this enzyme in a safe way.

G-score: A function to solve the puzzle of modeling the protonation states of β-secretase binding pocket

The population density concept has emerged as a proposal for the analysis of molecular dynamics results, the key characteristic of population density is the evaluation of the simultaneous occurrence of a set of relevant parameters for a system. However, despite its statistical strength, selection of the tolerance level for the comparison of different models may appear as arbitrary. This work introduces the G-score, a function which summarizes and categorizes the results of population density analysis. Additionally, it incorporates parameters based on rmsd and dihedral angles, besides the protein-protein and protein-ligand interatomic distances conventionally used, which complement each other to provide a better description of the behavior of the system. These newly-proposed tools were applied to determine the most probable protonation state of the aspartic dyad of BACE1, Asp93 and Asp289, in the presence of three types of transition state inhibitors namely: reduced amides, tertiary carbinamines and hydroxyethylamines. The results show a full agreement between G-score values and population density charts, with the advantage of allowing a quick and direct comparison among all the considered models. We anticipate that the simplicity of calculating the parameters employed in this study will permit the extensive use of population density and the G-score for other molecular systems.

Dissecting conformational changes in APP's transmembrane domain linked to ε-efficiency in familial Alzheimer's disease

The mechanism by which familial Alzheimer's disease (FAD) mutations within the transmembrane domain (TMD) of the Amyloid Precursor Protein (APP) affect ε-endoproteolysis is only poorly understood. Thereby, mutations in the cleavage domain reduce ε-efficiency of γ-secretase cleavage and some even shift entry into production lines. Since cleavage occurs within the TMD, a relationship between processing and TMD structure and dynamics seems obvious. Using molecular dynamic simulations, we dissect the dynamic features of wild-type and seven FAD-mutants into local and global components. Mutations consistently enhance hydrogen-bond fluctuations upstream of the ε-cleavage sites but maintain strong helicity there. Dynamic perturbation-response scanning reveals that FAD-mutants target backbone motions utilized in the bound state. Those motions, obscured by large-scale motions in the pre-bound state, provide (i) a dynamic mechanism underlying the proposed coupling between binding and ε-cleavage, (ii) key sites consistent with experimentally determined docking sites, and (iii) the distinction between mutants and wild-type.

Allosteric Modulation of Intact γ-Secretase Structural Dynamics

As a protease complex involved in the cleavage of amyloid precursor proteins that lead to the formation of amyloid β fibrils implicated in Alzheimer's disease, γ-secretase is an important target for developing therapeutics against Alzheimer's disease. γ-secretase is composed of four subunits: nicastrin (NCT) in the extracellular (EC) domain, presenilin-1 (PS1), anterior pharynx defective 1, and presenilin enhancer 2 in the transmembrane (TM) domain. NCT and PS1 play important roles in binding amyloid β precursor proteins and modulating PS1 catalytic activity. Yet, the molecular mechanisms of coupling between substrate/modulator binding and catalytic activity remain to be elucidated. Recent determination of intact human γ-secretase cryo-electron microscopy structure has opened the way for a detailed investigation of the structural dynamics of this complex. Our analysis, based on a membrane-coupled anisotropic network model, reveals two types of NCT motions, bending and twisting, with respect to PS1. These underlie the fluctuations between the "open" and "closed" states of the lid-like NCT with respect to a hydrophilic loop 1 (HL1) on PS1, thus allowing or blocking access of the substrate peptide (EC portion) to HL1 and to the neighboring helix TM2. In addition to this alternating access mechanism, fluctuations in the volume of the PS1 central cavity facilitate the exposure of the catalytic site for substrate cleavage. Druggability simulations show that γ-secretase presents several hot spots for either orthosteric or allosteric inhibition of catalytic activity, consistent with experimental data. In particular, a hinge region at the interface between the EC and TM domains, near the interlobe groove of NCT, emerges as an allo-targeting site that would impact the coupling between HL1/TM2 and the catalytic pocket, opening, to our knowledge, new avenues for structure-based design of novel allosteric modulators of γ-secretase protease activity.