DOI: 10.17344/acsi.2020.5863 Acta Chim. Slov. 2020, 67, - ©commons 1 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 generated 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 perpendicular 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 external medium, as well as, within the cell, between the lumens 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 proteins 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 curvature of the membrane,5-9 forcing the membrane to locally and/or globally adapt its shape.10-17 The phase separation of membrane inclusions (nanodomains) is an important 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 consuming forces are acting in the living cell.17 In experiments, the cells or lipid bilayer vesicles (as model systems) may also be deformed by extracellular forces due to pressure differences,21 stretching22 and fluid flow.23-25 An external force to the membrane surface can also be experimentally 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 response of the cell membrane shape to the force17 is important for different cell functions that range from adhesion Penič et al.: Active Forces of Myosin Motors May Control 2 Acta Chim. Slov. 2020, 67, - and migration to division, differentiation, and cell death.28 Therefore, the membrane mechanical properties are of great scientific interest. Consequently, new theoretical approaches for deeper understanding of the mechanisms that rule the cell's functions are constantly being developed. In red blood cells (RBC), besides the lipid bilayer also the membrane skeleton (spectrin-F-actin network attached 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 correspond 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 skeleton.30-32,35 The shear deformation of the membrane skeleton plays an essential role in the stability of echinocytic RBC shapes.30,31 Additionally, cytoskeleton induced protrusions in lipid membranes has been theoretically explained by the interplay between the elasticity of the membrane 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 layer, favors the concave local membrane shape, i.e. the inward 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 disruption 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 membrane 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 release 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 application of the local force.16,42,45,51,53,54 For isotropic membrane, 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 fluorescent 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 endovesicles with a unique ring-like toroidal shape joined with a central flat membrane segment, i.e. flat membrane structures with a bulby periphery called also torocyte endovesicles (Figure 1 A,B and E)57 with the shape very similar to the shape of Golgi bodies.58 It was proposed that the observed 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 mobile 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 vesicle 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 determination 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 local 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 Penič et al.: Active Forces of Myosin Motors May Control ... 2 Acta Chim. Slov. 2020, 67, - 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 microscopy 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 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 membrane inclusions (nanodomains) and the effect of the local active forces of NMIIA bipolar filament domains in a single 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 membranes 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 membrane into patches consisting of many molecules (Figure 2). A single patch is represented by a vertex in a triangulated 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 network can be implemented by ensuring that no vertex can penetrate through the triangular network. The maximal Penič et al.: Active Forces of Myosin Motors May Control ... 2 Acta Chim. Slov. 2020, 67, - 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 always 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 dj, d2, d3 and d4 are calculated and it is verified that they are between dmin and dmax. possible random displacement of the vertex in a single step (s) should be small enough so that the fourth vertex cannot move through the plane of the other three to the minimal 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 constraint dmax, see.64 The initial state of triangulated surface is a pentagonal 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 towards the exterior, negative curvature will force the membrane to bulge towards the interior compartment of the vesicle. The system is developed into the thermal equilibrium state. The evolution of the system is measured in Monte Carlo sweeps (mcs). One mcs consists of individual attempts 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 attempts 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 individual Monte Carlo step (vertex move or bond flip) is accepted with probability according to Metropolis Hastings algorithm, based on free energy change due to the Monte Carlo step. The energy W in simulation consists of three parts: where Wb is the bending energy of the membrane, W^ 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 membrane 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. where k is the bending stiffness of the membrane, c1, c2 and c0 are the two principal curvatures and the intrinsic curvature of the vesicle membrane at the point under consideration. 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 where w is a direct interaction constant. The energy is summed over all inclusion pairs with their in-plane distance rij , r0 is the range of direct interaction and H is a Heaviside step function with (r0-ri]) being function's argument. The energy contribution of the local protrusive forces due to the myosin motor domains/inclusions:55 where F is the magnitude of the force, the sum runs over all proteins, is the outwards facing normal to the membrane 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 inclusions (nanodomains) with negative intrinsic curvature. It can be seen in Figure 3 that the non-homogeneous lateral redistribution of the mobile inclusions with negative intrinsic 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 Penič et al.: Active Forces of Myosin Motors May Control ... 2 Acta Chim. Slov. 2020, 67, - Figure 3 depend on the inclusion concentration, the inclusion intrinsic curvature, the strength of the direct attractive interaction between inclusions and on the active forces exerted by the inclusions. Figure 3: Monte Carlo simulation of the RBC membrane transformation induced by mobile membrane inclusions with intrinsic curvature c0 = -1 imin-1. Two different membrane inclusion concentrations 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 kTHmn 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 triangulated 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: bending rigidity k = 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 protrusions (buds). Long undulated membrane protrusions may be further transformed into small independent spherical 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 detachment 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 perpendicular 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 observed in some experiments (see refs. 44,56 and the references 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 concentration of membrane inclusions exerting force on the membrane, the MC predicted RBC shape has one small and two large pancake-like torocyte membrane invaginations (Figure 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 complete invagination to the parents cell are supposed to be ruptured due to the frustrations in the orientational ordering of membrane components in the highly curved membrane necks.67 Note that in the torocytic membrane invaginations (Figure 3D), the myosin motor domains generated active forces that acted in the outward direction with respect to the torocyte invagination, i.e. in the direction towards 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 inclusions 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 additionally 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 invaginations and torocyte endovesicles in RBCs can be considered 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 (Figure 1) we have predicted in this work the invaginated stomatocyte RBC shapes having different shapes of invaginations, like torocytic, spherical, undulated necklace-like, etc. This is an extension of the previously theoretically predicted shape classes of the invaginated stomatocytic Penič et al.: Active Forces of Myosin Motors May Control ... 2 Acta Chim. Slov. 2020, 67, - shapes which were mostly limited to the simple (axisym-metric) stomatocytic shape with only one large invagination (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 membrane shapes and shape transformations is a widely used method to investigate physical models of various phenomena found experimentally in cellular systems. In the present 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 myosin domains generated forces on the RBC membrane favor 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 domains 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 Slovenian Research Agency (ARRS). Conflict of Interests The authors declare that there is no conflict of interests. 5. References 1. G. Cevc and D. Marsh, Phospholipid Bilayers: Physical Principles and Models, Wiley, 1987. 2. J. N. Israelachvili, Intermolecular and Surface Forces, Academic Press, 1997. 3. S. A. Safran, Statistical Thermodynamics of Surfaces, Interfaces, and Membranes, Addison-Wesley Publishing Company, Colorado, 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. Rezničkova, 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. Gozdz, 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. Mrowczynska, 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, Penič et al.: Active Forces of Myosin Motors May Control ... 2 Acta Chim. Slov. 2020, 67, - 53-60. D01:10.1016/j.ceb.2014.03.006 26. A. Boulbitch, Phys. Rev. E, 1998, 57, 2123. D0I: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. D0I: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. Iglic, J. Biomech., 1997, 30, 35-40. D0I:10.1016/S0021-9290(96)00100-5 31. A. Iglic, V. Kralj-Iglic and H. Hägerstrand, Eur. Biophys. J., 1998, 27, 335-339. D0I:10.1007/s002490050140 32. G. Lim H. W., M. Wortis and R. Mukhopadhyay, Proc. Natl. Acad. Sci. USA, 2002, 99, 16766-16769. D0I:10.1073/pnas.202617299 33. W. Helfrich, Z. Naturforsch. C, 1974, 29, 510-515. D0I:10.1515/znc-1974-9-1010 34. B. Stokke, A. Mikkelsen and A. Elgsaeter, Eur. Biophys. J., 1986, 13, 203-218. D0I:10.1007/BF00260368 35. R. Mukhopadhyay, H. G. Lim and M. Wortis, Biophys. J., 2002, 82, 1756-1772. D0I:10.1016/S0006-3495(02)75527-6 36. A. Iglic, S. Svetina and B. 2eks, Biophys. J., 1995, 69, 274-279. D0I:10.1016/S0006-3495(95)79899-X 37. H. Hägerstrand, V. Kralj-Iglic, M. Bobrowska- Hägerstrand and A. Iglic, Bull. Math. Biol., 1999, 61, 1019-1030. D0I: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. D0I:10.1103/PhysRevE.84.051906 39. E. A. Evans, Biophys. J., 1974, 14, 923-931. D0I: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. D0I:10.1103/PhysRevE.49.5389 41. M. P. Sheetz and S. Singer, Proc. Natl. Acad. Sci. USA, 1974, 71, 4457-4461. D0I:10.1073/pnas.71.11.4457 42. V. Kralj-Iglic, H. Hägerstrand, P. Veranic, K. Jezernik, B. Bab-nik, D. R. Gauger and A. Iglic, Eur. Biophys. J., 2005, 34, 1066-1070. D0I:10.1007/s00249-005-0481-0 43. B. Deuticke, Biochim. Biophys. Acta - Biomembranes, 1968, 163, 494-500. D0I:10.1016/0005-2736(68)90078-3 44. H. Hägerstrand and B. Isomaa, Biochim. Biophys. Acta - Biomembranes, 1992, 1109, 117-126. D0I:10.1016/0005-2736(92)90074-V 45. V. Kralj-Iglic, A. Iglic, H. Hägerstrand and P. Peterlin, Phys. Rev. E, 2000, 61, 4230. D0I:10.1103/PhysRevE.61.4230 46. T. M. Fischer, J. Phys. II, 1992, 2, 337-343. D0I:10.1051/jp2:1992137 47. T. M. Fischer, J. Phys. II, 1993, 3, 1795-1805. D0I:10.1051/jp2:1993230 48. J. Fournier, Phys. Rev. Lett., 1996, 76, 4436-4439. D01:10.1103/PhysRevLett.76.4436 49. C. Safinya, Coll. Surf. A, 1997, 128, 183-195. D0I:10.1016/S0927-7757(96)03914-3 50. J.-B. Fournier and P. Galatola, Braz. J. Phys., 1998, 28, 00-00. D01:10.1590/S0103-97331998000400008 51. V. Kralj-Iglic, Int. J. Nanomed., 2012, 7, 3579. D0I:10.2147/IJN.S29076 52. L. Mesarec, W. Gozdz, A. Iglic, V. Kralj-Iglic, E. Virga and S. Kralj, Sci. Rep., 2019, 9, 1-11. D0I:10.1038/s41598-019-56128-0 53. N. Bobrovska, W. Gozdz, V. Kralj-Iglic and A. Iglic, PloS One, 2013, 8, e73941. D0I:10.1371/journal.pone.0073941 54. L. Mesarec, W. Gozdz, S. Kralj, M. Fosnaric, S. Penic, V. Kralj-Iglic and A. Iglic, Eur. Biophys. J., 2017, 46, 705-718. D0I:10.1007/s00249-017-1212-z 55. M. Fosnaric, S. Penic, A. Iglic, V. Kralj-Iglic, 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 - Biomembranes, 1989, 982, 179-186. D01:10.1016/0005-2736(89)90053-9 57. M. Bobrowska-Hägerstrand, V. Kralj-Iglic, A. Iglic, K. Bi-alkowska, B. Isomaa and H. Hägerstrand, Biophys. J., 1999, 77, 3356-3362. D0I:10.1016/S0006-3495(99)77167-5 58. A. Iglic, M. Fosnaric, H. Hägerstrand and V. Kralj-Iglic, FEBS Lett., 2004, 574, 9-12. D0I:10.1016/j.febslet.2004.07.085 59. H. Hägerstrand, V. Kralj-Iglic, M. Fosnaric, M. Bobrows-ka-Hägerstrand, A. Wrobel, L. Mrowczynska, T. Söderström and A. Iglic, Biochim. Biophys. Acta - Biomembranes, 2004, 1665, 191-200. D0I:10.1016/j.bbamem.2004.08.010 60. A. S. Smith, R. B. Nowak, S. Zhou, M. Giannetto, D. S. Gokhin, J. Papoin, I. C. Ghiran, L. Blanc, J. Wan and V. M. Fowler, Proc. Natl. Acad. Sci. USA, 2018, 115, E4377-E4385. D0I:10.1073/pnas.1718285115 61. H. Noguchi, J. Phys. Soc. Jpn., 2009, 78, 041007. D01:10.1143/JPSJ.78.041007 62. G. Gompper and D. M. Kroll, in Statistical Mechanics of Membranes and Surfaces, eds. D. Nelson, T. Piran and S. Weinberg, World Scientific, Singapore, 2004, pp. 359-426. D0I:10.1142/9789812565518_0012 63. G. Gompper and D. M. Kroll, J. Phys. I, 1996, 6, 1305-1320. D0I:10.1051/jp1:1996246 64. M. Fosnaric, S. Penic, A. Iglic and I. Bivas, in Advances in Planar Lipid Bilayers and Liposomes, ed. A. Iglic & J. Genova, Elsevier, 2013, vol. 17, pp. 331-357. D0I:10.1016/B978-0-12-411516-3.00012-7 65. S. Penic, A. Iglic, I. Bivas and M. Fosnaric, Soft Matter, 2015, 11, 5004-5009. D0I: 10.1039/C5SM00431D. 66. W. Helfrich, Z. Naturforsch. C, 1973, 28, 693-703. D0I:10.1515/znc-1973-11-1209 67. D. Jesenek, S. Perutkova, W. Gozdz, V. Kralj-Iglic, A. Iglic and S. Kralj, Int. J. Nanomed., 2013, 8, 677-687. D0I:10.2147/IJN.S38314 68. M. Fosnaric, A. Iglic and S. May, Phys. Rev. E, 2006, 74, 051503. D0I:10.1103/PhysRevE.74.051503 69. A. Iglic, T. Slivnik and V. Kralj-Iglic, J. Biomech., 2007, 40, 2492-2500. D01:10.1016/j.jbiomech.2006.11.005 70. J. Käs and E. Sackmann, Biophys. J., 1991, 60, 825-844. D0I:10.1016/S0006-3495(91)82117-8 Penič et al.: Active Forces of Myosin Motors May Control ... 2 Acta Chim. Slov. 2020, 67, - 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 eritrocita 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. © (D Except when otherwise noted, articles in this journal are published under the terms and conditions of the Creative Commons Attribution 4.0 International License Penič et al.: Active Forces of Myosin Motors May Control ...