Ligand-receptor-mediated attachment of lipid vesicles to a supported lipid bilayer

The interaction of exosomes (cell-secreted \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim $$\end{document}∼100 nm-sized extracellular vesicles) or membrane-enveloped virions with cellular lipid membranes is often mediated by relatively weak ligand-receptor bonds. Interactions of this type can be studied using vesicles and observing their attachment to receptors located in a lipid bilayer formed at a solid surface. The contact region between a vesicle and the supported lipid bilayer and accordingly the number of ligand-receptor pairs there can be increased by deforming a vesicle. Herein, I (i) estimate theoretically the corresponding deformation energy assuming a disk-like or elongated shape of vesicles, (ii) present the equations allowing one to track such deformations by employing total internal reflection fluorescence microscopy and surface plasmon resonance, and (iii) briefly discuss some related experimental studies.


Introduction
Lipid vesicles can often directly attach to solid surfaces and be intact or ruptured and form a supported lipid bilayer [SLB; see e.g. the articles by Cho et al. (2010), Mapar et al. (2018) and references therein]. Such SLBs can also be fabricated by using other approaches including e.g. the solvent-assisted lipid bilayer and bicelle methods (Jackman and Cho 2020). In applications, an SLB can contain receptors, and this platform allows one to study various aspects of the ligand-receptor-mediated interaction of vesicles, membrane-enveloped virions, virus-like particles, or lipid nanoparticles (of ∼ 100 nm size) with lipid membranes as schematically shown in Fig. 1. In nature, such interactions between exosomes (cell-secreted ∼100 nm-sized extracellular vesicles) or membrane-enveloped virions and cellular lipid membranes are extremely important from various perspectives as reviewed by Kalluri and LeBleu (2020) and Mercer et al. (2010). The interaction of lipid nanoparticles with lipid membranes is also of high current interest from the perspective of the development of new drugs (Tibbitt et al. 2016;Zhdanov 2019).
If the attachment of a vesicle is mediated by one ligandreceptor pair, a vesicle can keep its shape, e.g., remain spherical (Fig. 1a). The formation of two or more ligand-receptor pairs may require appreciable vesicle deformation provided the distance between the pairs is comparable to or larger than the vesicle radius, R (Fig. 1b,c). If there are a few ( n > 1 ) or many ( n ≫ 1 ) ligand-receptor pairs, the gain in the energy due to the formation of additional ligand-receptor bonds (i.e., those formed after the formation of the first bond) is −(n − 1)I , where I > 0 is the bond energy, while the loss of energy is equal to the increase of the lipid-bilayer bending energy, ΔE b (provided the osmotic pressure is negligible). Thus, the formation of extra ligand-receptor pairs is favourable if This simple condition neglecting the entropic factors shows that the shape of a vesicle is determined by the interplay between the ligand-receptor-bond formation and membrane bending.
Experimentally, the deformation of vesicles during attachment to solid surfaces in general and in the case of ligandreceptor mediation in particular can be tracked optically e.g. by using total internal reflection fluorescence microscopy [TIRFM; this technique is reviewed by Boukobza et al.
1 3 (2001) and Olsson et al. (2015)], surface plasmon resonance [SPR; reviewed by Jung et al. (1998) and Rupert et al. (2016)], or localized surface plasmon resonance [LSPR; reviewed by Jackman et al. (2017)]. Vice versa, one can say that the vesicle deformation influences the TIRFM, SPR, and LSPR signals, and this effect can be important from the perspective of related applications of these techniques. Other techniques, e.g. crystal microbalance-dissipation [QCM-D; reviewed by Cho et al. (2010)], can be used here as well. In particular, some information about the deformation of vesicles directly attached to the support was already obtained by employing SPR [with emphasis on determining the size and concentration of sub-populations of extracellular vesicles (Rupert et al. 2016)], LSPR [with emphasis on the kinetics of vesicle attachment (Jackman et al. 2014), role of temperature (Oh et al. 2015) and divalent cations (Dacic et al. 2016), and deformation itself , and QCM-D [with emphasis on the role of osmotic pressure (Jackman et al. 2013)]. Despite these advances, one needs additional theoretical input for the use of these techniques.
Herein, I estimate ΔE b for n ≥ 2 and present the equations allowing one to calculate the scale of the vesicledeformation-related change of the TIRFM and SPR signals. The corresponding results are of intrinsic interest in the context of biophysics of multivalent interactions and also may potentially be useful in the context of development of sensors allowing the measurement of low concentrations of biological molecules via their attachment to a solid surface and playing the role of receptors by tracking their association with ligand-containing vesicles which can be employed as signal amplifiers.

Change of the bending energy
In the generic case without spontaneous curvature of a lipid bilayer, the vesicle bending energy can be represented as (Seifert 1997) where is the bending rigidity, c 1 = 1∕r 1 and c 2 = 1∕r 2 are the principal curvatures, r 1 and r 2 are the corresponding radii, and the integration is performed over the vesicle surface (ds is an element of the surface).
For spherically shaped vesicles with c 1 = c 2 = 1∕R , the bending energy is given by This conventional expression for the bending energy can be used to describe vesicles attached to the surface by one ligand-receptor pair (Fig. 1a).
In the case of a few ( n > 1 ) or many ( n ≫ 1 ) ligand-receptor pairs, the shape and exact value of the vesicle bending energy can be calculated via cumbersome numerical calculations (Seifert and Lipowsky 1990;Irajizad and Agrawal 2018). The results of such calculations are neither transparent nor convenient in applications. To get transparent results, I use simple physically reasonable approximations for the shape of deformed vesicles.
In particular, the vesicles attached to the surface via two ligand-receptor pairs ( n = 2 ) are represented as an elongated pellet ( Fig. 1b) with hemispherical left and right parts (with radius r) and cylindrically shaped central part (with radius r and length h; this length can be identified with the distance between the two ligand-receptor pairs). In this approximation, the curvatures of the hemispherical parts are c 1 = c 2 = 1∕r , the curvatures of the cylindrical part are c 1 = 0 and c 2 = 1∕r , and the vesicle bending energy is given by The change of the bending energy is obtained by subtracting the initial energy (3), i.e., The two sizes used in (4) and (5) are not independent, because the lipid-bilayer areas of a vesicle before and after deformation, 4 R 2 and 4 r 2 + 2 rh , should be equal. This condition yields According to this equation, r can be represented as a function of h, or h can be represented as a function of r, In addition, it is convenient to introduce the maximum size of a deformed vesicle, This relation can be rewritten as h = L − 2r . Substituting the latter relation into (6) yields Substituting this relation into (8), we obtain.
The vesicles attached to the surface via three or more ligand-receptor pairs ( n ≥ 3 ) can be viewed as a rounded disklike pellet with the peripheral radius r and cylindrically shaped central part of diameter h (Fig. 1c). Adopting this approximation, I employ the disk cross-section perpendicular to the surface and the maximum cross-section parallel to the surface to characterize the shape of the peripheral area. The corresponding curvatures are 1/r and 1∕(r + h∕2) , respectively. With these curvatures, the vesicle bending energy is represented as After subtracting the initial energy (3), the deformationrelated change of the bending energy is given by To relate h and r, I again take into account that the lipidbilayer areas of a vesicle before and after deformation, 4 R 2 and 4 r 2 + 2 hr + h 2 ∕2 , should be equal, or This equation yields (17) 4r 2 + hr + h 2 ∕2 = 4R 2 . Fig. 2 Increase of the vesicle bending energy as a function of h/R according to (12) and (16) The maximum size of a deformed vesicle is defined by analogy with (9), Using (18), (19), or (20), the change of the bending energy given by (16) can be expressed as a function of h, r, or L (see e.g. Fig. 2).

TIRFM, SPR, and vesicle deformation
In TIRFM or SPR experiments with vesicles immobilized at the interface, the signal measured is induced by the evanescent field. In particular, the signal can in general be described as [see, e.g., the articles by Boukobza et al. (2001) and Olsson et al. (2015) for TIRFM and by Jung et al. (1998) and Rupert et al. (2016)

for SPR, respectively]
where A is a constant proportional to the concentration of vesicles at the interface, J 0 is the intensity (square of the amplitude) of the incident light, is the decay length of the light intensity, z is the coordinate perpendicular to the interface, and the integration is performed over the vesicle surface.
For spherically-shaped vesicles ( Fig. 1b; Eq. 3), expression (21) yields (Olsson et al. 2015;Rupert et al. 2016) If vesicles are represented as an elongated pellet (Fig. 1b;, expression (21) results in where is the modified Bessel function of the first kind and order zero. Expression (23) can be rewritten as (23) is the dimensionless factor ( ≤ 1 ) describing the effect of the vesicle deformation on the TIRFM or SPR signal. If vesicles are represented as a rounded disk-like pellet ( Fig. 1c; Eqs.

Discussion and conclusions
The results presented (Fig. 2) indicate that the vesicle bending energy becomes comparable with 2 already for modest vesicle deformations, at h∕R ≃ 0.5 and 0.25 for the elongated and disk-like shapes, respectively. Such deformations can be tracked by using TIRFM or SPR (Fig. 3). My analysis forms a basis for the corresponding experiments. The first part of the analysis ("Change of the bending energy") can also be used for the interpretation of the data obtained by employing other techniques (e.g., QCM-D as briefly discussed below). (27) (b) Fig. 3 Factor describing the effect of the vesicle deformation on the TIRFM or SPR signal as a function of h/R at ∕R = 0.5 , 1, and 2: a for elongated vesicles according to (26), and b for disk-like shape of vesicles according (27) In particular, it allows one to discuss the interplay between the bending energy and the free energy of the formation of the ligand-receptor pairs. This interplay is complex, and usually one set of measurements does not make it possible to accurately estimate all the desirable parameters. Some limitations on or relation between the parameters can, however, be obtained. In the theoretical literature, one can often read that ≃ 25 k B T (Smith et al. 2004). The recent experiments with ∼ 100 nm-sized vesicles indicate that may be appreciably larger, up to ≃ 10 3 k B T , even in conventional lipid bilayers [see the LSPR study by Jackman et al. (2017) and references therein]. For biomimetic lipid compositions containing e.g. cholesterol [reviewed by Cebecauer et al. (2018) and Enkavi et al. (2019)], can be appreciably larger as well. Even if is relatively small, ≃ 25 k B T , the above-mentioned increase of the bending energy ( ≃ 2 ≃ 50 k B T ) is appreciable and can hardly be fully compensated by a few relatively weak bonds (smaller or about 5 kcal/mol or ≃ 8 k B T ) which are typical for the multivalent ligand-receptor interactions. With increasing vesicle size and/or length of ligand-receptor pairs, the number of ligand-receptor pairs located in the contact area can be large, and they can induce appreciable deformation of a vesicle.
To link the general conclusions above to real systems, one can look through a recent QCM-D study of attachment of biotinylated small (SUVs, ≃ 100 nm in diameter) and giant (GUVs, ≃ 15-20 m in diameter) unilamellar DOPC vesicles to a biotinylated SLB functionalized with streptavidin (Di Iorio et al. 2020). One of the advantages of this choice of species is that the biotin-streptavidin interaction has long been used for biomolecule immobilization, and its strength is well established [ K d ≃ 10 −14 M (Di Iorio et al. 2020) or ΔG = 18.3 kcal/mol (Weber et al. 1992)]. In the experiment, the fraction of biotinylated lipids, f, in an SLB was varied from 0.002 to 0.02, whereas in vesicles it was from 0.001 to 0.02. Analyzing the QCM-D data obtained for SUVs, the authors have concluded that the deformation of SUVs was modest (with the SUV surface-contact area ≤ 10% ) provided f ≤ 0.006 (in an SUV), and under such conditions the number of biotins (on the SUV side) involved in the contact was ≤ 130 . The formation of ligand-receptor pairs in the SUV-SLB contact area is accompanied by the loss of entropy of ligands and receptors. In the context under consideration, this loss per pair can be estimated as ≃ −2k B T ln(f ) or ≃ 6 kcal/mol (provided f = 0.006 ). The bond-formation gain in the free energy can be estimated by subtracting this value from ΔG , i.e., about 18-6=12 kcal/mol per pair. For 130 pairs, this gain is ≃ 1500 kcal/mole or ≃ 2500 k B T . Identifying the latter value with 2 (as suggested in the first paragraph of this section), I obtain ≃ 1200 k B T . Thus, this analysis also indicates that can be much larger than 25 k B T . Concerning GUVs, their deformation of attached was observed explicitly by using fluorescence microscopy (Di Iorio et al. 2020), but the discussion of the corresponding results is beyond our goals.
In another study, the deformation of SUVs ( ≃ 100 nm in diameter and of complex composition) adhered to poly-L-lysine coated glass was studied by using AFM (Vorselen et al. 2017). Referring to the conventional bending rigidity, ≃ 10-50 k B T , and accepting this scale of in the calculations, the authors have concluded that the bilayer bending alone cannot account for the high stiffness observed experimentally and refer to relatively high osmotic pressure (in their analysis, this pressure was used as a fitting parameter and estimated to be ≃ 0.15 MPa) to explain their observations. The discussion above shows that one cannot exclude that the bending rigidity might in reality be much higher than 10-50 k B T.
Finally, I can add that, although my work was motivated by referring to the ligand-receptor-mediated interaction of vesicles with an SLB, the results reported can be used in a more general context. In my analysis, the role of an SLB is in fact reduced to prevention of a direct contact of vesicles with a support, and accordingly the bilayer can be replaced by any other layer preventing such a contact (Kim et al. 2017). If a support itself does not interact with vesicles, it can be directly (without an SLB) functionalized by receptors (Oliverio et al. 2017;Liu et al. 2020). For example, a support can be isotropic or anisotropic and contain parallel protrusions so that it can induce the formation of elongated vesicles (Fig. 1b) after their attachment. In such situations, all the results presented can be valid. Some of the results may be valid even if vesicles interact not only with receptors but also directly with a functionalized support.