Next Article in Journal
Glucose Regulates m6A Methylation of RNA in Pancreatic Islets
Previous Article in Journal
Golgi Metal Ion Homeostasis in Human Health and Diseases
Previous Article in Special Issue
Plasma Membrane Fluidity: An Environment Thermal Detector in Plants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Electrostatic Basis of Diacylglycerol Pyrophosphate—Protein Interaction

1
Department of Chemistry & Biochemistry, Kent State University, 800 E. Summit St., Kent, OH 44242, USA
2
School of Natural and Social Sciences, Mount Vernon Nazarene University, 800 Martinsburg Road, Mount Vernon, OH 43050, USA
3
Department of Biological Sciences, Kent State University, 800 E. Summit St., Kent, OH 44242, USA
*
Authors to whom correspondence should be addressed.
Cells 2022, 11(2), 290; https://doi.org/10.3390/cells11020290
Submission received: 14 November 2021 / Revised: 4 January 2022 / Accepted: 7 January 2022 / Published: 15 January 2022

Abstract

:
Diacylglycerol pyrophosphate (DGPP) is an anionic phospholipid formed in plants, yeast, and parasites under multiple stress stimuli. It is synthesized by the phosphorylation action of phosphatidic acid (PA) kinase on phosphatidic acid, a signaling lipid with multifunctional properties. PA functions in the membrane through the interaction of its negatively charged phosphomonoester headgroup with positively charged proteins and ions. DGPP, like PA, can interact electrostatically via the electrostatic-hydrogen bond switch mechanism but differs from PA in its overall charge and shape. The formation of DGPP from PA alters the physicochemical properties as well as the structural dynamics of the membrane. This potentially impacts the molecular and ionic binding of cationic proteins and ions with the DGPP enriched membrane. However, the results of these important interactions in the stress response and in DGPP’s overall intracellular function is unknown. Here, using 31P MAS NMR, we analyze the effect of the interaction of low DGPP concentrations in model membranes with the peptides KALP23 and WALP23, which are flanked by positively charged Lysine and neutral Tryptophan residues, respectively. Our results show a significant effect of KALP23 on the charge of DGPP as compared to WALP23. There was, however, no significant effect on the charge of the phosphomonoester of DGPP due to the interaction with positively charged lipids, dioleoyl trimethylammonium propane (DOTAP) and dioleoyl ethyl-phosphatidylcholine (EtPC). Divalent calcium and magnesium cations induce deprotonation of the DGPP headgroup but showed no noticeable differences on DGPP’s charge. Our results lead to a novel model for DGPP—protein interaction.

Graphical Abstract

1. Introduction

Anionic signaling lipids such as phosphatidic acid (PA) and the phosphoinositides (PIPs) play crucial roles in the life of a cell [1,2,3]. Their diverse functions are mediated directly by interaction with membrane proteins and/or indirectly by modulation of key membrane properties, such as membrane curvature [4,5,6,7]. Here we focus on the PA metabolite, diacylglycerol pyrophosphate (DGPP), which is formed in plants and parasites upon stress signaling but has not been found in mammalian cells [8,9,10].
PA is present in small (at most a few percent) amounts in biomembranes but increases rapidly in plants after a stress stimulus [11,12]. This is followed by a decline in PA concentration and a subsequent increase in DGPP production as PA is converted to DGPP by phosphatidic acid kinase (PAK) [10]. PAK phosphorylates the phosphomonoester of PA to a pyrophosphate and has been reported in plants such as Arabidopsis thaliana, Craterostigma plantagineum and Vicia sativa following their treatment with ABA, hyperosmotic stress and Nod factor, respectively [8,13,14]. DGPP formation is purported to not only reduce the effect of PA (levels) but also actively function in the stress response pathway [15]. The signaling function of DGPP may involve either specific interactions with proteins, where protein function is modulated (such as activation and inactivation), or through changes in membrane physical properties such as membrane curvature, bending rigidity, and charge. Additionally, the charge on the pyrophosphate headgroup of DGPP is likely to influence protein–lipid interactions [16].
The ionization properties of DGPP, as elucidated by Strawn et al. in 2012 in model membranes of phosphatidylcholine (PC) and phosphatidylethanolamine (PE), differ from those of PA [17]. The pKa2 of the phosphomonoester of DGPP is lower than that of PA due to its location in the lipid headgroup–aqueous interface (with DGPP sitting higher than PA). Additionally, PE lowered the pKa2 of this phosphate following the electrostatic-hydrogen bond switch model, which we introduced in 2007 [17,18]. Structurally, the cylindrical shape of DGPP, unlike the cone shape of PA, does not induce a negative membrane curvature (stress) [17]. Differing from PA in both shape and charge, increased DGPP production may alter the dynamics of the membrane. After conversion of PA to DGPP [10], a DGPP enriched biomembrane will have less negative curvature stress and an increased negative surface potential, which is likely to influence the interaction with (cationic) ions and proteins [17]. However, the interaction of DGPP with cationic amino acids and/or divalent cations has not been fully explored and is the main focus of this work [19,20].
Positively charged amino acid side chains, such as Arg, His and Lys, mediate the interaction and anchoring of proteins and peptides to anionic phospholipid headgroups in biomembranes [21,22]. Such electrostatic interactions are vital in signaling pathways and enable the selective targeting of proteins to different lipid domains in biomembranes [23]. NMR and MD simulations on PA and model peptides (KALP and RALP) revealed a hydrogen bond formation between the phosphomonoester of PA and peptides containing Lys and Arg sidechains, and led to the formulation of the electrostatic-hydrogen bond switch model [18]. Additionally, mean field calculations led to further support for this model [24,25]. While numerous PA binding proteins are known [6,26,27], no specific PA binding domain has been identified [28]. This is in contrast to anionic lipids such as the PIPs, for which numerous lipid binding domains are known [29], and which allow for bioinformatic studies to identify novel lipid binding proteins [30]. In the case of PA, the specific binding of proteins is likely driven by the unique physical chemical properties of this lipid, namely the ionization properties and spontaneous curvature. In the case of DGPP, only cytosolic GAPDH from Hordeum vulgare has been identified as a potential binding protein [31]. The Lys residues in the NAD+ binding site of GAPDH have been reported to be very active in numerous activities including the stable bonds formed with the phosphate group of phosphatidylserine (PS) [32], and may be responsible for the identification of GAPDH in the DGPP binding screen. The implication of GAPDH-DGPP binding has not been further explored. Further work on the identification of unique DGPP binding proteins is thus an important area of future research.
Calcium (Ca2+) and magnesium (Mg2+) are divalent cations with prominent intracellular functions [33,34]. Calcium’s role as a secondary messenger is well recognized [35], and Mg2+ is the most abundant divalent cation in living cells with a vital role in plant signaling [36]. One of their mechanistic pathways is via interaction with zwitterionic and anionic lipids in the cell to mediate important physiological activities [37,38]. Changes in cytosolic calcium concentrations have been reported in stress responses [39,40] that can affect production of minor lipids in the membrane. Additionally, changes in lipid composition can alter calcium concentration [41]. Calcium can bind simultaneously to multiple anionic binding sites to affect lipid density and packing, change the conformation, and dehydrate the lipid headgroup region of the affected membrane [42,43]. The interaction of calcium and magnesium with membrane lipids is determined, in part, by the difference in their hydration energy [44]. The calcium ion more easily releases its associated water molecules and can therefore interact more strongly with anionic membrane lipids [44]. Magnesium is known to promote the fusion of lipid vesicles, regulate ion channels and signaling pathways [45,46].
The ionization properties of PI(4,5)P2 were differentially affected by the presence of these two divalent cations [47]. Compared to magnesium, calcium had a stronger interaction with the more accessible 5-phosphate of PI(4,5)P2 leading to clustering (domain formation) of this critical membrane lipid [44]. In comparison, magnesium induced clustering of PI(4,5)P2 only occurs at much higher concentrations [44]. As seen in PI(4,5)P2, phosphate position and orientation together with the inherent properties of these divalent cations affect their binding. Margutti et al. showed that calcium affected the surface pressure of DGPP (containing) monolayers differentially at pH 5 and 8 using Langmuir monolayers [19]. No literature exists on the interaction of Mg2+ with DGPP and how this compares with Ca2+.
DGPP is expected to interact with positively charged intracellular molecules to exert its function. To gain insight into these interactions, we use solid state 31P NMR to examine the impact of KALP23 and WALP23, transmembrane α-helical peptides [48,49,50] that differ in the presence of flanking lysine and tryptophan residues respectively, on the charge of DGPP. In addition, the effect of the important divalent Ca2+ and Mg2+ cations on the charge of DGPP together with the cationic lipid dioleoyl trimethylammonium-propane (DOTAP) and ethyl-phosphatidylcholine (EtPC) were also investigated. Here we find that KALP23 altered the charge of DGPP significantly as compared to WALP23. Surprisingly, the effect of DOTAP and EtPC on the charge of DGPP were very small. This observation can be rationalized by the deprotonation of the inner phosphate of the pyrophosphate of DGPP. While Ca2+ has a significant effect on the ionization of DGPP, we do not observe a differential interaction when comparing Ca2+ and Mg2+ ions, opposite to what we previously observed for PI(4,5)P2 [44,47].

2. Materials and Methods

2.1. Materials

1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC), 1,2-dioleoyl-3-trimethylammonium-propane (DOTAP), 1,2-dioleoyl-sn-glycero-3-ethylphosphocholine (EtPC), and 1,2-dioleoyl-sn-glycero-3-pyrophosphate (DGPP) were purchased from Avanti Polar Lipids (Birmingham, AL, USA). Lipids were dissolved in a 2:1 mixture of chloroform and methanol. Lipid stocks were prepared gravimetrically and their concentrations were confirmed by phosphate assay (according to Rouser [51]). Thin layer chromatography (TLC) was used to confirm the purity of prepared lipid stocks and only stocks that showed a single spot on the TLC plate were used. CaCl2, MgCl2, and A23187 ionophore were purchased in powder form from Sigma Aldrich (St Louis, MO, USA). A23187 ionophore was dissolved in pure ethanol. KALP23 and WALP23 peptides were a generous gift from Antoinette Killian. HPLC-grade water was purchased from Fisher Scientific.

2.2. Sample Preparation

NMR samples were prepared as described in Graber and Kooijman [52]. Dry lipid films of 10–15 μmol were prepared by mixing appropriate volumes of lipid stock solution (in 2:1 chloroform:methanol) in borosilicate glass tubes. For samples containing KALP23 or WALP23, the peptides were dissolved in 2,2,2-Trifluoroethanol (TFE) and added to the mixed lipid solution at a 100:2 or 100:4 lipid to peptide ratio. TFE is used here to maintain peptide secondary structure in our organic solvent mixtures. Chloroform, methanol, and TFE (where appropriate), in lipid or lipid–peptide mixtures were evaporated under nitrogen and the resulting lipid or lipid–peptide samples were left in the borosilicate glass and further dried overnight in a high vacuum at 45 deg Celsius to remove residual traces of organic solvent. Lipid or lipid–peptide films were hydrated with 2 mL of buffer and vortexed to form multilamellar vesicle dispersions. The following buffers were used for the indicated pH ranges; 20 mM citric acid, 30 mM MES for pH 4–6.5, 50 mM HEPES for pH 6.5–8.5, and 50 mM glycine for pH 8.5–10. Buffers also contained 100 mM NaCl and 2 mM EDTA to complex divalent cations, except divalent cation containing buffers which lacked EDTA. Experiments performed at a constant pH value of 7.20 ± 0.05 used a more concentrated 100 mM HEPES buffer to minimize pH variation. For samples that contained divalent cations, the cations were added to the buffer from concentrated aqueous stock solutions. In addition, A23187 ionophore was added at a 1:1000 ionophore-to-lipid ratio. Lipid dispersions were exposed to two freeze–thaw cycles to remove metastable states. The freeze–thaw procedure consisted of quickly freezing the lipid suspension in a mixture of ethanol and dry ice, and then thawing in warm water while occasionally vortexing the sample. The pH of each NMR sample was carefully measured using a Sentron Intelli pH probe and standard pH meter, and this measured bulk pH was used to prepare the pH titration curves. Lipid samples were spun down at 15,000 rpm at 4 °C for 45–60 min to create a concentrated lipid pellet. The lipid pellets were loaded into 4 mm zirconium MAS NMR sample tubes for NMR spectroscopy.

2.3. NMR Spectroscopy

31P NMR spectra were recorded as previously described [52]. An 85% H3PO4 standard solution was run before each sample solution to act as an external reference. Sample spectra were recorded under stable spinning conditions (5 kHz) at 21 ± 1.0 °C. Generally, 15,000 to 50,000 scans were collected. We only determined the upper part (pKa2) of the full titration curve for the systems explored here. We chose this to save precious peptide and lipid, and we previously showed that determining only the ionization behavior in the physiological pH range does not significantly affect the pKa2 value [53]. After the MAS spectra was recorded, the phase behavior of the lipid sample was determined via static experiment for select samples. The static experiments were recorded on the same probe and aided by low-power proton decoupling (the spinal 64 pulse program was used for proton decoupling). For each static experiment, 20,000 to 150,000 scans were collected.

2.4. Charge Determination from Titration Curve Data

The charge of the phosphomonoester of DGPP was determined by assuming that the chemical shift value is a weighted average of the chemical shifts of the singly and completely deprotonated state of the respective phosphate group [54]. The chemical shift can then be expressed as:
δ i o b s = f i , p   δ i , p + ( 1 f i , p   ) δ i , d
where δ i o b s is the pH dependent chemical shift of the phosphomonoester group i, while δi,p and δi,d are chemical shifts of the singly protonated and completely deprotonated phosphate group. Equation (1) can be re-arranged to
f i , p = δ i o b s δ i , d δ i , p   δ i , d  
where fi,p is the degree of protonation of the phosphate group i. The pKa2 value for the phosphomonoester of DGPP phosphate was determined based on a Henderson–Hasselbalch derived fit, as in Kooijman et al. [53], using the relation:
δ = δ i , p 10 p K a p H + δ i , d 1 + 10 p K a p H

2.5. Statistical Analysis

Experiments for DGPP’s interaction with KALP23, WALP23 and the control were performed at least three times. Those for DGPP’s interaction with divalent cations together with its control were performed at least 7 times. These were then used to calculate the means, standard deviations (SD), and standard errors of means (SEM). SigmaPlot (Systat Software Inc, San Jose, CA, USA) was used to plot the graphs using the means and SD. We performed one-way analysis of variance (ANOVA) followed by Dunnett’s test. p values < 0.001 were considered statistically significant.

3. Results and Discussion

The pyrophosphate headgroup of DGPP is unusual in membrane phospholipids, and we previously showed that in model membranes, only the phosphomonoester of DGPP carries charge, while the inner phosphate remains largely protonated [17]. We also showed that the phosphomonoester of DGPP follows the electrostatic hydrogen bond switch model, as PE increases the charge (decreases the pKa). However, how the pyrophosphate headgroup of DGPP interacts with cationic moieties in proteins or with divalent cations has only received limited attention. Here we investigate the ionization of DGPP in the presence of transmembrane peptides that are anchored in the membrane by specific flanking residues (see Figure 1A). Additionally, we study the effect of two crucial divalent cations on the charge (i.e., chemical shift) of DGPP.

3.1. Lysine Residues Interact Strongly with the Pyrophosphate Headgroup of DGPP and Increase Its Charge

Many peripheral membrane proteins (and some integral proteins as well) rely on cationic domains containing lysine and/or arginine residues to interact with anionic phospholipids. These interactions are often considered to be electrostatic in nature, but we have extensively shown that hydrogen bond formation is also an important driver of membrane recruitment and localization, especially in the case of PA [18,24].
The electrostatic-hydrogen bond switch model sums up these interactions for PA (and any phosphomonoester group for that matter) [18]. The anionic membrane potential recruits cationic proteins to the membrane. Depending on the potential, and other interactions (hydrogen bonding, hydrophobic, etc.), the protein then either binds or samples a different patch of the membrane. When a PA-binding protein binds to PA, it is able to hydrogen bond to the phosphomonoester headgroup. This leads to an increase in the negative charge of PA from −1 to −2, further strengthening the electrostatic interaction and thus docking of the protein onto the membrane. Both PA and DGPP contain a phosphomonoester in their headgroup, and DGPP may therefore use a similar mechanism to interact with membrane proteins.
To investigate this, we used the transmembrane peptide KALP23 as a model for cationic protein domain interaction with DGPP. The membrane spanning α-helix of KALP23 [48,49,50] consists of a repeating alanine-leucine motif and is flanked on either side by two lysine residues (see Figure 1A). The hydrophobic alanine-leucine repeat and positively charged flanking residues ensure the transmembrane orientation of the peptide and position the lysine residues in the right location to interact with the pyrophosphate of DGPP [48,55]. As a control, we used WALP23, which is flanked by two tryptophan residues on either side of the peptide (see Figure 1A). While the tryptophan residues anchor the peptide in a transmembrane orientation [48], it is not charged, and hence WALP23 cannot interact with DGPP through electrostatic interactions. The chemical shift (CS) of the DGPP phosphates was monitored for two concentrations of peptide (see Figure 1B,C). There is no significant change in CS for the mixtures containing WALP23, indicating that the presence of the peptide in the membrane does not affect the ionization of DGPP. When KALP23 is added to the membrane, we observe a concentration dependent increase (to more downfield values) for the CS of the DGPP phosphomonoester, indicating increased deprotonation. Additionally, we observe a significant broadening of the phosphomonoester peak of DGPP.
To further investigate this effect, the ionization of DGPP in the presence of KALP23 was observed over the pH range from 4 to 10 (see Figure 2A). The peaks were assigned as in our previous work, with the PC peak furthest downfield followed by the phosphomonoester peak and finally the phosphodiester peak. Figure 2B shows static 31P NMR spectra for selected pH values that show our model systems form lipid bilayers, characterized by a high field peak, and low field shoulder, under our experimental conditions (also see Figure 1D). As the pH increases the DGPP phosphomonoester peak shifts downfield as the phosphomonoester is deprotonated. As observed in the constant pH measurements, we observe that around pH 6.1–6.5, the phosphomonoester peak significantly broadens but does not disappear completely, allowing for a pH titration curve to be constructed. This significant broadening of the phosphomonoester peak of DGPP in the presence of KALP23 (Figure 2A) is not observed in the case of PA. The largest contribution to NMR peak line width is from T1 and T2 relaxation, which are strongly dependent on the dynamics of the nuclei [56,57,58]. In this case, the increased linewidth suggests that the lysine residues of KALP23 bind tightly to the pyrophosphate of DGPP resulting in altered relaxation dynamics. This tight binding is not observed in the case of PA as indicated by a lack of line broadening, and is likely to have important implications for the interaction of DGPP with peripheral membrane proteins.
This observation suggests that PA may have weaker electrostatic and hydrogen bond interactions with membrane proteins compared to DGPP (see further discussion below).
By plotting the CS values from Figure 2A with respect to pH, we obtain an ionization curve. This ionization curve is compared against our previous data for 5 mol% of DGPP in 95 mol% DOPC (Figure 2C). In the presence of KALP23 there is a significant shift of the ionization curve to the left, indicating a lower pKa2 value and hence higher DGPP charge (at constant pH). Based on a non-linear fit using Equation (3) (derived from the Henderson–Hasselbalch equation and the assumption that the observed chemical shift values are weighted averages of the chemical shifts of the singly dissociated and doubly dissociated states) we determine a pKa2 value of 6.38 ± 0.04 (see Table 1). This is a shift of 1.06 ± 0.06 pH values (logarithmic scale) with respect to the mixture in the absence of KALP23 (pKa 7.44 ± 0.02) [17]. This shift is larger compared to what we observed for PA in the presence of KALP23 (0.78 ± 0.09 pH values) [18], and suggests that DGPP interacts with KALP23 through a similar electrostatic-hydrogen bond switch mechanism. Interestingly, the interaction between the primary amine in the lysine residues of KALP23 is significantly “stronger” than those of PE. For a more than 2-fold larger concentration of PE compared to the lysines in KALP23, we find a pKa2 shift of only 0.73 ± 0.04 [17] (similar to what we observe for the PA-KALP23 interaction), compared to 1.06 ± 0.06 for KALP23. This observation shows that both hydrogen bond formation, and overall cationic charge (PE is zwitterionic) are important for this interaction, and that lysine residues interact “more strongly” with a pyrophosphate (as in DGPP) compared to a phosphomonoester (as in PA).

3.2. Positive Charge Alone Does Not Significantly Affect the Charge of the Phosphomonoester of DGPP

In the previous section we looked at the interaction of DGPP with the cationic peptide KALP23. While the increased deprotonation of DGPP in the presence of KALP23 suggests the formation of a hydrogen bond, it is also possible that KALP23 exerts its effect on DGPP ionization simply via its effect on the membrane potential. The cationic charge contributed by the lysine residues reduces the negative membrane potential and is likely to have a significant impact on ionizable groups at the membrane. The negative surface potential (due to the presence of DGPP) of the membrane lowers the pH at the membrane interface compared to the bulk pH. The presence of positively charged proteins or peptides at the membrane will decrease this negative membrane surface potential. To examine the effect of altering the membrane potential, the ionization of DGPP was examined in the presence of the cationic lipids DOTAP and EtPC.
Figure 3A shows the titration curve for the phosphomonoester of DGPP in a 95:5:16 mixture of PC:DGPP:DOTAP compared to same mixture without DOTAP (previous data, [17]). Since KALP23 contains four positive charges per peptide, while DOTAP has only one, we used a lipid-to-cation (DOTAP) ratio of 100:16 to be able to directly compare the KALP23 and cationic lipid data. The full 31P NMR titration curves are shown in supplementary data Figure S1. To consistently evaluate the significance of changes in pKa, we also calculate the degree of protonation and the charge at pH 7.2 for each of our titration data, see Figure S4, and Table 1.
Figure 3B shows the titration curves for the phosphomonoester of DGPP in a 95:5:16 and 95:5:32 mixture of PC:DGPP:EtPC. Again, the 31P NMR spectra as a function of pH are shown in the supplementary data, Figures S2 and S3. The protonation curves are shown in Figure S4. Interestingly, we now observe a small effect on the pKa2 and charge of the phosphomonoester of DGPP. For the 95:5:16 mixture the pKa2 is reduced to 7.25 compared to 7.44 for DGPP in a pure PC matrix, and doubling the amount of EtPC results in a further shift of the pKa2 to 7.14 (see Table 1). The effect on charge of the phosphomonoester is not significant for the 95:5:16 lipid mixture but becomes significant (difference of 0.15) for the mixture containing twice as much EtPC. The overall impact on the pKa2 of DGPP is thus surprisingly small compared to the large impact observed in the presence of KALP23 (where the pKa decreased by 1.06 pH units).

3.3. Proposed Model for the Effect of Positive Charge on the Ionization of DGPP

How do we reconcile the observed ionization behavior of the pyrophosphate of DGPP with basic electrostatics? Physics tells us that addition of positive charge to a surface should impact the surface potential and thus the ionization behavior of ionizable groups at this interface. Exactly as we have previously observed for PA [18]. Here, no such effect on the ionization of the DGPP phosphomonoester is observed. The most plausible explanation for this unusual behavior is the following: In model membranes containing only zwitterionic and anionic membrane lipids, the inner phosphate of the pyrophosphate headgroup remains protonated, as discussed previously [17]. Addition of a positive charge, right next to this pyrophosphate due to DOTAP and EtPC, now removes this proton from the inner phosphate resulting in an additional negative charge of DGPP. However, this negative charge is balanced by the positive charge on the DOTAP or EtPC molecule, resulting in no, or only very small local changes in the electrostatic membrane potential. A cartoon of this model is shown in Figure 4. This also explains the strong interaction of the lysine residues in KALP23 with the pyrophosphate headgroup. Not only is the charge in the phosphomonoester increased (as in the case of PA), but the inner phosphate now also becomes deprotonated, leading to a further increase in the total charge of DGPP (see Table 1). The almost −3 charge of DGPP (coupled to the hydrogen bonding with lysine residues) compared to the charge of −2 for PA also explains why we observe partial immobilization of the pyrophosphate headgroup as evidenced by the line broadening (see Figure 2).

3.4. Implications for DGPP—Protein Interactions

Our results have important implications for how DGPP may interact with membrane proteins and other membrane lipids. In the following discussion we will focus primarily on protein–lipid interactions for DGPP. Our observation of 31P NMR line width broadening coupled with the large shift in the titration curve due to KALP23 suggests that lysine residues have a much stronger interaction with the pyrophosphate of DGPP than the single phosphomonoester of PA. The primary amine of lysine likely deprotonates the entire pyrophosphate resulting in a total charge of DGPP of −3 (see Table 1), unlike the total charge for PA of −2. Couple this much larger charge with hydrogen bond formation (as summarized in the electrostatic-hydrogen bond switch model) and you get a much stronger binding interaction (higher binding constant) for lysine and presumably arginine residues with DGPP than PA [18]. It is interesting to note that PA binding proteins do not bind more strongly to DGPP as a result of the increased interaction between cationic and hydrogen bond forming amino acids and the pyrophosphate headgroup. We showed this previously for a random selection of PA binding proteins where we found either similar or significantly less binding to DGPP compared to PA [60]. This makes sense when we consider another critical lipid physical chemical property, namely spontaneous curvature. We have extensively shown that PA binding proteins require negative membrane curvature (or in other words lipids with negative spontaneous curvature such as PE or PA itself) to efficiently bind to a lipid membrane. DGPP, however, is a lipid with essentially zero spontaneous curvature [17], and hence PA binding proteins, while having a stronger electrostatic and hydrogen bond interaction with its headgroup, do not bind more efficiently to DGPP over PA. In the case of PA signaling in plants or parasites this means that conversion of PA, formed due to a biotic or abiotic signal, to DGPP may result in an attenuation of the PA signal. Fine tuning the amount of PA and DGPP in a membrane could thus be a sensitive regulator for protein membrane localization. One important avenue for future research is to carefully investigate whether there are indeed specific DGPP binding proteins, and then to explore in detail the interaction of those proteins with both PA and DGPP. Only one published study has looked at DGPP binding proteins to date [31].

3.5. The interaction of DGPP with Divalent Cations

Divalent cations have the potential to interact with anionic phospholipids and alter phospholipid ionization due to changes in membrane potential. Previous work with the anionic phospholipid PI(4,5)P2 has shown that calcium interacts strongly with the phosphomonoesters in its headgroup (specifically the 5-phosphate), while magnesium interacts much more weakly [47]. The significant difference in affinity was suggested to be due to differences in hydration energy and ion size [44].
Here we investigated the ionization of DGPP in the presence of varying concentrations of calcium and magnesium ions. The 31P chemical shift of the DGPP phosphomonoester was determined in DOPC/DGPP (95%/5%) MLVs at pH 7.2 ± 0.1 and 100 mM NaCl (see Figure 5). As the calcium concentration increases, the chemical shift for the DGPP phosphomonoester increases, indicating increased deprotonation. For calcium levels below 0.1 mM there is no significant effect, but for 1 mM calcium there is a change of more than 2 ppm, suggesting almost complete deprotonation of the phosphomonoester group. The relative ionization of DGPP can be estimated based on the change in chemical shift. The estimated DGPP charge is shown in Table 2. Please note that this charge assumes that the presence of Ca2+ or Mg2+ does not by itself impact the DGPP chemical shift (this assumption proved valid for PIP2 [47]). The DGPP phosphodiester is unaffected by the calcium content. Interestingly, with 1 mM Mg2+ there is as large an effect as for 1 mM Ca2+. This suggests that Ca2+ and Mg2+ bind to DGPP with equally high affinity. This is notably different from cation binding to PIP2, where a strong preference was shown for Ca2+ and indicates that the DGPP binding site is unique compared to PIP2 and can accept both ions without preference [44].

4. Conclusions

In this work, we have characterized the interaction between DGPP and several cationic membrane components. We found that the charge of DGPP is affected by the cellular cations Ca2+ and Mg2+, as expected. Unlike for PI(4,5)P2, however, DGPP is not specifically sensitive to Ca2+, but interacts with Ca2+ or Mg2+ to a similar degree. At least as far as our charge analysis can differentiate. These data also do not show a specific binding site for Ca2+, but additional experiments (e.g., X-ray fluorescence from Langmuir monolayers) [61] are needed to fully answer this question.
Based upon our investigation of the DGPP-KALP23 interaction, we suggest that DGPP interacts with membrane proteins in a similar manner to PA, i.e., via the electrostatic-hydrogen bond switch model. However, the interaction of Lys residues with DGPP is significantly stronger than that with PA. This may suggest that PA and DGPP compete for the same protein binding sites. However, PA binding proteins do not bind better to DGPP as we have shown and hence DGPP formation is likely to result in an attenuation of PA signaling. An unanswered question is how PA and DGPP interact when present in the same bilayer leaflet. Monolayer work suggests that PA and DGPP may form domains [19] but this needs to be corroborated in model lipid bilayers. It is tempting to speculate however that fine tuning of the concentration of PA and DGPP in a given membrane may thus be one mechanism by which to sensitively control PA—Protein interactions and lipid signaling. Taken together, our results give us important insights into how DGPP interacts with critical cellular components.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cells11020290/s1, Figure S1: 31P MAS NMR spectra as a function of pH for 5 mol% DGPP in DOPC/DOTAP (95:16) vesicles.; Figure S2: 31P MAS NMR spectra as a function of pH for 5 mol% DGPP in DOPC/EtPC (95:16) vesicles; Figure S3: 31P MAS NMR spectra as a function of pH for 5 mol% DGPP in DOPC/EtPC (95:32) vesicles; Figure S4: (A) Degree of protonation (fp) as a function of pH for DGPP in the presence of KALP23 (green) and DOTAP (blue) is compared for 5mol % DGPP in 95%DOPC vesicles in grey. (B) Degree of protonation (fp) as a function of pH for DGPP in the presence EtPC 16 (light red) and EtPC 32 (deep red) is compared for 5 mol% DGPP in 95%DOPC vesicles in grey.

Author Contributions

Conceptualization, D.O.K., Z.G. and E.E.K.; methodology, D.O.K., Z.G. and E.E.K.; formal analysis, S.M.L., Y.K., H.K., J.W.M., D.O.K. and Z.G.; resources, E.E.K.; data curation, D.O.K. and Z.G.; writing—original draft preparation, D.O.K., Z.G. and E.E.K.; writing—review and editing D.O.K., Z.G. and E.E.K.; visualization, D.O.K., Z.G. and E.E.K.; supervision, E.E.K.; project administration, E.E.K.; funding acquisition, E.E.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by the National Science Foundation under Grant No. CHE-1412920. D.O.K. and E.E.K. acknowledge financial support from the Department of Biological Sciences at Kent State University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the article and its Supplementary Data.

Acknowledgments

The authors would like to thank Mahinda Gangoda for technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, X.; Devaiah, S.P.; Zhang, W.; Welti, R. Signaling functions of phosphatidic acid. Prog. Lipid Res. 2006, 45, 250–278. [Google Scholar] [CrossRef]
  2. De Craene, J.O.; Bertazzi, D.L.; Bär, S.; Friant, S. Phosphoinositides, major actors in membrane trafficking and lipid signaling pathways. Int. J. Mol. Sci. 2017, 18, 634. [Google Scholar] [CrossRef]
  3. Balla, T. Phosphoinositides: Tiny lipids with giant impact on cell regulation. Physiol. Rev. 2013, 93, 1019–1137. [Google Scholar] [CrossRef]
  4. Kooijman, E.E.; Chupin, V.; de Kruijff, B.; Burger, K.N.J. Modulation of Membrane curvature by phosphatidic acid and lysophosphatidic acid. Traffic 2003, 4, 162–174. [Google Scholar] [CrossRef] [PubMed]
  5. Schink, K.O.; Tan, K.W.; Stenmark, H. Phosphoinositides in Control of Membrane Dynamics. Annu. Rev. Cell Dev. Biol. 2016, 32, 143–171. [Google Scholar] [CrossRef]
  6. Stace, C.L.; Ktistakis, N.T. Phosphatidic acid- and phosphatidylserine-binding proteins. Biochim. Biophys. Acta-Mol. Cell Biol. Lipids 2006, 1761, 913–926. [Google Scholar] [CrossRef] [PubMed]
  7. Shukla, S.; Jin, R.; Robustelli, J.; Zimmerman, Z.E.; Baumgart, T. PIP2 Reshapes Membranes through Asymmetric Desorption. Biophys. J. 2019, 117, 962–974. [Google Scholar] [CrossRef]
  8. Den Hartog, M.; Musgrave, A.; Munnik, T. Nod factor-induced phosphatidic acid and diacylglycerol pyrophosphate formation: A role for phospholipase C and D in root hair deformation. Plant J. 2001, 25, 55–65. [Google Scholar] [CrossRef] [PubMed]
  9. Marchesini, N.; Santander, V.; Machado-Domenech, E. Diacylglycerol pyrophosphate: A novel metabolite in the Trypanosoma cruzi phosphatidic acid metabolism. FEBS Lett. 1998, 436, 377–381. [Google Scholar] [CrossRef]
  10. Munnik, T.; De Vrije, T.; Irvine, R.F.; Musgrave, A. Identification of diacylglycerol pyrophosphate as a novel metabolic product of phosphatidic acid during G-protein activation in plants. J. Biol. Chem. 1996, 271, 15708–15715. [Google Scholar] [CrossRef]
  11. Testerink, C.; Munnik, T. Phosphatidic acid: A multifunctional stress signaling lipid in plants. Trends Plant Sci. 2005, 10, 368–375. [Google Scholar] [CrossRef]
  12. Vilchez, A.C.; Peppino Margutti, M.; Reyna, M.; Wilke, N.; Villasuso, A.L. Recovery from chilling modulates the acyl-editing of phosphatidic acid molecular species in barley roots (Hordeum vulgare L.). Plant Physiol. Biochem. 2021, 167, 862–873. [Google Scholar] [CrossRef] [PubMed]
  13. Zalejski, C.; Zhang, Z.; Quettier, A.L.; Maldiney, R.; Bonnet, M.; Brault, M.; Demandre, C.; Miginiac, E.; Rona, J.P.; Sotta, B.; et al. Diacylglycerol pyrophosphate is a second messenger of abscisic acid signaling in Arabidopsis thaliana suspension cells. Plant J. 2005, 42, 145–152. [Google Scholar] [CrossRef] [PubMed]
  14. Munnik, T.; Meijer, H.J.G.; Ter Riet, B.; Hirt, H.; Frank, W.; Bartels, D.; Musgrave, A. Hyperosmotic stress stimulates phospholipase D activity and elevates the levels of phosphatidic acid and diacylglycerol pyrophosphate. Plant J. 2000, 22, 147–154. [Google Scholar] [CrossRef]
  15. van Schooten, B.; Testerink, C.; Munnik, T. Signalling diacylglycerol pyrophosphate, a new phosphatidic acid metabolite. Biochim. Biophys. Acta-Mol. Cell Biol. Lipids 2006, 1761, 151–159. [Google Scholar] [CrossRef]
  16. Putta, P.; Creque, E.; Piontkivska, H.; Kooijman, E.E. Lipid-protein interactions for ECA1 an N-ANTH domain protein involved in stress signaling in plants. Chem. Phys. Lipids 2020, 231, 104919. [Google Scholar] [CrossRef]
  17. Strawn, L.; Babb, A.; Testerink, C.; Kooijman, E.E. The physical chemistry of the enigmatic phospholipid diacylglycerol pyrophosphate. Front. Plant Sci. 2012, 3, 40. [Google Scholar] [CrossRef]
  18. Kooijman, E.E.; Tieleman, D.P.; Testerink, C.; Munnik, T.; Rijkers, D.T.S.; Burger, K.N.J.; De Kruijff, B. An electrostatic/hydrogen bond switch as the basis for the specific interaction of phosphatidic acid with proteins. J. Biol. Chem. 2007, 282, 11356–11364. [Google Scholar] [CrossRef] [PubMed]
  19. Peppino Margutti, M.; Wilke, N.; Villasuso, A.L. Influence of Ca2+ on the surface behavior of phosphatidic acid and its mixture with diacylglycerol pyrophosphate at different pHs. Chem. Phys. Lipids 2020, 228, 104887. [Google Scholar] [CrossRef]
  20. Villasuso, A.L.; Wilke, N.; Maggio, B.; Machado, E. Zn2+-dependent surface behavior of diacylglycerol pyrophosphate and its mixtures with phosphatidic acid at different pHs. Front. Plant Sci. 2014, 5, 104887. [Google Scholar] [CrossRef]
  21. Liu, A.; Wenzel, N.; Qi, X. Role of lysine residues in membrane anchoring of saposin C. Arch. Biochem. Biophys. 2005, 443, 101–112. [Google Scholar] [CrossRef]
  22. Li, L.; Vorobyov, I.; Allen, T.W. The different interactions of lysine and arginine side chains with lipid membranes. J. Phys. Chem. B 2013, 117, 11906–11920. [Google Scholar] [CrossRef]
  23. Stahelin, R.V. Lipid binding domains: More than simple lipid effectors. J. Lipid Res. 2009, 50, S299–S304. [Google Scholar] [CrossRef] [PubMed]
  24. Loew, S.; Kooijman, E.E.; May, S. Increased pH-sensitivity of protein binding to lipid membranes through the electrostatic-hydrogen bond switch. Chem. Phys. Lipids 2013, 169, 9–18. [Google Scholar] [CrossRef] [PubMed]
  25. Mengistu, D.H.; Kooijman, E.E.; May, S. Ionization properties of mixed lipid membranes: A Gouy-Chapman model of the electrostatic-hydrogen bond switch. Biochim. Biophys. Acta-Biomembr. 2011, 1808, 1985–1992. [Google Scholar] [CrossRef] [PubMed]
  26. Jang, J.H.; Lee, C.S.; Hwang, D.; Ryu, S.H. Understanding of the roles of phospholipase D and phosphatidic acid through their binding partners. Prog. Lipid Res. 2012, 51, 71–81. [Google Scholar] [CrossRef]
  27. Zegarlińska, J.; Piaścik, M.; Sikorski, A.F.; Czogalla, A. Phosphatidic acid-A simple phospholipid with multiple faces. Acta Biochim. Pol. 2018, 65, 163–171. [Google Scholar] [CrossRef]
  28. Tanguy, E.; Kassas, N.; Vitale, N. Protein–phospholipid interaction motifs: A focus on phosphatidic acid. Biomolecules 2018, 8, 20. [Google Scholar] [CrossRef]
  29. Stahelin, R.V.; Scott, J.L.; Frick, C.T. Cellular and molecular interactions of phosphoinositides and peripheral proteins. Chem. Phys. Lipids 2014, 182, 3–18. [Google Scholar] [CrossRef]
  30. Dowler, S.; Currie, R.A.; Campbell, D.G.; Deak, M.; Kular, G.; Downes, C.P.; Alessi, D.R. Identification of pleckstrin-homology-domain-containing proteins with novel phosphoinositide-binding specificities. Biochem. J. 2000, 351, 19–31. [Google Scholar] [CrossRef]
  31. Astorquiza, P.L.; Usorach, J.; Racagni, G.; Villasuso, A.L. Diacylglycerol pyrophosphate binds and inhibits the glyceraldehyde-3-phosphate dehydrogenase in barley aleurone. Plant Physiol. Biochem. 2016, 101, 88–95. [Google Scholar] [CrossRef]
  32. Sirover, M.A. Structural analysis of glyceraldehyde-3-phosphate dehydrogenase functional diversity. Int. J. Biochem. Cell Biol. 2014, 57, 20–26. [Google Scholar] [CrossRef] [PubMed]
  33. Takahashi, A.; Camacho, P.; Lechleiter, J.D.; Herman, B. Measurement of intracellular calcium. Physiol. Rev. 1999, 79, 1089–1125. [Google Scholar] [CrossRef]
  34. Shaul, O. Magnesium transport and function in plants: The tip of the iceberg. BioMetals 2002, 15, 309–323. [Google Scholar] [CrossRef] [PubMed]
  35. Dodd, A.N.; Kudla, J.; Sanders, D. The language of calcium signaling. Annu. Rev. Plant Biol. 2010, 61, 593–620. [Google Scholar] [CrossRef] [PubMed]
  36. Kleczkowski, L.A.; Igamberdiev, A.U. Magnesium signaling in plants. Int. J. Mol. Sci. 2021, 22, 1159. [Google Scholar] [CrossRef]
  37. Papahadjopoulos, D.; Vail, W.J.; Newton, C.; Nir, S.; Jacobson, K.; Poste, G.; Lazo, R. Studies on membrane fusion. III. The role of calcium-induced phase changes. Biochim. et Biophys. Acta (BBA)-Biomembr. 1977, 465, 579–598. [Google Scholar] [CrossRef]
  38. Binder, H.; Zschörnig, O. The effect of metal cations on the phase behavior and hydration characteristics of phospholipid membranes. Chem. Phys. Lipids 2002, 115, 39–61. [Google Scholar] [CrossRef]
  39. Stael, S.; Wurzinger, B.; Mair, A.; Mehlmer, N.; Vothknecht, U.C.; Teige, M. Plant organellar calcium signalling: An emerging field. J. Exp. Bot. 2012, 63, 1525–1542. [Google Scholar] [CrossRef]
  40. Kader, M.A.; Lindberg, S. Cytosolic calcium and pH signaling in plants under salinity stress. Plant Signal. Behav. 2010, 5, 233–238. [Google Scholar] [CrossRef]
  41. Tang, D.; Dean, W.L.; Borchman, D.; Paterson, C.A. The influence of membrane lipid structure on plasma membrane Ca2+-ATPase activity. Cell Calcium 2006, 39, 209–216. [Google Scholar] [CrossRef]
  42. Garidel, P.; Blume, A.; Hübner, W. A Fourier transform infrared spectroscopic study of the interaction of alkaline earth cations with the negatively charged phospholipid 1,2-dimyristoyl-sn-glycero-3-phosphoglycerol. Biochim. Biophys. Acta (BBA)-Biomembr. 2000, 1466, 245–259. [Google Scholar] [CrossRef]
  43. Tsai, H.H.G.; Lai, W.X.; Lin, H.D.; Lee, J.B.; Juang, W.F.; Tseng, W.H. Molecular dynamics simulation of cation-phospholipid clustering in phospholipid bilayers: Possible role in stalk formation during membrane fusion. Biochim. Biophys. Acta-Biomembr. 2012, 1818, 2742–2755. [Google Scholar] [CrossRef]
  44. Wang, Y.H.; Collins, A.; Guo, L.; Smith-Dupont, K.B.; Gai, F.; Svitkina, T.; Janmey, P.A. Divalent cation-induced cluster formation by polyphosphoinositides in model membranes. J. Am. Chem. Soc. 2012, 134, 3387–3395. [Google Scholar] [CrossRef] [PubMed]
  45. Düzgüne, N.; Nir, S.; Wilschut, J.; Bentz, J.; Newton, C.; Portis, A.; Papahadjopoulos, D. Calcium- and magnesium-induced fusion of mixed phosphatidylserine/phosphatidylcholine vesicles: Effect of ion binding. J. Membr. Biol. 1981, 59, 115–125. [Google Scholar] [CrossRef]
  46. Martín-Molina, A.; Rodríguez-Beas, C.; Faraudo, J. Effect of calcium and magnesium on phosphatidylserine membranes: Experiments and all-atomic simulations. Biophys. J. 2012, 102, 2095–2103. [Google Scholar] [CrossRef]
  47. Graber, Z.T.; Gericke, A.; Kooijman, E.E. Phosphatidylinositol-4,5-bisphosphate ionization in the presence of cholesterol, calcium or magnesium ions. Chem. Phys. Lipids 2014, 182, 62–72. [Google Scholar] [CrossRef] [PubMed]
  48. De Planque, M.R.R.; Kruijtzer, J.A.W.; Liskamp, R.M.J.; Marsh, D.; Greathouse, D.V.; Koeppe, R.E.; De Kruijff, B.; Killian, J.A. Different membrane anchoring positions of tryptophan and lysine in synthetic transmembrane α-helical peptides. J. Biol. Chem. 1999, 274, 20839–20846. [Google Scholar] [CrossRef]
  49. De Planque, M.R.R.; Goormaghtigh, E.; Greathouse, D.V.; Koeppe, R.E.; Kruijtzer, J.A.W.; Liskamp, R.M.J.; De Kruijff, B.; Killian, J.A. Sensitivity of single membrane-spanning α-helical peptides to hydrophobic mismatch with a lipid bilayer: Effects on backbone structure, orientation, and extent of membrane incorporation. Biochemistry 2001, 40, 5000–5010. [Google Scholar] [CrossRef]
  50. Strandberg, E.; Morein, S.; Rijkers, D.T.S.; Liskamp, R.M.J.; Van der Wel, P.C.A.; Antoinette Killian, J. Lipid dependence of membrane anchoring properties and snorkeling behavior of aromatic and charged residues in transmembrane peptides. Biochemistry 2002, 41, 7190–7198. [Google Scholar] [CrossRef]
  51. Rouser, G.; Fleischer, S.; Yamamoto, A. Two dimensional thin layer chromatographic separation of polar lipids and determination of phospholipids by phosphorus analysis of spots. Lipids 1970, 5, 494–496. [Google Scholar] [CrossRef] [PubMed]
  52. Graber, Z.T.; Kooijman, E.E. Ionization behavior of polyphosphoinositides determined via the preparation of pH titration curves using solid-state 31P NMR. Methods Mol. Biol. 2013, 1009, 129–142. [Google Scholar] [CrossRef]
  53. Kooijman, E.E.; Carter, K.M.; Van Laar, E.G.; Chupin, V.; Burger, K.N.J.; De Kruijff, B. What makes the bioactive lipids phosphatidic acid and lysophosphatidic acid so special? Biochemistry 2005, 44, 17007–17015. [Google Scholar] [CrossRef]
  54. Appleton, T.G.; Hall, J.R.; Ralph, S.F.; Thompson, C.S.M. NMR study of acid-base equilibria and other reactions of ammineplatinum complexes with aqua and hydroxo ligands. Inorg. Chem. 1989, 28, 1989–1993. [Google Scholar] [CrossRef]
  55. De Planque, M.R.R.; Killian, J.A. Protein-lipid interactions studied with designed transmembrane peptides: Role of hydrophobic matching and interfacial anchoring (Review). Mol. Membr. Biol. 2003, 20, 271–284. [Google Scholar] [CrossRef]
  56. Lu, J.X.; Damodaraiv, K.; Blazyk, J.; Lorigan, G.A. Solid-state nuclear magnetic resonance relaxation studies of the interaction mechanism of antimicrobial peptides with phospholipid bilayer membranes. Biochemistry 2005, 44, 10208–10217. [Google Scholar] [CrossRef]
  57. Dufourc, E.J.; Mayer, C.; Stohrer, J.; Althoff, G.; Kothe, G. Dynamics of phosphate head groups in biomembranes. Comprehensive analysis using phosphorus-31 nuclear magnetic resonance lineshape and relaxation time measurements. Biophys. J. 1992, 61, 42–57. [Google Scholar] [CrossRef]
  58. Cullis, P.R.; De Kruyff, B.; Richards, R.E. Factors affecting the motion of the polar headgroup in phospholipid bilayers. A 31P NMR study of unsonicated phosphatidylcholine liposomes. BBA-Biomembr. 1976, 426, 433–446. [Google Scholar] [CrossRef]
  59. Kooijman, E.E.; King, K.E.; Gangoda, M.; Gericke, A. Ionization Properties of Phosphatidylinositol Polyphosphates in Mixed Model Membranes. Biochemistry 2009, 48, 9360–9371. [Google Scholar] [CrossRef] [PubMed]
  60. Putta, P.; Rankenberg, J.; Korver, R.A.; van Wijk, R.; Munnik, T.; Testerink, C.; Kooijman, E.E. Phosphatidic acid binding proteins display differential binding as a function of membrane curvature stress and chemical properties. Biochim. Biophys. Acta-Biomembr. 2016, 1858, 2709–2716. [Google Scholar] [CrossRef]
  61. Graber, Z.T.; Wang, W.; Singh, G.; Kuzmenko, I.; Vaknin, D.; Kooijman, E.E. Competitive cation binding to phosphatidylinositol-4,5-bisphosphate domains revealed by X-ray fluorescence. RSC Adv. 2015, 5, 106536–106542. [Google Scholar] [CrossRef]
Figure 1. Effect of transmembrane peptides on the chemical shift of DGPP. (A) Ribbon models for the transmembrane Alanine and Leucine (AL) peptides KALP23 and WALP23 with flanking lysine and tryptophan residues. (B) Representative solid-state 31P NMR spectra for control (DOPC:DGPP (95 mol%:5 mol%)), WALP23 and KALP23 at a 4:100 peptide–lipid ratio. (C) Quantified changes in the chemical shift of the phosphomonoester of DGPP in the presence of KALP23 and WALP23, compared to control, as a function of peptide–lipid ratio. Error bars show the standard deviation determined from at least three unique samples (* p < 0.001). (D) Representative static 31P NMR spectra for NMR samples from B that show samples formed lipid bilayers. Chemical shift values are relative to an 85% H3PO4 external reference.
Figure 1. Effect of transmembrane peptides on the chemical shift of DGPP. (A) Ribbon models for the transmembrane Alanine and Leucine (AL) peptides KALP23 and WALP23 with flanking lysine and tryptophan residues. (B) Representative solid-state 31P NMR spectra for control (DOPC:DGPP (95 mol%:5 mol%)), WALP23 and KALP23 at a 4:100 peptide–lipid ratio. (C) Quantified changes in the chemical shift of the phosphomonoester of DGPP in the presence of KALP23 and WALP23, compared to control, as a function of peptide–lipid ratio. Error bars show the standard deviation determined from at least three unique samples (* p < 0.001). (D) Representative static 31P NMR spectra for NMR samples from B that show samples formed lipid bilayers. Chemical shift values are relative to an 85% H3PO4 external reference.
Cells 11 00290 g001aCells 11 00290 g001b
Figure 2. Ionization curve for DGPP in the presence of KALP23. (A) 31P MAS NMR spectra as a function of pH for 5 mol% DGPP in 95% DOPC vesicles. (B) Static 31P NMR spectra for low, medium, and high pH values recorded for samples from A. (C) Ionization curve for the phosphomonoester of DGPP as a function of pH. The titration curve for DGPP in the presence of KALP23 (green) is compared against the control curve in the absence of KALP23 (grey), taken from [17]. The solid lines are from a Henderson–Hasselbalch fit of the ionization data (see Methods).
Figure 2. Ionization curve for DGPP in the presence of KALP23. (A) 31P MAS NMR spectra as a function of pH for 5 mol% DGPP in 95% DOPC vesicles. (B) Static 31P NMR spectra for low, medium, and high pH values recorded for samples from A. (C) Ionization curve for the phosphomonoester of DGPP as a function of pH. The titration curve for DGPP in the presence of KALP23 (green) is compared against the control curve in the absence of KALP23 (grey), taken from [17]. The solid lines are from a Henderson–Hasselbalch fit of the ionization data (see Methods).
Cells 11 00290 g002
Figure 3. Solid-state MAS 31P NMR-based titration curves for DGPP in the presence of cationic lipids. (A) DOTAP, blue diamonds and lines. Also shown are the KALP23 data from Figure 2. (B) EtPC 16 and EtPC 32 are shown in light red and dark red squares and lines, respectively. The DOTAP data in Figure 3A and Figure S4A are compared with our previous results for KALP23 and the 95:5 DOPC/DGPP control. Surprisingly, we observe no effect on the ionization behavior and charge of the phosphomonoester of DGPP due to the presence of the cationic DOTAP (the titration and especially the protonation curve for DOTAP completely overlaps with that of the PC/DGPP control). This is opposite to what we previously observed for PA. In the presence of DOTAP the ionization curve for PA moved to lower pH values, albeit less than in the presence of KALP23, indicating a lower pKa2, and hence a higher negative charge for the phosphomonoester as expected based on simple electrostatics. The effect of positive charge in the quaternary amine of DOTAP accounted for ~40% of the change in charge of PA [18] compared to the change in charge induced by KALP23 (due to both positive charge and hydrogen bonding). Since DOTAP has no effect on the ionization of the phosphomonoester of DGPP, something novel must be happening at the lipid headgroup–aqueous interface. To investigate this further, we used a cationic lipid that has a positive charge that is located further away from the hydrophobic interior of the membrane. In DOTAP, the positive charge is located in the quaternary amine, which is likely close to the pyrophosphate of DGPP. We thus chose a cationic PC where the phosphate is conjugated with an ethyl group, hence ethyl-PC (EtPC). The quaternary amine of this EtPC is located further away from the hydrophobic interior of the membrane and closer to the aqueous interface.
Figure 3. Solid-state MAS 31P NMR-based titration curves for DGPP in the presence of cationic lipids. (A) DOTAP, blue diamonds and lines. Also shown are the KALP23 data from Figure 2. (B) EtPC 16 and EtPC 32 are shown in light red and dark red squares and lines, respectively. The DOTAP data in Figure 3A and Figure S4A are compared with our previous results for KALP23 and the 95:5 DOPC/DGPP control. Surprisingly, we observe no effect on the ionization behavior and charge of the phosphomonoester of DGPP due to the presence of the cationic DOTAP (the titration and especially the protonation curve for DOTAP completely overlaps with that of the PC/DGPP control). This is opposite to what we previously observed for PA. In the presence of DOTAP the ionization curve for PA moved to lower pH values, albeit less than in the presence of KALP23, indicating a lower pKa2, and hence a higher negative charge for the phosphomonoester as expected based on simple electrostatics. The effect of positive charge in the quaternary amine of DOTAP accounted for ~40% of the change in charge of PA [18] compared to the change in charge induced by KALP23 (due to both positive charge and hydrogen bonding). Since DOTAP has no effect on the ionization of the phosphomonoester of DGPP, something novel must be happening at the lipid headgroup–aqueous interface. To investigate this further, we used a cationic lipid that has a positive charge that is located further away from the hydrophobic interior of the membrane. In DOTAP, the positive charge is located in the quaternary amine, which is likely close to the pyrophosphate of DGPP. We thus chose a cationic PC where the phosphate is conjugated with an ethyl group, hence ethyl-PC (EtPC). The quaternary amine of this EtPC is located further away from the hydrophobic interior of the membrane and closer to the aqueous interface.
Cells 11 00290 g003
Figure 4. (A) Schematic showing DOPC and DGPP in the membrane with no direct electrostatic interaction between DGPP’s inner phosphate and the phosphate of DOPC. (B) Schematic showing the electrostatic interaction between the positively charged amine of DOTAP and the inner phosphate of DGPP leading to its ionization (removal of the proton).
Figure 4. (A) Schematic showing DOPC and DGPP in the membrane with no direct electrostatic interaction between DGPP’s inner phosphate and the phosphate of DOPC. (B) Schematic showing the electrostatic interaction between the positively charged amine of DOTAP and the inner phosphate of DGPP leading to its ionization (removal of the proton).
Cells 11 00290 g004
Figure 5. DGPP Ionization in the presence of Divalent Cations. The DGPP phosphomonoester peak chemical shift value is plotted for varying cation concentrations. Chemical shift values are obtained from 31P MAS NMR spectra of DOPC:DGPP (95%/5%) MLVs in buffer at pH 7.2 ± 0.1 with 100 mM NaCl. Samples contained 0–1 mM Ca2+ or 1 mM Mg2+. Error bars show the standard deviation measured from at least seven samples. * p < 0.001.
Figure 5. DGPP Ionization in the presence of Divalent Cations. The DGPP phosphomonoester peak chemical shift value is plotted for varying cation concentrations. Chemical shift values are obtained from 31P MAS NMR spectra of DOPC:DGPP (95%/5%) MLVs in buffer at pH 7.2 ± 0.1 with 100 mM NaCl. Samples contained 0–1 mM Ca2+ or 1 mM Mg2+. Error bars show the standard deviation measured from at least seven samples. * p < 0.001.
Cells 11 00290 g005
Table 1. pKa2 values and charges for DGPP in DOPC, DOPC/KALP23, DOPC/DOTAP, DOPC/EtPC bilayers a.
Table 1. pKa2 values and charges for DGPP in DOPC, DOPC/KALP23, DOPC/DOTAP, DOPC/EtPC bilayers a.
Lipid CompositionpKa2 aPhosphomonoester Charge at pH 7.2Total Charge
at pH 7.2
(±0.1) b(±0.1) b
DOPC/DGPP 95:57.44 ± 0.02−1.38−1.38
DOPC/DGPP/KALP23 95:5:46.38 ± 0.04−1.84−2.84
DOPC/DGPP/DOTAP 95:5:167.43 ± 0.03−1.37−2.37
DOPC/DGPP/EtPC 95:5:167.25 ± 0.01−1.47−2.47
DOPC/DGPP/EtPC 95:5:327.14 ± 0.02−1.53−2.53
a The error indicated for the pKa2 values comes directly from the nonlinear curve fit (and is a lower limit). b The error in the total charge is an upper estimate which follows from the reproducibility in the calculated degree of protonation f taken from [59]. This error in the total charge is a conservative estimate, i.e., the actual error is likely to be lower.
Table 2. Calculated charge for DGPP in the presence of divalent cations. Estimated charge values for the DGPP phosphomonoester in the presence of varying amounts of divalent cations at pH 7.20 ± 0.05 and 100 mM NaCl. Charge values are calculated using Equation (2) and assuming that the chemical shift of the protonated and deprotonated states of the phosphomonoester are not influenced by the addition of divalent cations.
Table 2. Calculated charge for DGPP in the presence of divalent cations. Estimated charge values for the DGPP phosphomonoester in the presence of varying amounts of divalent cations at pH 7.20 ± 0.05 and 100 mM NaCl. Charge values are calculated using Equation (2) and assuming that the chemical shift of the protonated and deprotonated states of the phosphomonoester are not influenced by the addition of divalent cations.
Cation ContentPhosphomonoester Charge ± 0.1
at pH 7.20
0 mM−1.40
0.01 mM Ca2+−1.41
0.1 mM Ca2+−1.55
1 mM Ca2+−1.88
1 mM Mg2+−1.91
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Graber, Z.; Owusu Kwarteng, D.; Lange, S.M.; Koukounas, Y.; Khalifa, H.; Mutambuze, J.W.; Kooijman, E.E. The Electrostatic Basis of Diacylglycerol Pyrophosphate—Protein Interaction. Cells 2022, 11, 290. https://doi.org/10.3390/cells11020290

AMA Style

Graber Z, Owusu Kwarteng D, Lange SM, Koukounas Y, Khalifa H, Mutambuze JW, Kooijman EE. The Electrostatic Basis of Diacylglycerol Pyrophosphate—Protein Interaction. Cells. 2022; 11(2):290. https://doi.org/10.3390/cells11020290

Chicago/Turabian Style

Graber, Zachary, Desmond Owusu Kwarteng, Shannon M. Lange, Yannis Koukounas, Hady Khalifa, Jean W. Mutambuze, and Edgar E. Kooijman. 2022. "The Electrostatic Basis of Diacylglycerol Pyrophosphate—Protein Interaction" Cells 11, no. 2: 290. https://doi.org/10.3390/cells11020290

APA Style

Graber, Z., Owusu Kwarteng, D., Lange, S. M., Koukounas, Y., Khalifa, H., Mutambuze, J. W., & Kooijman, E. E. (2022). The Electrostatic Basis of Diacylglycerol Pyrophosphate—Protein Interaction. Cells, 11(2), 290. https://doi.org/10.3390/cells11020290

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