Single-particle Cryo-EM and molecular dynamics simulations: A perfect match

Knowledge of the structure and dynamics of biomolecules is key to understanding the mechanisms underlying their biological functions. Single-particle cryo-electron microscopy (cryo-EM) is a powerful structural biology technique to characterize complex biomolecular systems. Here, we review recent advances of how Molecular Dynamics (MD) simulations are being used to increase and enhance the information extracted from cryo-EM experiments. We will particularly focus on the physics underlying these experiments, how MD facilitates structure refinement, in particular for heterogeneous and non-isotropic resolution, and how thermodynamic and kinetic information can be extracted from cryo-EM data.

Towards Interpretable Cryo-EM: Disentangling Latent Spaces of Molecular Conformations

Molecules are essential building blocks of life and their different conformations (i.e., shapes) crucially determine the functional role that they play in living organisms. Cryogenic Electron Microscopy (cryo-EM) allows for acquisition of large image datasets of individual molecules. Recent advances in computational cryo-EM have made it possible to learn latent variable models of conformation landscapes. However, interpreting these latent spaces remains a challenge as their individual dimensions are often arbitrary. The key message of our work is that this interpretation challenge can be viewed as an Independent Component Analysis (ICA) problem where we seek models that have the property of identifiability. That means, they have an essentially unique solution, representing a conformational latent space that separates the different degrees of freedom a molecule is equipped with in nature. Thus, we aim to advance the computational field of cryo-EM beyond visualizations as we connect it with the theoretical framework of (nonlinear) ICA and discuss the need for identifiable models, improved metrics, and benchmarks. Moving forward, we propose future directions for enhancing the disentanglement of latent spaces in cryo-EM, refining evaluation metrics and exploring techniques that leverage physics-based decoders of biomolecular systems. Moreover, we discuss how future technological developments in time-resolved single particle imaging may enable the application of nonlinear ICA models that can discover the true conformation changes of molecules in nature. The pursuit of interpretable conformational latent spaces will empower researchers to unravel complex biological processes and facilitate targeted interventions. This has significant implications for drug discovery and structural biology more broadly. More generally, latent variable models are deployed widely across many scientific disciplines. Thus, the argument we present in this work has much broader applications in AI for science if we want to move from impressive nonlinear neural network models to mathematically grounded methods that can help us learn something new about nature.

In Silico Exploration of Potential Phytoconstituents from the Bark Extract of Boswellia serrata for Hemorrhoidal Disease: Molecular Docking and Molecular Dynamics Analysis

Boswellia serrata Roxb. Ex Colebr is a popular medicinal plant used traditionally in herbal medicinal preparations to treat a variety of diseases. The purpose of the present investigation was to investigate the anti-hemorrhoidal property of the bark extract of B. serrata (BS). For this, the sequential Soxhlet extraction method was carried out by using different solvents such as hexane, chloroform, and methanol. After the extraction, the obtained dry extracts were tested for quantitative determinations such as total alkaloid content (TAC), total flavonoid content (TFC), total phenol content (TPC), and total tannin content (TTC) for all the extracts. Moreover, in vitro antioxidant activity was measured using 2,2-diphenyl-1-picrylhydrazyl (DPPH) free radical scavenging activity and scavenging activity against 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid (ABTS). Methanolic bark extract showed the highest TPC (67.10±1.83), TFC (372.73±4.45), TAC (9.732±1.06), and TTC (48.932±1.82), as well as the antioxidant assays DPPH (IC50=9.88 μg/ml) and ABTS (IC50=15.09 μg/ml). In this study, both LC-MS and GC-MS were performed to identify the chemical composition of all the extracts. Consequently, 19 compounds were identified by GC-MS and 27 compounds were identified by LC-MS analysis. The identified phytoconstituent(s) that could potentially interact with the target protein cyclooxygenase-2 (COX-2) (PDB: 4RRW) using molecular dynamics simulation and in silico docking were studied. Three compounds that have passed in drug-likeness and ADME-Tox properties are having more docking score than the standard. In this study, camptothecin, justicidin B, and taxiphyllin are identified as potential lead compounds with anti-hemorrhoidal properties and may be helpful in the process of drug development and discovery of novel drugs. Hence, these results demonstrate that BS is a good source of pharmacologically active components with potential applications against hemorrhoidal disease.

Predicting virus Fitness: Towards a structure-based computational model

Predicting the impact of new emerging virus mutations is of major interest in surveillance and for understanding the evolutionary forces of the pathogens. The SARS-CoV-2 surface spike-protein (S-protein) binds to human ACE2 receptors as a critical step in host cell infection. At the same time, S-protein binding to human antibodies neutralizes the virus and prevents interaction with ACE2. Here we combine these two binding properties in a simple virus fitness model, using structure-based computation of all possible mutation effects averaged over 10 ACE2 complexes and 10 antibody complexes of the S-protein (∼380,000 computed mutations), and validated the approach against diverse experimental binding/escape data of ACE2 and antibodies. The ACE2-antibody selectivity change caused by mutation (i.e., the differential change in binding to ACE2 vs. immunity-inducing antibodies) is proposed to be a key metric of fitness model, enabling systematic error cancelation when evaluated. In this model, new mutations become fixated if they increase the selective binding to ACE2 relative to circulating antibodies, assuming that both are present in the host in a competitive binding situation. We use this model to categorize viral mutations that may best reach ACE2 before being captured by antibodies. Our model may aid the understanding of variant-specific vaccines and molecular mechanisms of viral evolution in the context of a human host.

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.

Modelling SARS-CoV-2 spike-protein mutation effects on ACE2 binding

The binding affinity of the SARS-CoV-2 spike (S)-protein to the human membrane protein ACE2 is critical for virus function. Computational structure-based screening of new S-protein mutations for ACE2 binding lends promise to rationalize virus function directly from protein structure and ideally aid early detection of potentially concerning variants. We used a computational protocol based on cryo-electron microscopy structures of the S-protein to estimate the change in ACE2-affinity due to S-protein mutation (ΔΔG_bind) in good trend agreement with experimental ACE2 affinities. We then expanded predictions to all possible S-protein mutations in 21 different S-protein-ACE2 complexes (400,000 ΔΔG_bind data points in total), using mutation group comparisons to reduce systematic errors. The results suggest that mutations that have arisen in major variants as a group maintain ACE2 affinity significantly more than random mutations in the total protein, at the interface, and at evolvable sites. Omicron mutations as a group had a modest change in binding affinity compared to mutations in other major variants. The single-mutation effects seem consistent with ACE2 binding being optimized and maintained in omicron, despite increased importance of other selection pressures (antigenic drift), however, epistasis, glycosylation and in vivo conditions will modulate these effects. Computational prediction of SARS-CoV-2 evolution remains far from achieved, but the feasibility of large-scale computation is substantially aided by using many structures and mutation groups rather than single mutation effects, which are very uncertain. Our results demonstrate substantial challenges but indicate ways forward to improve the quality of computer models for assessing SARS-CoV-2 mutation effects.

Stability and expression of SARS-CoV-2 spike-protein mutations

Protein fold stability likely plays a role in SARS-CoV-2 S-protein evolution, together with ACE2 binding and antibody evasion. While few thermodynamic stability data are available for S-protein mutants, many systematic experimental data exist for their expression. In this paper, we explore whether such expression levels relate to the thermodynamic stability of the mutants. We studied mutation-induced SARS-CoV-2 S-protein fold stability, as computed by three very distinct methods and eight different protein structures to account for method- and structure-dependencies. For all methods and structures used (24 comparisons), computed stability changes correlate significantly (99% confidence level) with experimental yeast expression from the literature, such that higher expression is associated with relatively higher fold stability. Also significant, albeit weaker, correlations were seen between stability and ACE2 binding effects. The effect of thermodynamic fold stability may be direct or a correlate of amino acid or site properties, notably the solvent exposure of the site. Correlation between computed stability and experimental expression and ACE2 binding suggests that functional properties of the SARS-CoV-2 S-protein mutant space are largely determined by a few simple features, due to underlying correlations. Our study lends promise to the development of computational tools that may ideally aid in understanding and predicting SARS-CoV-2 S-protein evolution.

Structural heterogeneity and precision of implications drawn from cryo-electron microscopy structures: SARS-CoV-2 spike-protein mutations as a test case

Protein structures may be used to draw functional implications at the residue level, but how sensitive are these implications to the exact structure used? Calculation of the effects of SARS-CoV-2 S-protein mutations based on experimental cryo-electron microscopy structures have been abundant during the pandemic. To understand the precision of such estimates, we studied three distinct methods to estimate stability changes for all possible mutations in 23 different S-protein structures (3.69 million ΔΔG values in total) and explored how random and systematic errors can be remedied by structure-averaged mutation group comparisons. We show that computational estimates have low precision, due to method and structure heterogeneity making results for single mutations uninformative. However, structure-averaged differences in mean effects for groups of substitutions can yield significant results. Illustrating this protocol, functionally important natural mutations, despite individual variations, average to a smaller stability impact compared to other possible mutations, independent of conformational state (open, closed). In summary, we document substantial issues with precision in structure-based protein modeling and recommend sensitivity tests to quantify these effects, but also suggest partial solutions to the problem in the form of structure-averaged “ensemble” estimates for groups of residues when multiple structures are available.

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.

Construction of the small intestine on molecular dynamics simulation and preliminary exploration of drug intestinal absorption prediction

In this study, molecular dynamics simulation was applied to the construction of the small intestinal epithelial cell membrane and prediction of drug absorption. First, we constructed a system of a small intestinal epithelial cell membrane that was close to the real proportion and investigated the effects of temperature, water layer thickness, and ionic strength on membrane properties to optimize environmental parameters. Next, three drugs with different absorptivity, including Ephedrine (EPH), Quercetin (QUE), and Baicalin (BAI), were selected as model drugs to study the ability of drugs through the membrane by the free diffusion and umbrella sampling simulation, and the drug permeation ability was characterized by the free diffusion coefficient D and free energy barrier (△G) in the processes. The results showed that the free diffusion coefficient D and △G orders of the three drugs were consistent with the classical experimental absorption order, indicating that these two parameters could be used to jointly characterize the membrane permeability of the drugs.

Microsecond melting and revitrification of cryo samples: protein structure and beam-induced motion

A novel approach to time-resolved cryo-electron microscopy (cryo-EM) has recently been introduced that involves melting a cryo sample with a laser beam to allow protein dynamics to briefly occur in the liquid, before trapping the particles in their transient configurations by rapidly revitrifying the sample. With a time resolution of just a few microseconds, this approach is notably fast enough to study the domain motions that are typically associated with the activity of proteins but which have previously remained inaccessible. Here, crucial details are added to the characterization of the method. It is shown that single-particle reconstructions of apoferritin and Cowpea chlorotic mottle virus from revitrified samples are indistinguishable from those from conventional samples, demonstrating that melting and revitrification leaves the particles intact and that they do not undergo structural changes within the spatial resolution afforded by the instrument. How rapid revitrification affects the properties of the ice is also characterized, showing that revitrified samples exhibit comparable amounts of beam-induced motion. The results pave the way for microsecond time-resolved studies of the conformational dynamics of proteins and open up new avenues to study the vitrification process and to address beam-induced specimen movement.

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

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

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

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

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

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

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

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

Effects of cryo-EM cooling on structural ensembles

Structure determination by cryo electron microscopy (cryo-EM) provides information on structural heterogeneity and ensembles at atomic resolution. To obtain cryo-EM images of macromolecules, the samples are first rapidly cooled down to cryogenic temperatures. To what extent the structural ensemble is perturbed during cooling is currently unknown. Here, to quantify the effects of cooling, we combined continuum model calculations of the temperature drop, molecular dynamics simulations of a ribosome complex before and during cooling with kinetic models. Our results suggest that three effects markedly contribute to the narrowing of the structural ensembles: thermal contraction, reduced thermal motion within local potential wells, and the equilibration into lower free-energy conformations by overcoming separating free-energy barriers. During cooling, barrier heights below 10 kJ/mol were found to be overcome, which is expected to reduce B-factors in ensembles imaged by cryo-EM. Our approach now enables the quantification of the heterogeneity of room-temperature ensembles from cryo-EM structures.

The rapid temperature drop during plunge-freezing affects the structural ensembles obtained by cryo-EM. To quantify the extent of perturbation, Bock and Grubmüller combined continuum calculations, MD simulations, and kinetic models.

Structure and Mutations of SARS-CoV-2 Spike Protein: A Focused Overview

The spike protein (S-protein) of SARS-CoV-2, the protein that enables the virus to infect human cells, is the basis for many vaccines and a hotspot of concerning virus evolution. Here, we discuss the outstanding progress in structural characterization of the S-protein and how these structures facilitate analysis of virus function and evolution. We emphasize the differences in reported structures and that analysis of structure-function relationships is sensitive to the structure used. We show that the average residue solvent exposure in nearly complete structures is a good descriptor of open vs closed conformation states. Because of structural heterogeneity of functionally important surface-exposed residues, we recommend using averages of a group of high-quality protein structures rather than a single structure before reaching conclusions on specific structure-function relationships. To illustrate these points, we analyze some significant chemical tendencies of prominent S-protein mutations in the context of the available structures. In the discussion of new variants, we emphasize the selectivity of binding to ACE2 vs prominent antibodies rather than simply the antibody escape or ACE2 affinity separately. We note that larger chemical changes, in particular increased electrostatic charge or side-chain volume of exposed surface residues, are recurring in mutations of concern, plausibly related to adaptation to the negative surface potential of human ACE2. We also find indications that the fixated mutations of the S-protein in the main variants are less destabilizing than would be expected on average, possibly pointing toward a selection pressure on the S-protein. The richness of available structures for all of these situations provides an enormously valuable basis for future research into these structure-function relationships.

Energetics and J-coupling constants for Ala, Gly, and Val peptides demonstrated using ABEEM polarizable force field in vacuo and an aqueous solution

The development of an atom-bond electronegativity equalisation method at the σπ-level (ABEEM) polarisable force field (PFF) for peptides is presented. ABEEM PFF utilises a fluctuating charge model to explicitly describe...

Design of Transmembrane Mimetic Structural Probes to Trap Different Stages of γ-Secretase-Substrate Interaction

The transmembrane domain (TMD) of the amyloid precursor protein of Alzheimer's disease is cut processively by γ-secretase through endoproteolysis and tricarboxypeptidase "trimming". We recently developed a prototype substrate TMD mimetic for structural analysis-composed of a helical peptide inhibitor linked to a transition-state analogue-that simultaneously engages a substrate exosite and the active site and is pre-organized to trap the carboxypeptidase transition state. Here, we developed variants of this prototype designed to allow visualization of transition states for endoproteolysis, TMD helix unwinding, and lateral gating of the substrate, identifying potent inhibitors for each class. These TMD mimetics exhibited non-competitive inhibition and occupy both the exosite and the active site, as demonstrated by inhibitor cross-competition experiments and photoaffinity probe binding assays. The new probes should be important structural tools for trapping different stages of substrate recognition and processing via ongoing cryo-electron microscopy with γ-secretase, ultimately aiding rational drug design.

Chemical stabilization of γ-polyglutamate by chitosan and the effect of co-solvents on the stability

In protein-based formulations, conformational distortions and attractive interactions may cause insoluble and undesired aggregates. In the case of ionic peptides, including cationic or anionic, commonly electrostatic interactions are the main factors that control structure assembling. In this study, it was proposed that grafting of chitosan (CS) to γ-polyglutamic acid (γ-PGA) might exhibit much strong inhibiting effect on the formation of protein aggregates due to multiple amino groups and hydrophilic properties. To guarantee stable and safe biopharmaceutical formulation, the potency of a variety of stabilizers including sugars (glucose, sucrose), polyols (sorbitol, glycerol), surfactant (Tween 20), salting-out salt (PBS), and also different pH values have been evaluated on stabilizing or destabilizing the native state of CS-g-PGA copolymer using FTIR, CD, DLS, and SDS-PAGE. The comparable analysis revealed that the stability of CS-g-PGA was strongly dependent on pH owing to the polyelectrolyte characteristics of the polymers. Altogether these results implied that CS at optimized conditions might be an important precursor for the pharmaceutical industry and function as a new polymer for aggregation suppression and protein stabilization.

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.

Behavior of Hydrated Lipid Bilayers at Cryogenic Temperatures

Neutron diffraction was used to study the behavior of water present in phospholipid multilamellar stacks from 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) at cryogenic temperatures. Evidence was found for the existence of a highly viscous phase of water that exists between 180 and 220 K based on the observation that water can leave the intermembrane space at these low temperatures. Similar measurements are described in the literature for purple membrane (PM) samples. From a comparison with results from this natural membrane by using the same flash-cooling protocol, it is found that in the case of pure lipid samples, less water is trapped and the water flows out at lower temperatures. This suggests that the water is less hindered in its movements than in the PM case. It is shown that at least the Lβ′-phase of DMPC can be trapped likely by flash cooling; upon heating to about 260 K, it transforms to another phase that was not fully characterized.

Identification of Structural Calcium Binding Sites in Membrane-Bound Presenilin 1 and 2

Variants of presenilin (PS1 and PS2) are the main genetic risk factors of familial Alzheimer's disease and thus central to the disease etiology. Although mostly studied as catalytic units of γ-secretase controlling Aβ production, presenilins also affect calcium levels, which are disturbed in Alzheimer's disease. We investigated the interaction of calcium with both PS1 and PS2 using all-atom molecular dynamics (MD) simulations in realistic membrane models, with the specific aim to identify any Ca²⁺ sites. We did not observe any complete Ca²⁺ leak event, but we identified four persistent Ca²⁺ sites in membrane-bound PS1 and PS2: One in HL2 near the C-terminal of TM6, one in HL2 toward the N-terminal of TM7, a site at the catalytic aspartate on TM7, and a site at the PALP motif on TM9. The sites feature negatively charged glutamates and aspartates typical of calcium binding. Structural homology to diaspartate calcium transport proteins and mutation studies of calcium efflux support our identified calcium sites. Calcium consistently dampens HL2 motions in all comparisons (PS1, protonated PS1, PS2, protonated PS2). Due to their location in HL2 and the active site, we propose that the calcium sites control autoproteolytic maturation of presenilin by a pH-dependent conformational restriction of the HL2 recognition loop, which also regulates calcium transport proteins such as inositol 1,4,5-triphosphate receptor and ryanodine receptor. Our structural dynamics could provide a possible molecular basis for the need of both calcium and presenilin for lysosome proteolytic function, perhaps relevant also to other protein misfolding diseases.