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

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.

Conformational Changes and Unfolding of β-Amyloid Substrates in the Active Site of γ-Secretase

Alzheimer's disease (AD) is the leading cause of dementia and is characterized by a presence of amyloid plaques, composed mostly of the amyloid-β (Aβ) peptides, in the brains of AD patients. The peptides are generated from the amyloid precursor protein (APP), which undergoes a sequence of cleavages, referred as trimming, performed by γ-secretase. Here, we investigated conformational changes in a series of β-amyloid substrates (from less and more amyloidogenic pathways) in the active site of presenilin-1, the catalytic subunit of γ-secretase. The substrates are trimmed every three residues, finally leading to Aβ40 and Aβ42, which are the major components of amyloid plaques. To study conformational changes, we employed all-atom molecular dynamics simulations, while for unfolding, we used steered molecular dynamics simulations in an implicit membrane-water environment to accelerate changes. We have found substantial differences in the flexibility of extended C-terminal parts between more and less amyloidogenic pathway substrates. We also propose that the positively charged residues of presenilin-1 may facilitate the stretching and unfolding of substrates. The calculated forces and work/energy of pulling were exceptionally high for Aβ40, indicating why trimming of this substrate is so infrequent.

Cholesterol-dependent amyloid β production: space for multifarious interactions between amyloid precursor protein, secretases, and cholesterol

Amyloid β is considered a key player in the development and progression of Alzheimer's disease (AD). Many studies investigating the effect of statins on lowering cholesterol suggest that there may be a link between cholesterol levels and AD pathology. Since cholesterol is one of the most abundant lipid molecules, especially in brain tissue, it affects most membrane-related processes, including the formation of the most dangerous form of amyloid β, Aβ42. The entire Aβ production system, which includes the amyloid precursor protein (APP), β-secretase, and the complex of γ-secretase, is highly dependent on membrane cholesterol content. Moreover, cholesterol can affect amyloidogenesis in many ways. Cholesterol influences the stability and activity of secretases, but also dictates their partitioning into specific cellular compartments and cholesterol-enriched lipid rafts, where the amyloidogenic machinery is predominantly localized. The most complicated relationships have been found in the interaction between cholesterol and APP, where cholesterol affects not only APP localization but also the precise character of APP dimerization and APP processing by γ-secretase, which is important for the production of Aβ of different lengths. In this review, we describe the intricate web of interdependence between cellular cholesterol levels, cholesterol membrane distribution, and cholesterol-dependent production of Aβ, the major player in AD.

Docking cholesterol to integral membrane proteins with Rosetta

Lipid molecules such as cholesterol interact with the surface of integral membrane proteins (IMP) in a mode different from drug-like molecules in a protein binding pocket. These differences are due to the lipid molecule's shape, the membrane's hydrophobic environment, and the lipid's orientation in the membrane. We can use the recent increase in experimental structures in complex with cholesterol to understand protein-cholesterol interactions. We developed the RosettaCholesterol protocol consisting of (1) a prediction phase using an energy grid to sample and score native-like binding poses and (2) a specificity filter to calculate the likelihood that a cholesterol interaction site may be specific. We used a multi-pronged benchmark (self-dock, flip-dock, cross-dock, and global-dock) of protein-cholesterol complexes to validate our method. RosettaCholesterol improved sampling and scoring of native poses over the standard RosettaLigand baseline method in 91% of cases and performs better regardless of benchmark complexity. On the β2AR, our method found one likely-specific site, which is described in the literature. The RosettaCholesterol protocol quantifies cholesterol binding site specificity. Our approach provides a starting point for high-throughput modeling and prediction of cholesterol binding sites for further experimental validation.

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.

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.