- Research article
- Open access
- Published:
Lipid bilayer properties potentially contributed to the evolutionary disappearance of betaine lipids in seed plants
BMC Biology volume 21, Article number: 275 (2023)
Abstract
Background
Many organisms rely on mineral nutrients taken directly from the soil or aquatic environment, and therefore, developed mechanisms to cope with the limitation of a given essential nutrient. For example, photosynthetic cells have well-defined responses to phosphate limitation, including the replacement of cellular membrane phospholipids with non-phosphorous lipids. Under phosphate starvation, phospholipids in extraplastidial membranes are replaced by betaine lipids in microalgae. In higher plants, the synthesis of betaine lipid is lost, driving plants to other strategies to cope with phosphate starvation where they replace their phospholipids by glycolipids.
Results
The aim of this work was to evaluate to what extent betaine lipids and PC lipids share physicochemical properties and could substitute for each other. By neutron diffraction experiments and dynamic molecular simulation of two synthetic lipids, the dipalmitoylphosphatidylcholine (DPPC) and the dipalmitoyl-diacylglyceryl-N,N,N-trimethylhomoserine (DP-DGTS), we found that DP-DGTS bilayers are thicker than DPPC bilayers and therefore are more rigid. Furthermore, DP-DGTS bilayers are more repulsive, especially at long range, maybe due to unexpected unscreened electrostatic contribution. Finally, DP-DGTS bilayers could coexist in the gel and fluid phases.
Conclusion
The different properties and hydration responses of PC and DGTS provide an explanation for the diversity of betaine lipids observed in marine organisms and for their disappearance in seed plants.
Background
Many organisms, such as bacteria, plants, and fungi, rely on mineral nutrients taken up directly from the environment, and therefore, specialize their metabolism to cope with their limited availability. For example, most organisms have well-defined responses to phosphate limitation, including the replacement of cellular membrane phospholipids with non-phosphorous lipids. This phenomenon has been well documented in plants [1, 2], where glycolipids replace phospholipids; whereas in bacteria, fungi or microalgae, membrane phospholipids could also be replaced by another kind of non-phosphorous lipids called betaine lipids [3,4,5,6]. Presently, three basic types of betaine lipids are known in photosynthetic organisms; DGTS (diacylglyceryl-N,N,N-trimethylhomoserine), DGTA (diacylglycerohydroxymethyl-N,N,N-trimethyl-b-alanine), and DGCC (diacylglycerylcarboxy-N-hydroxymethylcholine) [7,8,9]. They are diversely represented in organisms, DGTS being the most frequent whereas DGTA and DGCC have been solely observed in marine organisms [10]. In algae, under phosphate starvation, a situation commonly met in the environment, betaine lipids replace phospholipids in extraplastidic membranes. Because betaine lipids are localized in these membranes [11, 12] and share a common structural fragment with the main extraplastidic phospholipid phosphatidylcholine (PC) (Fig. 1A and B), it can be speculated that these two lipid classes are interchangeable, but this was never demonstrated. Betaine lipid synthesis is located in the endoplasmic reticulum (ER) [13, 14] and betaine lipids are expected to be absent in photosynthetic membranes [12]. Therefore, this PC-betaine lipid replacement is not expected to occur in photosynthetic membranes. However, it might occur at the surface of the chloroplast envelope where PC might be present [15,16,17]. Nothing is known about the composition of mitochondrial membranes in algae but because PC is a major lipid component in plant and fungal mitochondria, this replacement might also occur in mitochondria.
The identification and description of betaine lipids in cells are so far only based on biochemical analyses and therefore poorly cover a wide range of organisms [7, 10, 18,19,20]. In a simplified view, betaine lipids are present in non-vascular plants, algae, and some fungi, but even these groups show several exceptions, and they are absent from seed plants, i.e., gymnosperms and angiosperms [21]. The presence of betaine lipids is not linked to the synthesis of betaine, a soluble compound present in almost every organism including most animals, plants, and microorganisms, acting as a protectant against osmotic stress [22]. All these studies were accomplished before the amendment of the classification and the use of genetic tools to improve the phylogeny. Because the biosynthetic pathway to DGTS is now known, new phylogenetic inference on the appearance and disappearance of DGTS biosynthesis enzymes could be achieved.
Biosynthesis of DGTS was first discovered in bacteria and is realized in two steps by the two enzymes BtaA and BtaB, whereas in Chlamydomonas and in fungi, it is catalyzed by one bifunctional enzyme Bta1 containing the BtaA and BtaB domains. BtaA uses diacylglycerol and S-adenosylmethionine (SAM) as substrates to produce diacylglycerol-O-homoserine that will then be used by BtaB to produce DGTS by adding successively three methyl groups on the nitrogen from the SAM donor [5, 14]. DGTA is assumed to be synthesized by decarboxylation/recarboxylation of DGTS [23] but the enzyme and the mechanisms are still unknown and nothing is known so far about DGCC synthesis. The reconstruction of phylogenetic relationships could help to decipher the evolution of the lipid metabolism from algae to higher plants, the adaptation of plant and algae strategy to nutrient starvation, and the disappearance of betaine lipid synthesis pathways through evolution.
To unravel the functions of betaine lipids as well as their ability to act as phosphorous-free PC lipid substitutes, we investigated their physicochemical properties. The only physicochemical study was carried out by Sato and Murata in 1991 [24]. They showed that the gel-to-fluid phase transition temperature of DGTS bilayers is only slightly higher than that of PC lipids, suggesting similarities in the headgroups’ mutual interactions. To which extent betaine lipids are suited to replace PC lipids also depends on other factors such as structural aspects of lipid packing and their influence on membrane fluidity, bending rigidity, and on membrane-membrane interactions. Here, we have investigated these properties by neutron diffraction and computer simulations.
Results
Phase transition for DP-DGTS bilayers is broader than for DPPC bilayers
Lipids with the same diacylglycerol backbone were chosen to study molecules that differ only by their polar head. Because DPPC has been widely studied and DP-DGTS is commercially available, our work was achieved with these two molecules with saturated C16 tails (Fig. 1A and B). To determine the chain melting temperature for DPPC and DP-DGTS, differential scanning calorimetry (DSC) thermograms were collected from monodisperse unilamellar vesicles obtained by extrusion (Fig. 1C).
The DSC data show a sharp phase transition at 40.2 ± 0.1 °C for DPPC corresponding to the transition between the ripple phase and the fluid phase, which is consistent with earlier reports on DPPC large unilamellar vesicles [25]. For DP-DGTS, the phase transition is much broader and occurs within a range of 4 °C around 43.8 °C, which confirms the previous report [24]. With this in mind and to favor the fluid phase state, all following experiments were carried out at 50 °C, at least 5 °C above the phase transition in excess water.
We then reconstituted multilayer stacks of membranes from pure DPPC or DP-DGTS at the surface of silicon wafers and collected neutron diffraction patterns (Additional file 1) that contained three or four diffraction orders, indicating that the lattice disorder of the sample is small [26]. In the studied range of relative humidities (i.e., covering high to low osmotic pressures), the coexistence of gel and fluid phase was observed with two sets of Bragg peaks (Additional file 1E and E’, Additional file 2), each set corresponding to a lamellar phase, the gel phase to the peaks at wider angles and the fluid phase to the peaks at smaller angles.
As described in previous studies [27,28,29], we extracted from these data the pressure-distance curves, i.e., the dehydrating osmotic pressure Π as a function of the repeat distance d. As seen in Fig. 2A, by increasing the humidity, DPPC molecules transit from the gel to the fluid phase via a ripple phase through a narrow window of osmotic pressures as previously reported [30, 31]. In contrast, DP-DGTS bilayers show a phase coexistence that can be observed over a wide Π-range and without the appearance of a third phase that could be attributed to a distinct ripple phase (Fig. 2B) before forming a single fluid phase at high humidity (i.e., at low Π). Based on DSC and neutron diffraction as two independent techniques, we can safely conclude that the phase transition for DP-DGTS is broad. This observation indicates that the free energy difference between the two phases is very small over a wide osmotic pressure range and may be connected to the shapes of the pressure-distance relations in the two phases, which are discussed further below.
Hydration repulsion is stronger between DP-DGTS bilayer than between DPPC bilayers
Hydration forces are repulsive and of short range [32]; they are important for diverse processes such as preventing membrane fusion and adhesion [33] or protein adsorption [34]. To evaluate hydration forces, using membrane diffraction data at 8% D2O contrast as described in the material and method and Additional file 3, we extracted the water layer thickness dw for each osmotic pressure (Fig. 3). The exponential fit to the data, \(\Pi ={\Pi }_{0}{e}^{-{d}_{w}/{\lambda }_{w}}\) in Fig. 3A and B, yields the decay length λw, representing the short-range repulsion between the bilayers, which is dominated by hydration forces [35].
For the gel phase of DPPC, because we only have two data points, we could not evaluate the decay length. In the fluid phase, we found that DPPC λw(fluid) = 5.2 Å, in acceptable agreement with the 4 Å reported by Kowalik et al. [35]. For DP-DGTS, because the gel and fluid phases coexist in a wide range of pressures, we could calculate the decay length of both. We obtain λw(gel) = 6.9 Å and λw(fluid) = 7.7 Å. λw(gel) is about 3 times higher for DP-DGTS bilayers than what is reported in the literature for DPPC (λw(gel) = 2.1 Å [35] whereas λw(fluid) is only 1.4 times higher for DP-DGTS bilayers than for DPPC. Moreover, there is no significant difference in the decay length between gel and fluid phases in DP-DGTS bilayers, which is in contrast to DPPC. These differences suggest that the hydration repulsion between the surfaces of these two bilayers is significantly different.
Short-range interactions between adjacent bilayers in a stack could be represented by the inter-membrane compression modulus B [29]. As explained in the “Methods” section and supported by Additional file 4, we could extract B = 29 MPa for DPPC at RH = 99%. This is consistent with an earlier report of B = 46 MPa measured at 60 °C and 95% RH [36], because B decreases systematically with increasing hydration due to the associated decay of the membrane interactions [37]. For DP-DGTS, B = 1.5 MPa, a much smaller value than the one for DPPC. The difference can be attributed to the water layer, which is thicker for DP-DGTS at this humidity (Fig. 3B). Therefore, the compression modulus of DP-DGTS is lower than that of DPPC, reflecting weaker interactions between DP-DGTS bilayers than between DPPC bilayers.
To evaluate the long-range interaction between lipid bilayers, we determine the lamellar period of the bilayer at full hydration, by studying multilamellar vesicles (MLV) of DPPC and DP-DGTS in excess water by Small Angle Neutron Scattering (SANS) (Fig. 4). For DPPC, the periodic lamellar structure indicates a repeat distance of 65.8 Å in the gel phase at T = 20 °C and 67.8 Å in the fluid phase at T = 50 °C, confirming the literature values [38]. However, for DP-DGTS a Bragg peak could hardly be identified, indicating an only very weak correlation between membranes. A fit to the weak intensity modulations yields the repeat distance of 190 Å in the gel phase at 20 °C and 210 Å in the fluid phase at 50 °C.
The behavior of DP-DGTS indicates the presence of a long-range repulsive contribution to the interaction between DP-DGTS bilayers that is not present for DPPC bilayers. On top of the hydration repulsion, the repulsion may be strengthened due to the presence of charges. Namely, one cannot exclude that a small fraction of lipids deviates from the zwitterionic state, resulting in a tiny effective surface charge, which might be more pronounced for DP-DGTS because the pKa values of phosphates and carboxyl groups are not identical. Altogether, these data indicate that the interaction between DP-DGTS bilayers is more repulsive than that between DPPC bilayers, both at short and long-range.
DP-DGTS bilayers are thicker and more rigid than DPPC bilayers
Because the bilayer thickness affects membrane protein insertion as well as membrane mechanical properties, we also analyzed the bilayer thickness db as a function of Π (Fig. 3C and D). The db values obtained for DPPC are 42.5 ± 0.6 Å and 36.1 ± 0.7 Å respectively for the gel and the fluid phase, with a slight decrease in db when Π increased. These values are slightly smaller than those described in the literature (44.2 Å for the gel phase [29] and 37.2 Å for the fluid phase [39], consistent with our bilayer thickness definition, measured between the middle of the two opposite polar heads and therefore not comprising the complete polar heads as the references cited above. Our SLD analysis from Additional file 3 and Fig. 3C and D reproduces the known dehydration-thickening of fluid PC lipid multi-bilayers. Indeed, the slight increase of the bilayer thickness upon osmotic dehydration is expected because the hydration favors disordering of the lipid components of the bilayer as previously reported for dioleoylPC [40]. This phenomenon explains the transition between the gel and the fluid phase by simply decreasing Π at constant temperature. Altogether, by comparison with literature, the results obtained on DPPC bilayers validate the methodology employed here.
By following the same experimental procedure as for DPPC bilayers, we observed for DP-DGTS a decrease of bilayer thickness between the gel and fluid phase, as in DPPC bilayers and as expected upon fatty acid chain melting. However, whatever the phase, the thickness of DP-DGTS bilayers (45.2 ± 0.9 Å in the gel phase and 40.9 ± 3.8 Å in the fluid phase) is significantly higher than that of DPPC. Furthermore, the DP-DGTS bilayer thickness in the fluid phase seems to respond differently to hydration. In fact, whereas DPPC bilayers follow a dehydration-thickening trend, DP-DGTS bilayers seem to follow a dehydration-thinning trend.
The membrane bending modulus κ is a measure of the energy required to bend a membrane from planarity to a defined curvature. This value, as the compression modulus B mentioned earlier, is accessible by Bragg sheet analysis (Additional file 4 and “Methods” section). We obtained κ ≈ 14 ± 2 kBT and κ ≈ 22 ± 2 kBT for DPPC and DP-DGTS, respectively. The value found for DPPC is in agreement with the literature (18.3 ± 3.1 kBT or 15.0 ± 1.6 kBT) [41]. The higher bending modulus for DP-DGTS can be explained by the higher bilayer thickness. In fact, according to classical beam theory the bending rigidity of a homogeneous planar object scales with the third power of the thickness. Application of this relation to the bilayers while simplifying them as homogeneous objects yields the estimate \({\kappa }_{DGTS}={\kappa }_{DPPC}\bullet {\left({d}_{b}^{DGTS}/{d}_{b}^{DPPC}\right)}^{3}\approx 21 k\mathrm{B}T\), which is in good agreement with our data on bilayer thicknesses.
Molecular dynamics simulations confirm stronger hydration repulsion of DGTS and reveal differences in the water distribution
To further investigate the differences observed between DPPC and DP-DGTS bilayers, we employed atomistic molecular dynamics (MD) simulations. The thermodynamic extrapolation method [42, 43] allows us to perform simulations at prescribed water chemical potential and thus to obtain the interaction pressures between membranes with high precision (Additional file 5A and B). Snapshots of our simulations in fluid and gel states are presented in Additional file 6.
We extracted from the simulations the Π dependence of the bilayer thickness (Additional file 5C and D) and of the water layer thickness (Additional file 5E and F) and compared the experimental data (circles) to the simulations (triangles). While the simulations largely reproduce the experimental data for DPPC, indicating that the simulation model represents well the hydration forces between DPPC bilayers, the hydration dependence of the bilayer thickness of DP-DGTS is not so well reproduced. However, the simulations clearly confirm that DP-DGTS bilayers are thicker and interact more repulsively than DPPC bilayers, although the difference obtained by the simulations is not quite as pronounced as observed in the experiments. This discrepancy may be attributed to the fact that any deviation from charge-neutrality, which may occur in the experimental system due to lipid (de-) protonation effects, is not captured by the simulations, in which the molecules are represented as strictly zwitterionic.
However, to investigate the difference observed in the simulations between DPPC and DP-DGTS bilayers, in the spirit of our earlier work [28, 43], we analyzed the polar headgroup orientation with respect to the membrane normal (Fig. 5A). For DPPC, we find that the headgroup is predominantly aligned parallel to the surface, especially at low hydration as previously reported [28]. In contrast, for DP-DGTS the headgroup vector is more perpendicular to the surface, which suggests a more repulsive dipole configuration in line with the more repulsive pressure-distance curve of DP-DGTS. However, a more quantitative analysis of the headgroup dipole moment \({M}_{z}\) in the direction normal to the membranes (see inset of Fig. 5A and the “Methods” section) reveals that Mz for DGTS is only about half as strong as for PC. Indeed, the electric dipole of DGTS is smaller than that of PC (1 carbon linking the nitrogen to the carboxylate group in DGTS versus 2 carbons and 1 oxygen linking the nitrogen to the phosphate in PC) and the intramolecular charge distribution is different (Additional file 7). Thus, the more long-ranged nature of the repulsion observed between DP-DGTS bilayers is not due to the dipole strength or orientation.
We then investigated the water orientation profiles (Fig. 5B). For DPPC, water is strongly polarized and therefore oriented close to the membrane surfaces and this polarization decays smoothly towards the center of the water layer [28], which is expected from simple continuum modeling [44]. In contrast, for DP-DGTS the water polarization is strongly influenced by the negatively charged carboxylate group inside the polar head region, indicating a strong interaction with water molecules, not present with the phosphate group in DPPC. We therefore extracted the component density profiles of the simulated bilayers to analyze the interconnection between water molecules and polar heads (Fig. 5C and D). Interestingly, a water layer strongly bound to the polar head of DP-DGTS appears as a shoulder in the profile in Fig. 5C, whereas this shoulder is absent in Fig. 5D for DPPC bilayers. Therefore, the simulation indicates different interactions with water molecules for DP-DGTS bilayers that could explain partly the observed stronger repulsion between DP-DGTS bilayers. However, the long-ranged distance repulsion not entirely reproduced by molecular dynamics could also be attributed to incomplete charge-neutrality if a small fraction of DP-DGTS molecules deviates from the zwitterionic state.
Evolution of betaine lipid is influenced by environmental conditions
Our results indicate that the polar heads of DGTS and PC behave differently due to their different chemical structures, which in turn leads to different properties of the formed bilayers. This raises the question of whether the disappearance of betaine lipids through evolution is linked to these differences. To address this question, phylogenetic analyses were done on the gene involved in DGTS synthesis. DGTS is biosynthesized in a two-step reaction catalyzed either by two enzymes BtaA and BtaB in prokaryotes, or by a bifunctional enzyme Bta1 in eukaryotes enclosing both BtaA and BtaB domains [5, 14]. By domain sequence homology, we retrieved Bta1 amino acid sequences from at least 14 families of eukaryotes, bacteria, and archaea and performed separate Bayesian phylogenetic analyses on each domain (Fig. 6 and Additional file 8) because in eukaryotes BtaA and BtaB are arranged in two different configurations [45].
The phylogenetic tree shows a robust separation between marine and freshwater organisms (posterior probability 1.00) and does not corroborate taxonomy [46,47,48]. Such clustering overlaps with the domain arrangement, i.e., BtaA-BtaB in marine organisms vs BtaB-BtaA in freshwater ones. Interestingly, the planctomycete sequences, a clade of marine bacteria [49], cluster robustly (PP 0.92) at the base of the clade containing all marine protists and algae, whereas Archea and other bacteria form a clade in basal position. Also intriguing is the position of the marine Mamiellales (Chlorophyta) such as Ostreococcus taurii clustering together with the haptophyte Emiliania huwleyi (PP 0.99), and diatoms (PP 0.51) such as Phaeodactylum tricornutum, while other terrestrial chlorophytes cluster together in a separate clade.
The habitat was taken into consideration as a character in a bi-partitioned Bayesian analysis and results suggest that either two distinct lateral gene transfers occurred for marine and freshwater organisms, or the two BtaA domains split very early in the evolution (around 1.5 billion years ago) and then evolved separately in marine and freshwater organisms. In this case, the probability that the Last Universal Common Ancestor harbored a BtaA domain similar to the marine or the freshwater organisms is not decisive (0.42 ± 0.14 and 0.49 ± 0.14, respectively). Contrariwise, the Most Recent Common Ancestor of freshwater plants highly likely presented a freshwater BtaA domain (0.83 ± 0.09). Along evolution, habitat switches occur and are detectable as homoplasies (Chloropicon primus, a marine species, clusters with liverwort and non-vascularized plants) or synapomorphies (the terrestrial fungi). As a result, the phylogenetic analysis indicates that the evolution of the Bta1 gene might be driven by the physicochemical properties of the environment, such as the osmotic pressure.
Discussion
Betaine lipids have been reported in the literature for a long time to be a good substitute for phospholipids in extraplastidial membranes, especially under phosphate starvation conditions. We show that DGTS biosynthesis evolution splits between marine and freshwater organisms and does not follow the taxonomy (Fig. 6 and Additional file 8), suggesting a different selection pressure acting on marine and freshwater environments. Noticeably, marine organisms are naturally rich in very long-chain polyunsaturated fatty acids (VLC-PUFAs) — i.e., 20:5 and 22:6 — whereas freshwater organisms contain only medium-chain (16 or 18 carbon) PUFAs. A correlation between the marine habitat and the synthesis of VLC-PUFA was suggested [50, 51].
We show that DP-DGTS bilayers in the fluid phase are on average 6 Å thicker than DPPC bilayers (Fig. 3), have higher bending rigidity, and a greater tendency than DPPC bilayers to coexist in the gel and fluid phase (Figs. 1 and 2). Including PUFAs in phospholipid membranes affects the physical properties of the membranes by promoting the fluid phase, increasing disorder, and decreasing thickness [52]. In vivo, in marine organisms, DGTS and PC do not share the same fatty acid composition. For instance, in Microchloropsis gaditana, a eustigmatophyte (stramenopile) marine microalgae, DGTS is richer in 20:5 than PC that contains mainly 16 and 18 carbon fatty acids with 0, 1, or 2 unsaturations [53]. Thus, depleting the PC of VLC-PUFAs could increase the membrane thickness as well as the energy required for deformation, leading to a higher bending rigidity [54]. Therefore, the variation of fatty acid composition observed between PC and DGTS in vivo could contribute to have matching physicochemical properties of these two lipids and favor the replacement of one lipid by another depending on the availability of phosphorus in the marine environment [55, 56]. This adaptation mechanism could allow through the evolution of marine organisms to produce other forms of betaine lipids such as DGTA and DGCC.
However, freshwater organisms and plants do not produce VLC-PUFA, they mainly contain fatty acids with 16 or 18 carbons. Therefore, there is no possible adaptation of the membrane thickness by changing the fatty acid chain length [50]. Freshwater organisms with a significant fraction of betaine lipids in their composition usually have little or no PC [57]. Even some algae contain only betaine lipids and no PC such as Chlamydomonas reinhardtii whereas some others do not contain betaine lipids but only PC such as Chlorella sorokiniana [58, 59]. Nonetheless, it was found very recently that in phosphate starvation Chlorella kessleri replaced PC by DGTS completely although it was previously thought to be an alga that produces no betaine lipid [60]. Therefore, we could suppose that the inability of freshwater organisms to make PC and DGTS “compatible” through fatty acid chain adjustments forces these organisms to avoid PC/DGTS mixed bilayers. However, these organisms conserve their pathway to synthesize DGTS to maintain the possibility to use betaine lipids to face phosphate starvation by erasing PC.
In land plants, the loss of DGTS synthesis is concomitant to the loss of desiccation tolerance [61]. DGTS was detected in some mosses and ferns but is degraded during the dry season [21, 62, 63] and DGTS synthesis genes are absent in gymnosperms and angiosperms. We showed that DP-DGTS bilayers have more repulsive hydration interactions than DPPC bilayers at short and long ranges, possibly due to polar head strong interaction with water molecules and a stronger deviation from the zwitterionic state. In land plants during phosphate starvation, PC is partly replaced by DGDG which has a much less repulsive hydration interaction than PC [28]. Interaction modalities between water and DGTS are possibly deleterious in terrestrial environments. Further investigations are needed to address this question; for instance by producing Arabidopsis thaliana containing DGTS and by investigating its properties and membrane architecture during drought.
Conclusions
In marine organisms where VLC-PUFAs are present, DGTS and PC could coexist within the same bilayer because bilayer thickness could be kept homogenous by adapting fatty acid composition. The environment pressure on betaine lipid evolution was therefore different and allowed marine organisms to produce other forms of betaine lipids such as DGTA and DGCC. Betaine lipids in these conditions are a good substitute for PC during phosphate starvation. In freshwater organisms, the coexistence of PC and DGTS is scarce usually with one dominating the other in terms of abundance in the membrane. The hypothesis of the absence of VLC-PUFA that would allow DGTS and PC bilayer to have similar thickness needs to be experimentally validated. The correlation between the loss of desiccation tolerance and the loss of DGTS is intriguing. Is the difference of bilayer thickness between DGTS and PC involved in this process? This property might be an explanation for the progressive disappearance of betaine lipids in land plants. Further investigations are needed to address this question; one would be to produce an Arabidopsis thaliana plant containing DGTS and investigate its properties and architecture during drought.
Methods
Lipids
1,2-Dipalmitoyl-sn-glycero-3-phosphocholine (16:0/16:0 PC) and 1,2-dipalmitoyl-sn-glycero-3-O-4′-(N,N,N-trimethyl)-homoserine (16:0/16:0 DGTS) were purchased at Avanti Polar Lipids (Alabaster, AL, USA).
Differential scanning calorimetry
Differential scanning calorimetry (DSC) experiments were performed using a Micro DSC III differential scanning calorimeter (Setaram, France). The aluminum capsules containing the samples were filled with 500 µL of 0.1% and 0.2% (w/v) solutions of DPPC and DP-DGTS, respectively. The lipid suspensions were extruded with a thermostated extrusion kit from Avanti Polar Lipids (Alabaster, AL, USA) by passing the suspensions 11 times through 0.1 µm pore diameter polycarbonate membranes. The two samples were heated and cooled between − 4° and 55 °C, with a heating rate of 0.5 °C/min. Only the last three of four cycles were identical and used in the data analysis, the first one being always different from the following cycles due to the annealing of the vesicles after the extrusion step. The phase transition temperature was defined as the maximum of the endo- and exothermic peaks.
Neutron experiments
Neutron diffraction experiments
All neutron diffraction data were collected on the D16 diffractometer of the Institut Laue-Langevin (ILL, Grenoble, France), according to previous works [27]. The wavelength of the neutron beam was λ = 4.47 Å. The experiments were conducted using the BerILL humidity chambers developed at ILL [64] to control in situ the sample’s hydration and temperature. For both lipids, the temperature in the chamber was maintained at 50 °C during the measurements. The range of humidities investigated was 30% to 100% RH by changing the temperature of the water reservoir generating the water vapor. The Ω-scans were collected by steps of 0.05 deg in an Ω range of − 1 to 15 deg. The lamellar periodicity d was calculated using Bragg’s law as explained in the previous article.
Neutron scattering length density profiles (NSLD) were calculated from the integrated intensities of Bragg peaks corrected for the beam geometry, the neutron absorption (Cabs), and the Lorentz factor correction (CLor) according to [65], resulting in the corrected discrete structure factor of order n:
where In is the intensity of the Bragg peak at the order n.
The corrections are given by:
where µ is the absorption coefficient (5 cm−1), t is the sample thickness calculated from the deposited amount of dry lipid (0.5 mg) and the sample area (10 cm2). It equals 50 µm before hydration, for a deposited amount of 0.5 mg.
The NSLD were constructed with the Fourier transform of the structure factor as follows [66]:
where z is the direction perpendicular to the bilayer planes and vn corresponds to the phase of the structure factor of order n. According to the literature, at the 8% D2O contrast, we tested different hypotheses for the discrete structure factor signs to obtain a centrosymmetric SLD profile with a minimum at z = 0, the bilayer mid-plane, where hydrogen-rich methyl groups are located and yield the lowest SLD, and the highest SLD for the polar head regions. The assigned phases that give the best agreement with these constraints were − , − , + , and − in agreement with [67] and our molecular dynamics analysis.
We consider here as a definition of the bilayer thickness (db) the center-to-center distance between lipid polar heads as obtained from a fit to the two headgroup layer positions in the NSLD profile (as illustrated in Additional file 3) [29]. The error associated to db is given by the standard deviation between the fit and the calculated NSLD. Finally, the water layer thickness (dw) is calculated from the known d-spacing d and the bilayer thickness according to:
and the error is defined as \(\delta {d}_{w}= \sqrt{{\delta }^{2}d+{\delta }^{2}{d}_{b}}\).
Bragg sheet analysis
According to the Discrete Smectic Hamiltonian description of interacting multilayers [68] the membrane fluctuation self and cross correlations that give rise to characteristic off-specular scattering are governed by the mechanical properties of the interacting membranes in terms of the membrane bending modulus κ and the inter-membrane compression modulus B. As we have shown earlier, the experimentally obtained reciprocal space maps within this framework can be satisfactorily modeled solely based on the underlying mechanical parameters κ and B, and on an empirical cut-off parameter termed R [69]. In practice, this procedure relies on the kinematic approximation (KA) of wave scattering, because application of the more accurate distorted-wave Born approximation [70] would require detailed additional knowledge of the sample structure, which is unavailable. As a consequence, our KA-based treatment, which is only valid wherever the intensity is weak compared to the incident beam, does not correctly capture the specular maximum of the first Bragg sheet, where this condition is typically violated. In the past, we therefore ignored the first Bragg sheet and relied on the second one [36, 69]. In line with our more recent work [71], we combine here information from the first two Bragg sheets (see Additional file 4): while the Caillé parameter
is obtained from the specular/diffuse scattering intensity ratio in the second Bragg sheet, the de Gennes parameter
is obtained from the decay of the off-specular intensity in the first Bragg sheet along q||, excluding the specular intensity which violates the KA. The best-matching values of η, Λ, and R are then determined by their systematic variation in the model until the best agreement with the experimental data is achieved. Finally, the mechanical parameters are obtained by solving Eqs. 6 and 7 for κ and B.
Small angle neutron scattering (SANS)
For the SANS experiments, aqueous solutions of DPPC and DGTS were prepared in pure D2O at 10% (w/v) and 3% (w/v), respectively, and then vortexed before transfer into 1-mm path length quartz cells (Hellma, Germany) for the data acquisition. For DPPC, this is well above the swelling limit commonly observed at 40% water for phosphatidylcholines.
SANS data were collected at the D16 cold neutron diffractometer of the Institut Laue-Langevin (ILL, Grenoble, France). D16 uses a monochromator made of nine highly oriented pyrolytic graphite (HOPG) crystals whose orientation can be set to focus the beam vertically in a continuous manner from unfocussed to detector or sample focusing. For the SANS experiments, the D16 instrument was set to pinhole geometry by collimating the beam both horizontally and vertically to produce a symmetrical angular resolution and beam size in the two directions. As described above, the samples were bulk suspensions of multilamellar vesicles prepared in a large excess of pure D2O. The sample temperature was changed using a water bath connected to the sample changer. In these conditions, the d-spacing measured from the Debye-Sherrer ring observed on the SANS patterns yields the d-spacing at maximum swelling of the lamellar phase with excess water in equilibrium with the membrane stacks.
Set up of the computer model
The computer model of the hydrated bilayers employs atomistic representations of lipid and water molecules (for simulation snapshots, see Additional file 6). Following our previous work [35], we use the assisted freezing method [72] for the construction of fully hydrated membranes in the Lβ (gel) phase at an initial temperature of T = 270 K. Initially, 16 lipid molecules in each leaflet are arranged on a hexagonal lattice with random orientation but avoiding strong overlap of the atomic positions. 960 water molecules are inserted randomly in a 3-nm slab, corresponding to \({n}_{w}=30\) waters per lipid. The system is then carefully relaxed during a 1 ns molecular dynamics simulation forcing the dihedral angles of the lipid tails in gauche transformation as described elsewhere [72]. To enforce a Lβ configuration as shown in Additional file 6A, one of the hydrated leaflets including the hydrating waters is then selected and rotated by 180° around the z-axis, after which the system is replicated in the x–y plane resulting in a total of 128 lipid and 3840 water molecules. In all simulations, periodic boundary conditions are employed using triclinic unit cells with a lattice angle γ = 60° and the initial area per lipid is fixed to 0.58 nm2. This system is then relaxed for 100 ns at 270 K and finally heated up to 300 K using a heating rate of 0.1 K/ns which is sufficiently low to avoid spontaneous melting of the computational model [73]. During the relaxation semi-isotropic pressure coupling is employed using the Berendsen barostat [74] which in the heating step is replaced by anisotropic coupling to allow for full relaxation of the lattice modes.
To construct the corresponding bilayers in the Lα (fluid) phase, the system is constructed equivalently, but instead of anisotropic pressure coupling a semi-isotropic coupling scheme is employed in the last step at a temperature of 330 K. All scripts employed to set up these systems as well as equilibrated structures and simulation input files are available in https://doi.org/10.18419/darus-2360.
For DPPC, we use the well-established united-atom Berger force-field [75]—which is based on the Optimized Potentials for Liquid Simulations (OPLS) force field with refined parameters for the hydrocarbon tail interactions [76]—and the simple point charge/extended (SPC/E) water model [77]. Although this forcefield is known to be problematic when comparing structural properties from NMR experiments, it reproduces experimental pressure-distance data (Additional file 5) [35, 43] and the chain melting thermodynamics [73]. Furthermore, comparison with simulations results that use different force fields for lipids and water have shown that the hydration thermodynamics are robust with respect to force field variations [35, 78]. Contrary to the original parametrization [75] electrostatic interactions are calculated by the Particle − Mesh − Ewald (PME) method [79, 80] using a relative accuracy of 10–5 and we checked in Additional file 9 systematically the influence of the cutoff for Lennard–Jones interactions by using published experimental values [81,82,83]. We found that rc = 1.4 nm yields lipid molecular areas and water densities consistent both with the LJ-PME method [84] corresponding to an infinite cutoff and also in agreement with experimental data.
An atomistic representation of DP-DGTS is obtained by adapting the Berger DPPC atomistic topology with refined OPLS parameters for the carboxyl group [85]. In detail, the nitrogen group and the lipid tails are identical and the additional methylene groups are taken from the OPLS force field which results in total charge neutrality on the betaine headgroup. All headgroup bonded parameters stem from OPLS. An initial structure of a DP-DGTS molecule was created using the geometry optimization of the Avogadro simulation package [86]. All the simulation input files are available at https://doi.org/10.18419/darus-2360 and corresponding molecular topologies are displayed in Additional file 7.
Computer simulations
All atomistic Molecular Dynamic (MD) simulations are performed using versions 2020 and 2021 of the GROMACS simulation package [87]. All simulations are performed with an integration time step of 1 fs in the canonical constant pressure ensemble as described above. Temperature was maintained at T = 300 K for the gel phase and T = 330 K for the fluid phase using the canonical velocity rescaling thermostat [88] with a characteristic time of 0.5 ps. The time constant for the pressure coupling is set to 1 ps. Lennard–Jones interactions are truncated and shifted to zero at \({r}_{c}=1.4 {\text{nm}}\). Water molecules are kept rigid using the SETTLE algorithm [89]. Analysis of the simulations is performed using our freely available MAICoS package (https://www.maicos-analysis.org/) and chemical potentials are evaluated using the multistate Bennet acceptance ratio method (MBAR) [90] and the toolchain of the Alchemistry project (https://alchemistry.org/).
In the computer simulations, the dehydrating (osmotic) pressure is evaluated using the Gibbs–Duhem relation,
where \({\mu }_{0}\) and \({v}_{w}^{0}\) denote the chemical potential and the partial molecular volume, respectively, of pure water in bulk and \(\mu\) is the chemical potential of water between the bilayers at a given hydration level [43]. Equation (8) indicates that the chemical potential of water in atomistic MD simulations needs to be evaluated with a precision as high as \(\delta \mu \sim 0.01 {k}_{B}T\). We independently measure the excess and ideal contributions, \({\mu }^{ex}\) and \({\mu }^{id}={k}_{B}T\mathrm{ln}({\rho }_{w}{\Lambda }^{3}/{m}_{w})\), where \({\rho }_{w}\) is the water density, \(\Lambda\) its thermal wave length, \({N}_{A}\) the Avogadro number and \({m}_{w}\) the water molecular mass. \({\mu }^{ex}\) is evaluated using a Free Energy Perturbation ansatz where a test particle is brought from the vacuum to the water using 38 discrete \(\lambda\) states for the MBAR analysis, where first the Lennard Jones interactions are turned on using a soft-core potential approach and then the electrostatic interactions are added. While by definition \(\mu\) is constant over the simulation volume in thermal equilibrium, its contributions \({\mu }^{ex}\) and \({\mu }^{id}\) are not due to the inhomogeneous water distribution perpendicular and depend on the coordinate perpendicular to the membrane surface. We thus choose to measure the chemical potential at the center of the water slab between the bilayers and fix the test particle in the plane parallel to the membrane surface using a harmonic potential acting on its center of mass. To overcome potential sampling issues due to the slow dynamics of lipids in the gel phase five independent systems are constructed as described above at \({n}_{w}=30\) water per lipid and then step-wise dehydrated by a 100 ns equilibration run. The five equilibrated systems at each \({n}_{w}\) are then sampled at each \(\lambda\) state for 100 ns, corresponding to \(5\times 38\times 100 \mathrm{ns}=19 \mathrm{\mu s}\) sampling time per data point shown in Fig. 5 and Additional file 5. For simulations in the fluid phase, one system is melted at \(T=330 {\text{K}}\) and dehydrated equivalently, where for consistency each \(\lambda\) state is sampled for 500 ns and analyzed in blocks of 100 ns.
NSLD profiles can be obtained readily from the simulation data by multiplying the united atom density probability with the corresponding scattering length densities, which we took from NIST (https://www.ncnr.nist.gov/resources/activation/). However, this approach corresponds to taking infinitely many terms into account in Eq. (4), which in the experimental analysis is impossible since only a few Bragg peaks are present. We thus first calculate the corresponding structure factors \({\left|{F}_{2}\right|}^{2}\) via discrete cosine transform and limit to three terms in the reciprocal space, which is the typical number of Bragg peaks present in the experimental analysis. After back transform, this yields NSLD profiles from which the bilayer thickness \({d}_{b}\) is evaluated. Note that the phases \({\nu }_{n}\) obtained from the simulations perfectly agree with the experimental analysis. We observe a constant shift in the \(\Pi -d\) curves in Additional file 5 when comparing our simulation results to the experimental data of Δdfluid=3Å for the fluid phase and Δdgel=8Å for the gel phase, which we attribute to the united-atom representation of the lipid tails and consequent underestimation of the chain-chain repulsion between lipid monolayers. This explanation agrees well with the observed larger shift for the gel phase where the tails are more ordered and stretched. Consequently, simulation data for \(d\) and \({d}_{b}\) have been shifted by \(\Delta d\) in Additional file 5.
Density profiles for Fig. 5 as well as for the NSLD profiles have been extracted from the simulation trajectory using our open-source analysis framework MAICoS (https://www.maicos-analysis.org/), where in Fig. 5C and D, we additionally split into the zwitterionic headgroup parts, the water and the lipid densities. The water polarization profiles in Fig. 5B are obtained from the water charge density profiles \({\rho }_{\mathrm{w}}^{(c)}\) as \({m}_{\mathrm{w}}\left(z\right)= {\int }_{-d/2 }^{z}{\rho }_{\mathrm{w}}^{(c)}\left(z\right)dz{\prime}\), where the lamellar repeat distance \(d\) equals the size of the simulation box in \(z\)-direction. Similarly, the dipole moment of the polar heads in the direction perpendicular to the surface is obtained from the polarization densities \({m}_{\mathrm{l}}\) of a single membrane leaflet according to \({M}_{z}= {A}_{\mathrm{l}}{\int }_{-d/2}^{d/2}{m}_{\mathrm{l}}\left(z\right)dz\), where the subscript \(\mathrm{l}\) denotes that only the lipid charges are considered in the double integration and \({A}_{\mathrm{l}}\) denotes the lateral area per lipid. In detail, for DP-DGTS the tetramethylammonium-carboxylate vector is used, whereas for DPPC the phosphate-tetramethylammonium vector is used. Since by construction of the models, the lipids carry no net charge, this is exactly the dipole moment per unit volume of the membrane.
Phylogeny
The Microchloropsis gaditana locus Naga_100016g36 (annotated as Protein of unknown function DUF3419) was identified as putative betaine lipid synthase by homology with Nannochloropsis oceanica BTA1 [45]. The Naga_100016g36 CDS was translated and the amino acid sequence used as a query for a series of BLAST searches in all the publicly available databases. Both tBLASTn and PSI-BLAST [91, 92] were performed together with a human-curated process in order to obtain the widest dataset possible and the most robust one. The same was done with all the sequences used by [45]. The final dataset contained amino acid sequences from at least 14 families of eukaryotes, bacteria, and archaea. Because in most eukaryotes the two domains of the BTA protein (the 3-amino-3-carboxypropyl-transferase and the S-adenosyl-L-methionine-dependent methyltransferase (SAM) domains) are arranged in two different configurations, two separate alignments were produced. The amino acid sequences of the two domains were compiled in two separate fasta files and aligned using MUSCLE [93] then curated using BMGE (Block Mapping and Gathering with Entropy) software [94] to select the phylogenetic informative regions. For each alignment, the substitution model that best fits the data was selected by running the “Find Best DNA/Protein Model” utility implemented in MEGA X (Molecular Evolutionary Genetics Analyses) software [95]. The Le-Gascuel substitution model for Maximum Likelihood (ML) phylogenetic inference method was chosen as the best fitting for both datasets [96]. Non-uniformity of evolutionary rates among sites was modeled using a discrete Gamma distribution (+ G). For the SAM domain alignment, a certain fraction of sites were considered to be evolutionarily invariable (+ I). The ML phylogenetic analyses were supported by 5000 (bootstrap) pseudoreplicates. Initial trees for the heuristic search were obtained automatically by applying Neighbor Joining (NJ) and BioNJ algorithms to a matrix of pairwise distances estimated using a JTT model, and then selecting the topology with superior log likelihood value. All positions with less than 95% site coverage were eliminated. Trees were drawn to scale with the branch length measured in the number of substitutions per site.
In order to infer both the fusion history and the habitat of the organisms possessing a BTA gene, a Bayesian inference was carried out using MrBayes v3.2.7 [97, 98] using a partitioned model [99]. After some trial runs, the conditions for the Bayesian analysis were set up to ensure that the Average standard deviation of split frequencies reached stationarity over the course of the sampling. For each analysis, a total of 1,600,000 generations was implemented, with successive samples separated by 100 generations after an initial “burn in” period of 25% of the number of samples. The Bayesian posterior probabilities (BPP) were estimated by two independent runs of four Metropolis Coupled chains (MCMCMC). For each dataset, the model selection was done during the analysis by estimating the posterior probabilities of the different models together with their parameters.
Availability of data and materials
All data generated or analyzed during this study are included in this published article, its supplementary information files and publicly available repositories. The neutron data that support the findings of this study are available upon request at the ILL with this link https://doi.ill.fr/10.5291/ILL-DATA.8-02-820, https://doi.ill.fr/10.5291/ILL-DATA.TEST-3119 and the data that support the MD simulation are publicly available at https://doi.org/10.18419/darus-2360.
Abbreviations
- BPP:
-
Bayesian Posterior Probability
- DP:
-
Dipalmitoyl
- DGCC:
-
Diacylglycerylcarboxy-N-hydroxymethylcholine
- DGTA:
-
Diacylglycerohydroxymethyl-N,N,N-trimethyl-b-alanine
- DGTS:
-
Diacylglyceryl-N,N,N-trimethylhomoserine
- DSC:
-
Differential scanning calorimetry
- ER:
-
Endoplasmic reticulum
- NSLD:
-
Neutron Scattering Length Density
- MD:
-
Molecular dynamics
- PC:
-
Phosphatidylcholine
- SAM:
-
S-adenosylmethionine
- SANS:
-
Small Angle Neutron Scattering
References
Andersson MX, Larsson KE, Tjellström H, Liljenberg C, Sandelius AS. Phosphate-limited oat. The plasma membrane and the tonoplast as major targets for phospholipid-to-glycolipid replacement and stimulation of phospholipases in the plasma membrane. J Biol Chem. 2005;280:27578–86.
Jouhet J, Maréchal E, Baldan B, Bligny R, Joyard J, Block MA. Phosphate deprivation induces transfer of DGDG galactolipid from chloroplast to mitochondria. J Cell Biol. 2004;167:863–74.
Abida H, Dolch L-J, Meï C, Villanova V, Conte M, Block MA, et al. Membrane glycerolipid remodeling triggered by nitrogen and phosphorus starvation in Phaeodactylum tricornutum. Plant Physiol. 2015;167:118–36.
Mühlroth A, Winge P, El Assimi A, Jouhet J, Maréchal E, Hohmann-Marriott MF, et al. Mechanisms of Phosphorus Acquisition and Lipid Class Remodeling under P Limitation in a Marine Microalga. Plant Physiol. 2017;175:1543–59.
Riekhof WR, Andre C, Benning C. Two enzymes, BtaA and BtaB, are sufficient for betaine lipid biosynthesis in bacteria. Arch Biochem Biophys. 2005;441:96–105.
Senik SV, Maloshenok LG, Kotlova ER, Shavarda AL, Moiseenko KV, Bruskin SA, et al. Diacylglyceryltrimethylhomoserine content and gene expression changes triggered by phosphate deprivation in the mycelium of the basidiomycete Flammulina velutipes. Phytochemistry. 2015;117:34–42.
Dembitsky VM. Betaine ether-linked glycerolipids: chemistry and biology. Prog Lipid Res. 1996;35:1–51.
Kato M, Sakai M, Adachi K, Ikemoto H, Sano H. Distribution of betaine lipids in marine algae. Phytochemistry. 1996;42:1341–5.
Sato N. Betaine lipids Bot Mag Tokyo. 1992;105:185–97.
Cañavate JP, Armada I, Ríos JL, Hachero-Cruzado I. Exploring occurrence and molecular diversity of betaine lipids across taxonomy of marine microalgae. Phytochemistry. 2016;124:68–78.
Eichenberger W, Araki S, Müller DG. Betaine lipids and phospholipids in brown algae. Phytochemistry. 1993;34:1323–33.
Künzler K, Eichenberger W, Radunz A. Intracellular localization of two betaine lipids by cell fractionation and immunomicroscopy. Z Naturforsch C J Biosci. 1997;52:487–95.
Moore TS, Du Z, Chen Z. Membrane lipid biosynthesis in Chlamydomonas reinhardtii. In vitro biosynthesis of diacylglyceryltrimethylhomoserine. Plant Physiol. 2001;125:423–9.
Riekhof WR, Sears BB, Benning C. Annotation of genes involved in glycerolipid biosynthesis in Chlamydomonas reinhardtii: discovery of the betaine lipid synthase BTA1Cr. Eukaryot Cell. 2005;4:242–52.
Botella C, Jouhet J, Block MA. Importance of phosphatidylcholine on the chloroplast surface. Prog Lipid Res. 2017;65:12–23.
Karki N, Johnson BS, Bates PD. Metabolically distinct pools of phosphatidylcholine are involved in trafficking of fatty acids out of and into the chloroplast for membrane production. Plant Cell. 2019;31:2768–88.
Flori S, Jouneau P-H, Finazzi G, Maréchal E, Falconet D. Ultrastructure of the periplastidial compartment of the diatom phaeodactylum tricornutum. Protist. 2016;167:254–67.
Künzler K, Eichenberger W. Betaine lipids and zwitterionic phospholipids in plants and fungi. Phytochemistry. 1997;46:883–92.
Makewicz A, Gribi C, Eichenberger W. Lipids of Ectocarpus fasciculatus (Phaeophyceae). Incorporation of [l-14C]Oleate and the role of TAG and MGDG in lipid metabolism. Plant Cell Physiol. 1997;38:952–62.
Vaskovsky VE, Khotimchenko SV, Boolugh EM. Distribution of diacylglycerotrimethylhomoserine and phosphatidylcholine in mushrooms. Phytochemistry. 1998;47:755–60.
Rozentsvet OA. Comparative examination of distribution of phospholipids and a betaine lipid DGTS in tropical fern species. Biochem Syst Ecol. 2004;32:303–11.
Chen THH, Murata N. Enhancement of tolerance of abiotic stress by metabolic engineering of betaines and other compatible solutes. Curr Opin Plant Biol. 2002;5:250–7.
Sato N. Lipids in Cryptomonas CR-1. II. Biosynthesis of betaine lipids and galactolipids. Plant Cell Physiol. 1991;32:845–51.
Naoki Sato, Norio Murata. Transition of lipid phase in aqueous dispersions of diacylglyceryltrimethylhomoserine. Biochim Biophys Acta. 1991;1082:108–11.
Drazenovic J, Wang H, Roth K, Zhang J, Ahmed S, Chen Y, et al. Effect of lamellarity and size on calorimetric phase transitions in single component phosphatidylcholine vesicles. Biochim Biophys Acta. 2015;1848:532–43.
Wiener MC, White SH. Fluid bilayer structure determination by the combined use of x-ray and neutron diffraction. I. Fluid bilayer models and the limits of resolution. Biophys J. 1991;59:162–73.
Bolik S, Albrieux C, Schneck E, Demé B, Jouhet J. Sulfoquinovosyldiacylglycerol and phosphatidylglycerol bilayers share biophysical properties and are good mutual substitutes in photosynthetic membranes. Biochim Biophys Acta. 2022;1864:184037.
Kanduč M, Schlaich A, de Vries AH, Jouhet J, Maréchal E, Demé B, et al. Tight cohesion between glycolipid membranes results from balanced water-headgroup interactions. Nat Commun. 2017;8:14899.
Nagle JF, Tristram-Nagle S. Structure of lipid bilayers. Biochim Biophys Acta. 2000;1469:159–95.
Braganza LF, Worcester DL. Hydrostatic pressure induces hydrocarbon chain interdigitation in single-component phospholipid bilayers. Biochemistry. 1986;25:2591–6.
Matsuki H, Goto M, Tada K, Tamai N. Thermotropic and barotropic phase behavior of phosphatidylcholine bilayers. Int J Mol Sci. 2013;14:2282–302.
Lis LJ, Mcalister M, Fuller N, Rand RP, Parsegian VA. Interactions between neutral phospholipid bilayer membranes. Biophys J. 1982;37:657–65.
Leikin SL, Kozlov MM, Chernomordik LV, Markin VS, Chizmadzhev YA. Membrane fusion: overcoming of the hydration barrier and local restructuring. J Theor Biol. 1987;129:411–25.
Lu DR, Lee SJ, Park K. Calculation of solvation interaction energies for protein adsorption on polymer surfaces. J Biomater Sci Polym Ed. 1992;3:127–47.
Kowalik B, Schlaich A, Kanduč M, Schneck E, Netz RR. Repulsion difference between ordered and disordered membranes due to cancellation of membrane-membrane and water-mediated interactions. J Phys Chem Lett. 2017;8:2869–74.
Schneck E, Rehfeldt F, Oliveira RG, Gege C, Demé B, Tanaka M. Modulation of intermembrane interaction and bending rigidity of biomembrane models via carbohydrates investigated by specular and off-specular neutron scattering. Phys Rev E. 2008;78: 061924.
Mennicke U, Constantin D, Salditt T. Structure and interaction potentials in solid-supported lipid membranes studied by X-ray reflectivity at varied osmotic pressure. Eur Phys J E. 2006;20:221–30.
Zhang R, Tristram-Nagle S, Sun W, Headrick RL, Irving TC, Suter RM, et al. Small-angle x-ray scattering from lipid bilayers is well described by modified Caillé theory but not by paracrystalline theory. Biophys J. 1996;70:349–57.
Kučerka N, Tristram-Nagle S, Nagle JF. Closer look at structure of fully hydrated fluid phase DPPC bilayers. Biophys J. 2006;90:L83–5.
Stachura SS, Malajczuk CJ, Kuprusevicius E, Mancera RL. Influence of bilayer size and number in multi-bilayer DOPC simulations at full and low hydration. Langmuir. 2019;35:2399–411.
Nagle JF. Experimentally determined tilt and bending moduli of single-component lipid bilayers. Chem Phys Lipid. 2017;205:18–24.
Schlaich A, Kowalik B, Kanduč M, Schneck E, Netz RR. Physical mechanisms of the interaction between lipid membranes in the aqueous environment. Physica A. 2015;418:105–25.
Schneck E, Sedlmeier F, Netz RR. Hydration repulsion between biomembranes results from an interplay of dehydration and depolarization. Proc Natl Acad Sci. 2012;109:14405–9.
Marčelja S, Radić N. Repulsion of interfaces due to boundary water. Chem Phys Lett. 1976;42:129–30.
Murakami H, Nobusawa T, Hori K, Shimojima M, Ohta H. Betaine lipid is crucial for adapting to low temperature and phosphate deficiency in nannochloropsis. Plant Physiol. 2018;177:181–93.
Burki F, Kudryavtsev A, Matz MV, Aglyamova GV, Bulman S, Fiers M, et al. Evolution of Rhizaria: new insights from phylogenomic analysis of uncultivated protists. BMC Evol Biol. 2010;10:377.
Ochoa de Alda JAG, Esteban R, Diago ML, Houmard J. The plastid ancestor originated among one of the major cyanobacterial lineages. Nat Commun. 2014;5:4937.
Stiller JW, Schreiber J, Yue J, Guo H, Ding Q, Huang J. The evolution of photosynthesis in chromist algae through serial endosymbioses. Nat Commun. 2014;5:5764.
Fuerst JA, Sagulenko E. Beyond the bacterium: planctomycetes challenge our concepts of microbial structure and function. Nat Rev Microbiol. 2011;9:403–13.
Li-Beisson Y, Thelen JJ, Fedosejevs E, Harwood JL. The lipid biochemistry of eukaryotic algae. Prog Lipid Res. 2019;74:31–68.
Valentine RC, Valentine DL. Omega-3 fatty acids in cellular membranes: a unified concept. Prog Lipid Res. 2004;43:383–402.
Fernandes MX, Castanho MARB, Garcı́a de la Torre J. Brownian dynamics simulation of the unsaturated lipidic molecules oleic and docosahexaenoic acid confined in a cellular membrane. Biochim Biophys Acta. 2002;1565:29–35.
Billey E, Magneschi L, Leterme S, Bedhomme M, Andres-Robin A, Poulet L, et al. Characterization of the bubblegum acyl-CoA synthetase of microchloropsis gaditana. Plant Physiol. 2021;185:815–35.
Vásquez V, Krieg M, Lockhead D, Goodman MB. Phospholipids that contain polyunsaturated fatty acids enhance neuronal cell mechanics and touch sensation. Cell Rep. 2014;6:70–80.
Van Mooy BAS, Fredricks HF, Pedler BE, Dyhrman ST, Karl DM, Koblízek M, et al. Phytoplankton in the ocean use non-phosphorus lipids in response to phosphorus scarcity. Nature. 2009;458:69–72.
Cañavate JP, Armada I, Hachero-Cruzado I. Interspecific variability in phosphorus-induced lipid remodelling among marine eukaryotic phytoplankton. New Phytol. 2017;213:700–13.
Sato N, Furuya M. Distribution of diacylglyceryltrimethylhomoserine and phosphatidylcholine in non-vascular green plants. Plant Sci. 1985;38:81–5.
Giroud C, Gerber A, Eichenberger W. Lipids of chlamydomonas reinhardtii. Analysis of molecular species and intracellular site(s) of biosynthesis. Plant Cell Physiol. 1988;29:587–95.
Widzgowski J, Vogel A, Altrogge L, Pfaff J, Schoof H, Usadel B, et al. High light induces species specific changes in the membrane lipid composition of Chlorella. Biochem J. 2020;477:2543–59.
Oishi Y, Otaki R, Iijima Y, Kumagai E, Aoki M, Tsuzuki M, et al. Diacylglyceryl-N, N, N-trimethylhomoserine-dependent lipid remodeling in a green alga. Chlorella kessleri Commun Biol. 2022;5:19.
Oliver MJ, Tuba Z, Mishler BD. The evolution of vegetative desiccation tolerance in land plants. Plant Ecol. 2000;151:85–100.
Guschina IA, Harwood JL, Smith M, Beckett RP. Abscisic acid modifies the changes in lipids brought about by water stress in the moss Atrichum androgynum. New Phytol. 2002;156:255–64.
Rozenstvet OA, Saksonov SV, Filin VR, Dembitsky VM. Seasonal changes of lipid content in the leaves of some ferns. Physiol Plant. 2001;113:59–63.
Gonthier J, Barrett MA, Aguettaz O, Baudoin S, Bourgeat-Lami E, Demé B, et al. BerILL: The ultimate humidity chamber for neutron scattering. JNR. 2019;21:65–76.
Harroun TA, Katsaras J, Wassall SR. Cholesterol is found to reside in the center of a polyunsaturated lipid membrane. Biochemistry. 2008;47:7090–6.
Worcester DL, Franks NP. Structural analysis of hydrated egg lecithin and cholesterol bilayers II. Neutron diffraction. J Mol Biol. 1976;100:359–78.
Pan J, Heberle FA, Carmichael JR, Ankner JF, Katsaras J. Time-of-flight Bragg scattering from aligned stacks of lipid bilayers using the Liquids Reflectometer at the Spallation Neutron Source. J Appl Cryst. 2012;45:1219–27.
Lei N, Safinya CR, Bruinsma RF. Discrete harmonic model for stacked membranes: theory and experiment. J Phys II France. 1995;5:1155–63.
Schneck E, Demé B, Gege C, Tanaka M. Membrane adhesion via homophilic saccharide-saccharide interactions investigated by neutron scattering. Biophys J. 2011;100:2151–9.
Sinha SK, Sirota EB, Garoff S, Stanley HB. X-ray and neutron scattering from rough surfaces. Phys Rev B. 1988;38:2297–311.
Latza VM, Demé B, Schneck E. Membrane adhesion via glycolipids occurs for abundant saccharide chemistries. Biophys J. 2020;118:1602–11.
Schubert T, Schneck E, Tanaka M. First order melting transitions of highly ordered dipalmitoyl phosphatidylcholine gel phase membranes in molecular dynamics simulations with atomistic detail. J Chem Phys. 2011;135: 055105.
Kowalik B, Schubert T, Wada H, Tanaka M, Netz RR, Schneck E. Combination of MD simulations with two-state kinetic rate modeling elucidates the chain melting transition of phospholipid bilayers for different hydration levels. J Phys Chem B. 2015;119:14157–67.
Berendsen HJC, Postma JPM, van Gunsteren WF, DiNola A, Haak JR. Molecular dynamics with coupling to an external bath. J Chem Phys. 1984;81:3684–90.
Berger O, Edholm O, Jähnig F. Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature. Biophys J. 1997;72:2002–13.
Jorgensen WL, Tirado-Rives J. The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. J Am Chem Soc. 1988;110:1657–66.
Berendsen HJC, Grigera JR, Straatsma TP. The missing term in effective pair potentials. J Phys Chem. 1987;91:6269–71.
Wolde-Kidan A, Dat Pham Q, Schlaich A, Loche P, Sparr E, R. Netz R, et al. Influence of polar co-solutes and salt on the hydration of lipid membranes. Phys Chem Chem Phys. 2019;21:16989–7000.
Darden T, York D, Pedersen L. Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems. J Chem Phys. 1993;98:10089–92.
Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG. A smooth particle mesh Ewald method. J Chem Phys. 1995;103:8577–93.
Lide DR. Abundance of elements in the Earth’s crust and in the sea. CRC handbook of chemistry and physics, Internet Version. 2005;:14–7.
Lis LJ, Mcalister M, Fuller N, Rand RP, Parsegian VA. Measurement of the lateral compressibility of several phospholipid bilayers. Biophys J. 1982;37:667–72.
Petrache HI, Dodd SW, Brown MF. Area per lipid and acyl length distributions in fluid phosphatidylcholines determined by 2H NMR spectroscopy. Biophys J. 2000;79:3172–92.
Wennberg CL, Murtola T, Hess B, Lindahl E. Lennard-jones lattice summation in bilayer simulations has critical effects on surface tension and lipid properties. J Chem Theory Comput. 2013;9:3527–37.
Maciejewski A, Pasenkiewicz-Gierula M, Cramariuc O, Vattulainen I, Rog T. Refined OPLS all-atom force field for saturated phosphatidylcholine bilayers at full hydration. J Phys Chem B. 2014;118:4571–81.
Hanwell MD, Curtis DE, Lonie DC, Vandermeersch T, Zurek E, Hutchison GR. Avogadro: an advanced semantic chemical editor, visualization, and analysis platform. Journal of cheminformatics. 2012;4:1–17.
Abraham MJ, Murtola T, Schulz R, Páll S, Smith JC, Hess B, et al. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX. 2015;1–2:19–25.
Bussi G, Donadio D, Parrinello M. Canonical sampling through velocity rescaling. J Chem Phys. 2007;126: 014101.
Miyamoto S, Kollman PA. Settle: an analytical version of the SHAKE and RATTLE algorithm for rigid water models. J Comput Chem. 1992;13:952–62.
Shirts MR, Chodera JD. Statistically optimal analysis of samples from multiple equilibrium states. J Chem Phys. 2008;129: 124105.
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.
Altschul SF, Wootton JC, Gertz EM, Agarwala R, Morgulis A, Schäffer AA, et al. Protein database searches using compositionally adjusted substitution matrices. Febs J. 2005;272:5101–9.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.
Criscuolo A, Gribaldo S. BMGE (Block Mapping and Gathering with Entropy): a new software for selection of phylogenetic informative regions from multiple sequence alignments. BMC Evol Biol. 2010;10:210.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35:1547–9.
Le SQ, Gascuel O. An improved general amino acid replacement matrix. Mol Biol Evol. 2008;25:1307–20.
Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient bayesian phylogenetic inference and model choice across a large model space. System Biol. 2012;61:539–42.
Ronquist F, Klopfstein S, Vilhelmsen L, Schulmeister S, Murray DL, Rasnitsyn AP. A total-evidence approach to dating with fossils, applied to the early radiation of the hymenoptera. Syst Biol. 2012;61:973–99.
Acknowledgements
Not applicable.
Funding
This work was supported by ANR-DFG ANR-18-CE92-0015 and SCHN 1396/2. SB was supported by a joint funding Glyco@Alps / ILL PhD Program, by the French National Research Agency in the framework of the "Investissements d’avenir” program Glyco@Alps (ANR-15-IDEX-02), “Origin Of Life” (ANR-17-EURE-0003), “Oceanomics” (ANR-11-BTBR-0008) and the Labex GRAL, financed within the University Grenoble Alpes graduate school (Ecoles Universitaires de Recherche) CBH-EUR-GS (ANR-17-EURE-0003). AS acknowledges funding from the DFG under Germany's Excellence Strategy – EXC 2075 – 390740016. Computing time was provided by the state of Baden-Württemberg through bwHPC.
Author information
Authors and Affiliations
Contributions
JJ, BD, and ES conceived and designed the research; SB and BD performed the DSC and neutron diffraction experiments, TM performed the SAXS experiment, AS performed the molecular dynamic simulations, OB and AA did the phylogeny analysis. SB, BD, JJ, AS, and ES analyzed and interpreted the data. All authors were involved in drafting and revising the manuscript. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:
Figure S1. Diffraction scan I(2q, Ω) showing the Bragg peak positions of DPPC (a-c and’ a’-c’) and DP-DGTS (d-f and d’-f’) at three different humidities: low (a, a’ and d, d’) intermediate (b, b’ and e, e’) and high (c, c’ and f, f’). The measurements were done at 50°C. a to f: Solvent composition 8% D2O. Relative humidity values: DPPC 30.3, 85.0, 95.8%. DP-DGTS 32.2, 87.2, 97.9%. a’ to f’: Solvent composition 100 % D2O. Relative humidity values: DPPC 30.3, 85.0, 95.8%. DP-DGTS 32.2, 87.2, 97.9%.
Additional file 2: Figure S2.
Intensity of the Bragg peaks at 3 humidities. The gel and fluid phases in the DP-DGTS sample at RH=87.1% are noted 1 and 2. Phase 1 is probably in the gel state and phase 2 the fluid one. The temperature of the sample was 50°C.
Additional file 3:
Figure S3. Neutron Scattering Length Density (NSLD) of DPPC (a and c) and DP-DGTS (b and d) measured at 50°C as described in [27]. The gel and fluid phase are separated in two graphs for the two lipids, a and b for the gel phase, c and d for the fluid phase. The bilayer thickness db corresponds to the center to center distance between headgroups as described in [29] and the water layer thickness dw corresponds to d-db with d being the lamellar repeat distance or d-spacing and db the bilayer thickness. These are represented in panel a. The solid lines are the NSLD profile calculated using the minimum of three Bragg peaks and the dash lines are calculated using 4 or 5 peaks when visible at the highest humidities for DP-DGTS. The dotted black lines show the shift upon hydration of the maxima in the SLD profile, showing the increase of the water layer thickness.
Additional file 4:
Figure S4. 2q-integrated Bragg sheet intensities of DPPC membranes (top) and DP-DGTS membranes (bottom) at RH=99 % of the first (left) and second (right) Bragg sheets as a function of Ω, featuring the respective central specular maxima symmetrically flanked by the slowly decaying diffuse scattering intensity. The latter is locally decorated with minima at conditions of high absorption (Ω ≈ 0 and Ω ≈ 2q) and peaks arising from multiple scattering effects [69, 71]. The solid lines (in red) superimposed to the experimental data points represent simulated Bragg sheet intensities corresponding to the best-matching parameters in the continuum-mechanical model simultaneously describing the first and second Bragg sheets as explained in the material and method section. Absorption close to Ω ≈ 0 and Ω ≈ 2q was modeled as described in [71]. The best-matching model parameters for both systems are summarized with the empirical cut-off parameter termed R and the mechanical parameters k and B derived thereof according to Eqs. 6 and 7 as described in the material and method section [69].
Additional file 5:
Figure S5. a and b. Simulated lamellar period of DPPC (a) and DP-DGTS (b) bilayer stacks as a function of the humidity, measured at 50°C (circle) or computed (triangle). c to f. Comparison of the bilayer thickness (c and d) and of the water layer thickness (e and f) as a function of the hydration pressure obtained by simulation (triangles) or by experiment (circles) in DPPC bilayers (c and e) and in DP-DGTS bilayers (d and f).
Additional file 6:
Figure S6. a,b. Snapshots of the simulations in the Lβ(gel) and Lα (fluid) phases, respectively. Blue lines show the periodic simulation box. c. Top-view of a simulation in the gel phase where the head-group atoms are highlighted using spheres. The tilt angle is free to adjust according to the lipid lattice structure with a typical value of . d. For the simulations in the fluid phase the tilt angle is fixed at . Note that due to the minimum image convention in the periodic boundary conditions employed, i.e., the box angle in C is allowed to flip its direction.
Additional file 7:
Figure S7. Molecular representation of DP-DGTS (a) and DPPC (b). (i) Shows the chemical structure for both lipids and (ii) the united atom representation, where hydrogens are implicitly included on the heavy carbons (green spheres). Red spheres denote oxygen atoms, blue nitrogen atoms and gold phosphor atoms, respectively. (iii) Partial charge distribution with color coding on the right. The expected molecular dipole of the zwitterionic headgroups gets largely compensated due to the polarity of the covalent bonds.
Additional file 8:
Figure S8. Bayesian phylogenetic tree of the S-adenosyl-L-methionine-dependent methyltransferase domain of the BTA1 protein. Bayesian Posterior Probability values are reported next to the branches. The color of the node circles represents BPP. Color code is reported on the figure. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. On the right-hand side the configuration of the protein is reported; red: BtaA domain: 3-amino-3-carboxypropyl-transferase domain (Domain of Unknown Function DUF3419); blue: BtaB domain: S-adenosyl-L-methionine-dependent methyltransferase (SAM) domain.
Additional file 9:
Figure S9. a. Influence of the Lennard-Jones cutoff on the bulk water density at 300 K and at 330 K. Vertical dashed lines denote the experimental values [81] b. Area per lipid in a DPPC membrane for the different treatment of the Lennard Jones interactions. The gray are denotes the typical values obtained experimentally; the dashed line is obtained by Lis et al. [82], solid line by Petrache et al. [83].
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Bolik, S., Schlaich, A., Mukhina, T. et al. Lipid bilayer properties potentially contributed to the evolutionary disappearance of betaine lipids in seed plants. BMC Biol 21, 275 (2023). https://doi.org/10.1186/s12915-023-01775-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12915-023-01775-z