Next Article in Journal
Interactions between Host Immunity and Skin-Colonizing Staphylococci: No Two Siblings Are Alike
Next Article in Special Issue
Macroglossia, Motor Neuron Pathology, and Airway Malacia Contribute to Respiratory Insufficiency in Pompe Disease: A Commentary on Molecular Pathways and Respiratory Involvement in Lysosomal Storage Diseases
Previous Article in Journal
In Vitro Activity of the Bacteriophage Endolysin HY-133 against Staphylococcus aureus Small-Colony Variants and Their Corresponding Wild Types
Previous Article in Special Issue
Molecular Pathways and Respiratory Involvement in Lysosomal Storage Diseases
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 3
10.3390/ijms20030717
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lipid–Protein Interactions in Niemann–Pick Type C Disease: Insights from Molecular Modeling

by
Simon Wheeler
1,
Ralf Schmid
2,3 and
Dan J Sillence
1,*
1
School of Pharmacy, De Montfort University, The Gateway, Leicester LE1 9BH, UK
2
Leicester Institute of Structural and Chemical Biology, Henry Wellcome Building, University of Leicester, Lancaster Road, Leicester LE1 7RH, UK
3
Department of Molecular and Cell Biology, Henry Wellcome Building, University of Leicester, Lancaster Road, Leicester LE1 7RH, UK
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(3), 717; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20030717
Submission received: 16 January 2019 / Revised: 31 January 2019 / Accepted: 3 February 2019 / Published: 7 February 2019
(This article belongs to the Special Issue Molecular Features of Lysosomal Storage Disorders)
Browse Figures
Versions Notes

Abstract

:
The accumulation of lipids in the late endosomes and lysosomes of Niemann–Pick type C disease (NPCD) cells is a consequence of the dysfunction of one protein (usually NPC1) but induces dysfunction in many proteins. We used molecular docking to propose (a) that NPC1 exports not just cholesterol, but also sphingosine, (b) that the cholesterol sensitivity of big potassium channel (BK) can be traced to a previously unappreciated site on the channel’s voltage sensor, (c) that transient receptor potential mucolipin 1 (TRPML1) inhibition by sphingomyelin is likely an indirect effect, and (d) that phosphoinositides are responsible for both the mislocalization of annexin A2 (AnxA2) and a soluble NSF (N-ethylmaleimide Sensitive Fusion) protein attachment receptor (SNARE) recycling defect. These results are set in the context of existing knowledge of NPCD to sketch an account of the endolysosomal pathology key to this disease.

1. Introduction

Niemann–Pick type C disease (NPCD) is a devastating neurological disorder that mostly affects children. Its root cause is dysfunction in the cholesterol-binding proteins NPC1 or NPC2. Low-density lipoprotein (LDL)-derived cholesterol is delivered to the cell in the form of cholesteryl esters. These are carried by the endocytic pathway to the lysosomes where they are hydrolyzed by lysosomal acid lipase (LAL). Unesterified cholesterol (sometimes termed “free cholesterol”) is collected from intraluminal membranes by NPC2 which binds to loop 2 of NPC1 [1,2] and transfers the lipid to the N-terminal domain (NTD) of NPC1 [3,4] (see Figure S1, for the domains of NPC1). This domain then relays the cholesterol to the NPC1 sterol-sensing domain (SSD) [3,5,6] from where the lipid transfers to the endoplasmic reticulum (ER) via a partner protein [7]. The most up-to-date version of this account [2] proposes that cholesterol diffuses from the NPC1 SSD into the luminal side of the lysosomal membrane, spontaneously flips, and is then collected by other organelles. However, cholesterol flipping is rapid [8]; thus, it is unclear why NPC1 deletion would lead to cholesterol accumulation in the lysosome or a cholesterol deficit in the ER, which are in fact cellular hallmarks of NPCD [9,10]. Thus, we asked if there was a cholesterol-binding site on NPC1 aligned with the cytosolic face of the late endosomal/lysosomal (LEL) membrane. If such a site existed, it would explain more fully NPC1′s role in cholesterol export and, thus, the accumulation of cholesterol that results from NPC1′s deletion or mutation.
This accumulation of cholesterol triggers a secondary build-up of other fats including sphingolipids and phosphoinositides (Figure 1) [11]. This pathology prompted us to postulate that the accumulating lipids could bind to LEL proteins and lead to changes in their function. For example (Figure 1), cholesterol accumulation will likely affect the big potassium channel (BK, also known as Slo1, MaxiK and KCa1.1), which was recently found on lysosomes [12,13]. It is long recognized as being sensitive to its lipid environment, with cholesterol suggested as the key species [14,15,16,17], and there is evidence, in the form of structure–activity studies, to suggest a specific interaction rather than an influence on bilayer thickness [18,19]. Thus, we would expect cholesterol to bind to the BK protein. The lysosomal role of the channel is unclear, though it is probably significant that it acts to regulate Transient receptor potential mucolipin 1 (TRPML1), a Ca2+-efflux channel required for membrane fusion and thus endocytosis [12,20]. BK is composed of four identical subunits, each of which contains a voltage sensor and two regulation of conductance for potassium (RCK) domains. Each subunit also contributes two helices to the channel pore (Figure S2). A 2012 paper [21] traced cholesterol sensitivity to a motif (Val509–Lys518) positioned in RCK1; this sensitivity was lost when the protein was cleaved before Val509. This cleavage was quite drastic, removing all of RCK1 and -2, thus significantly weakening the link between Ca2+ binding and channel opening, and questioning the physiological relevance of these findings. More importantly still, in the light of more recent structural work on the protein [22,23], it appears that the region identified as binding cholesterol is present in the cytosol where levels of cholesterol will be lower than in the membrane. Given that the membrane region has four cholesterol-binding motifs (Figure S2), this seems questionable. Thus, we wished to revisit the question of where cholesterol binds to BK.
The accumulation of sphingomyelin (SM) in the LELs of NPCD cells is triggered by increased levels of cholesterol. SM was reported to inhibit Ca2+ efflux through the TRPML1 channel [24] (Figure 1), something that is necessary for the correct completion of endocytosis [20]. SM accumulation, therefore, will contribute to endocytic failure; thus, we wished to identify where SM binds to TRPML1. A recent structural study [25] revealed that increased endolysosomal pH also compromises TRPML1 function. This channel exists as a homotetramer whose assembly is driven by interactions between the loop 1 regions on adjacent subunits (Figure S3). An aspartate-rich region on each loop, thus, forms a pore which confers pH-sensitivity on the channel. At low pH, the aspartate residues are all protonated; thus, Ca2+ ions can pass freely through. As pH rises, the aspartates become deprotonated, and the resulting anionic carboxylate moieties bind Ca2+, thereby inhibiting its own conductance. In agreement with previous reports [26,27], we have unpublished data [28] that the pH of the LEL compartment in NPCD is increased, and that the increase stems from defective activation of the vATPase proton pump by glucosylceramide (Figure 1). This is another reason for expecting reduced TRPML1 function in NPCD.
Whilst the assembly of cholesterol and sphingolipids in ordered domains termed lipid rafts is widely appreciated, it is less well known that cholesterol also mediates clustering of phosphoinositides (Figure 1, reviewed in References [29,30]). Therefore, it seems likely that cholesterol accumulation in NPCD will lead to trapping of phosphoinositides in LEL vesicles; indeed, phosphoinositides do increase in the LEL compartment of NPCD cells [11]. These negatively charged lipids are likely to bind to particular proteins; thus, the accumulation of phosphoinositides may lead to the accumulation of such proteins and, thus, impaired protein function.
It is known for some time that the peripheral membrane protein annexin A2 (AnxA2) mislocalizes in NPCD, associating with late, not early, endosomes [31,32]. Multiple studies reported the association of AnxA2 with lipid rafts [33,34,35], although phosphatidylinositol-4,5-bisphosphate (PI(4,5)P2) was also reported as the key binding lipid [36,37]. The clustering of phosphoinositides by cholesterol offers a way of unifying these disparate observations. Therefore, we asked if phosphoinositides could bind to AnxA2 (Figure 1).
Endocytosis and autophagy, both of which are dysfunctional in NPCD [38,39], rely on membrane fusion events, which in turn are dependent on soluble NSF (N-ethylmaleimide Sensitive Fusion) protein attachment receptor (SNARE) proteins such as syntaxins and vesicle-associated membrane proteins (VAMPs). The subsequent re-use of these proteins, therefore, requires their release from the assembled SNARE complex, and this process is retarded in cholesterol-rich membranes [40]. Significant work was done on the part of syntaxin 1 (Stx1) adjacent to the membrane (the juxta-membrane region) which is polybasic and thus polycationic in vivo [41,42,43,44]. Many other syntaxins have a similar region, as do their VAMP partners (see multi-sequence alignments, Figure S4). Thus, we wished to examine the question of whether phosphoinositides, not cholesterol, are directly responsible for retarded SNARE recycling (Figure 1).
To examine lipid–protein interactions, it is necessary to identify the region of the protein where the lipid associates or binds. A protein’s sequence can be searched for established lipid-binding motifs using computer algorithms (e.g., FuzzPro, bioinformatics.nl/cgi-bin/emboss/fuzzpro), but the predictive power of this technique is variable. For instance, the most common cholesterol-binding motif, termed CRAC [45,46] (cholesterol recognition amino-acid consensus) is defined as L/V-X1-5-Y-X1-5-R/K and, thus, is of variable length and strictly defined at only one position. Thus, algorithmic searching for this motif typically generates numerous false positives for any given protein sequence. The same is true of the reverse motif CARC [47], which also binds cholesterol. In contrast, a sphingolipid-binding pattern [48] was postulated as a post hoc rationalization of experimental observations in a small number of proteins, though this was not extrapolated to a general motif. More recently, the protein p24 was found to bind to sphingomyelin via a defined region [49]. Based on this finding, a set of sphingolipid-binding motifs, e.g., VX2V2X2LF, was elucidated [50]. However, to our knowledge, no interactions between sphingolipids and specific residues within these motifs were identified, nor is it clear where the required terminal aromatic residue of the motif should be positioned relative to the membrane in which the protein is embedded. Thus, algorithmic searching of sequences for binding motifs offers only an approximate guide to where lipid–protein interactions may occur, and gives only limited information of three-dimensional (3D) features of such interactions. We set out to use molecular docking, which not only allows the exploration and rationalization of lipid binding in structural detail, but also allows ruling out some sites if binding proves impossible or implausible. This technique is, hence, much more reliable than merely identifying motifs.
Such docking requires a structural model of the protein in question, which can be gained by experiment (X-ray, cryoEM, or NMR) or derived from such a structure (by homology modeling or threading). Lipid binding to specific sties of the protein can then be estimated with the protein mostly rigid, but giving flexibility to some designated side chains, while the lipid is allowed to vary its position, orientation, and bond angles. We started our calculations with the bond angles in an energy-minimized state; the docking software randomly selects a starting position and orientation. All these parameters are then systematically altered until an energetically favorable binding “pose”—assessed by an energy scoring function—is discovered, if possible. This procedure is repeated multiple times. and the results are ranked to show more or less energetically favorable binding poses. These poses are grouped into clusters defined by root-mean-square (RMS) difference of atomic positions with the cut-off for membership of a cluster set at >2 Å. The calculations examine only a sample of all possible combinations of parameters and are limited by the accuracy of the scoring function. Thus, the results are hypotheses that need to be inspected and tested—a binding pose with a favorable energy score might be biologically impossible.
With membrane-bound and -associated proteins, such as those considered here, the orientation of the protein in the membrane (assessed here using the Orientation of Proteins in Membranes (OPM) database, opm.phar.umich.edu/ [51]) offers a significant clue as to the biological plausibility of binding poses. Thus, a lipid, illustrated by cholesterol (orange) in Figure 2, may bind to a trans-membrane protein (light grey) parallel to the lipophilic section of the membrane and with the hydroxy group aligned with other lipid head groups (Figure 2A). This is the most energetically favorable situation. However, the lipid may also be parallel to the membrane, but vertically displaced from it (Figure 2B). Whilst possible [52], this requires a deformation in the membrane which attracts an energy penalty; such a pose would, therefore, be energetically less favored. Burying the hydroxy group inside the lipophilic portion of the membrane by orienting it in a perpendicular (Figure 2C) or antiparallel (Figure 2D) fashion is still possible, but energetically disfavored. In this work, we considered binding poses similar to those in Figure 2B–D only when calculations did not generate poses corresponding to the more favorable situation shown in Figure 2A.
We made our calculations more computationally tractable by neglecting the membrane; thus, we assessed only lipid–protein interactions. Lipid–protein interactions considered in isolation may be energetically favorable with the lipid in perpendicular or antiparallel situations. As the lipid starting position and orientation are randomized, some of the iterations of the docking calculations will start with the lipid in a position perpendicular or antiparallel to the lipophilic section of the membrane. Thus, some docking runs will start with the lipid in an orientation that is energetically favorable when considered in isolation but energetically disfavorable when considered in the context of the membrane as a whole. Such an orientation is biologically implausible. Energetic refinement of this orientation will not necessarily result in the lipid adopting a parallel orientation. Therefore, some docking runs may result in binding poses that appear energetically favorable because they only consider lipid–protein interactions but would be energetically disfavorable overall and, thus, biologically implausible. This is a consequence of the randomization of the initial position of the lipid, and not a flaw in the docking experiments.
Whilst this discussion applies mostly to cholesterol, similar considerations pertain to other lipids. Aliphatic lipid tails generally make non-specific interactions with trans-membrane proteins; thus, we routinely perform docking calculations with only the lipid head group to focus attention on the distinctive part of the molecule. These head groups include the first carbon atoms of the aliphatic chains, and it is quite possible that the randomized initial positioning for the docking calculation will orient these atoms and, thus, the alkyl chains, away from the membrane and toward the cytosol. It is also possible that the positioning will point the alkyl chains in opposite directions, one toward the membrane and one toward solvent (referred to as a splayed arrangement in Table S1). It is equally possible that the energy minimization, because done in isolation from the membrane and solvent, will retain either of these orientations if adopted by the initial random positioning. Thus, once again, docking calculations may output binding poses that are energetically favorable in isolation but disfavorable overall and, thus, biologically implausible. This does not indicate a flawed method.
For this reason, it is necessary that the output of such docking be firstly inspected for biological plausibility. Our figures show representative binding poses that result from the docking runs—these are both energetically favorable and biologically plausible. Table S1 presents our extended data including number of clusters, cluster size, assessment of favorability, and interacting residues. Once filtered for energetic favorability and biological plausibility, the results of docking experiments should next be tested by “wet” experiments. Whilst we present no new lab data here, our results are fully consistent with existing findings from other labs which are cited extensively in the text. The main proposal of this work is a set of lipid–protein interactions that explain and connect disparate data from the literature and are testable in vitro. We propose that errors in lipid distribution in NPCD are the proximate cause of protein dysfunction, and that this failure leads eventually to defects at the cellular level. We depict this in graphical form in Figure 1.

2. Results

2.1. What Does the NPC System Do?

As discussed above, the current view of cholesterol export via NPC1 [2] does not fully explain the lipid misdistribution resulting from mutation or deletion of this protein. This cholesterol export narrative can, however, be made consistent with current knowledge if it is refined to include a binding site on the cytosolic face of the LEL membrane which delivers cholesterol to the ER. Thus, we performed cholesterol docking to the entire SSD of NPC1 using the recently published X-ray structure (Protein Data Bank (PDB) 5U74 [6]). We found lipid-binding poses in two energetically favorable clusters aligned to both the luminal and cytosolic faces of the membrane. Representative examples are shown in Figure 3A; extended results are given in Table S1. In both cases, binding is anchored by an H-bond between the cholesterol hydroxy group and the side chain of an Asp residue, whilst the remainder of the sterol is encased in a cavity formed from lipophilic side chains (Figure 3A; Table S1). These side chains make van der Waals contact with the cholesterol molecule. Crucially, Pro691 is among the residues predicted to interact with cholesterol (pink in Figure 3). Mutation of this residue was shown to reduce cholesterol export [53] and to reduce labeling of NPC1 by a photoactivatable cholesterol analog [54], implying that Pro691 plays an important role in cholesterol binding. The emergence of this residue from docking calculations forms an important agreement between in vitro and in silico approaches. The suggestion of two binding pockets means that NPC1 can potentially be added to the list of proteins featuring cholesterol binding in both leaflets of a lipid bilayer [55,56], a phenomenon termed the “mirror code” [57]. The path between these two binding sites remains to be uncovered, but we propose that NPC1 interacts with the ER to deliver cholesterol from the flipped pocket. This modified narrative takes account of the rapid flipping of cholesterol [8], and the cholesterol deficit at the ER [10] in NPCD, and is, thus, more consistent with available data than that proposed previously [2].
Experiments with lipid profiling and photoactivatable probes demonstrated increased levels of sphingosine (Sph) in NPCD, and that this accumulation can be traced to the LEL [58,59,60], supporting the idea that NPC1 can act as a sphingosine export protein. To investigate this further, we performed our docking experiments with Sph and the SSD of NPC1. Consistent with the studies just cited, we found plausible binding poses for sphingosine aligned to both sides of the LEL membrane. Reminiscent of cholesterol, the polar head group of Sph forms H-bonds or ionic interactions with the side chains of Asp residues, while the lipophilic tail makes van der Waals contact with the same set of lipophilic side chains as cholesterol (Figure 3B; Table S1). We also docked Sph to NPC2 and the NTD of NPC1 (Figure S5). Our calculations gave energetically favorable binding of Sph with the polar head group oriented the same way as cholesterol. This is known to be critical for the correct functioning of the cholesterol export machinery [3,5] and, thus, further supports the idea that Sph may use the same pathway. The consequences of impaired Sph export for LEL function are unclear, though it is perhaps significant that this lipid was recently reported as an inhibitor of glucosylceramidase 2 [61], a target of miglustat [62], which is the only clinically used treatment of NPCD. Alternatively, sphingosine (pKa 6.6 [63]) is likely to be protonated in the endolysosome (pH 4–5) but deprotonated in the cytosol (pH 7.4) and, thus, sphingosine export may act to move H+ ions from one to the other. NPC1 deficiency will, therefore, lead to changes in pH in both the LEL and cytosol. This hypothesis could be tested by looking for changes in pH homeostasis in both compartments.

2.2. Revising the BK–Cholesterol Interaction Site

As stated in the introduction, we wished to reconsider cholesterol binding to BK in the membrane region. To do this, we built a model of the membrane section of the human BK channel based on a recent cryoEM structure (PDB 5TJI) [22,23]. We attempted docking of cholesterol at each of each of the intra-membrane cholesterol binding motifs. For CARC, CRAC1, or CRAC3 motifs (defined as in [21]), no poses could be found that were both energetically favorable and biologically plausible (assessed as described in the introduction). We then attempted to dock cholesterol in the space near the CRAC2 cholesterol-binding motif (Val258–Arg266, located on the voltage sensor), allowing various combinations of side chains to be flexible. Using this procedure, we found a number of plausible binding poses (representative examples are shown in Figure 4; extended results are given in Table S1). The cholesterol hydroxy group can form H-bonds with Lys211 and/or Asn265 (in some poses, these residues also H-bonded to each other, data not shown). The lipophilic portion of cholesterol makes van der Waals contacts with various lipophilic side chains (Table S1), in particular Leu226, which has favorable interactions with the terminus of cholesterol. In some poses, the side chain of Tyr263 makes CH–π interactions with the smooth α-face of cholesterol (Figure 4B). Thus, these poses feature the classical elements of a CRAC motif (basic, aromatic, and lipophilic residues) [46], but are not derived from a contiguous section of protein sequence, something previously observed in the β2 adrenoceptor [64]. Whilst our results cannot disprove cholesterol binding in the cytosolic domains, we believe that cholesterol binding in the membrane region is more likely. Our docking is fully consistent with previous studies demonstrating BK’s cholesterol sensitivity [14,15,16,17,18]. In particular, both of the representative poses shown in Figure 4 rely on the correct cholesterol stereochemistry, which was shown to be key to its effects on the BK channel [18,19]. Our hypothesis could be tested by mutagenesis experiments involving key hydrogen-bonding residues Lys211 and Asn265.
The implications of cholesterol-induced BK dysfunction for NPCD are multiple (Figure 1). Firstly, this channel as reported as a regulator of LEL–ER membrane contact sites (MCSs) and, thus, to be important for refilling lysosomal Ca2+ stores [13] which are reported deficient in NPCD [58,65]. It is also significant that blocking lysosomal Ca2+ refilling gave a cellelar phenotype resembling that of lysosomal storage disorders LSDs [66], a class of which NPCD is a member. BK was also reported as mutually regulating lysosomal Ca2+ channel TRPML1 [12] (see next section). Whatever the details, the importance of BK for NPCD is attested by the repair of the endocytic defect by either channel agonism [67] or over-expression [12].

2.3. LEL Ca2+ and TRPML1

Impaired formation of LEL–ER contact sites is likely to result in reduced concentration of lysosomal Ca2+, which was indeed found in NPCD [58,65], though not without some disagreement [24]. Reduced LEL Ca2+ will in turn result in impaired calcium release via the TRPML1 channel (Figure 1), something that is necessary for the correct completion of endocytosis [20]. As discussed above, sphingomyelin (SM) accumulates in the LELs of NPCD cells and is reported to inhibit conductance through this channel [24]. There are at least three possibilities for the site of this interaction: (1) SM binds at the same site as synthetic ligand ML-SA1, (2) SM competes with the natural ligand PI(3,5)P2, or (3) SM inhibits assembly of the four proteins into the channel (Figure S3A). The first of these is precluded, as the ML-SA1 binding site lies in the intra-membrane region [68,69] and will, thus, be inaccessible to the SM head group. The second was investigated by docking the relevant lipid head groups to a region of the human protein (PDB 5WJ9, [68]) recently identified in a marmoset TRPML3 channel [70] (Figure S3C). Use of just the head groups focused attention on the distinctive part of the lipid (hydrophobic tails are common to most classes), as well as rendering the calculation more tractable. Both PI(3,5)P2 and natural inhibitor PI(4,5)P2 were found to bind to this polybasic region, thus supporting mutagenesis experiments [70] locating binding here (poses not shown, extended data in Table S1). SM was also successfully docked to this region (poses not shown; extended data in Table S1) but the energy of interaction was much weaker, suggesting that SM will not compete with phosphoinositides for this binding site. Inhibiting channel assembly (the third possibility) was previously shown to lead to reduced TRPML1 function [25] and two sites of inter-chain binding were identified in structural studies: Arg146–Val175 (H-bonded) in the linker region [25] and Arg486–Glu276 (zwitterionic) in the juxta-membrane region [68] (Figure S3A). The Arg146–Val175 site is too far from the membrane to be relevant (Figure S3). We carried out docking of the SM head group to the Arg486–Glu276 region. Very few energetically favorable binding poses positioned the SM head group in a way that would align correctly with the LEL membrane. Those that did feature interactions between the lipid and the protein, but also the protein side chains, preserve interactions between them rather than become separated by SM. Channel assembly is, thus, maintained rather than disrupted in the presence of this polar lipid (Figure 5C,D). Considering this, direct effects of SM on the TRPML1 channel appear unlikely and we believe that the inhibitory properties of this lipid may well derive from indirect effects, perhaps mediated through membrane organization.
TRPML1 plays roles in endocytosis [20,71] and autophagy [72]; thus, the failure of these processes in NPCD [38,39,73,74] is unsurprising (Figure 1), although there may be some redundancy with the two-pore channels (TPCs) [75,76], P2X4 [77] and P/Q-type voltage-gated calcium channel [78,79].

2.4. Cholesterol Clusters Phosphoinositides—Implications for Annexin Localization

In the introduction, we raised the possibility that cholesterol-mediated clustering of phosphoinositides [29,30] might be invoked to explain the mislocalization of Annexin A2 (AnxA2) in NPCD [31,32]. AnxA2 is a peripheral membrane protein, suggesting that it is likely to associate with membranes via interactions between surface residues with side chains that extend from the protein and head groups that protrude from the membrane. Indeed, previous site-directed mutagenesis identified Lys279 and Lys281 as important for PI(4,5)P2 binding [36].
With this in mind, we attempted to dock AnxA2 (PDB 2HYW [80]) with PI(4,5)P2 using just the lipid head group for the reasons discussed above. The search was concentrated in a space including Lys279, Lys281, and a nearby basic residue Arg284. These residues are expected to be protonated and, thus, cationic in vivo. Consistent with the previously reported results [36,37], our docking experiments show these cationic residues making zwitterionic interactions with the two anionic phosphate groups from the lipids (example binding pose shown in Figure 6, ionic interactions in magenta; extended data shown in Table S1). The two lipid side chains point away from the protein and remain in the membrane. The start of the lipophilic tails is indicated in Figure 6 (one is obscured behind Lys281). Docking of SM was also performed, and the results were much less energetically favorable than for PI(4,5)P2 (Table S1) suggesting that the phosphoinositide is genuinely responsible for AnxA2 binding. Thus, in NPCD, AnxA2 will be expected to mislocalize away from early endosomes and to cholesterol-rich late endosomes as observed [31,32] and contribute to the documented endocytic defect [38,73]

2.5. Cholesterol Clusters Phosphoinositides—Implications for SNARE Complex Disassembly

In the introduction, we noted that SNARE recycling from cholesterol-rich membranes is retarded [40]. It was postulated that the cholesterol mediates the clustering of phosphoinositides [29,30] and that phosphoinositides might bind to SNARE components and, thus, trap them in membranes. To test this idea, we built a model of a SNARE complex featuring Stx7 and VAMP8—a combination reported as necessary for endocytosis [81]—based on an X-ray structure (PDB 3HD7 [82]) of a related bundle. The juxta-membrane basic residues (VAMP8: Arg67, Lys68, Lys72; Syntaxin 7: Arg232, Lys233, Arg235), as with AnxA2, will be expected to be protonated and, thus, cationic. These extend from the protein in different directions, and it was, therefore, necessary to select appropriate combinations of these to be flexible in the docking experiments. When we attempted docking of PI(3,5)P2 to this protein complex, we found numerous plausible poses. As expected, the anionic phosphates on the lipid form zwitterionic interactions with the cationic residues on the protein (magenta, Figure 7A,B; extended data in Table S1). The first carbon atoms of the lipophilic tails are indicated and are oriented toward the membrane (Figure 7A,B). Importantly, the phosphates interact simultaneously with residues on both SNARE partners. This would be expected to promote SNARE interaction and, thus, inhibit complex disassembly. Our results are, therefore, consistent with both the cholesterol-mediated clustering of PI(3,5)P2 [29,30] and the recycling defect observed in cholesterol-rich membranes [40], and they add a previously unappreciated subtlety to dysfunctional endocytosis in NPCD. Repeated attempts to dock PI(4,5)P2 to our model consistently failed to give plausible binding poses, suggesting that clustering is dependent on only the 3,5-isomer. The situation with autophagy is more complicated, as this was reported to rely on Stx17 in both mammals and Drosophila [83,84]; Stx7 does not appear to be involved [85]. Stx17 is unusual in its possession of two trans-membrane helices, neither of which has a polybasic juxta-membrane region (Figure S4).
The widespread effects of phosphoinositides on ion channels [86], not least as agonists of TRPML1 [87,88] and the TPCs [89], suggest that cholesterol-mediated clustering of these lipids will be a fruitful area for future research in lysosomal storage disorders.

3. Discussion

We argued here that the accumulation of cholesterol in the LEL compartment results from dysfunctional NPC1 and affects some proteins directly (e.g., BK). At the same time, it triggers the secondary build-up of other lipids which affect other proteins (e.g., AnxA2). Thus, the misdistribution of lipids in the LEL membrane results in widespread protein dysfunction. Ultimately, these defects manifest as errors at the cellular level (Figure 1). The account presented here concerns proteins of the LEL membrane and, even on that restricted basis, is a simplification as, to say nothing of other organelles and processes, there is evidence that AnxA6 [90], CLC-6 [26,79], mTOR [55], rab9 [91,92], and TPC1 [60] are also involved. The absence of high-quality structural data and/or inconclusive in vitro results meant a modeling approach was inappropriate for these proteins at this time.

4. Materials and Methods

Multi-sequence alignments were performed in Clustal Omega (ebi.ac.uk/Tools/msa/clustalo/) [93] and visualized in JalView 2.10.5 (http://www.jalview.org/) (version, manufacture, city, if any state, country) [94]. Lipid-binding motifs were located using Fuzzpro (bioinformatics.nl/cgi-bin/emboss/fuzzpro).
Protein structures were either downloaded from the PDB (5U74 for NPC1, 5WJ9 for TRPML1, 2HYW for AnxA2) or models were built using SwissModel (https://swissmodel.expasy.org/interactive) (version, manufacture, city, if any state, country) [95] (5TJI as a template for BK, 3HDY as a template for SNARE bundle). Quality was assessed using QMEANBrane (https://swissmodel.expasy.org/qmean/) (version, manufacture, city, if any state, country) [96,97]; for details, see the Figure S5. Approximate positions in the membrane were found using the OPM database (http://opm.phar.umich.edu/) [51]. The AnxA2 structure features calcium ions which by default are set to zero charge by AutoDock. Thus, atomic charges for these and spatially proximate atoms were calculated using the Atomic Charge Calculator (webchem.ncbr.muni.cz/Platform/ChargeCalculator) [98], and the relevant AutoDock files were manually edited. Ligand structures were prepared in Avogadro (avogadro.cc/) and minimized using the MMF94 force field with at least 5000 steps; other settings were defaults. The lipids considered are shown in Figure S6.
Docking of lipids to proteins was performed using AutoDock 4.2.6 (http://autodock.scripps.edu/) (version, manufacture, city, if any state, country) [99,100] using default settings. Search spaces, the residues allowed to be flexible, and the number of algorithm runs are given in Figure S7. AutoDock clusters binding poses by RMS distance (cut-off >2 Å). Docking scores were used as a preliminary assessment, followed by manual inspection for biological plausibility as discussed in the introduction.
Protein structures, including docking results, were visualized in UCSF Chimera (https://www.cgl.ucsf.edu/chimera/) (version, manufacture, city, if any state, country) [101].

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/20/3/717/s1.

Author Contributions

Conceptualization, S.W. and D.S.; investigation, S.W.; writing—original draft preparation, S.W.; writing—review and editing, S.W., R.S., and D.S.; supervision, R.S. and D.S.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

AnxAnnexin
BKBig potassium
CRACCholesterol recognition amino-acid consensus
EREndoplasmic reticulum
LDLLow-density lipoprotein
LELLate endolysosome
LSDLysosomal Storage Disorder
MCSMembrane contact site
mTORmechanistic Target of Rapamycin
NPCNiemann–Pick type C
NPCDNiemann–Pick type C disease
NSFN-ethylmaleimide Sensitive Fusion
NTDN-terminal domain
PI(3,5)P2Phosphatidylinositol-3,5-bisphosphate
PI(4,5)P2Phosphatidylinositol-4,5-bisphosphate
RCKRegulation of conductance by potassium
RMSRoot mean square
SMSphingomyelin
SNARESoluble NSF protein attachment receptor
SphSphingosine
SSDSterol-sensing domain
StxSyntaxin
TPCTwo-pore channel
TRPMLTransient receptor potential mucolipin
VAMPVesicle-associated membrane protein

References

  1. Deffieu, M.S.; Pfeffer, S.R. Niemann–Pick Type C 1 Function Requires Lumenal Domain Residues That Mediate Cholesterol-Dependent NPC2 Binding. Proc. Natl. Acad. Sci. USA 2011, 108, 18932–18936. [Google Scholar] [CrossRef] [PubMed]
  2. Li, X.; Saha, P.; Li, J.; Blobel, G.; Pfeffer, S.R. Clues to the Mechanism of Cholesterol Transfer from the Structure of NPC1 Middle Lumenal Domain Bound to NPC2. Proc. Natl. Acad. Sci. USA 2016, 113, 10079–10084. [Google Scholar] [CrossRef] [PubMed]
  3. Infante, R.E.; Wang, M.L.; Radhakrishnan, A.; Kwon, H.J.; Brown, M.S.; Goldstein, J.L. NPC2 Facilitates Bidirectional Transfer of Cholesterol between NPC1 and Lipid Bilayers, a Step in Cholesterol Egress from Lysosomes. Proc. Natl. Acad. Sci. USA 2008, 105, 15287–15292. [Google Scholar] [CrossRef] [PubMed]
  4. Wang, M.L.; Motamed, M.; Infante, R.E.; Abi-Mosleh, L.; Kwon, H.J.; Brown, M.S.; Goldstein, J.L. Identification of Surface Residues on Niemann-Pick C2 Essential for Hydrophobic Handoff of Cholesterol to NPC1 in Lysosomes. Cell Metab. 2010, 12, 166–173. [Google Scholar] [CrossRef] [PubMed]
  5. Kwon, H.J.; Abi-Mosleh, L.; Wang, M.L.; Deisenhofer, J.; Goldstein, J.L.; Brown, M.S.; Infante, R.E. Structure of N-Terminal Domain of NPC1 Reveals Distinct Subdomains for Binding and Transfer of Cholesterol. Cell 2009, 137, 1213–1224. [Google Scholar] [CrossRef] [PubMed]
  6. Li, X.; Lu, F.; Trinh, M.N.; Schmiege, P.; Seemann, J.; Wang, J.; Blobel, G. 3.3 Å Structure of Niemann–Pick C1 Protein Reveals Insights into the Function of the C-Terminal Luminal Domain in Cholesterol Transport. Proc. Natl. Acad. Sci. USA 2017, 114, 9116–9121. [Google Scholar] [CrossRef] [PubMed]
  7. Pfisterer, S.G.; Peränen, J.; Ikonen, E. LDL-cholesterol Transport to the Endoplasmic Reticulum: Current Concepts. Curr. Opin. Lipidol. 2016, 27. [Google Scholar] [CrossRef]
  8. Steck, T.L.; Lange, Y. How Slow Is the Transbilayer Diffusion (Flip-Flop) of Cholesterol? Biophys. J. 2012, 102, 945–946. [Google Scholar] [CrossRef] [PubMed]
  9. Abi-Mosleh, L.; Infante, R.E.; Radhakrishnan, A.; Goldstein, J.L.; Brown, M.S. Cyclodextrin Overcomes Deficient Lysosome-to-Endoplasmic Reticulum Transport of Cholesterol in Niemann-Pick Type C Cells. Proc. Natl. Acad. Sci. USA 2009, 106, 19316–19321. [Google Scholar] [CrossRef] [PubMed]
  10. Frolov, A.; Zielinski, S.E.; Crowley, J.R.; Dudley-Rucker, N.; Schaffer, J.E.; Ory, D.S. NPC1 and NPC2 Regulate Cellular Cholesterol Homeostasis through Generation of Low Density Lipoprotein Cholesterol-Derived Oxysterols. J. Biol. Chem. 2003, 278, 25517–25525. [Google Scholar] [CrossRef]
  11. Tharkeshwar, A.K.; Trekker, J.; Vermeire, W.; Pauwels, J.; Sannerud, R.; Priestman, D.A.; Te Vruchte, D.; Vints, K.; Baatsen, P.; Decuypere, J.-P.; et al. A Novel Approach to Analyze Lysosomal Dysfunctions through Subcellular Proteomics and Lipidomics: The Case of NPC1 Deficiency. Sci. Rep. 2017, 7, 41408. [Google Scholar] [CrossRef] [PubMed]
  12. Cao, Q.; Zhong, X.Z.; Zou, Y.; Zhang, Z.; Toro, L.; Dong, X.-P. BK Channels Alleviate Lysosomal Storage Diseases by Providing Positive Feedback Regulation of Lysosomal Ca2+ Release. Dev. Cell 2015, 33, 427–441. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, W.; Zhang, X.; Gao, Q.; Lawas, M.; Yu, L.; Cheng, X.; Gu, M.; Sahoo, N.; Li, X.; Li, P.; et al. A Voltage-Dependent K+ Channel in the Lysosome Is Required for Refilling Lysosomal Ca2+ Stores. J. Cell Biol. 2017, 216, 1715–1730. [Google Scholar] [CrossRef] [PubMed]
  14. Yuan, C.; O’Connell, R.J.; Feinberg-Zadek, P.L.; Johnston, L.J.; Treistman, S.N. Bilayer Thickness Modulates the Conductance of the BK Channel in Model Membranes. Biophys. J. 2004, 86, 3620–3633. [Google Scholar] [CrossRef] [PubMed]
  15. Lam, R.S.; Shaw, A.R.; Duszyk, M. Membrane Cholesterol Content Modulates Activation of BK Channels in Colonic Epithelia. Biochim. Biophys. Acta-Biomembr. 2004, 1667, 241–248. [Google Scholar] [CrossRef] [PubMed]
  16. Yuan, C.; O’Connell, R.J.; Jacob, R.F.; Mason, R.P.; Treistman, S.N. Regulation of the Gating of BKCa Channel by Lipid Bilayer Thickness. J. Biol. Chem. 2007, 282, 7276–7286. [Google Scholar] [CrossRef] [PubMed]
  17. Purcell, E.K.; Liu, L.; Thomas, P.V.; Duncan, R.K. Cholesterol Influences Voltage-Gated Calcium Channels and BK-Type Potassium Channels in Auditory Hair Cells. PLoS ONE 2011, 6, e26289. [Google Scholar] [CrossRef]
  18. Bukiya, A.N.; Belani, J.D.; Rychnovsky, S.; Dopico, A.M. Specificity of Cholesterol and Analogs to Modulate BK Channels Points to Direct Sterol–channel Protein Interactions. J. Gen. Physiol. 2010, 137, 93–110. [Google Scholar] [CrossRef]
  19. Yuan, C.; Chen, M.; Covey, D.F.; Johnston, L.J.; Treistman, S.N. Cholesterol Tuning of BK Ethanol Response Is Enantioselective, and Is a Function of Accompanying Lipids. PLoS ONE 2011, 6, e27572. [Google Scholar] [CrossRef]
  20. Pryor, P.R.; Reimann, F.; Gribble, F.M.; Luzio, J.P. Mucolipin-1 Is a Lysosomal Membrane Protein Required for Intracellular Lactosylceramide Traffic. Traffic 2006, 7, 1388–1398. [Google Scholar] [CrossRef]
  21. Singh, A.K.; McMillan, J.; Bukiya, A.N.; Burton, B.; Parrill, A.L.; Dopico, A.M. Multiple Cholesterol Recognition/Interaction Amino Acid Consensus (CRAC) Motifs in Cytosolic C Tail of Slo1 Subunit Determine Cholesterol Sensitivity of Ca2+- and Voltage-Gated K+ (BK) Channels. J. Biol. Chem. 2012, 287, 20509–20521. [Google Scholar] [CrossRef] [PubMed]
  22. Tao, X.; Hite, R.K.; MacKinnon, R. Cryo-EM Structure of the Open High-Conductance Ca2+-Activated K+ Channel. Nature 2017, 541, 46–51. [Google Scholar] [CrossRef] [PubMed]
  23. Hite, R.K.; Tao, X.; MacKinnon, R. Structural Basis for Gating the High-Conductance Ca2+-Activated K+ Channel. Nature 2017, 541, 52–57. [Google Scholar] [CrossRef] [PubMed]
  24. Shen, D.; Wang, X.; Li, X.; Zhang, X.; Yao, Z.; Dibble, S.; Dong, X.; Yu, T.; Lieberman, A.P.; Showalter, H.D.; et al. Lipid Storage Disorders Block Lysosomal Trafficking by Inhibiting a TRP Channel and Lysosomal Calcium Release. Nat. Commun. 2012, 3, 731. [Google Scholar] [CrossRef] [PubMed]
  25. Li, M.; Zhang, W.K.; Benvin, N.M.; Zhou, X.; Su, D.; Li, H.; Wang, S.; Michailidis, I.E.; Tong, L.; Li, X.; et al. Structural Basis of Dual Ca2+/PH Regulation of the Endolysosomal TRPML1 Channel. Nat. Struct. Mol. Biol. 2017, 24. [Google Scholar] [CrossRef] [PubMed]
  26. Chakraborty, K.; Leung, K.; Krishnan, Y. High Lumenal Chloride in the Lysosome Is Critical for Lysosome Function. Elife 2017, 6, e28862. [Google Scholar] [CrossRef] [PubMed]
  27. Lafourcade, C.; Sobo, K.; Kieffer-Jaquinod, S.; Garin, J.; van der Goot, F.G. Regulation of the V-ATPase along the Endocytic Pathway Occurs through Reversible Subunit Association and Membrane Localization. PLoS ONE 2008, 3, e2758. [Google Scholar] [CrossRef]
  28. Wheeler, S.; Haberkant, P.; Sullo, N.; Ferraz, M.J.; Halter, D.; Sprong, H.; Aerts, J.M.F.G.; Sillence, D.J. Cytoplasmic Glucosylceramide Modulates Endolysosomal Function in NPC Disease. Neurobiol. Dis. 2019. In press. [Google Scholar]
  29. Jiang, Z.; Redfern, R.E.; Isler, Y.; Ross, A.H.; Gericke, A. Cholesterol Stabilizes Fluid Phosphoinositide Domains. Chem. Phys. Lipids 2014, 182, 52–61. [Google Scholar] [CrossRef]
  30. Picas, L.; Gaits-Iacovoni, F.; Goud, B. The Emerging Role of Phosphoinositide Clustering in Intracellular Trafficking and Signal Transduction [Version 1; Referees: 4 Approved]. F1000Research 2016, 5. [Google Scholar] [CrossRef]
  31. Mayran, N.; Parton, R.G.; Gruenberg, J. Annexin II Regulates Multivesicular Endosome Biogenesis in the Degradation Pathway of Animal Cells. EMBO J. 2003, 22, 3242–3253. [Google Scholar] [CrossRef] [PubMed]
  32. te Vruchte, D.; Lloyd-Evans, E.; Veldman, R.J.; Neville, D.C.A.; Dwek, R.A.; Platt, F.M.; van Blitterswijk, W.J.; Sillence, D.J. Accumulation of Glycosphingolipids in Niemann-Pick C Disease Disrupts Endosomal Transport. J. Biol. Chem. 2004, 279, 26167–26175. [Google Scholar] [CrossRef] [PubMed]
  33. Sargiacomo, M.; Sudol, M.; Tang, Z.; Lisanti, M.P. Signal Transducing Molecules and Glycosyl-Phosphatidylinositol-Linked Proteins Form a Caveolin-Rich Insoluble Complex in MDCK Cells. J. Cell Biol. 1993, 122, 789–807. [Google Scholar] [CrossRef]
  34. Babiychuk, E.B.; Draeger, A. Annexins in Cell Membrane Dynamics. J. Cell Biol. 2000, 150, 1113–1124. [Google Scholar] [CrossRef] [PubMed]
  35. Liu, Y.; Casey, L.; Pike, L.J. Compartmentalization of Phosphatidylinositol 4,5-Bisphosphate in Low-Density Membrane Domains in the Absence of Caveolin. Biochem. Biophys. Res. Commun. 1998, 245, 684–690. [Google Scholar] [CrossRef] [PubMed]
  36. Gokhale, N.A.; Abraham, A.; Digman, M.A.; Gratton, E.; Cho, W. Phosphoinositide Specificity of and Mechanism of Lipid Domain Formation by Annexin A2-P11 Heterotetramer. J. Biol. Chem. 2005, 280, 42831–42840. [Google Scholar] [CrossRef] [PubMed]
  37. Drücker, P.; Pejic, M.; Grill, D.; Galla, H.J.; Gerke, V. Cooperative Binding of Annexin A2 to Cholesterol- and Phosphatidylinositol-4,5-Bisphosphate-Containing Bilayers. Biophys. J. 2014, 107, 2070–2081. [Google Scholar] [CrossRef] [PubMed]
  38. Pipalia, N.H.; Hao, M.; Mukherjee, S.; Maxfield, F.R. Sterol, Protein and Lipid Trafficking in Chinese Hamster Ovary Cells with Niemann-Pick Type C1 Defect. Traffic 2007, 8, 130–141. [Google Scholar] [CrossRef] [PubMed]
  39. Ward, C.; Martinez-Lopez, N.; Otten, E.G.; Carroll, B.; Maetzel, D.; Singh, R.; Sarkar, S.; Korolchuk, V.I. Autophagy, Lipophagy and Lysosomal Lipid Storage Disorders. Biochim. Biophys. Acta-Mol. Cell Biol. Lipids 2016, 1861, 269–284. [Google Scholar] [CrossRef] [PubMed]
  40. Fraldi, A.; Annunziata, F.; Lombardi, A.; Kaiser, H.-J.; Medina, D.L.; Spampanato, C.; Fedele, A.O.; Polishchuk, R.; Sorrentino, N.C.; Simons, K.; et al. Lysosomal Fusion and SNARE Function Are Impaired by Cholesterol Accumulation in Lysosomal Storage Disorders. EMBO J. 2010, 29, 3607–3620. [Google Scholar] [CrossRef] [PubMed]
  41. Murray, D.H.; Tamm, L.K. Molecular Mechanism of Cholesterol- and Polyphosphoinositide-Mediated Syntaxin Clustering. Biochemistry 2011, 50, 9014–9022. [Google Scholar] [CrossRef]
  42. Murray, D.H.; Tamm, L.K. Clustering of Syntaxin-1A in Model Membranes Is Modulated by Phosphatidylinositol 4,5-Bisphosphate and Cholesterol. Biochemistry 2009, 48, 4617–4625. [Google Scholar] [CrossRef] [PubMed]
  43. Khelashvili, G.; Galli, A.; Weinstein, H. Phosphatidylinositol 4,5-Biphosphate (PIP(2)) Lipids Regulate the Phosphorylation of Syntaxin N-Terminus by Modulating Both Its Position and Local Structure. Biochemistry 2012, 51, 7685–7698. [Google Scholar] [CrossRef]
  44. Singer-Lahat, D.; Barak-Broner, N.; Sheinin, A.; Greitzer-Antes, D.; Michaelevski, I.; Lotan, I. The Dual Function of the Polybasic Juxtamembrane Region of Syntaxin 1A in Clamping Spontaneous Release and Stimulating Ca2+-Triggered Release in Neuroendocrine Cells. J. Neurosci. 2018, 38, 220–231. [Google Scholar] [CrossRef] [PubMed]
  45. Li, H.; Papadopoulos, V. Peripheral-Type Benzodiazepine Receptor Function in Cholesterol Transport. Identification of a Putative Cholesterol Recognition/Interaction Amino Acid Sequence and Consensus Pattern. Endocrinology 1998, 139, 4991–4997. [Google Scholar] [CrossRef] [PubMed]
  46. Fantini, J.; Barrantes, F.J. How Cholesterol Interacts with Membrane Proteins: An Exploration of Cholesterol-Binding Sites Including CRAC, CARC and Tilted Domains. Front. Physiol. 2013, 4. [Google Scholar] [CrossRef]
  47. Baier, C.J.; Fantini, J.; Barrantes, F.J. Disclosure of Cholesterol Recognition Motifs in Transmembrane Domains of the Human Nicotinic Acetylcholine Receptor. Sci. Rep. 2011, 1, 69. [Google Scholar] [CrossRef]
  48. Mahfoud, R.; Garmy, N.; Maresca, M.; Yahi, N.; Puigserver, A.; Fantini, J. Identification of a Common Sphingolipid-Binding Domain in Alzheimer, Prion, and HIV-1 Proteins. J. Biol. Chem. 2002, 277, 11292–11296. [Google Scholar] [CrossRef]
  49. Contreras, F.-X.; Ernst, A.M.; Haberkant, P.; Bjorkholm, P.; Lindahl, E.; Gonen, B.; Tischer, C.; Elofsson, A.; von Heijne, G.; Thiele, C.; et al. Molecular Recognition of a Single Sphingolipid Species by a Protein/’s Transmembrane Domain. Nature 2012, 481, 525–529. [Google Scholar] [CrossRef]
  50. Björkholm, P.; Ernst, A.M.; Hacke, M.; Wieland, F.; Brügger, B.; von Heijne, G. Identification of Novel Sphingolipid-Binding Motifs in Mammalian Membrane Proteins. Biochim. Biophys. Acta-Biomembr. 2014, 1838, 2066–2070. [Google Scholar] [CrossRef]
  51. Lomize, M.A.; Pogozheva, I.D.; Joo, H.; Mosberg, H.I.; Lomize, A.L. OPM Database and PPM Web Server: Resources for Positioning of Proteins in Membranes. Nucleic Acids Res. 2012, 40, 370–376. [Google Scholar] [CrossRef]
  52. Khelashvili, G.; Pabst, G.; Harries, D. Cholesterol Orientation and Tilt Modulus in DMPC Bilayers. J. Phys. Chem. B 2010, 114, 7524–7534. [Google Scholar] [CrossRef]
  53. Millard, E.E.; Gale, S.E.; Dudley, N.; Zhang, J.; Schaffer, J.E.; Ory, D.S. The Sterol-Sensing Domain of the Niemann-Pick C1 (NPC1) Protein Regulates Trafficking of Low Density Lipoprotein Cholesterol. J. Biol. Chem. 2005, 280, 28581–28590. [Google Scholar] [CrossRef]
  54. Ohgami, N.; Ko, D.C.; Thomas, M.; Scott, M.P.; Chang, C.C.Y.; Chang, T.-Y. Binding between the Niemann–Pick C1 Protein and a Photoactivatable Cholesterol Analog Requires a Functional Sterol-Sensing Domain. Proc. Natl. Acad. Sci. USA 2004, 101, 12473–12478. [Google Scholar] [CrossRef]
  55. Castellano, B.M.; Thelen, A.M.; Moldavski, O.; Feltes, M.; van der Welle, R.E.N.; Mydock-McGrane, L.; Jiang, X.; van Eijkeren, R.J.; Davis, O.B.; Louie, S.M.; et al. Lysosomal Cholesterol Activates MTORC1 via an SLC38A9-Niemann-Pick C1 Signaling Complex. Science 2017, 355, 1306–1311. [Google Scholar] [CrossRef]
  56. Di Scala, C.; Fantini, J.; Yahi, N.; Barrantes, F.J.; Chahinian, H. Anandamide Revisited: How Cholesterol and Ceramides Control Receptor-Dependent and Receptor-Independent Signal Transmission Pathways of a Lipid Neurotransmitter. Biomolecules 2018, 8. [Google Scholar] [CrossRef]
  57. Fantini, J.; Di Scala, C.; Evans, L.S.; Williamson, P.T.F.; Barrantes, F.J. A Mirror Code for Protein-Cholesterol Interactions in the Two Leaflets of Biological Membranes. Sci. Rep. 2016, 6, 21907. [Google Scholar] [CrossRef]
  58. Lloyd-Evans, E.; Morgan, A.J.; He, X.; Smith, D.A.; Elliot-Smith, E.; Sillence, D.J.; Churchill, G.C.; Schuchman, E.H.; Galione, A.; Platt, F.M. Niemann-Pick Disease Type C1 Is a Sphingosine Storage Disease That Causes Deregulation of Lysosomal Calcium. Nat. Med. 2008, 14, 1247–1255. [Google Scholar] [CrossRef]
  59. Höglinger, D.; Nadler, A.; Haberkant, P.; Kirkpatrick, J.; Schifferer, M.; Stein, F.; Hauke, S.; Porter, F.D.; Schultz, C. Trifunctional Lipid Probes for Comprehensive Studies of Single Lipid Species in Living Cells. Proc. Natl. Acad. Sci. USA 2017, 114, 1566–1571. [Google Scholar] [CrossRef]
  60. Höglinger, D.; Haberkant, P.; Aguilera-Romero, A.; Riezman, H.; Porter, F.D.; Platt, F.M.; Galione, A.; Schultz, C. Intracellular Sphingosine Releases Calcium from Lysosomes. Elife 2015, 4, e10616. [Google Scholar] [CrossRef]
  61. Schonauer, S.; Körschen, H.G.; Penno, A.; Rennhack, A.; Breiden, B.; Sandhoff, K.; Gutbrod, K.; Dörmann, P.; Raju, D.N.; Haberkant, P.; et al. Identification of a Feedback Loop Involving Beta-Glucosidase 2 and Its Product Sphingosine Sheds Light on the Molecular Mechanisms in Gaucher Disease. J. Biol. Chem. 2017. [Google Scholar] [CrossRef]
  62. Ridley, C.M.; Thur, K.E.; Shanahan, J.; Thillaiappan, N.B.; Shen, A.; Uhl, K.; Walden, C.M.; Rahim, A.A.; Waddington, S.N.; Platt, F.M.; et al. Beta-Glucosidase 2 (GBA2) Activity and Imino Sugar Pharmacology. J. Biol. Chem. 2013. [Google Scholar] [CrossRef]
  63. Sasaki, H.; Arai, H.; Cocco, M.J.; White, S.H. PH Dependence of Sphingosine Aggregation. Biophys. J. 2009, 96, 2727–2733. [Google Scholar] [CrossRef]
  64. Hanson, M.A.; Cherezov, V.; Griffith, M.T.; Roth, C.B.; Jaakola, V.-P.; Chien, E.Y.T.; Velasquez, J.; Kuhn, P.; Stevens, R.C. A Specific Cholesterol Binding Site Is Established by the 2.8 Å Structure of the Human Β2-Adrenergic Receptor. Structure 2008, 16, 897–905. [Google Scholar] [CrossRef]
  65. Xu, M.; Liu, K.; Swaroop, M.; Porter, F.D.; Sidhu, R.; Finkes, S.; Ory, D.S.; Marugan, J.J.; Xiao, J.; Southall, N.; et al. D-Tocopherol Reduces Lipid Accumulation in Niemann-Pick Type C1 and Wolman Cholesterol Storage Disorders. J. Biol. Chem. 2012, 287, 39349–39360. [Google Scholar] [CrossRef]
  66. Garrity, A.G.; Wang, W.; Collier, C.M.D.; Levey, S.A.; Gao, Q.; Xu, H. The Endoplasmic Reticulum, Not the PH Gradient, Drives Calcium Refilling of Lysosomes. Elife 2016, 5, e15887. [Google Scholar] [CrossRef]
  67. Zhong, X.Z.; Sun, X.; Cao, Q.; Dong, G.; Schiffmann, R.; Dong, X.P. BK Channel Agonist Represents a Potential Therapeutic Approach for Lysosomal Storage Diseases. Sci. Rep. 2016, 6, 33684. [Google Scholar] [CrossRef]
  68. Schmiege, P.; Fine, M.; Blobel, G.; Li, X. Human TRPML1 Channel Structures in Open and Closed Conformations. Nature 2017, 550, 366. [Google Scholar] [CrossRef]
  69. Fine, M.; Schmiege, P.; Li, X. Structural Basis for PtdInsP2-Mediated Human TRPML1 Regulation. Nat. Commun. 2018, 9, 4192. [Google Scholar] [CrossRef]
  70. Hirschi, M.; Herzik, M.A., Jr.; Wie, J.; Suo, Y.; Borschel, W.F.; Ren, D.; Lander, G.C.; Lee, S.-Y. Cryo-Electron Microscopy Structure of the Lysosomal Calcium-Permeable Channel TRPML3. Nature 2017, 550, 411. [Google Scholar] [CrossRef]
  71. Venkatachalam, K.; Kiselyov, K. Chapter 25—TRPML1-Dependent Processes as Therapeutic Targets. In TRP Channels as Therapeutic Targets; Szallasi, A., Ed.; Academic Press: Boston, MA, USA, 2015; pp. 469–482. [Google Scholar]
  72. Vergarajauregui, S.; Connelly, P.S.; Daniels, M.P.; Puertollano, R. Autophagic Dysfunction in Mucolipidosis Type IV Patients. Hum. Mol. Genet. 2008, 17, 2723–2737. [Google Scholar] [CrossRef]
  73. Choudhury, A.; Dominguez, M.; Puri, V.; Sharma, D.K.; Narita, K.; Wheatley, C.L.; Marks, D.L.; Pagano, R.E. Rab Proteins Mediate Golgi Transport of Caveola-Internalized Glycosphingolipids and Correct Lipid Trafficking in Niemann-Pick C Cells. J. Clin. Investig. 2002, 109, 1541–1550. [Google Scholar] [CrossRef]
  74. Harvald, E.B.; Olsen, A.S.B.; Faergeman, N.J. Autophagy in the Light of Sphingolipid Metabolism. Apoptosis 2015, 20, 658–670. [Google Scholar] [CrossRef]
  75. Ruas, M.; Rietdorf, K.; Arredouani, A.; Davis, L.C.; Lloyd-Evans, E.; Koegel, H.; Funnell, T.M.; Morgan, A.J.; Ward, J.A.; Watanabe, K.; et al. Purified TPC Isoforms Form NAADP Receptors with Distinct Roles for Ca2+ Signaling and Endolysosomal Trafficking. Curr. Biol. 2010, 20, 703–709. [Google Scholar] [CrossRef]
  76. Lin, P.-H.; Duann, P.; Komazaki, S.; Park, K.H.; Li, H.; Sun, M.; Sermersheim, M.; Gumpper, K.; Parrington, J.; Galione, A.; et al. Lysosomal Two-Pore Channel Subtype 2 (TPC2) Regulates Skeletal Muscle Autophagic Signaling. J. Biol. Chem. 2015, 290, 3377–3389. [Google Scholar] [CrossRef]
  77. Cao, Q.; Zhong, X.Z.; Zou, Y.; Murrell-Lagnado, R.; Zhu, M.X.; Dong, X.-P. Calcium Release through P2X4 Activates Calmodulin to Promote Endolysosomal Membrane Fusion. J. Cell Biol. 2015, 209, 879–894. [Google Scholar] [CrossRef]
  78. Tian, X.; Gala, U.; Zhang, Y.; Shang, W.; Nagarkar Jaiswal, S.; di Ronza, A.; Jaiswal, M.; Yamamoto, S.; Sandoval, H.; Duraine, L.; et al. A Voltage-Gated Calcium Channel Regulates Lysosomal Fusion with Endosomes and Autophagosomes and Is Required for Neuronal Homeostasis. PLOS Biol. 2015, 13, e1002103. [Google Scholar] [CrossRef]
  79. Reddy, J.V.; Ganley, I.G.; Pfeffer, S.R. Clues to Neuro-Degeneration in Niemann-Pick Type C Disease from Global Gene Expression Profiling. PLoS ONE 2006, 1, e19. [Google Scholar] [CrossRef] [PubMed]
  80. Shao, C.; Zhang, F.; Kemp, M.M.; Linhardt, R.J.; Waisman, D.M.; Head, J.F.; Seaton, B.A. Crystallographic Analysis of Calcium-Dependent Heparin Binding to Annexin A2. J. Biol. Chem. 2006, 281, 31689–31695. [Google Scholar] [CrossRef]
  81. Mullock, B.M.; Smith, C.W.; Ihrke, G.; Bright, N.A.; Lindsay, M.; Parkinson, E.J.; Brooks, D.A.; Parton, R.G.; James, D.E.; Luzio, J.P.; et al. Syntaxin 7 Is Localized to Late Endosome Compartments, Associates with Vamp 8, and Is Required for Late Endosome–Lysosome Fusion. Mol. Biol. Cell 2000, 11, 3137–3153. [Google Scholar] [CrossRef] [PubMed]
  82. Stein, A.; Weber, G.; Wahl, M.C.; Jahn, R. Helical Extension of the Neuronal SNARE Complex into the Membrane. Nature 2009, 460, 525–528. [Google Scholar] [CrossRef] [PubMed]
  83. Itakura, E.; Kishi-Itakura, C.; Mizushima, N. The Hairpin-Type Tail-Anchored SNARE Syntaxin 17 Targets to Autophagosomes for Fusion with Endosomes/Lysosomes. Cell 2012, 151, 1256–1269. [Google Scholar] [CrossRef] [PubMed]
  84. Takáts, S.; Nagy, P.; Varga, Á.; Pircs, K.; Kárpáti, M.; Varga, K.; Kovács, A.L.; Hegedűs, K.; Juhász, G. Autophagosomal Syntaxin17-Dependent Lysosomal Degradation Maintains Neuronal Function in Drosophila. J. Cell Biol. 2013, 201, 531–539. [Google Scholar] [CrossRef] [PubMed]
  85. Furuta, N.; Fujita, N.; Noda, T.; Yoshimori, T.; Amano, A. Combinational Soluble N-Ethylmaleimide-Sensitive Factor Attachment Protein Receptor Proteins VAMP8 and Vti1b Mediate Fusion of Antimicrobial and Canonical Autophagosomes with Lysosomes. Mol. Biol. Cell 2010, 21, 1001–1010. [Google Scholar] [CrossRef] [PubMed]
  86. Rosenhouse-Dantsker, A.; Mehta, D.; Levitan, I. Regulation of Ion Channels by Membrane Lipids. In Comprehensive Physiology; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2011. [Google Scholar]
  87. Dong, X.; Shen, D.; Wang, X.; Dawson, T.; Li, X.; Zhang, Q.; Cheng, X.; Zhang, Y.; Weisman, L.S.; Delling, M.; et al. PI(3,5)P2 Controls Membrane Trafficking by Direct Activation of Mucolipin Ca2+ Release Channels in the Endolysosome. Nat. Commun. 2010, 1, 38. [Google Scholar] [CrossRef] [PubMed]
  88. Zhang, X.; Li, X.; Xu, H. Phosphoinositide Isoforms Determine Compartment-Specific Ion Channel Activity. Proc. Natl. Acad. Sci. USA 2012, 109, 11384–11389. [Google Scholar] [CrossRef] [PubMed]
  89. Wang, X.; Zhang, X.; Dong, X.; Samie, M.; Li, X.; Cheng, X.; Goschka, A.; Shen, D.; Zhou, Y.; Harlow, J.; et al. TPC Proteins Are Phosphoinositide- Activated Sodium-Selective Ion Channels in Endosomes and Lysosomes. Cell 2012, 151, 372–383. [Google Scholar] [CrossRef]
  90. Domon, M.M.; Besson, F.; Bandorowicz-Pikula, J.; Pikula, S. Annexin A6 Is Recruited into Lipid Rafts of Niemann–Pick Type C Disease Fibroblasts in a Ca2+-Dependent Manner. Biochem. Biophys. Res. Commun. 2011, 405, 192–196. [Google Scholar] [CrossRef]
  91. Ganley, I.G.; Pfeffer, S.R. Cholesterol Accumulation Sequesters Rab9 and Disrupts Late Endosome Function in NPC1-Deficient Cells. J. Biol. Chem. 2006, 281, 17890–17899. [Google Scholar] [CrossRef]
  92. Kaptzan, T.; West, S.A.; Holicky, E.L.; Wheatley, C.L.; Marks, D.L.; Wang, T.; Peake, K.B.; Vance, J.; Walkley, S.U.; Pagano, R.E. Development of a Rab9 Transgenic Mouse and Its Ability to Increase the Lifespan of a Murine Model of Niemann-Pick Type C Disease. Am. J. Pathol. 2009, 174, 14–20. [Google Scholar] [CrossRef]
  93. Sievers, F.; Wilm, A.; Dineen, D.; Gibson, T.J.; Karplus, K.; Li, W.; Lopez, R.; McWilliam, H.; Remmert, M.; Söding, J.; et al. Fast, Scalable Generation of High-quality Protein Multiple Sequence Alignments Using Clustal Omega. Mol. Syst. Biol. 2011, 7. [Google Scholar] [CrossRef] [PubMed]
  94. Waterhouse, A.M.; Procter, J.B.; Martin, D.M.A.; Clamp, M.; Barton, G.J. Jalview Version 2—A Multiple Sequence Alignment Editor and Analysis Workbench. Bioinformatics 2009, 25, 1189–1191. [Google Scholar] [CrossRef] [PubMed]
  95. Arnold, K.; Bordoli, L.; Kopp, J.; Schwede, T. The SWISS-MODEL Workspace: A Web-Based Environment for Protein Structure Homology Modelling. Bioinformatics 2006, 22, 195–201. [Google Scholar] [CrossRef] [PubMed]
  96. Benkert, P.; Biasini, M.; Schwede, T. Toward the Estimation of the Absolute Quality of Individual Protein Structure Models. Bioinformatics 2011, 27, 343–350. [Google Scholar] [CrossRef] [PubMed]
  97. Studer, G.; Biasini, M.; Schwede, T. Assessing the Local Structural Quality of Transmembrane Protein Models Using Statistical Potentials (QMEANBrane). Bioinformatics 2014, 30, i505–i511. [Google Scholar] [CrossRef]
  98. Ionescu, C.-M.; Sehnal, D.; Falginella, F.L.; Pant, P.; Pravda, L.; Bouchal, T.; Svobodová Vařeková, R.; Geidl, S.; Koča, J. AtomicChargeCalculator: Interactive Web-Based Calculation of Atomic Charges in Large Biomolecular Complexes and Drug-like Molecules. J. Cheminform. 2015, 7, 50. [Google Scholar] [CrossRef]
  99. Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. AutoDock4 and AutoDockTools4: Automated Docking with Selective Receptor Flexibility. J. Comput. Chem. 2009, 30, 2785–2791. [Google Scholar] [CrossRef]
  100. Sanner, M.F. Python: A Programming Language for Software Integration and Development. J. Mol. Graph. Model. 1999, 17, 57–61. [Google Scholar]
  101. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera—A Visualization System for Exploratory Research and Analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef]
Ijms 20 00717 g001 550
Figure 1. Pathological cascade in Niemann–Pick type C disease (NPCD). All steps are supported by at least one published report referenced in the text; steps shown with bold arrows are the subject of this publication.
Figure 1. Pathological cascade in Niemann–Pick type C disease (NPCD). All steps are supported by at least one published report referenced in the text; steps shown with bold arrows are the subject of this publication.
Ijms 20 00717 g001
Ijms 20 00717 g002 550
Figure 2. Four possibilities for cholesterol binding to a trans-membrane protein displayed in cartoon form. Trans-membrane proteins (light grey) are embedded in the lipid bilayer (dark grey). Cholesterol (orange) may bind parallel to the membrane (A, energetically favored), parallel to the membrane but vertically displaced from it (B, energetically less favored), perpendicular to the membrane (C, energetically disfavored), or antiparallel (D, energetically disfavored).
Figure 2. Four possibilities for cholesterol binding to a trans-membrane protein displayed in cartoon form. Trans-membrane proteins (light grey) are embedded in the lipid bilayer (dark grey). Cholesterol (orange) may bind parallel to the membrane (A, energetically favored), parallel to the membrane but vertically displaced from it (B, energetically less favored), perpendicular to the membrane (C, energetically disfavored), or antiparallel (D, energetically disfavored).
Ijms 20 00717 g002
Ijms 20 00717 g003 550
Figure 3. Docking identifies two putative binding pockets on NPC1 aligned to each side of the lysosomal membrane. Sample poses for both cholesterol (orange, A) and sphingosine (purple, B) show binding in these pockets consistent with evidence of NPC1’s role as an exporter of both lipids. Selected interacting residues are labeled; H-bonds are shown in green (Sph and cholesterol), and ionic interactions are shown in magenta (Sph only). Pro691 is shown in pink; mutation of this residue results in reduced cholesterol binding [53,54]. Lines indicate the approximate position of the membrane (derived from the Orientation of Proteins in Membranes OPM database [51]).
Figure 3. Docking identifies two putative binding pockets on NPC1 aligned to each side of the lysosomal membrane. Sample poses for both cholesterol (orange, A) and sphingosine (purple, B) show binding in these pockets consistent with evidence of NPC1’s role as an exporter of both lipids. Selected interacting residues are labeled; H-bonds are shown in green (Sph and cholesterol), and ionic interactions are shown in magenta (Sph only). Pro691 is shown in pink; mutation of this residue results in reduced cholesterol binding [53,54]. Lines indicate the approximate position of the membrane (derived from the Orientation of Proteins in Membranes OPM database [51]).
Ijms 20 00717 g003
Ijms 20 00717 g004 550
Figure 4. Sample cholesterol-binding poses at the putative binding site on the voltage sensor of big potassium channel (BK). Cholesterol (orange) is able to form H-bonds (green) with Lys211 (A,B) and Asn265 (A), a CH–π interaction with Tyr263 (B), and a hydrophobic interaction with Leu226 (A,B).
Figure 4. Sample cholesterol-binding poses at the putative binding site on the voltage sensor of big potassium channel (BK). Cholesterol (orange) is able to form H-bonds (green) with Lys211 (A,B) and Asn265 (A), a CH–π interaction with Tyr263 (B), and a hydrophobic interaction with Leu226 (A,B).
Ijms 20 00717 g004
Ijms 20 00717 g005 550
Figure 5. Putative binding sites of lipids at Transient receptor potential mucolipin 1 (TRPML1). (A,B) Modeling supports mutagenesis experiments identifying a polybasic region as key the binding site for phosphatidylinositol-3,5-bisphosphate (PI(3,5)P2) (green, A) and PI(4,5)P2 (green, B). (C) Glu276 from one chain (blue) forms a charged interaction with Arg486 from another (grey); interatomic distance: 3.77 Å. (D) Modeling suggests plausible binding poses for SM (purple, typical example shown), but these do not disturb the inter-residue interactions present in the apo state; interatomic distance: 2.97 Å. Protein oriented as in Figure S3; distance markers are shown in yellow.
Figure 5. Putative binding sites of lipids at Transient receptor potential mucolipin 1 (TRPML1). (A,B) Modeling supports mutagenesis experiments identifying a polybasic region as key the binding site for phosphatidylinositol-3,5-bisphosphate (PI(3,5)P2) (green, A) and PI(4,5)P2 (green, B). (C) Glu276 from one chain (blue) forms a charged interaction with Arg486 from another (grey); interatomic distance: 3.77 Å. (D) Modeling suggests plausible binding poses for SM (purple, typical example shown), but these do not disturb the inter-residue interactions present in the apo state; interatomic distance: 2.97 Å. Protein oriented as in Figure S3; distance markers are shown in yellow.
Ijms 20 00717 g005
Ijms 20 00717 g006 550
Figure 6. PI(4,5)P2 can bind to annexin A2 (AnxA2). PI(4,5)P2 (turquoise) can bind to the convex face of AnxA2 (grey) driven by zwitterionic interactions (magenta) with Lys279, Lys281, and Arg284.
Figure 6. PI(4,5)P2 can bind to annexin A2 (AnxA2). PI(4,5)P2 (turquoise) can bind to the convex face of AnxA2 (grey) driven by zwitterionic interactions (magenta) with Lys279, Lys281, and Arg284.
Ijms 20 00717 g006
Ijms 20 00717 g007 550
Figure 7. PI(4,5)P2 can bind to two members of a SNARE bundle simultaneously. Sample poses of PI(3,5)P2 (turquoise) bridging vesicle-associated membrane protein 8( VAMP8) (blue) and syntaxin 7 (grey) in the assembled Soluble NSF (N-ethylmaleimide Sensitive Fusion) protein attachment receptor (SNARE) bundle with zwitterionic interactions (magenta) (A,B).
Figure 7. PI(4,5)P2 can bind to two members of a SNARE bundle simultaneously. Sample poses of PI(3,5)P2 (turquoise) bridging vesicle-associated membrane protein 8( VAMP8) (blue) and syntaxin 7 (grey) in the assembled Soluble NSF (N-ethylmaleimide Sensitive Fusion) protein attachment receptor (SNARE) bundle with zwitterionic interactions (magenta) (A,B).
Ijms 20 00717 g007

Share and Cite

MDPI and ACS Style

Wheeler, S.; Schmid, R.; Sillence, D.J. Lipid–Protein Interactions in Niemann–Pick Type C Disease: Insights from Molecular Modeling. Int. J. Mol. Sci. 2019, 20, 717. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20030717

AMA Style

Wheeler S, Schmid R, Sillence DJ. Lipid–Protein Interactions in Niemann–Pick Type C Disease: Insights from Molecular Modeling. International Journal of Molecular Sciences. 2019; 20(3):717. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20030717

Chicago/Turabian Style

Wheeler, Simon, Ralf Schmid, and Dan J Sillence. 2019. "Lipid–Protein Interactions in Niemann–Pick Type C Disease: Insights from Molecular Modeling" International Journal of Molecular Sciences 20, no. 3: 717. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20030717

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop