Elastic deformations of bolalipid membranes

Archaeal membranes have unique mechanical properties that enable these organisms to survive under extremely aggressive environmental conditions. The so-called bolalipids contribute to this exceptional stability. They have two polar heads joined by two hydrocarbon chains. The two headgroups can face different sides of the membrane (O-shape conformation) or the same side (U-shape conformation). We have developed an elasticity theory for bolalipid membranes and show that the energetic contributions of (i) tilt deformations, (ii) area compression/stretching deformations, (iii) as well as those of Gaussian splay from the two membrane surfaces are additive, while splay deformations yield a cross-term. The presence of a small fraction of U-shaped molecules resulted in spontaneous membrane curvature. We estimated the tilt modulus to be approximately equal to that of membranes in eukaryotic cells. In contrast to conventional lipids, the bolalipid membrane possesses two splay moduli, one of which is estimated to be an order of magnitude larger than that of conventional lipids. The projected values of elastic moduli act to hamper pore formation and to decelerate membrane fusion and fission.


Introduction
Much of the membrane's mechanical behavior has been attributed to the protein scaffold that is anchored to the bilayer core. 1 However, the lipid bilayer dominates the total cellular area compliance and bending stiffness. 2 Moreover, the lipid bilayer's mechanical properties are determinants of fusion and fission, 3 and the lipid core's resistance to tilt and bending makes significant contributions to the energetics of membrane protein conformational transitions. 4 Transmembrane signal transduction is thought to invoke the assembly of liquid-ordered domains. 5 In turn, registration of ordered domains from the two monolayers is governed by line tension 6 and, probably to a minor degree, by the mutual attraction of stiff regions in both monolayers, since their registration minimizes spatial restraints on membrane undulations, i.e. maximizes entropy. 7 The process of phase separation into domains, as well as tether formation, cell shape changes, and budding all require that the lipids in the two monolayers slide against each other. This sliding is dominated by the interleaflet drag that appears to be much larger than the in-layer viscosity. 8 However, it is totally unclear what governs lipid membrane reshaping in the absence of interleaflet drag. This is the case in archaea's cell membrane, which is formed by unique components, so-called bolalipids (bipolar lipids), as opposed to membranes of bacteria and eukaryotes. Bolalipids are believed to be responsible for the phenomenal stability of archaeal membranes, i.e. they allow archaea to exist under extreme conditions, such as high pressure (~400 atm), high temperatures (~ 100 °C), high methane concentrations and very low or high environmental acidity. 9 In contrast to "conventional" lipids, bolalipids consist of two polar heads and two hydrocarbon chains. In a hydrophilic environment, these lipids self-assemble into monolayers. 10 In addition, bolalipid membranes are considered promising materials for various scientific and engineering applications, 11,12 which emphasizes the necessity of investigating their distinctive mechanical properties.
Theoretical investigations of conventional lipids' mechanics have been carried out in the framework of microscopic and macroscopic models. Microscopic models are represented by various molecular dynamic models 13 and analytical solutions of statistical mechanics equations. 14 Macroscopic models use the elasticity theory to treat membranes as a continuum elastic medium. Here we have focused on the lipid membrane elasticity theories. The first elasticity theory for conventional lipid membranes was developed by Helfrich. 15 Despite the simplicity of Helfrich's model, it was successfully utilized for theoretical investigations of membrane structures and membrane-associated phenomena. [16][17][18][19] Another big advance towards complete elasticity theory was work accomplished by Hamm and Kozlov,20 in which the authors accounted for the bilayer's intrinsic structure within the framework of so-called tilt deformation. This theory is still widely used for the investigation of various membrane processes and phenomena, such as poration, fission, fusion, and domain formation. 6,[21][22][23][24][25][26][27][28][29] These theory-based models enable the systematization of available experimental data and possess substantial predictive power. However' the afore-mentioned elasticity theory still requires an adaptation for bolalipid membranes.
Bolalipids have been experimentally investigated for a long time. 11 However, not much theoretical research has been carried out' and all of it was completed in the framework of microscopic models: by means of molecular dynamics 13,30 and analytical solutions of equations of statistical mechanics. 14,31 A macroscopic elasticity theory for bolalipid membranes has not yet been developed.
Bolalipids have two conformations: (1) the so-called, O-shapes, in which polar heads are located on different sides of the membrane (Fig. 1a); (2) the so-called U-shapes, in which both polar heads are located on the same side of the membrane (Fig. 1b and c). content depends on the particular experimental setup and may reach up to 60%. A membrane that mainly consists of O-shapes is likely to differ in its mechanical properties from a membrane in which two monolayers interact at the membrane midplane. The elasticity formalism has not yet been developed for the bolalipid membranes. The main aim of the present work is to fill that gap.
Firstly, we derive a general expression for the energy surface density of bolalipid membranes that exclusively consist of O-shaped lipids. As a starting point, we use the general elasticity theory of lipid membranes. 20 Secondly, we consider U-shapes' contribution to the elastic energy. Thirdly, we suggest possible experiments and theoretical estimations for defining elasticity moduli and others parameters of the model.

Statement of the problem
We treat the membrane as a continuous medium which can be subjected to elastic deformations. We aim for the development of an elasticity theory for bolalipid membranes and assume that all deformations are small, so we calculate their energy up to the second order. First, we will only consider bolalipid membranes that consist of O-shapes which thus possess mirror symmetry with respect to the surface that passes through the middle of the membrane in its undeformed state, the so-called "midplane".
We abide by the previously established algorithm for conventional lipid monolayers. 20 For convenience, we reproduced the basic equations without excessive mathematical details. Eqn (1) is the general expression for the elastic energy F of a laterally liquid medium, written up to the second-order term: 20 dF = dV σ L ε + 1 2 λ L ε 2 + 1 2 4λ T u zβ u zβ , (1) where u i denotes the components of the displacement vector u = r -r 0 in local tangential Cartesian coordinates Oxyz. r 0 and r are the radius-vectors of the volume element in the non-deformed state and deformed states, respectively. u ij is the deformation tensor components: u ij = 1 2 ∇ i u j + ∇ j u i + ∇ u T ∇ u ji ; σ L , λ T , and λ L are elastic moduli. u's components are related to the lateral expansion ɛ of a volume element by the volumetric incompressibility condition: (1 + ɛ)(1 + Δ z u z ) = 1. With the second order terms of ɛ it reads: Δ z u z = -ɛ + ɛ 2 +… The deformations are further written in terms of ɛ rather than through u, for convenience.
The final expression for F is written in terms of splay and tilt deformations. 20 Tilt deformations are characterized by the tilt-vector t. It describes the deviation of the average direction n (also called "director", the unit vector) of lipids from the normal N to the membrane surface: t = n/(nN) -N. Splay deformations are characterized by the mean curvature J and the Gaussian curvature k of the pivotal surface of lipid monolayers. By definition, a surface is called pivotal when it does not stretch upon splay. Experimental evidence locates the pivotal surface of "normal" phospholipids in the region of the carbonyl groups. 18 The curvatures are also expressed through n: J = -div(n), k = ∂n x ∂x ∂n y ∂y . Bolalipid membranes' splay deformations (Fig. 2e) would locate the pivotal surface in the vicinity of the membrane midplane. However, for symmetric barrel-like deformations (Fig.  2d), two pivotal surfaces are required as this deformation resembles the symmetric splay of conventional lipid bilayers. They should be located near the head-group regions. It thus appears more convenient to abandon the pivotal surface and to define all deformations with respect to the surface at the membrane midplane. The drawback of such an approach is that we can no longer consider splay and compressing/stretching to be independent of each other.
A vector field of unit normal N to the midplane defines the shape of the midplane. Characterizing membrane deformations requires a pair of unit field vectors n. Otherwise the membrane would be reduced to an infinitesimally thin film with some independent internal structure, merely defined by the bolalipid's tilt. Such a description would neither capture highly curved membranes nor the influence of local disturbances, like those displayed upon protein insertion. Moreover, it would fail to describe the case of asymmetric content of Ushape molecules (see Fig. 1b). With a pair of n, the average orientation of bolalipids in the upper and the bottom parts of the membrane can be described relative to the midplane (see Fig. 2). All parameters corresponding to the upper and lower membrane halves will be denoted by indices 1 and 2 (see Fig. 2). In the unstrained symmetric membrane, the midplane is flat and N, n 1 and n 2 are collinear. Bolalipid membranes are considered both laterally liquid and locally volumetrically incompressible 6,20,[24][25][26][27][28][29] which are similar to membranes made from conventional lipids.

Solution to the problem
We introduced an additional local Cartesian coordinate system Oxyζ. Its origin O coincides with the local tangential Cartesian coordinate system Oxyz that has already been introduced. It is located at the midplane. For both systems N is directed along the z axis and Oxyz forms the local tangential basis. However, in the Oxyz system, the scale of the Oz axis does not change upon deformation of a membrane element; in the Oxyζ system the Oζ axis linearly scales with the deformation along this axis (see Fig. 2a and b).

Tilt deformation of bolalipid membranes
In tilt deformation, both directors deviate from N ( Fig. 2c). Membrane thickness remains unchanged, since the membrane is volumetrically incompressible (compare Fig. 2a and c).
Tilt can be described by the following dependence of u on z 20 : For small deformations: t ≈ n -N. The only nonzero deformation tensor components u z β take the following form: u zβ = 1 2 t 1, β · θ z + t 2, β · θ −z β = x, y , where θ(z) is the Heaviside step function, defined as θ(z) = 0 for z < 0 and θ(z) = 1 for z > 0. By inserting these deformation tensor components into eqn (1) to find the integral over the membrane thickness, we obtained the contribution of tilt deformation to F: where h is half of the membrane thickness. Eqn (3) does not contain any cross-terms on tilts t 1 t 2 from the upper and bottom parts. The reason is that the cross-terms would be linked to contributions from the average curvature of the lipid hydrocarbon tails. They are negligible in comparison to F, as indicated by the observation of a substantial population of U-shapes in bolalipid membranes. 10,13 That is, the energy of such a significantly curved hydrocarbon chain is comparable with the characteristic energy of thermal fluctuations, k B T ~ 4 × 10 −21 J. Consequently, it is safe to assume that chain bending occurs at a negligible energetic expense.

Splay and compressing/stretching of bolalipid membranes
The local curvature J of a lipid monolayer is given as 21 J = -div(n). Splay does not lead to the shearing of volume elements so that u z β = 0. Splay contributions to F are due to the stretching of the hydrocarbon region (ε ≠ 0). Treating a small deformed patch of conventional lipids in terms of a curvilinear trapezium 20 yields proportionality between e and three variables: J, k, and α (the relative area change of the whole membrane): ε = α + ζJ + ζ 2 k. ζ is the distance between the midplane and the volume element ( Fig. 2d and e). α adopts values different from zero if the membrane is subjected to lateral tension, σ. For further calculations, we note that J and k are of the first and second orders of smallness, respectively. 20 In contrast, deformations of bolalipid membranes are parameterized by two pairs of curvatures: J 1 , k 1 and J 2 , k 2 . Thus, stretching a volume element located at distance ζ from the midplane takes the following form: ε = 1 + ε m + α 1 + ζJ 1 + ζ 2 k 1 θ ζ + 1 + ε m + α 1 − ζJ 2 + ζ 2 k 2 θ − ζ − 1, (4) The equilibrium midplane stretching value, ɛ m , upon splay deformation (see Fig. 2d) is found by minimizing F.
For further calculations, we switch from the Oxyζ coordinate system to the Oxyz coordinate system using the volumetric incompressibility condition for a membrane patch of area A 0 : A 0 z = A 0 ∫ 0 ζ 1 + ε ζ′ dζ′. z may be expressed via ζ as: Substituting eqn (5) into eqn (4), we obtain ɛ(z): z 2 + k 1 z 2 , z < 0 . (6) We obtain the splay contribution to the free energy, F j by inserting eqn (6) into eqn (1) and zeroing u zβ (pure splay). Assuming laterally uniform tension allows for the minimization of F J with respect to ɛ m : where the elastic moduli are defined as follows: for symmetric membranes. B d and B s respectively characterize the splay of the whole membrane (the curvatures of the membrane parts above and below the midplane with equal absolute values and opposite signs) and the intrinsic membrane splay that acts to preserve a flat membrane on average (the curvatures of the membrane parts are equal, both in absolute value and sign). K G is the Gaussian curvature modulus. With these definitions, F J of an arbitrarily deformed small patch of a bolalipid membrane can be expressed as: where J ss is determined by B s J ss = τ 1 = τ 2 . J ss is similar to the spontaneous curvatures of conventional lipid membranes.
We have disregarded mixed deformations, such as simultaneously occurring splay and tilt, since they are energetically decoupled as has been derived for conventional lipids. The reasoning is that a linear vector term cannot be part of an energy expression. However, the second-order cross-term of the scalar quantity splay (div(n)) and the vector quantity tilt (tiltvector t) does not obey that requirement. Thus, tilt and splay must be considered independently.
In contrast to the expressions for F of conventional lipid bilayers, the cross-term for the curvatures of opposing membrane parts exists in the corresponding expressions for bolalipids. The cross-terms for the tilts of opposing membrane parts, for α and J 1 or J 2 are absent. That does not transform the midplane surface into a neutral one, since midplane stretching ɛ m still depends on curvature (eqn (7)). It only means that the energy contributions from the deformations induced by the lateral tension s and by the applied torques are independent of each other. The cross-term in J originates from the fact that the Europe PMC Funders Author Manuscripts upper and lower halves share a common midplane. Model accuracy allows us to neglect the Gaussian curvature cross-terms.

Limiting case
O-shape membranes topologically differ from conventional lipid bilayers. Nevertheless, symmetric deformations of both bolalipid and conventional membranes should be describable via the same equations, because their midplanes are similarly deformed. For conventional lipids, the pivotal surface is located in the region of the carbonyl groups. 18 In our model, the pivotality (inextensibility) of this surface is equivalent to the infinitely large stretching modulus in the carbonyl group region, λ L (z = h) → ∞. We thus obtain ℎ, which matches the corresponding value in the model for conventional monolayers. 20 Upon substitution of λ L (z = h) → ∞ eqn (9) reduces to In the penultimate equation we substituted h -z for z′ to match the previously used notations. 20 We obtained J ss by assuming a mechanically stable membrane: ∫ 0 ℎ σ L z dz = 0.
Thus, for defined conditions eqn (10) coincides with the corresponding equations for conventional monolayers. 20

Spontaneous curvature
On the one hand, the spontaneous curvature arises in an expression for F for a bolalipid layer composed only of O-shape bolalipids (eqn (9)). On the other hand, it is evident that the Ushape lipids would make a stronger contribution to the spontaneous curvature. Combining the contributions of O-shapes and U-shapes is a non-trivial task since bilayers from U-shape molecules do not necessarily obey eqn (9). They are best described by the equations for conventional lipid membranes (eqn (1)). However, if we consider that the concentration of U-shapes is usually below 10%, 10,13 we can assume that the U-shape concentration is of the same order of smallness as the elastic deformations. Thus, the effect of U-shapes on the elastic moduli must be negligible, since accounting for it would involve third or fourth order corrections to F. Under such conditions, the presence of U-shapes affects only those terms that are linear on deformations, i.e. τ 1 , τ 2 . Thus, we may limit all further considerations to the linear terms x 1 and x 2 of U-shape concentrations, since the limited accuracy of the model does not allow us to discern the contribution of higher order terms. Assuming that the headgroups of the U-shapes localize toward the external surface of the bolalipid membrane (Fig. 1b), we must account for different areas a 1 and a 2 of the two head-group regions.
Repeating the derivation of eqn (8) for asymmetric membranes yields eqn (11) for its parameters: The coefficient τ u reflects the additional stress induced by the presence of U-shapes. τ 1 and τ 2 exclusively act to affect J. Since the per surface areas occupied by one U-shape molecule and by one O-shape molecule differ twofold, the ratio (a 2 -a 1 )/a 1 can be estimated as (a 2 -a 1 )/a 1 ≈ x 2 -x 1 . Thus, the effect of up to 10% U-shapes on F J of bolalipid membranes may be taken into account by inserting into eqn (9) J's dependence on x 1 and x 2 : where J s0 = τ 0 /B s , J su = τ u /B d + τ A /B d . The spontaneous curvatures J su and J du reflect the contribution of U-shapes. J su is an additive to J s0 , the spontaneous curvature of O-shapes. Jsu and J du are linear in composition since the concentration of U-shapes is small, and there is no reason to suspect non-ideal mixing (e.g. phase separation). Even if nonideal mixing would occur, it is unlikely, that its effect would significantly alter the model prediction, since the model is of limited accuracy.
Eqn (12) ignores the entropic contribution of mixing U-shapes with O-shapes, which should be encountered when the deformational energy is comparable to or smaller than the thermal energy k B T. Any lateral inhomogeneity of U-shapes may favor membrane deformations that are laterally non-uniform. For this reason we estimate the spontaneous curvature of Ushapes below.

Elasticity modulus of tilt
The elasticity moduli B s , B d , K t , K G , K A (compare eqn (12)) should be measured experimentally, calculated from microscopic models, or otherwise assessed. Simple calculations show that K t of conventional lipids should be close to the surface tension of the oil-water interface, 20 which was experimentally confirmed. 32 Extending the same considerations to bolalipid membranes, we estimate its K t to be equal to that of conventional lipids, i.e. K t ~ 50 dyn cm -1 .
All other elasticity moduli depend on lipid structures and properties, thereby precluding this type of simple estimation. They should be experimentally measured. However, assessing K G is very difficult even in the case of conventional lipids. At the same time, K G G only needs to be accounted for in a narrow and peculiar set of problems, in which membrane topology changes. Below we focus on how to estimate B s and B d .

1.
B d . K a and B d measurements are commonly based on monitoring the increment in vesicular membrane area upon application of hydrostatic pressure. The change in surface area of a giant unilamellar vesicle (GUV) is associated with undulations and expansion of the area per lipid molecule. [33][34][35][36] GUVs with a diameter of about 10 μm are well suited for this purpose because the average curvature is small. Since J 1 and J 2 have different signs, J 1 + J 2 is much smaller than J 1 -J 2 . This means that B d and K A can be determined in such experiments. [33][34][35][36] The energetic contribution of the Gaussian curvature is constant because the system's topology does not change during the experiment (Gauss-Bonnet theorem).
The indices "1" and "2" correspond to external and internal parts of the membrane that form the nanotube, respectively. We define tube radius R = 1/J at the membrane midplane; h is equal to half of the membrane thickness, J ss and J sd are spontaneous curvatures (eqn (10)).
The energy density is multiplied by the area of the non-deformed state, 40 which with sufficient accuracy may be assumed to be equal to the area of the nanotube midplane. F given by eqn (13) should be minimized with respect to J and the concentration of U-shapes, x 1 and x 2 . As a result we find R as a function of σ. B s is then obtained by varying σ via the application of transmembrane voltage. 37 For conventional lipids, the elastic moduli are much greater than the characteristic energy of thermal fluctuations, k B T. For instance, B s is about 36 10k B T. Similarly, we may thus assume that the lateral distribution of U-shapes is only governed by F. The formation of nanotubes occurs much faster than the lateral redistribution of membrane components with non-zero spon-taneous curvature 39 (U-shapes). Consequently the U-shape concentrations in the internal and the external parts of the nanotube membrane immediately after formation are equal to the U-shape concentration in a flat membrane and J sd = 0. Minimizing F given by eqn (13) with respect to nanotube curvature, we obtain B d : − B s 1 + 3 ℎ/R 2 1 − ℎ/R 2 ℎJ ss − 6 ℎ/R 2 1 + ℎ/R 2 3 ℎ/R 4 + 8 ℎ/R 2 + 1 (14) R' is the derivative of R with respect to σ (R' < 0 for real systems); R' could be measured experimentally. The expression can be simplified for large R (h/R ≪ 1): The nanotube radius subsequently relaxes due to the lateral redistribution of U-shapes in both monolayers. The relaxation is governed by the independent minimization of elastic energy in each monolayer. As a result, a local effective transversal asymmetry of U-shaped molecules may emerge. The characteristic time amounts to about 1 s for conventional lipids (dioleoyl-phosphatidylethanolamine, DOPE). 39 The resulting equilibrium curvature of the nanotube as well as x 1 and x 2 can be obtained by minimizing F (eqn (13)) with respect to R, x 1 , and x 2 . Energy minimization demands the absence of U-shapes (possessing positive J s ) on the internal part of the nanotube having negative curvature. This yields the expression for B d in the limit of large R: Thus, B d can be obtained from experiments with GUVs, while B s and J ss are attainable by measuring nanotube radii in and out of equilibrium.

Discussion
We have obtained a general expression for the surface energy density of elastic deformations for bolalipid membranes which consist of two types of molecules: O-shapes and U-shapes. F includes cross-terms for (i) curvatures of opposing membrane parts and (ii) curvatures with U-shape concentrations. Tilt cross-terms are absent because they are determined by the average bending of hydrocarbon chains which make a negligibly small contribution to F. For specific cases where transmembrane peptides are embedded in a bilayer, the absence of the tilt cross-term would no longer apply. In the framework of our model, however, this should be taken into account by means of coupling boundary conditions set on deformations rather than the free energy functional.
Our theoretical considerations should help to ascertain the major differences in the roles archaeal and conventional lipids play during cellular processes that involve membrane reshaping. Archaea possess at least three distinct membrane remodelling systems. 41 The first uses an archaeal actin-related protein, the "cell division A" CdvA protein. The second is comprised of the bacterial-type system FtsZ. The third alternative cell division apparatus is homologous to the eukaryotic ESCRT-III (endosomal sorting complex required for transport). Remarkably enough, membrane scission by the yeast ESCRT-III complex does not require a special type of lipid in addition to certain amounts of anionic lipids to preserve a negative net charge 42 -a requirement that can easily be met by bolalipids. However, the energetics of scission should be fundamentally different when bolalipids are involved because their elastic moduli are different. In contrast to scission, early stages of fusion, i.e.
hemifusion and fusion pore formation depend on lipid curvature, 43 and it would be interesting to see how bolalipids may meet these requirements. To get a first impression about the energetics involved in membrane remodelling, we will estimate the basic elastic parameters of bolalipids from our theory. The spontaneous curvature of a monolayer composed of U-shapes can be estimated using a toy-model. Symmetrical insertion of U-shaped bolalipid molecules into a membrane that consists of O-shapes (Fig. 3a) enlarges the membrane surface, S h , more than it increases the midplane area, S t . In the limiting case of a pure U-shape monolayer, the spontaneous curvature is positive, since all polar headgroups are located at the same side of the membrane ( Fig. 3a and b).
If a h is the area per lipid headgroup of both O-shapes and U-shapes, and if a t is the area of an O-shaped or U-shaped molecule at the membrane midplane, we attain the following expressions for the membrane surface areas: where h 0 is the equilibrium thickness of the monolayer from U-shapes, which can be taken as equal to h-half of the thickness of membranes made from O-shapes; N is the number of total lipid molecules. eqn (17) assumes linear dependence of the headgroup and tail region areas S h and S t on x. J 1,2 can be found as: For vesicles made from conventional lipids, Eves does also not depend on radius, but with only 500k B T it is tenfold smaller. However, E ves of bolalipid vesicles decreases significantly when more U-shapes face the outer surface than the vesicle lumen. To account for lipid asymmetry, we fix R and minimize F J (eqn 12) with respect to the U-shape concentrations in both membrane halves, x 1 and x 2 . This allows us to obtain the equilibrium concentrations x 1 = 4h/R, x 2 = 0, which in turn yield E ves ≈ 4πR 2 · 4B s R 2 = 16πB s ≈ 500k B T . Thus lipid asymmetry brings E ves down to the value obtained for conventional lipids. However, the reduction of E ves comes at a substantial cost: for vesicles with an R of 40 or 100 nm, U-shapes have to be enriched in the outer half by x 1 = 20% or x 1 = 8%, respectively.
Based on the estimates for E ves , we expect that the fusion of bolalipid membranes requires an asymmetrical U-shape distribution. Components with non-zero spontaneous curvature substantially alter the rate of membrane fusion even in the case of conventional lipid membranes. 3,44 It decelerates if the contacting (proximal) leaflets have positive spontaneous curvature, and accelerates if the positive spontaneous curvature is acquired by the distal monolayers. Thus, enrichment of U-shapes in the distal halves of the membrane should facilitate membrane fusion. The asymmetry has to be locally restricted to the fusion zone. While lysolipids, which play that role in conventional membranes, may be produced at little cost by phospholipases and selectively enriched by protein imposed curvature, the corresponding mechanisms in bolalipid membranes are not known. Both translocation of a charged bolalipid head-group and protein-induced bending of the bolalipid membrane are certainly energetically much more costly than in the case of conventional lipid bilayers.
We conclude that bolalipids' unique chemical structure sustains the unique stability of archaeal membranes. Their self-assembly into a monolayer, instead of into a bilayer as is the case with conventional membranes, should significantly hinder membrane reshaping by fusion and fission. From our estimations, bolalipids possess a splay modulus that is an order of magnitude larger than that of conventional mammalian lipids. It acts to further hamper fusion and inhibit pore formation, thus allowing archaea to maintain the membrane barrier to ions and other molecules even in extremely aggressive environments. However, the projected price for this stability is rate deceleration in cell division or endocytotic uptake.