1. Introduction
Myo-inositol polyphosphates (InsPs) are a family of messenger molecules that control a wide array of biological processes in eukaryotic cells (
Figure 1a). In these molecules, the phosphoryl groups are arranged in different numbers and patterns around the inositol scaffold, creating great structural variety that ranges from species with only one phosphoryl group up to the inositol pyrophosphates (PP-InsPs), which carry seven or eight phosphates [
1]. The PP-InsPs contain one or two high-energy diphosphate groups in addition to monophosphoryl groups, thereby accommodating an extraordinary negative charge density.
The most highly phosphorylated canonical PP-InsP in mammals is 1,5(PP)
2-InsP
4, (from here on abbreviated as InsP
8), which is generated
via the phosphorylation of 5PP-InsP
5 (also called InsP
7 or 5-IP
7) by PPIP5Ks (diphosphoinositol-pentakisphosphate kinases) or phosphorylation of 1PP-InsP
5 by IP6Ks (inositol hexakisphosphate kinases). In recent years, InsP
8 has emerged as a sensor and regulator of inorganic phosphate in various organisms. The enzymatic activity of PPIP5Ks is directly regulated by phosphate concentrations, and cellular levels of InsP
8 correlate with phosphate availability [
2]. In fission yeast
S.pombe, InsP
8 activates the vacuolar VTC complex by binding to an SPX domain, driving polyphosphate synthesis to store excess phosphate [
3]. In plants, it was demonstrated that the activity of the InsP
8-producing enzymes VIH1/2 negatively regulates phosphate starvation responses [
4]. Mechanistically, InsP
8 binds to a standalone SPX domain (SPX1) in
arabidopsis thaliana, which enables dimerization with the transcription factor PHR1 and suppresses downstream targets [
5]. In cultured human cell lines, InsP
8 was found to activate phosphate efflux by binding to the SPX domain of XPR1, removing excess phosphate and preventing tissue calcifications [
6,
7,
8]. Moreover, at an organismal level, mutations in PPIP5Ks are associated with hearing loss [
9,
10].
The function of 5PP-InsP
5, has been investigated more closely, especially its role in cellular energy signaling, where its concentration reacts to ATP levels and regulates glycolytic flux accordingly [
11,
12]. Notably, 5PP-InsP
5 stimulates exocytosis of insulin-containing granules from pancreatic β-cells, ostensibly by competing with the plasma membrane lipid phosphatidylinositol-4,5-bisphosphate (PIP
2) for binding of synaptotagmin 7 [
13,
14]. 5PP-InsP
5 also inhibits the activity of Akt kinase, a central metabolic regulator, by binding to the PH domain of Akt, releasing it from the plasma membrane, preventing its activation [
15,
16,
17]. As a result, IP6K1 knockout mice display abnormally high Akt activity and are resistant to weight-gain when fed a high-fat diets [
15]. IP6K1 knockout in mice also increases the activity of AMPK, inhibiting fat accumulation in favor of thermogenesis, and further supporting the lean phenotype [
18].
In many cases it is difficult to attribute an observed phenotype to an individual messenger, because genetic deletion of IP6Ks inevitably reduces the levels of both 5PP-InsP
5 and InsP
8 [
19]. Even when biochemical insight into the mechanisms of action is available, it sometimes cannot explain why PP-InsPs possess distinct biological functions in spite of their very similar structure and similarly high charge density. For example, phosphate efflux
via XPR1 is activated almost exclusively by InsP
8, even though InsP
6 and 5PP-InsP
5 bind with relatively similar dissociation constants (K
d) to the SPX domain [
6]. Binding affinities are also similar for InsP
8 or 5PP-InsP
5 towards the SPX1 protein in rice, but genetic perturbation experiments showed that only InsP
8 is the physiologically relevant ligand [
5]. Similarly, in
in vitro experiments, InsP
8 is 20-fold more effective than 5PP-InsP
5 at activating polyphosphate synthesis by the VTC complex from budding yeast
S.cerevisiae, despite a charge difference of only one. It was also more effective than other InsP
8 isomers [
20], suggesting that affinity is determined by the exact shape of the molecule, rather than sheer charge density. Due to its higher concentrations, 5PP-InsP
5 is currently assumed to be the relevant ligand. In another example from fission yeast
S.pombe, polyphosphate synthesis by the VTC complex clearly depended on InsP
8 by the PPIP5K orthologue Asp1 [
3]. Possible explanations how substrate specificity of PP-InsP-binding proteins may be regulated include ternary interactions with more than one binding partner, local enrichment of individual PP-InsPs and localized differences in solution conditions, which change protonation and metal complexation.
A property unique to highly phosphorylated InsPs is their ability to adopt an alternate 1-axial / 5-equatorial conformation conformation at elevated pH, where the substituent at the 2-position is equatorial, and all others are axial, which separates phosphoryl groups further in space and reduces charge repulsion between them (
Figure 1b,c) [
21,
22,
23]. The PP-InsPs and their non-hydrolyzable methylene bisphosphonate analogs (PCP-InsPs,
Figure 1b) also display this behavior, and appear to have a higher propensity to adopt an axial conformation relative to the InsPs, especially in the presence of magnesium cations [
23].
Figure 1.
a) Structure of InsP
6, 5PP-InsP
5, InsP
8 and approximate concentrations in mammalian cells. [
24,
25,
26,
27]
b) pH-dependent equilibrium between axial (ax.) and equatorial (eq.) conformation of the non-hydrolyzable 5PP-InsP
5 analogue 5PCP-InsP
5 [
23]
c) Electrostatic potential mapped on an isodensity surface for fully deprotonated, axial 5PCP-InsP
5 (L
13-) and monoprotonated, equatorial 5PCP-InsP
5 (HL
12-) (B3LYP/3-21+G* geometries; isodensity value = 0.004, scale: -1.3 V (red) to -1.0 V (blue)). Atom color code: C (grey), H (white), O (red), P (orange).
Figure 1.
a) Structure of InsP
6, 5PP-InsP
5, InsP
8 and approximate concentrations in mammalian cells. [
24,
25,
26,
27]
b) pH-dependent equilibrium between axial (ax.) and equatorial (eq.) conformation of the non-hydrolyzable 5PP-InsP
5 analogue 5PCP-InsP
5 [
23]
c) Electrostatic potential mapped on an isodensity surface for fully deprotonated, axial 5PCP-InsP
5 (L
13-) and monoprotonated, equatorial 5PCP-InsP
5 (HL
12-) (B3LYP/3-21+G* geometries; isodensity value = 0.004, scale: -1.3 V (red) to -1.0 V (blue)). Atom color code: C (grey), H (white), O (red), P (orange).
PP-InsPs can potentially form a wide range of different species, depending on their conformation, protonation state and complexation of metal cations, each of which can be expected to interact in a different way with other biomolecules. We therefore wanted to obtain a more detailed understanding of PP-InsP speciation, specifically under conditions approximating a cytosolic setting. NMR spectroscopy is the method of choice for characterizing these equilibria due to the structural information it provides, but the lack of sensitivity has limited past investigations to experiments using non-physiological concentrations of InsPs and PP-InsPs (> 100 µM) [
23].
Here, we have conducted a systematic comparison of InsP6, 5PP-InsP5 and InsP8 using NMR spectroscopy of 13C-labeled InsP messengers, to better understand the intricacies of their speciation. Intriguingly, we found that InsP8 is able to adopt axial conformation under conditions very reminiscent of a cytosolic environment. We then used 31P-NMR to identify likely protonation states and K+- and mononuclear Mg2+-complexes of InsP8 across the pH range. Finally, ITC experiments revealed that addition of Mg2+, favoring axial conformation, influenced the binding parameters of the interaction between InsP8 and an SPX protein domain. Our results imply that conditions for biochemical and biophysical characterization of PP-InsP protein interactions should always be chosen with great care, and highlight InsP8 as a potential pH-, metal ion-, and temperature-dependent intracellular molecular switch.
2. Materials and Methods
Synthesis and BIRD-HMQC NMR analysis of InsPs
13C-labeled InsPs and PP-InsPs were synthesized chemo-enzymatically, as previously published, with slight adjustments [
25,
28]. Enzymatic synthesis of InsP
8 was carried out at pH 6.0 instead of 6.4, and in the presence of additional 150 mM (NH
4)
2SO
4.
Samples for HMQC-NMR studies contained 50 µM 13C6-labeled InsPs (InsP6 / 5PP-InsP5 / InsP8), 2 mM bis-tris-propane, 130 mM KCl, 10 mM NaCl and 0 / 50 / 250 µM MgCl2 in D2O.
BIRD-HMQC NMR spectra were recorded at 277 K and 600 MHz (
1H frequency) on a Bruker AV-III NMR spectrometer (Bruker Biospin, Rheinstetten, Germany) using cryogenically cooled 5 mm QCI-triple resonance probe equipped with one-axis self-shielded gradients. The spectrometer was operated using topspin 3.5 pl6 software. Instrument temperature was calibrated against d
4-methanol according to Findeisen
et al. [
29]. Acquisition parameters were SW(
13C): 60-90 ppm, NS: 128, TD(
13C): 64
NMR titrations
In order to create a system free of coordinating counterions, a batch of InsP8 was synthesized as described above, by enzymatic phosphorylation of 5PP-InsP5 and precipitation with Mg2+, followed by Mg2+ chelation on Amberlite® cation exchange resin. Unlike in the standard procedure, the resin was loaded with NMe4Cl instead of NH4CO3, resulting in a batch of InsP8 with non-coordinating NMe4+ counterions. 1 mM EDTA was added to all 31P-NMR titrations to chelate remaining traces of Mg2+.
NMR-samples contained 1 mM InsP8, 1 mM EDTA, pH 3.0 – 12.5 (in H2O, steps of 0.5 pH units) and a) for the non-coordinating condition: 150 mM NMe4Cl or b) 150 mM KCl or c) 150 mM KCl and 1 mM MgCl2. Sealed glass capillaries containing 50 mM phosphonoacetic acid in D2O were added into the sample tubes for locking and chemical shift calibration.
31P-NMR spectra were recorded at 295 K on a Bruker spectrometer (see above) operating at 600 MHz for protons and 244 MHz for phosphorous nuclei. SW: -40-20 ppm, NS: 1024.
Spectra at pH 3.0 were assigned by various 2D-NMR methods (
Figure S5) and chemical shift changes were tracked across the pH range. The data were analyzed using the HypNMR 2006 software [
31]. Different possible stoichiometries were tested, and the final chemical models were selected on the basis of the σ parameter (scaled sum of square differences between predicted and experimental chemical shift values), the model confidence level estimator (chi square), and the internal consistency of data reflected in standard deviations of the formation constants [
32]. Species distribution diagrams were produced using the HySS program [
33].
Isothermal titration calorimetry
The VTC2 SPX domain (residues 1-182) was expressed with a C-terminal His
6-tag and purified by Ni-affinity and size-exclusion chromatography, as previously published [
34].
Protein stocks were diluted to 300 µl with final buffer conditions: 25 mM HEPES pH 7.4, 150 mM KCl, 40 mM NaCl, 0.5 mM TCEP (ITC buffer). The exact protein concentration for each replicate was determined separately by Bradford assay. InsP8 was diluted to 500 µM in ITC buffer. For binding experiments in the presence of magnesium ions, both dilution buffers were supplemented with MgCl2 to give a final concentration of 1 mM after dilution.
ITC experiments were carried out at 25°C in a MicroCal PEAQ-ITC calorimeter (Malvern Panalytical GmbH, Germany), with ca. 50 µM protein in the cell and 500 µM ligand in the syringe. InsP8 was titrated into the solution in nineteen 2 µl-steps. Spacing between injections was 150 s.
The corresponding instrument software (MicroCal PEAQ-ITC Analysis) was used for baseline correction, peak integration, data fitting and determination of binding parameters.
DFT calculations
The input geometries were built employing the protonation and complexation patterns determined by
31P NMR in this report. Three water molecules were included in the first coordination sphere of the Mg
2+ cation. The initial geometries were pre-optimized by means of a molecular mechanic method (MM+), in order to explore the potential energy surface. Then, a Density Functional Theory (DFT) optimization protocol was performed in Gaussian 09 [
35], using the B3LYP functional, with an ultrafine integration grid and the 3-21+G* split valence basis set. The potassium ions were treated using the effective core potential LANL2DZ relativistic procedure [
36]. The solvent was modelled through the Truhlar and coworkers’ SMD solvation model [
37]. All final structures were minima in the potential energy surface, being the nature of the stationary points verified through vibrational analysis.
4. Discussion
We have conducted a systematic comparison of inositol poly- and pyrophosphates regarding their conformation at varying pH and ionic composition. Despite their structural similarities and comparable charge, clear differences between InsP
6, 5PP-InsP
5 and InsP
8 became apparent. Generally speaking, the more phosphoryl groups the molecules carry, the lower the pH range in which the conformational change takes place. It was previously demonstrated for InsP
6, that phosphoryl groups are spatially better separated in the axial conformation, minimizing steric and electrostatic repulsion between the anionic groups and providing one of the driving forces for the conformational change of higher InsPs (
Figure 1c) [
22,
23]. Additional phosphoryl groups lead to accumulation of more negative charge in the equatorial plane, which means a critical charge density is reached at lower pH, facilitating the conformational change.
Our results also confirmed previous studies, which demonstrated that metal cation coordination can promote the transition to the axial conformation at lower pH. Frost
et al. reported as early as 1979, that alkali cations stabilize the axial conformation of InsP
6 [
38]. There have since been multiple studies showing that complexation of various metals facilitates the transition to axial InsP
6 [
42,
45,
46,
47]. Hager
et al. could subsequently show that Mg
2+ ions help to stabilize 5PCP-InsP
5 in an axial conformation [
23,
38]. Notably, we now found that InsP
8 begins to change conformation at pH* 7.5 in the presence of 5 equiv. MgCl
2 - reminiscent of cytosolic conditions - which led us to focus our next efforts on this molecule.
Thermodynamic studies provided more information about the conformational equilibrium of InsP
8. The transition from the equatorial to axial conformation is an exothermic process, and the formation of the axial conformer is enthalpically favored but entropically hindered. Addition of Mg
2+ ions reduced both the enthalpic driving force and the entropic penalty, shifting the equilibrium more toward the axial side. As has been previously demonstrated for InsP
6, the phosphoryl groups are more exposed to the solvent in the axial conformation, leading to a more ordered hydration shell and overall loss of entropy compared to equatorial conformation, hence a negative ΔS [
22]. Positively charged Mg
2+ ions will coordinate to the phosphoryl groups, thereby shielding some of the negative charge, which reduces hydration, and could explain why there is less entropic penalty associated with the conformational change at higher Mg
2+ concentrations. Similarly, coordination of Mg
2+ cations might also reduce electrostatic repulsion between phosphoryl groups in the equatorial plane, thus reducing the enthalpic driving force towards the axial conformation. Of course, the forces governing the conformational equilibrium are likely more complex than the interpretation above. A plethora of different multinuclear complexes can be expected to coexist and interconvert, each with its own protonation / complexation / conformation equilibria. The parameters we determined experimentally should therefore be understood as the sum of all these processes.
Trying to make a prediction about InsP
8 in biological settings based on our results, we considered the following: Our assay conditions can only approximate cytosolic conditions, and we did not increase MgCl
2 content beyond 0.25 mM or five equivalents relative to 50 µM InsP
8, to avoid precipitation. Furthermore, the concentration of InsP
8 could not be reduced below 50 µM, without compromising the NMR detection. However, the endogenous concentration InsP
8 is thought to be around 1 µM, and cytosolic concentrations of free Mg
2+ are approximately 0.5 - 1 mM [
27,
48]. Therefore, there is a far greater excess of Mg
2+ in a cytosolic setting than in our assays. This excess should push the conformational equilibrium further toward the axial side than the ca. 30% we observed under our conditions. Nevertheless, five equivalents Mg
2+ should largely saturate InsP
8, which at pH 7.4 carries 10 – 11 negative charges (
Figure 5). For comparison, InsP
6 is also predicted to form pentamagnesium complexes under cytosolic conditions [
41,
49]. Taken together, these considerations strongly suggest that InsP
8 is actively interconverting between conformations under cytosolic conditions. Based on our experiments, we expect more than 30% of the cytosolic pool of InsP
8 to adopt axial conformation under physiological conditions.
Phosphorous NMR-studies revealed further details of the multifaceted speciation of InsP8 regarding protonation and complexation with potassium and magnesium ions. At physiological pH and in the presence of K+ and an equimolar amount of Mg2+, the most abundant species is the equatorial complex [MgK3(H3L)]6-, in which Mg2+ is coordinated to P2, P1α and P1β. Another species present at physiological pH is the axial complex [MgK3(H2L)]7- in which Mg2+ is coordinated by P5α and P1α. Overall, the pyrophosphate groups in positions 1 and 5 play an essential role in Mg2+ binding, as the ion is always coordinated to either the two pyrophosphate groups (in axial complexes) or P2 and PP1 (in equatorial complexes). Notably, all metal complexes of InsP8 are far more stable than the equivalent ones formed by InsP6 or 5PP-InsP5.
The Mg
2+ complexes of PP-InsPs were also found to be more stable than those of ATP and ADP, two other well-known cytosolic Mg
2+ chelators. Formation constants of complexes between Mg
2+ and ATP, previously reported at 150 mM NaCl and 37°C [
50] are: Mg
2+ + ATP = [MgATP]
2-: log(K) = 4.34, Mg
2+ + HATP = [Mg(HATP)]
-: log(K) = 2.39. In contrast, we measured log formation constants of 15.0 and 8.7 for [MgK
3(H
3L)] complexes of InsP
8 and 5PP-InsP
5. The question that inevitably comes to mind, is to what extent the significant stability of these complexes influences their interactions with other biomolecules. Do metal ions bridge binding interactions between PP-InsPs and proteins? And if so, is this effect more pronounced for InsP
8 than its less densely phosphorylated relatives? And how does this ultimately influence the strength and specificity of binding?
It is easy to envision signaling functions associated with the metal complexation and the conformational equilibrium of PP-InsPs. In light of their role as cellular ATP and phosphate sensors, other sensor functions seem plausible [
5,
43]. Given how sensitively PP-InsP speciation reacts to solution conditions
in vitro, it is highly likely to change upon local subcellular perturbations in pH or metal composition, which may promote selective engagement of signaling partners. Such a pH-sensing mechanism has been proven in the case of the yeast transcription factor Opi1, which is retained on the ER membrane by binding to phosphatidic acid (PA). Upon decrease of intracellular pH and protonation of the PA headgroup, Opi1 was released and activated its downstream target genes, involved in inositol biosynthesis [
51]. The authors proposed that phosphatidyl inositol lipids might sense pH through similar mechanisms, and we envision the same might be true for soluble InsPs. Interestingly, the enzymes KCS1 and PLC1, part of the inositol pyrophosphate synthesis pathway in yeast, were found in a genome-wide screen for proteins involved in intracellular pH sensing [
52].
In the large majority of cases, it is unclear by which mechanisms proteins manage to recognize and distinguish the different PP-InsPs [
6,
20]. Differential metal binding and / or a drastically altered molecular shape in the axial conformation might provide a convenient way to recognize the appropriate ligand. The idea that solution conditions, specifically the presence or absence of divalent cations, can have a tremendous influence on the outcome of
in vitro binding studies with InsPs, has already been proposed more than twenty years ago [
53]. Our ITC experiments now provide a first small hint that the axial conformation might play a role in this differentiation. We showed before that InsP
8 is partially present in axial conformation under the ITC assay conditions with Mg
2+. Dissociation constants for InsP
8 binding to the VTC2 SPX domain were almost cut in half by addition of 1 mM Mg
2+, and binding shifted towards a more enthalpically dominated interaction in the presence of Mg
2+. It is tempting to speculate that the decrease in ΔH might reflect a more specific fit of an axial magnesium complex into the binding site, compared to the free, equatorial molecule, but structural evidence would be a prerequisite to support this hypothesis and is currently unavailable.
Overall, our results highlight the immense complexity of PP-InsP speciation and how this speciation may influence their behavior. Solution conditions for PP-InsP-protein binding experiments in the literature differ widely regarding pH, salts, and other additives.
For future biochemical studies, it will be important to consider carefully, which molecular species are formed under the given assay conditions, how those conditions might affect the experimental outcome, and to what extent the conditions mirror cellular settings. An accurate understanding of this complexity will be indispensable in decoding the diverse roles of InsPs and PP-InsPs as cellular messengers.
Author Contributions
Conceptualization, D.F., L.K., N.V.; methodology, L.K., N.V., P.S.; validation, L.K., N.V.; formal analysis, N.V.; investigation, L.K.; resources, L.K., P.S.; data curation, L.K., P.S.; writing—original draft preparation, L.K., N.V., D.F.; writing—review and editing, L.K., P.S., N.V., D.F.; visualization, L.K., N.V.; supervision, D.F., N.V.; project administration, D.F., N.V.; funding acquisition, D.F., N.V. All authors have read and agreed to the published version of the manuscript.