674 Acta Chim. Slov. 2020, 67, 674–681 Penič et al.: Active Forces of Myosin Motors May Control ... DOI: 10.17344/acsi.2020.5863 Scientific paper Active Forces of Myosin Motors May Control Endovesiculation of Red Blood Cells Samo Penič,1,* Miha Fošnarič,2 Luka Mesarec,1 Aleš Iglič1,3,4 and Veronika Kralj-Iglič2,4 1 Faculty of Electrical Engineering, University of Ljubljana, Tržaška cesta 25, SI-1000 Ljubljana 2 Faculty of Health Sciences, University of Ljubljana, Zdravstvena pot 5, SI-1000 Ljubljana 3 Faculty of Medicine, University of Ljubljana, Zaloška 9, SI-1000 Ljubljana, Slovenia 4 Institute of Biosciences and Bioresources, National Research Council, Pietro Castellino 111, 80131 Napoli, Italy * Corresponding author: E-mail: samo.penic@fe.uni-lj.si, Tel: +386 14768 822 Received: 01-27-2020 Abstract By using Monte Carlo (MC) simulations, we have shown that the active forces generated by (NMIIA) motor domains bound to F-actin may partially control the endovesiculation of the red blood cell (RBC) membrane. The myosin gener- ated active forces favor pancake-like (torocyte) RBC endovesicles with a large flat central membrane region and a bulby periphery. We suggest that the myosin generated active forces acting on the RBC membrane in the direction perpendic- ular to the membrane surface towards the interior of the RBC may influence also other RBC shape transformations and the stability of different types of RBC shapes and should be therefore considered in the future theoretical studies of the RBC vesiculation and shape transformations. Keywords: Myosin generated active forces; Monte Carlo simulations; endovesiclulation; intrinsic curvature; red blood cell 1. Introduction A biological membrane forms a physical boundary between the inner volume of a biological cell and the ex- ternal medium, as well as, within the cell, between the lu- mens of intracellular organelles and cytosol. The main building block of the biological membrane lipid bilayer consists of two monolayers of phospholipid molecules.1 In the lipid bilayer, other constituents like membrane pro- teins are intercalated.2,3 The transmembrane proteins are the pining points for membrane skeleton or/and cytoskel- eton.4 Furthermore, some of the inclusions in the lipid bi- layer, like proteins or nanoparticles, can impose local cur- vature of the membrane,5–9 forcing the membrane to locally and/or globally adapt its shape.10–17 The phase sep- aration of membrane inclusions (nanodomains) is an im- portant mechanism that may induce the local changes of membrane curvature and is therefore the driving force for transformations of the cell shape.17–19 The spontaneous phase separation of membrane proteins, driven by the forces of actin polymerization and cell-substrate adhesion, was in a quantitative manner predicted for the first time by Veksler and Gov.20 Besides the membrane inclusions’ driven membrane shape changes, other mechanisms also determine the membrane cell shape. Among them, constant energy con- suming forces are acting in the living cell.17 In experi- ments, the cells or lipid bilayer vesicles (as model systems) may also be deformed by extracellular forces due to pres- sure differences,21 stretching22 and fluid flow.23–25 An ex- ternal force to the membrane surface can also be experi- mentally generated, for example by a cantilever of the atomic force microscope.26 The cell membrane resistance to the shape changes depends on its mechanical properties27 mainly governed by the composition of the membrane. The dynamic re- sponse of the cell membrane shape to the force17 is impor- tant for different cell functions that range from adhesion 675Acta Chim. Slov. 2020, 67, 674–681 Penič et al.: Active Forces of Myosin Motors May Control ... and migration to division, differentiation, and cell death.28 Therefore, the membrane mechanical properties are of great scientific interest. Consequently, new theoretical ap- proaches for deeper understanding of the mechanisms that rule the cell’s functions are constantly being devel- oped. In red blood cells (RBC), besides the lipid bilayer also the membrane skeleton (spectrin-F-actin network at- tached to the inner lipid layer)29 plays an important role in cell shape determination and cell transformations.30–32 It has been shown that equilibrium RBC shapes that corre- spond to the minimal membrane elastic energy, should take into account the local and non-local bending ener- gy27,33,34 and the elastic energy of the membrane skele- ton.30–32,35 The shear deformation of the membrane skele- ton plays an essential role in the stability of echinocytic RBC shapes.30,31 Additionally, cytoskeleton induced pro- trusions in lipid membranes has been theoretically ex- plained by the interplay between the elasticity of the mem- brane lipid bilayer and the membrane skeleton – firstly for axisymmetric shapes36,37 and later also for non-axisym- metric shapes.38 Non-local bending energy27,33,39,40 depends on the change in the relaxed areas of the two lipid layers of the membrane bilayer.41 Decreasing the difference between the relaxed areas of the outer layer and the inner lipid lay- er, favors the concave local membrane shape, i.e. the in- ward bending of the membrane.32,35,41 On the other hand, increasing the difference between the relaxed areas of the outer and inner membrane layers promotes the formation of convex local shape, i.e. the outward bending.31,35,42 It was shown that exogenously added amphiphiles predominantly bound in the outer lipid layer induce the transformation of the discoid RBC into the spiculated echinocytic RBC, while amphiphiles predominantly bound in the inner lipid layer induce the transformation into invaginated stomatocytic shapes.43 When RBCs are incubated with high sublytic concentrations of amphi- philes, the microexovesiculation starts.42,44,45 The amphi- phile induced RBC microexovesicles are highly depleted in the membrane skeleton37 suggesting that a local disrup- tion of the interactions between the membrane skeleton and the membrane bilayer occurred prior to micro- or na- no-exovesiculation.45 In RBCs, most types of amphiphiles induce spherical microexovesicles that are formed from sphere-like mem- brane buds and are free of the membrane skeleton. On the other hand, there are also anisotropic amphiphiles which induce the growing of tubular membrane buds and the re- lease of stable tubular microexovesicles.44,45 The tubular budding and vesiculation of the RBC membrane can be theoretically explained by deviatoric membrane properties due to the in-plane orientational ordering of the aniso- tropic membrane inclusions induced by intercalated ani- sotropic amphiphiles.44,45 The deviatoric properties of the membrane12,15,42,45–52 may explain the experimentally ob- served tubular membrane protrusion without the applica- tion of the local force.16,42,45,51,53,54 For isotropic mem- brane, the application of the local force is necessary to theoretically predict the tubular membrane protru- sions.16,54,55 It was reported44,56 that some amphiphilic molecules can induce large membrane invaginations (i.e. stomatocyte RBC shape) as shown in Figure 1C and endovesicles in RBCs. By means of transmission electron microscopy (TEM) and confocal laser scanning microscopy using fluo- rescent markers, it was also observed that many stomatocy- togenic amphiphiles (for example chlorpromazine hydro- chloride) can induce in RBCs with large concentrations of amphiphiles small spherical endovesicles (Figure 1D).44,56,57 On the other hand, the stomatocytogenic detergent poly- ethyleneglycol dodecylether (C12E8) induces large endove- sicles with a unique ring-like toroidal shape joined with a central flat membrane segment, i.e. flat membrane struc- tures with a bulby periphery called also torocyte endovesi- cles (Figure 1 A,B and E)57 with the shape very similar to the shape of Golgi bodies.58 It was proposed that the ob- served RBC torocyte endovesicles were formed in a process in which an initially stomatocytic RBC invagination loses volume while maintaining a large surface area.57,59 The results of theoretical modeling indicated that torocyte RBC endovesicles can be mechanically stabilized by non-homogeneous lateral distribution of laterally mo- bile anisotropic membrane inclusions, like for example by anisotropic detergent-membrane component complex- es.58,59 It was further shown in 2019 that the mechanical stability of torocyte (pancake-like) closed membrane vesi- cle shapes can be explained also by the coupling of the curved isotropic membrane inclusions and active (cy- toskeletal) forces.55 Until recently,60 it was believed that the active forces are absent in the mechanisms of the determi- nation of the RBC shape and vesiculation, as discussed above. It has been shown in 2018 that nonmuscle myosin IIA (NMIIA) motor domains may generate tension in spectrin-F-actin also in a 2-dimensional RBC membrane skeleton and in this way partially control the RBC shape.60 The role of NMIIA contractility in generating tension in the RBC network and partially controlling the RBC shape was thus in the past completely neglected, until recently.60 The length of NMIIA filament is around 200 nm.60 The influence of NMIIA bipolar filaments, associated with a 2-dimensional RBC membrane skeleton can be thus described by the local force acting on the RBC membrane in the direction to the interior of the RBC. In the present paper, we shall study the interplay between the laterally mobile membrane inclusions and the inward-oriented lo- cal forces of the NMIIA bipolar filaments in the formation of torocyte endovesicles (flat membrane structures with a bulby periphery). The present study was motivated by the results presented in references 44, 56, 57, 59, 60, but it is not entirely/directly connected to the experimental and theoretical results presented in these articles. For the sake 676 Acta Chim. Slov. 2020, 67, 674–681 Penič et al.: Active Forces of Myosin Motors May Control ... of simplicity, we shall consider the closed lipid bilayer membrane with one type of inclusions only. In our model, the single inclusion can induce local membrane bending due to its negative intrinsic curvature and also because of inward-oriented active forces. This means that we have joined the effect of the intrinsic curvature of the mem- brane inclusions (nanodomains) and the effect of the local active forces of NMIIA bipolar filament domains in a sin- gle type of membrane inclusions (nanodomains). 2. Monte Carlo Simulations Monte Carlo (MC) triangulated mesh was used to numerically model and investigate the vesicle shape and the lateral distribution of membrane inclusions by means of computer simulation.55 Phospholipid bilayer mem- branes can be treated due to their small thickness in the first approximation as a two dimensional surface, allowing the continuum approach in the theoretical description of membrane surfaces.61 In the model we discretize the mem- brane into patches consisting of many molecules (Figure 2). A single patch is represented by a vertex in a triangulat- ed surface model. The main model parameter that defines mechanical bilayer properties is bending stiffness. The vesicle is represented by a set of N vertices that are linked by bonds of flexible length d to form a closed, randomly triangulated, self-avoiding network.62,63 The lengths of the tethers can vary between a minimal (dmin) and a maximal (dmax) value. The self-avoidance of the net- work can be implemented by ensuring that no vertex can penetrate through the triangular network. The maximal Figure 1: TEM and SEM micrographs showing stomatocytic human RBC with plate-like torocytic invaginations/endovesicles with a bulby periphery (panels A and B) formed in a process in which initially stomatocytic RBC invaginations (panel C) loose volume while maintaining a large surface area. RBCs were incubated by amphiphilic molecules octaethyleneglycol dodecylether (C12E8). Panels D and E show confocal laser scanning micros- copy of RBCs incubated with amphiphilic molecules chlorpromazine inducing small spherical endovesicles (panel D) and with C12E8 inducing toroidal endovesicles (panel E) as presented also in panels A and B. Adapted from.57 677Acta Chim. Slov. 2020, 67, 674–681 Penič et al.: Active Forces of Myosin Motors May Control ... possible random displacement of the vertex in a single step (s) should be small enough so that the fourth vertex can- not move through the plane of the other three to the min- imal allowed distance, dmin, from the three vertices. In our scenario, we use s=0.15 dmin and dmax=1.7 dmin. For details about the expressions to calculate self-avoidance con- straint dmax, see.64 The initial state of triangulated surface is a pentago- nal bipyramid with all the edges divided into equilateral bonds so that the network is composed of 3(N-2) bonds forming 2(N-2) triangles. Nc randomly selected vertices are given non-zero isotropic intrinsic curvature of c0, thus they become the model of the membrane inclusion. The rest of the vertices have zero intrinsic curvature. Positive curvature means that the membrane will locally bulge to- wards the exterior, negative curvature will force the mem- brane to bulge towards the interior compartment of the vesicle. The system is developed into the thermal equilibri- um state. The evolution of the system is measured in Mon- te Carlo sweeps (mcs). One mcs consists of individual at- tempts to displace each of the N vertices by a random increment in the sphere with radius s – the action we will refer to as a vertex move. Membrane fluidity is maintained by flipping bonds within the triangulated network. In each mcs, the vertex move attempts are followed by a 3N at- tempts to flip a randomly chosen bond. A single bond flip involves the four vertices of two neighboring triangles. The tether between the two vertices is cut and re-established between the other two, previously unconnected vertices (a detailed description was published elsewhere).65 Each in- dividual Monte Carlo step (vertex move or bond flip) is accepted with probability according to Metropolis Hast- ings algorithm, based on free energy change due to the Monte Carlo step. Figure 2: Schematic figure of the triangulated membrane surface in our MC model of the RBC membrane. The links between vertices are bonds, forming triangulated mesh. Lengths of all bonds must satisfy the bond length constraint so that the bond lengths are al- ways in the range between dmin and dmax. The grayed out vertex is displaced by step size s to a new position. After the move, the lengths of bonds d1, d2, d3 and d4 are calculated and it is verified that they are between dmin and dmax. The energy W in simulation consists of three parts: (1) where Wb is the bending energy of the membrane, Wd is the energy of the direct interaction between vertices with intrinsic curvature and WF is the energy due to myosin forces acting on the membrane. For the bending energy Wb of the membrane, we use the standard Helfrich expression for a tensionless mem- brane with a term that represents intrinsic curvature.66 The membrane keeps fixed topology, thus the contribution of the Gaussian curvature to the change of bending energy is cancelled out. (2) where κ is the bending stiffness of the membrane, c1, c2 and c0 are the two principal curvatures and the intrinsic curva- ture of the vesicle membrane at the point under consider- ation. Note that only points where inclusions are located have non-zero intrinsic curvature c0, whereas all other points on the membrane have intrinsic curvature c0 set to zero. The integration is performed over the membrane area A. For modelling attraction force between the vertices with intrinsic curvature energy term:55 (3) where w is a direct interaction constant. The energy is summed over all inclusion pairs with their in-plane dis- tance rij , r0 is the range of direct interaction and H is a Heaviside step function with (r0-rij) being function’s argu- ment. The energy contribution of the local protrusive forc- es due to the myosin motor domains/inclusions:55 (4) where F is the magnitude of the force, the sum runs over all proteins, is the outwards facing normal to the mem- brane at the location of the vertex with the inclusion i and xi is the position vector of the vertex with the inclusion i. 3. Results and Discussion Figure 3 shows the MC simulations of RBC shape transformations induced by laterally mobile membrane in- clusions (nanodomains) with negative intrinsic curvature. It can be seen in Figure 3 that the non-homogeneous later- al redistribution of the mobile inclusions with negative in- trinsic shape causes the RBC membrane to locally change the curvature resulting in the global transformation of the RBC shape. The MC predicted RBC shapes presented in 678 Acta Chim. Slov. 2020, 67, 674–681 Penič et al.: Active Forces of Myosin Motors May Control ... Figure 3 depend on the inclusion concentration, the inclu- sion intrinsic curvature, the strength of the direct attrac- tive interaction between inclusions and on the active forc- es exerted by the inclusions. Figure 3: Monte Carlo simulation of the RBC membrane transfor- mation induced by mobile membrane inclusions with intrinsic cur- vature c0 = –1 lmin–1. Two different membrane inclusion concentra- tions p are considered. The upper figures A and B shows the MC results in the absence of myosin generated actin acting forces (F = 0), while in the lower panels C and D the myosin generated actin active forces F = –1 kT/lmin are taken into account. The negative sign of F denotes that the myosin generated forces points into the vesicle interior (perpendicular to the local membrane surface). The trian- gulated membrane surface is drawn semitransparent to uncover its interior shape. Red vertices on the mesh represents the locations of the membrane inclusions. In the panels C and D black and white outlines are added to better visualize the shape of endovesicles and membrane necks. The values of other model parameters are: bend- ing rigidity κ = 25 kT and direct interaction parameter w = 1.25 kT. Figures 3A and B show that the lateral accumulation of membrane inclusions (nanodomains) can induce the formation of long undulated thin inward membrane pro- trusions (buds). Long undulated membrane protrusions may be further transformed into small independent spher- ical endovesicles, as observed in Figure 1D, due to the frustrations in the orientational ordering of membrane components in the highly curved membrane necks.67 The same mechanism can also be responsible for the detach- ment of the complete inward membrane protrusion from the parent membrane67 and the consequent formation of the endovesicles as shown in experimental Figure 1. In calculations presented in Figures 3C and D, it is taken into account that the membrane inclusions (nano- domains) exert also active forces in the direction perpen- dicular to the membrane surface towards the interior of the RBC. As already discussed above, the active forces in the RBC membrane are generated by myosin (NMIIA) motor domains (inclusions) bound to F-actin of the RBC membrane skeleton.60 It can be seen in panels of Figure 3 that at a smaller concentration of membrane inclusions exerting force on the membrane, the MC predicted RBC shape has one large invagination (Fig. 3C) as can be ob- served in some experiments (see refs. 44,56 and the refer- ences therein). Large invaginations can be separated from the parent cell due to the frustrations in the orientational ordering of membrane components in the highly curved membrane necks connecting the invagination and the parent cell,67 resulting in the formation of a large endove- sicle. Furthermore, Figure 3 shows that a larger concentra- tion of membrane inclusions exerting force on the mem- brane, the MC predicted RBC shape has one small and two large pancake-like torocyte membrane invaginations (Fig- ure 3D) as can also be observed in the experiments (Figure 1, panels A, B and E). Again, the necks connecting the torocyte structures as well as the neck connecting the com- plete invagination to the parents cell are supposed to be ruptured due to the frustrations in the orientational order- ing of membrane components in the highly curved mem- brane necks.67 Note that in the torocytic membrane invag- inations (Figure 3D), the myosin motor domains generated active forces that acted in the outward direction with re- spect to the torocyte invagination, i.e. in the direction to- wards the inner RBC solution. It is also important to point out that the myosin motor domains/inclusions are mainly accumulated at the bulby rim of the torocyte membrane invaginations (Figure 3D). The stability of torocyte RBC endovesicles can be theoretically explained also by anisotropic membrane in- clusions which exhibit orientational ordering in a highly curved bulby periphery.57,59 In this work we have shown that the stability of torocyte endovesicles may be addition- ally favored also by active forces on the RBC membrane, generated by myosin motor nanodomains (i.e. membrane inclusions/nanodomains composed of myosin-actin-spec- trin-lipids complex, see also references 7, 19, 68, 69). Note that in the present work the total number of myosin motor nanodomains (inclusions) in Figures 3C and D is larger than the actual number of myosin motor nanodomains found in the RBC membrane.60 Therefore, the proposed mechanisms of myosin inclusions/nanodo- mains driven formation and stabilization of torocyte in- vaginations and torocyte endovesicles in RBCs can be con- sidered as an additional and complementary mechanism to the non-homogeneous lateral distribution and the ori- entational ordering of anisotropic membrane inclusions/ nanodomains in the RBCs membrane.57,59 In accordance with experimental observations (Fig- ure 1) we have predicted in this work the invaginated stomatocyte RBC shapes having different shapes of invag- inations, like torocytic, spherical, undulated necklace-like, etc. This is an extension of the previously theoretically pre- dicted shape classes of the invaginated stomatocytic 679Acta Chim. Slov. 2020, 67, 674–681 Penič et al.: Active Forces of Myosin Motors May Control ... shapes which were mostly limited to the simple (axisym- metric) stomatocytic shape with only one large invagina- tion (Figure 3C) (see for example ref. 32), experimentally observed also in a giant unilamellar lipid vesicle.70 4. Conclusions Numerical computer modelling of the cell mem- brane shapes and shape transformations is a widely used method to investigate physical models of various phenom- ena found experimentally in cellular systems. In the pres- ent work, we used MC simulations to theoretically study the influence of active forces on the red blood cell (RBC) shape and vesiculation. We have shown that the active forces, generated by myosin motor domains, may partially control endovesiculation of the RBC membrane and also the RBC shape changes in general. Among others the my- osin domains generated forces on the RBC membrane fa- vor experimentally observed pancake-like (torocyte) RBC endovesicles with a large flat central membrane and a bul- by periphery which were also experimentally observed.57 These theoretically predicted shapes (Figure 3D) are very similar to the shapes of Golgi bodies, so the theoretical study presented in this work may be relevant also for better understanding of the physical mechanisms determining the shape of Golgi body.58 In conclusion, by using MC simulations, it is shown in this work that the recently discovered myosin motor do- mains generated active forces on the RBC membrane60 may partially control the endovesiculation of the red blood cells and the RBC shape changes in general. We conclude that the myosin generated active forces acting on the RBC membrane should be therefore necessarily considered in the future relevant theoretical studies of the RBC vesicula- tion and shape transformations. Acknowledgments This project has received funding from the European Union’s Horizon 2020 research and innovation program under grant agreement No. 801338 (VES4US project). The authors also acknowledge the financial support from the grants No. P2-0232, P3-0388 and J2-8166 from the Slove- nian Research Agency (ARRS). Conflict of Interests The authors declare that there is no conflict of inter- ests. 5. References 1. G. Cevc and D. Marsh, Phospholipid Bilayers: Physical Princi- ples and Models, Wiley, 1987. 2. J. N. Israelachvili, Intermolecular and Surface Forces, Academ- ic Press, 1997. 3. S. A. Safran, Statistical Thermodynamics of Surfaces, Interfaces, and Membranes, Addison-Wesley Publishing Company, Col- orado, USA, 1994. 4. Z. Peng, R. J. Asaro and Q. Zhu, Phys. Rev. E, 2010, 81, 031904. DOI:10.1103/PhysRevE.81.031904 5. I. Szleifer, D. Kramer, A. Ben-Shaul, W. M. Gelbart and S. A. Safran, J. Chem. Phys., 1990, 92, 6800–6817. DOI:10.1063/1.458267 6. C. Nielsen, M. Goulian and O. S. Andersen, Biophys. J., 1998, 74, 1966–1983. DOI:10.1016/S0006-3495(98)77904-4 7. M. Fošnaric, K. Bohinc, D. R. Gauger, A. Iglic, V. Kralj-Iglic and S. May, J. Chem. Inf. Model., 2005, 45, 1652–1661. DOI:10.1021/ci050171t 8. S. May, Langmuir, 2002, 18, 6356–6364. DOI:10.1021/la025747c 9. M. Daniel, J. Řezníčková, M. Handl, A. Iglič and V. Kralj-Iglič, Sci. Rep., 2018, 8, 10810. DOI:10.1038/s41598-018-28965-y 10. V. Markin, Biophys. J., 1981, 36, 1–19. DOI:10.1016/S0006-3495(81)84713-3 11. S. Leibler, J. Phys., 1986, 47, 507–516. DOI:10.1051/jphys:01986004703050700 12. V. Kralj-Iglič, S. Svetina and B. Žekž, Eur. Biophys. J., 1996, 24, 311–321. DOI:10.1007/BF00180372 13. B. Božič, V. Kralj-Iglič and S. Svetina, Phys. Rev. E, 2006, 73, 041915. DOI:10.1103/PhysRevE.73.041915 14. A. Iglič, B. Babnik, K. Bohinc, M. Fošnarič, H. Hägerstrand and V. Kralj-Iglič, J. Biomech., 2007, 40, 579–585. DOI:10.1016/j.jbiomech.2006.02.006 15. N. Walani, J. Torres and A. Agrawal, Proc. Natl. Acad. Sci. USA, 2015, 112, E1423–E1432. DOI:10.1073/pnas.1418491112 16. L. Mesarec, W. Góźdź, V. K. Iglič, S. Kralj and A. Iglič, Coll. Surf. B, 2016, 141, 132–140. DOI:10.1016/j.colsurfb.2016.01.010 17. N. Gov, Philos. Trans. R. Soc. B Biol. Sci., 2018, 373, 20170115. DOI:10.1098/rstb.2017.0115 18. H. Hägerstrand, L. Mrówczyńska, U. Salzer, R. Prohaska, K. A. Michelsen, V. Kralj-Iglič and A. Iglič, Mol. Membr. Biol., 2006, 23, 277–288. DOI:10.1080/09687860600682536 19. A. Iglič, M. Lokar, B. Babnik, T. Slivnik, P. Veranič, H. Hägerstrand and V. Kralj-Iglič, Blood Cells. Mol. Dis., 2007, 39, 14–23. DOI:10.1016/j.bcmd.2007.02.013 20. A. Veksler and N. S. Gov, Biophys. J., 2007, 93, 3798–3810. DOI:10.1529/biophysj.107.113282 21. D. Hunter and B. Frisken, Biophys. J., 1998, 74, 2996–3002. DOI:10.1016/S0006-3495(98)78006-3 22. E. Karatekin, O. Sandre, H. Guitouni, N. Borghi, P.-H. Puech and F. Brochard-Wyart, Biophys. J., 2003, 84, 1734–1749. DOI:10.1016/S0006-3495(03)74981-9 23. R. Bruinsma, Physica A, 1996, 234, 249–270. DOI:10.1016/S0378-4371(96)00358-5 24. V. Vitkova, M. Mader and T. Podgorski, EPL Europhys. Lett., 2004, 68, 398. DOI:10.1209/epl/i2004-10211-9 25. M. M. Kozlov, F. Campelo, N. Liska, L. V. Chernomordik, S. J. Marrink and H. T. McMahon, Curr. Opin. Cell Biol., 2014, 29, 680 Acta Chim. Slov. 2020, 67, 674–681 Penič et al.: Active Forces of Myosin Motors May Control ... 53–60. DOI:10.1016/j.ceb.2014.03.006 26. A. Boulbitch, Phys. Rev. E, 1998, 57, 2123. DOI:10.1103/PhysRevE.57.2123 27. E. Evans and R. Skalak, Mechanics and Thermodynamics of Biomembranes., Boca Raton FL CRC, 1980. DOI:10.1115/1.3138234 28. D. E. Discher, in Physics of Biological Membranes, Springer, 2018, pp. 263–285. DOI:10.1007/978-3-030-00630-3_11 29. N. Mohandas and E. Evans, Annu. Rev. Biophys. Biomol. Struct., 1994, 23, 787–818. DOI:10.1146/annurev.bb.23.060194.004035 30. A. Iglič, J. Biomech., 1997, 30, 35–40. DOI:10.1016/S0021-9290(96)00100-5 31. A. Iglič, V. Kralj-Iglič and H. Hägerstrand, Eur. Biophys. J., 1998, 27, 335–339. DOI:10.1007/s002490050140 32. G. Lim H. W., M. Wortis and R. Mukhopadhyay, Proc. Natl. Acad. Sci. USA, 2002, 99, 16766–16769. DOI:10.1073/pnas.202617299 33. W. Helfrich, Z. Naturforsch. C, 1974, 29, 510–515. DOI:10.1515/znc-1974-9-1010 34. B. Stokke, A. Mikkelsen and A. Elgsaeter, Eur. Biophys. J., 1986, 13, 203–218. DOI:10.1007/BF00260368 35. R. Mukhopadhyay, H. G. Lim and M. Wortis, Biophys. J., 2002, 82, 1756–1772. DOI:10.1016/S0006-3495(02)75527-6 36. A. Iglič, S. Svetina and B. Žekš, Biophys. J., 1995, 69, 274–279. DOI:10.1016/S0006-3495(95)79899-X 37. H. Hägerstrand, V. Kralj-Iglič, M. Bobrowska- Hägerstrand and A. Iglič, Bull. Math. Biol., 1999, 61, 1019–1030. DOI:10.1103/PhysRevE.61.4230 38. E. J. Spangler, C. W. Harvey, J. D. Revalee, P. B. S. Kumar and M. Laradji, Phys. Rev. E, 2011, 84, 051906. DOI:10.1103/PhysRevE.84.051906 39. E. A. Evans, Biophys. J., 1974, 14, 923–931. DOI:10.1016/S0006-3495(74)85959-X 40. L. Miao, U. Seifert, M. Wortis and H.-G. Döbereiner, Phys. Rev. E, 1994, 49, 5389. DOI:10.1103/PhysRevE.49.5389 41. M. P. Sheetz and S. Singer, Proc. Natl. Acad. Sci. USA, 1974, 71, 4457–4461. DOI:10.1073/pnas.71.11.4457 42. V. Kralj-Iglič, H. Hägerstrand, P. Veranič, K. Jezernik, B. Bab- nik, D. R. Gauger and A. Iglič, Eur. Biophys. J., 2005, 34, 1066–1070. DOI:10.1007/s00249-005-0481-0 43. B. Deuticke, Biochim. Biophys. Acta - Biomembranes, 1968, 163, 494–500. DOI:10.1016/0005-2736(68)90078-3 44. H. Hägerstrand and B. Isomaa, Biochim. Biophys. Acta - Bi- omembranes, 1992, 1109, 117–126. DOI:10.1016/0005-2736(92)90074-V 45. V. Kralj-Iglič, A. Iglič, H. Hägerstrand and P. Peterlin, Phys. Rev. E, 2000, 61, 4230. DOI:10.1103/PhysRevE.61.4230 46. T. M. Fischer, J. Phys. II, 1992, 2, 337–343. DOI:10.1051/jp2:1992137 47. T. M. Fischer, J. Phys. II, 1993, 3, 1795–1805. DOI:10.1051/jp2:1993230 48. J. Fournier, Phys. Rev. Lett., 1996, 76, 4436–4439. DOI:10.1103/PhysRevLett.76.4436 49. C. Safinya, Coll. Surf. A, 1997, 128, 183–195. DOI:10.1016/S0927-7757(96)03914-3 50. J.-B. Fournier and P. Galatola, Braz. J. Phys., 1998, 28, 00–00. DOI:10.1590/S0103-97331998000400008 51. V. Kralj-Iglič, Int. J. Nanomed., 2012, 7, 3579. DOI:10.2147/IJN.S29076 52. L. Mesarec, W. Góźdź, A. Iglič, V. Kralj-Iglič, E. Virga and S. Kralj, Sci. Rep., 2019, 9, 1–11. DOI:10.1038/s41598-019-56128-0 53. N. Bobrovska, W. Góźdź, V. Kralj-Iglič and A. Iglič, PloS One, 2013, 8, e73941. DOI:10.1371/journal.pone.0073941 54. L. Mesarec, W. Góźdź, S. Kralj, M. Fošnarič, S. Penič, V. Kralj-Iglič and A. Iglič, Eur. Biophys. J., 2017, 46, 705–718. DOI:10.1007/s00249-017-1212-z 55. M. Fosnarič, S. Penič, A. Iglič, V. Kralj-Iglič, M. Drab and N. Gov, Soft Matter, 2019. 15, 5319–5330. DOI: 10.1039/C8SM02356E 56. H. Hägerstrand and B. Isomaa, Biochim. Biophys. Acta – Bi- omembranes, 1989, 982, 179–186. DOI:10.1016/0005-2736(89)90053-9 57. M. Bobrowska-Hägerstrand, V. Kralj-Iglič, A. Iglič, K. Bi- alkowska, B. Isomaa and H. Hägerstrand, Biophys. J., 1999, 77, 3356–3362. DOI:10.1016/S0006-3495(99)77167-5 58. A. Iglič, M. Fošnarič, H. Hägerstrand and V. Kralj-Iglič, FEBS Lett., 2004, 574, 9–12. DOI:10.1016/j.febslet.2004.07.085 59. H. Hägerstrand, V. Kralj-Iglič, M. Fošnarič, M. Bobrows- ka-Hägerstrand, A. Wróbel, L. Mrówczyńska, T. Söderström and A. Iglič, Biochim. Biophys. Acta - Biomembranes, 2004, 1665, 191–200. DOI:10.1016/j.bbamem.2004.08.010 60. A. S. Smith, R. B. Nowak, S. Zhou, M. Giannetto, D. S. Gokh- in, J. Papoin, I. C. Ghiran, L. Blanc, J. Wan and V. M. Fowler, Proc. Natl. Acad. Sci. USA, 2018, 115, E4377–E4385. DOI:10.1073/pnas.1718285115 61. H. Noguchi, J. Phys. Soc. Jpn., 2009, 78, 041007. DOI:10.1143/JPSJ.78.041007 62. G. Gompper and D. M. Kroll, in Statistical Mechanics of Mem- branes and Surfaces, eds. D. Nelson, T. Piran and S. Weinberg, World Scientific, Singapore, 2004, pp. 359–426. DOI:10.1142/9789812565518_0012 63. G. Gompper and D. M. Kroll, J. Phys. I, 1996, 6, 1305–1320. DOI:10.1051/jp1:1996246 64. M. Fošnarič, S. Penič, A. Iglič and I. Bivas, in Advances in Planar Lipid Bilayers and Liposomes, ed. A. Iglič & J. Genova, Elsevier, 2013, vol. 17, pp. 331–357. DOI:10.1016/B978-0-12-411516-3.00012-7 65. S. Penič, A. Iglič, I. Bivas and M. Fošnarič, Soft Matter, 2015, 11, 5004–5009. DOI: 10.1039/C5SM00431D. 66. W. Helfrich, Z. Naturforsch. C, 1973, 28, 693–703. DOI:10.1515/znc-1973-11-1209 67. D. Jesenek, S. Perutková, W. Góźdź, V. Kralj-Iglič, A. Iglič and S. Kralj, Int. J. Nanomed., 2013, 8, 677–687. DOI:10.2147/IJN.S38314 68. M. Fošnarič, A. Iglič and S. May, Phys. Rev. E, 2006, 74, 051503. DOI:10.1103/PhysRevE.74.051503 69. A. Iglič, T. Slivnik and V. Kralj-Iglič, J. Biomech., 2007, 40, 2492–2500. DOI:10.1016/j.jbiomech.2006.11.005 70. J. Käs and E. Sackmann, Biophys. J., 1991, 60, 825–844. DOI:10.1016/S0006-3495(91)82117-8 681Acta Chim. Slov. 2020, 67, 674–681 Penič et al.: Active Forces of Myosin Motors May Control ... Except when otherwise noted, articles in this journal are published under the terms and conditions of the  Creative Commons Attribution 4.0 International License Povzetek Z uporabo Monte Carlo (MC) simulacij stomatocitnih oblik eritrocita smo pokazali, da lahko aktivne sile miozinskih domen, pripetih na F-aktinske molekule na notranji strani membrane eritrocita, delno kontrolirajo endovezikulacijo membrane eritrocita. Miozinsko-aktinsko generirane aktivne sile, ki delujejo na membrano eritrocita, tako med drugim vplivajo tudi na nastanek ploščatih torocitnih endoveziklov eritrocitov. Na osnovi dobljenih rezultatov MC simulacij ve- zikulazije eritrocitov zaključujemo, da lahko miozinsko-aktinsko generirane aktivne sile, ki delujejo v smeri pravokotno na površino membrane eritrocita proti notranjosti eritrocita, pomembno vplivajo tudi na ostale transformacije oblik eri- trocita ter na stabilnost različnih oblik eritrocita in bi jih bilo zato potrebno upoštevati pri bodočih teoretičnih študijah transformacije oblik in vezikulacije eritrocitov.