The Scientific World JournalVolume 2012 (2012), Article ID 983138, 7 pagesdoi:10.1100/2012/983138

Research Article

Fusion Pore Diameter Regulation by Cations Modulating Local Membrane Anisotropy

Doron Kabaso,1,2 Ana I. Calejo,³ Jernej Jorgačevski,2,4 Marko Kreft,2,4,5 Robert Zorec,² and Aleš Iglič¹

¹Laboratory of Biophysics, Faculty of Electrical Engineering, University of Ljubljana, 1000 Ljubljana, Slovenia²Laboratory of Neuroendocrinology-Molecular Cell Physiology, Faculty of Medicine, University of Ljubljana, 1000 Ljubljana, Slovenia³Departamento de Biologia e CESAM, Universidade de Aveiro, 3810-193 Aveiro, Portugal⁴Celica Biomedical Center, Tehnoloski Park 24, 1000 Ljubljana, Slovenia⁵Department of Biology, Biotechnical Faculty, University of Ljubljana, Vecna pot 111, 1000 Ljubljana, Slovenia

Received 18 November 2011; Accepted 8 December 2011

Academic Editors: F. Benfenati and W. Voos

Copyright © 2012 Doron Kabaso et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The fusion pore is an aqueous channel that is formed upon the fusion of the vesicle membrane with the plasma membrane. Once the pore is open, it may close again (transient fusion) or widen completely (full fusion) to permit vesicle cargo discharge. While repetitive transient fusion pore openings of the vesicle with the plasma membrane have been observed in the absence of stimulation, their frequency can be further increased using a cAMP-increasing agent that drives the opening of nonspecific cation channels. Our model hypothesis is that the openings and closings of the fusion pore are driven by changes in the local concentration of cations in the connected vesicle. The proposed mechanism of fusion pore dynamics is considered as follows: when the fusion pore is closed or is extremely narrow, the accumulation of cations in the vesicle (increased cation concentration) likely leads to lipid demixing at the fusion pore. This process may affect local membrane anisotropy, which reduces the spontaneous curvature and thus leads to the opening of the fusion pore. Based on the theory of membrane elasticity, we used a continuum model to explain the rhythmic opening and closing of the fusion pore.

1. Introduction

During exocytosis the membrane-attached vesicle can release its cargo through full fusion with the plasma membrane, or it can form a stable membranous neck (i.e., fusion pore) (Figure 1) mediating the cargo release via the repetitive fusion pore openings and closures [1, 2]. The advantages of the latter mechanism are that the vesicle can be reused in another round of exocytosis [3] and that the release process can be further modulated by the interaction between the vesicle lumen and the physical properties of the membranous neck region. The fusion pore is initially narrow (appears closed), and then it widens up to such an extent that the pore can be electrically measured as a transient increase in membrane capacitance. The time interval for the fusion pore opening or closure can range from milliseconds to minutes [4–6]. Recently, it has been shown that the rhythmic opening and closure of the fusion pore [4] can be further increased using a cAMP-increasing agent, which drives the opening of nonspecific cation channels. Yet, the underlying mechanisms for the enhanced fusion pore periodicity remain unclear. In the present paper, we constructed a computational model investigating the contribution of cation fluxes through the vesicle membrane and the fusion pore to explain the observed periodic transient fusion pore events induced by the application of a cAMP-increasing agent.

Figure 1: Schematic illustration of the vesicle fusion with the plasma membrane. Note that the fusion pore is a narrow membranous neck connecting the vesicle to the plasma membrane.

The energetics of the formation and stability of the fusion pore is not trivial. The increase of membrane bending energy, due to the transformation of the fusion pore into a thin tubular structure, is substantial [7, 8]. In addition, the narrow fusion pore increases the electrostatic repulsive forces between negatively charged lipids enriched in the opposing cell membranes [9]. On the other hand, SNARE (soluble N-ethylmaleimide-sensitive factor attachment protein receptor) proteins were suggested to facilitate the fusion of the vesicle with the plasma membrane [10, 11]. In addition, cholesterol and PIP₂ lipids may form a nanodomain raft that stabilizes the highly curved fusion pore region [8, 12]. It was proposed that, without these stabilizing components, the fusion pore will spontaneously close or open irreversibly during complete exocytosis [13]. These stabilizing components may either increase the area difference between the two membrane leaflets or interact unequally with one of the membrane leaflets. Moreover, the increased stability could be due to the reduction in the curvature energy mismatch between the spontaneous curvature of the anisotropic membrane components and the tubular membrane curvature of the pore [8, 14–17]. The deviatoric theory of membrane elasticity has been previously employed to describe the stability of highly curved cell membrane geometries such as dimeric detergent-induced membrane nanotubes [18, 19], nanotubes of giant unilamellar vesicles [18, 19], cell membrane pores [20], narrow necks connecting daughter vesicles with the parent membrane [8, 21, 22], and inverted hexagonal phospholipid phases [23, 24]. The large difference between the two principal membrane curvatures of these structures increases the importance of anisotropic membrane components (either lipids, proteins, or their complexes) that have different spontaneous curvature along the two principal directions of the membrane [8, 22]. For example, a recent computational model has demonstrated that the narrow fusion pore can be reproduced only by the inclusion of anisotropic membrane components into the highly curved membrane nanodomain [8].

The aim of the present paper is to study the underlying mechanisms responsible for the enhanced periodicity of fusion pore events following the application of a cAMP-increasing agent (Figure 2). We constructed a computational model, in which the fusion pore is modelled as a single anisotropic membrane domain (nanodomain). The membrane force exerted to open and close the fusion pore was derived from the differentiation of the membrane bending energy. The cation concentration in the vesicle was dependent on the morphologic features of the vesicle as well as the membrane area fraction of cation fluxes. The cation flux into the vesicle through the area fraction of cation channels increased the cation concentration at the closed (extremely narrow) pore, while the cation flux out of the vesicle through the pore was possible only when the pore was sufficiently wide (open). The rhythmic opening and closure of the fusion pore were reproduced by coupling the spontaneous curvature of the anisotropic membrane nanodomain and the cation concentration at the fusion pore (Figure 3). The physiological consequences of this novel model are discussed.

(a)

(b)

Figure 2: The application of a cAMP increasing agent increases the rhythmic openings and closures of an individual fusion pore, recorded in resting lactotrophs. Membrane capacitance recordings demonstrate the effects of a cAMP increasing agent (dbcAMP, 10 mM) (b) in comparison to the control experiment (a). The capacitance trace demonstrates changes between low and high states, corresponding to effectively closed and open fusion pore states.

Figure 3: Schematic model for the rhythmic opening and closure of the fusion pore. The stable fusion pore of a vesicle can either be in open or closed states (a). The influx of cations via ion channels in the vesicle membrane is constant, while the efflux through the fusion pore out of the vesicle is allowed only when the fusion pore is open (see insets). The spontaneous curvature at the fusion pore is modulated by the cation concentration in the vesicle. The membrane force, derived from the differentiation of the system free energy, is exerted in the radial direction (see arrows) of the modelled axisymmetric fusion pore (see insets). The fusion pore is composed of a single anisotropic membrane domain (b). The membrane spontaneous curvature () of the fusion pore is plotted as a function of the relative cation concentration () (c). The proposed dynamics is that the cation flux through the vesicle increases the cation concentration in the vesicle, which then reduces the intrinsic membrane curvature () leading to the opening of the fusion pore. The opening of the pore causes the abrupt reduction in cation concentration, which increases the spontaneous membrane curvature of the pore leading to the abrupt closing of the fusion pore.

2. The Model

The aim of the present computational model is to explain the rhythmic opening and closing of the fusion pore (Figure 2). A vesicle is attached to the plasma membrane through the narrow fusion pore (Figure 1). It was reported previously that the diameter of the fusion pore may be <0.5 nm [5, 8]. In the model, the attached vesicle is composed of two compartments (Figure 3). The first compartment is a sphere, modelling the vesicle without the fusion pore. The second compartment is an axisymmetric tubular structure, modelling the fusion pore (neck) region. The dimensions of the vesicle and the fusion pore regions are according to previously published data [8]. It is assumed that there is a constant influx of cations through the vesicle membrane and that the efflux of cations occurs only through the fusion pore opening. The two membrane leaflets of the fusion pore are modelled as a single membrane entity. The diameter of the fusion pore is the only morphologic parameter being varied along the simulation. Previously, it has been shown that isotropic membrane components alone are not sufficient for the stability of the narrow fusion pore [8]. On the other hand, anisotropic membrane components are crucial for its stability [8, 21, 22]. For the sake of simplicity, in this work the membrane of the fusion pore includes only a single homogeneous anisotropic membrane domain. To summarize, the underlying assumptions in our model are as follows. (1)The spontaneous curvature of the anisotropic membrane domain (nanodomain) is affected by the cation concentration in the vesicle. The accumulation of cations at the entrance of the fusion pore leads to lipid demixing, which reduces the spontaneous curvature of the anisotropic domain thereby opening the fusion pore. (2) The flux through the pore reduces the cation concentration in the vesicle, enabling the reorganization of the anisotropic membrane domain in the fusion pore and the increase in its spontaneous curvature, which closes the fusion pore (Figure 3).

2.1. Derivation of the System’s Free Energy

The membrane bending (free) energy is considered within continuum approach [25–27]: where is the membrane bending modulus, is the mean membrane curvature, is the curvature deviator of the membrane, is the spontaneous membrane curvature and is the intrinsic (spontaneous) membrane curvature deviator, and are the principal membrane curvatures, is the surface membrane tension, and is the infinitesimal membrane area element. It can be seen from (1) that the membrane bending properties can be expressed in a simple way as a function of and . The difference between the membrane curvature and the membrane intrinsic (spontaneous) curvature determines the energy cost for bending the membrane away from its favourable curvature. The cation concentration in the vesicle () affects the intrinsic membrane curvatures, as follows: where and are the intrinsic membrane curvature and membrane curvature deviator, respectively, and is the maximum concentration of cations. Note that the spontaneous curvature of the anisotropic membrane nanodomain depends on the normalized cation concentration ().

The anisotropic membrane nanodomain has cylindrical intrinsic curvatures, where (see Figure 3(b)), and, consequently, . The radius of the fusion pore is where is a small deviation from the initial equilibrium radius of the tubular fusion pore. More details on the curvatures of the tubular fusion pore are outlined in the Supplementary Material available online at doi:10.1100/2012/983138.

2.2. Equilibrium Model

The equation of motion of the membrane height deflection (i.e., in the radial direction) along the cylindrical main axis is given by [28, 29] the following where is the friction coefficient describing the drag of the fluid surrounding the membrane. The forces (per unit area) derived from the differentiation of the free energy (1) are outlined in the Supplementary Material. The fusion pore and the vesicle share the same cation concentration. As a result, the effects of cations on the spontaneous curvature of the fusion pore is uniform throughout the pore (Figure 3).

By a small perturbation to the uniform initial state, we obtain where is a small deviation from the uniform value.

The linearization of (4) yields where the function describes the small undulation in the membrane force, and describes higher-order terms of the membrane undulation. The force acting on the membrane at equilibrium state is described by .

The equilibrium radius is obtained by assuming that the undulation equals zero. The following equation was solved for : The resulted equilibrium radius of the fusion pore is

2.3. Dynamic Model

The fusion pore is allowed to evolve only along the radial direction. We here consider the second layer of model complexity, which is the cation dynamics. The cation concentration is increased due to the influx through the vesicle membrane, while the efflux through the fusion pore reduces the cation concentration. The dynamics of cations is where is time, and are the flux rates through the fusion pore and vesicle, respectively, is the (dimensionless) area fraction of the cation fluxes, and and are the radii of the spherical vesicle and the fusion pore, respectively. Note that and are defined at the start of the simulation, while and are the two variables in the system. For the sake of simplicity, it is assumed that the cation concentration in the vesicle is equal to the cation concentration in the fusion pore.

3. Results

3.1. Equilibrium Shape Analysis

We analyzed the effects of varying the relative cation concentration () on the equilibrium radius (). The parameter values used for the analysis were as follows:  μm⁻¹,  g s⁻²,  nm, and . In Figure 4, we plot the relationship between and . In this relationship, the equilibrium radius was abruptly increased when approached 0.9. The increase in the equilibrium radius was due to the reduction in the spontaneous curvature of the anisotropic membrane domain at the fusion pore. On the other hand, the reduction in the cation concentration from 0.9 to 0.1 led to the abrupt decrease in (i.e., the abrupt closure of the fusion pore).

Figure 4: The relationship between the equilibrium radius of the pore () and the relative cation concentration in the vesicle (). The relationship between and is derived from the minimization of the membrane bending energy. Note that an increase in the cation concentration (bold line) leads to an abrupt increase in the fusion pore radius, (i.e., the opening of the fusion pore), while, in the opposite direction (dotted line), a reduction in the cation concentration leads to an abrupt closure of the fusion pore.

3.2. Numerical Simulations: Dynamic Model

In our dynamic model, the membrane force exerted along the perpendicular direction to the major axis of the fusion pore membrane was derived from the anisotropic membrane bending energy. To have a better understanding of how the cation fluxes through the vesicle membrane and fusion pore affected the opening and closure of the fusion pore, numerical simulations were employed to solve this nonlinear system of fusion pore dynamics. It was assumed that the forces exerted on the fusion pore membrane were the same in the open and closed states but in the opposite direction (Figure 4). In particular, the equilibrium pore radius diagram as a function of the cation concentration in the closed state of the pore was reversed to the one in the open state. Figure 5 shows the functional dependence between changes in the fusion pore radius (Figure 5(b)) and cation concentration in the vesicle (Figure 5(c)). To obtain the same time interval for the open and closed states of the fusion pore, the cation flux rate through the fusion pore was set to be twice as large as the flux rate through the vesicle membrane. Note that the cation flux through the fusion pore effectively occurred only after its abrupt opening, while the cation flux through the vesicle membrane was active throughout the simulation. In the pore opened state, the cation concentration in the vesicle was increased from 0.1 to 0.9 over a time interval of 0.3 sec. In the pore closed state, the concentration of the cations in the vesicle was reduced from 0.9 to 0.1 over the same time interval.

(a)

(b)

(c)

Figure 5: The effects of changes in cation concentration in the vesicle on the rhythmic opening and closure of the fusion pore. Cell-attached patch-clamp recordings (capacitance trace) following the application of a cAMP-increasing agent demonstrate fast rhythmic opening and closure of an individual fusion pore (a). Note that the capacitance change between closed and open fusion pore states is about 5 fF. Numerical simulations demonstrating the time dynamics of the pore radius () (b). Note the abrupt opening and closure of the fusion pore. The membrane force derived from the membrane bending energy is exerted along the perpendicular direction to the major axis of the fusion pore (see arrows). The cation concentration dynamics includes the efflux through the fusion pore and the influx through the vesicle membrane (c).

4. Discussion

In this study, we have constructed a model of fusion pore dynamics to explain the enhanced periodicity in the opening and closure of individual fusion pores observed following the application of a cAMP-increasing agent (Figure 2). The assumed axisymmetric tubular fusion pore was connected to a vesicle of spherical geometry (Figure 3). The membrane of the fusion pore was composed of a single anisotropic membrane domain. The bending energy of the fusion pore membrane took into account the curvature mismatch between the spontaneous curvature of the anisotropic membrane pore domain and the curvature of the pore tubular membrane (Figure 3). The spontaneous curvature of the anisotropic membrane pore was considered to depend on the local concentration of cations. The differentiation of the free energy with respect to the membrane pore deformation gave a membrane force in the radial direction of the axisymmetric fusion pore. The influx of cations was through the membrane of the vesicle, while the cation efflux was through the fusion pore. In equilibrium state, the force exerted on the fusion pore membrane was zero, and the fusion pore diameter was mainly a function of the coupling between the intrinsic (spontaneous) curvature of the anisotropic membrane pore domain and the cation concentration (Figure 4).

Pituitary cells express various cation channels, for instance HCN channels which are modulated by cAMP [30]. These channels are permeable to Na⁺ and K⁺ and are therefore likely candidates for modulating local membrane anisotropy. Following the application of a cAMP-increasing agent, the patch-clamp membrane capacitance data from hormone secreting pituitary cells revealed faster rhythmic openings and closures of the fusion pore (Figure 2). To reproduce the sharp changes in the fusion pore diameter, it was necessary that the spontaneous curvature of the anisotropic membrane pore domain was sufficiently large. The influx of cations through the vesicle membrane (i.e., through ion channels) increased the cation concentration in the vesicle, which reduced the spontaneous curvature of the anisotropic membrane pore domain leading to the fast opening of the fusion pore. On the other hand, the cation efflux through the fusion pore reduced the cation concentration in the vesicle, which increased the spontaneous curvature of the anisotropic membrane pore domain causing the rapid closure of the fusion pore (Figure 5). To conclude, the complete cycle of fusion pore opening and closure could be driven by cation fluxes through the vesicle membrane and the fusion pore.

The fusion pore is a membranous neck of an attached vesicle with the plasma membrane (Figure 1). The stability of the fusion pore is evident from the periodic opening and closure of the same fusion pore [4–6]. The alternative scenario is the new formation of the fusion pore before each cycle of opening and closure, but this is of low probability due to the accompanied high bending energy costs [9]. By the cell-attached patch-clamp recordings in resting lactotrophs, the conductance of the fusion pore was estimated as low as a few pS, accounting for an equilibrium radius below 1 nm [5, 8]. In the present study, we revealed that the coupling of the pore membrane domain spontaneous curvature and the cation concentration may cause fast transitions between the open and closed fusion pore states.

To date, the composition of the proposed anisotropic pore membrane domain remains unclear. In general, the specific shape of fusion pore membrane can be stabilized by either lipids, proteins, or their complexes. Among the many possibilities, the SNARE protein complex is one of the major candidates responsible for the stability of the fusion pore. The SNARE protein complex is reported to direct the attachment of the vesicle to the plasma membrane, thereby overcoming the large energy barrier [9]. The SNARE protein complex at the plasma membrane is bound to synaptobrevin on the membrane of the vesicle [10]. This attachment has been shown to precede the formation of the fusion pore. The availability of synaptobrevin is regulated by sphingosine, a releasable backbone of sphingolipid [31]. Possibly, the electrostatic repulsion between sphingosines and positively charged residues of synaptobrevin could account for a change in synaptobrevin’s structure enabling its contact with the SNARE protein complex [31]. However, the SNARE complex may change the fusion pore properties at the postfusion stage [8]. We here propose that the accumulation of cations at the entrance of the fusion pore (Figure 3(a)) may interact electrostatically with the positively charged sphingosines leading to the local dispersion (or demixing) of sphingosines in the anisotropic pore domain. The large anisotropic spontaneous curvature of the fusion pore could be stabilized also by negatively charged lipids such as PIP₂ (Phosphatidylinositol 4,5-bisphosphate) lipids [12]. The aggregation of cations in the vesicle could attract PIP₂ lipids from the pore region, thereby changing the local anisotropicity of the fusion pore.

Another possible mechanism for the stability and reversibility of the fusion pore is the dynamic reorganization of a cholesterol-sphingomyelin raft at the fusion pore. The raft is stabilized by hydrogen bonds and hydrophobic interactions between cholesterol and sphingomyelin lipids [32]. The anisotropic spontaneous curvature of the lipid raft might be, according to the molecular interactions, between the raft constituents. The narrow (closed) fusion pore could be stabilized by the gel-like properties such as increased rigidity of the cholesterol-sphingomyelin nanodomain raft. The increase in cation concentration may modulate the molecular interactions within the raft, thereby reorganizing the raft into a more disordered (liquid) state. The rapid and reversible phase transition between gel-like and disordered states may account for the rhythmic closing and opening of the fusion pore, respectively.

The proposed theory of fusion pore opening and closure dynamics awaits confirmation in future experiments. The usage of ion channel blockers against cAMP-channels may also modulate the rhythmic opening and closure of the fusion pore. The relationship between the fusion pore opening time interval (i.e., dwell time) and the number of cation channels can be investigated by the partial blocking of the ion channels. The origin of the anisotropic mechanical properties of the pore membrane can be identified by mutagenic assays and overexpression of specific lipids or proteins. Finally, we hope that this computational model would shed light on the intricate biophysical mechanisms underlying the rhythmic opening and closure of the fusion pore.

Acknowledgments

This work was supported by the Slovenian Research Agency, Grant nos: P3 310, J3-3632, J3-4146, and J3 4051, and the CipKeBip (to R. Zorec).

References

J. E. Heuser and T. S. Reese, “Evidence for recycling of synaptic vesicle membrane during transmitter release at the frog neuromuscular junction,” Journal of Cell Biology, vol. 57, no. 2, pp. 315–344, 1973.

G. A. de Toledo, R. Fernandez-Chacon, and J. M. Fernandez, “Release of secretory products during transient vesicle fusion,” Nature, vol. 363, no. 6429, pp. 554–558, 1993.

B. Ceccarelli, W. P. Hurlbut, and A. Mauro, “Depletion of vesicles from frog neuromuscular junctions by prolonged tetanic stimulation,” Journal of Cell Biology, vol. 54, no. 1, pp. 30–38, 1972.

M. Stenovec, M. Kreft, I. Poberaj, W. J. Betz, and R. Zorec, “Slow spontaneous secretion from single large dense-core vesicles monitored in neuroendocrine cells,” The FASEB Journal, vol. 18, no. 11, pp. 1270–1272, 2004.

N. Vardjan, M. Stenovec, J. Jorgačevski, M. Kreft, and R. Zorec, “Subnanometer fusion pores in spontaneous exocytosis of peptidergic vesicles,” Journal of Neuroscience, vol. 27, no. 17, pp. 4737–4746, 2007.

J. Jorgačevski, M. Potokar, S. Grilc et al., “Munc18-1 tuning of vesicle merger and fusion pore properties,” Journal of Neuroscience, vol. 31, no. 24, pp. 9055–9066, 2011.

Š. Perutková, V. Kralj-Iglič, M. Frank, and A. Iglič, “Mechanical stability of membrane nanotubular protrusions influenced by attachment of flexible rod-like proteins,” Journal of Biomechanics, vol. 43, no. 8, pp. 1612–1617, 2010.

J. Jorgačevski, M. Fošnarič, N. Vardjan et al., “Fusion pore stability of peptidergic vesicles,” Molecular Membrane Biology, vol. 27, no. 2-3, pp. 65–80, 2010.

M. M. Kozlov and V. S. Markin, “Possible mechanism of membrane fusion,” Biofizika, vol. 28, no. 2, pp. 242–247, 1983.

R. Jahn and R. H. Scheller, “SNAREs—engines for membrane fusion,” Nature Reviews Molecular Cell Biology, vol. 7, no. 9, pp. 631–643, 2006.

J. G. Duman and J. G. Forte, “What is the role of SNARE proteins in membrane fusion?” American Journal of Physiology, vol. 285, no. 2, pp. C237–C249, 2003.

B. Rituper, B. Davletov, and R. Zorec, “Lipid-protein interactions in exocytotic release of hormones and neurotransmitters,” Clinical Lipidology, vol. 5, no. 5, pp. 747–761, 2010.

M. B. Jackson and E. R. Chapman, “The fusion pores of Ca(2+)-triggered exocytosis,” Nature Structural and Molecular Biology, vol. 15, no. 7, pp. 684–689, 2008.

A. Iglič, H. Hägerstrand, P. Veranič, A. Plemenitaš, and V. Kralj-Iglič, “Curvature-induced accumulation of anisotropic membrane components and raft formation in cylindrical membrane protrusions,” Journal of Theoretical Biology, vol. 240, no. 3, pp. 368–373, 2006.

A. Tian and T. Baumgart, “Sorting of lipids and proteins in membrane curvature gradients,” Biophysical Journal, vol. 96, no. 7, pp. 2676–2688, 2009.

B. Sorre, A. Callan-Jones, J. B. Manneville et al., “Curvature-driven lipid sorting needs proximity to a demixing point and is aided by proteins,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 14, pp. 5622–5626, 2009.

D. Kabaso, E. Gongadze, P. Elter et al., “Attachment of rod-like (BAR) proteins and membrane shape,” Mini-Reviews in Medicinal Chemistry, vol. 11, no. 4, pp. 272–282, 2011.

J. B. Fournier and P. Galatola, “Bilayer membranes with 2D-nematic order of the surfactant polar heads,” Brazilian Journal of Physics, vol. 28, no. 4, pp. 329–338, 1998.

V. Kralj-Iglič, A. Iglič, H. Hägerstrand, and P. Peterlin, “Stable tubular microexovesicles of the erythrocyte membrane induced by dimeric amphiphiles,” Physical Review E, vol. 61, no. 4, pp. 4230–4234, 2000.

M. Fošnarič, V. Kralj-Iglič, K. Bohinc, A. Iglič, and S. May, “Stabilization of pores in lipid bilayers by anisotropic inclusions,” Journal of Physical Chemistry B, vol. 107, no. 45, pp. 12519–12526, 2003.

V. Kralj-Iglič, V. Heinrich, S. Svetina, and B. Žekš, “Free energy of closed membrane with anisotropic inclusions,” European Physical Journal B, vol. 10, no. 1, pp. 5–8, 1999.

A. Iglič, B. Babnik, K. Bohinc, M. Fošnarič, H. Hägerstrand, and V. Kralj-Iglič, “On the role of anisotropy of membrane constituents in formation of a membrane neck during budding of a multicomponent membrane,” Journal of Biomechanics, vol. 40, no. 3, pp. 579–585, 2007.

M. Rappolt, A. Hodzic, B. Sartori, M. Ollivon, and P. Laggner, “Conformational and hydrational properties during the Lβ- to Lα- and Lα- to HII-phase transition in phosphatidylethanolamine,” Chemistry and Physics of Lipids, vol. 154, no. 1, pp. 46–55, 2008.

A. Perutková, M. Daniel, M. Rappolt et al., “Elastic deformations in hexagonal phases studied by small-angle X-Ray diffraction and simulations,” Physical Chemistry Chemical Physics, vol. 13, no. 8, pp. 3100–3107, 2011.

A. Iglič, B. Babnik, U. Gimsa, and V. Kralj-Iglič, “On the role of membrane anisotropy in the beading transition of undulated tubular membrane structures,” Journal of Physics A, vol. 38, no. 40, pp. 8527–8536, 2005.

A. Iglič, M. Lokar, B. Babnik et al., “Possible role of flexible red blood cell membrane nanodomains in the growth and stability of membrane nanotubes,” Blood Cells, Molecules, and Diseases, vol. 39, no. 1, pp. 14–23, 2007.

D. Kabaso, E. Gongadze, J. Jorgačevski et al., “Exploring the binding dynamics of BAR proteins,” Cellular and Molecular Biology Letters, vol. 16, no. 3, pp. 398–411, 2011.

A. Veksler and N. S. Gov, “Phase transitions of the coupled membrane-cytoskeleton modify cellular shape,” Biophysical Journal, vol. 93, no. 11, pp. 3798–3810, 2007.

D. Kabaso, R. Shlomovitz, K. Schloen, T. Stradal, and N. S. Gov, “Theoretical model for cellular shapes driven by protrusive and adhesive forces,” PLoS Computational Biology, vol. 7, no. 5, Article ID e1001127, 2011.

K. Kretschmannova, M. Kucka, A. E. Gonzalez-Iglesias, and S. S. Stojilkovic, “The expression and role of hyperpolarization-activated and cyclic nucleotide-gated channels in endocrine anterior pituitary cells,” Molecular Endocrinology. In press

F. Darios, C. Wasser, A. Shakirzyanova et al., “Sphingosine facilitates SNARE complex assembly and activates synaptic vesicle exocytosis,” Neuron, vol. 62, no. 5, pp. 683–694, 2009.

D. Lingwood and K. Simons, “Lipid rafts as a membrane-organizing principle,” Science, vol. 327, no. 5961, pp. 46–50, 2010.