Effects of grafted polymers on the lipid membrane fluidity.

Since the membrane fluidity controls the cellular functions, it is important to identify the factors that determine the cell membrane viscosity. Cell membranes are composed of not only lipids and proteins but also polysaccharide chain-anchored molecules, such as glycolipids. To reveal the effects of grafted polymers on the membrane fluidity, in this study, we measured the membrane viscosity of polymer-grafted giant unilamellar vesicles (GUVs), which were prepared by introducing the poly (ethylene glycol) (PEG)-anchored lipids to the ternary GUVs composed of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), and cholesterol. The membrane viscosity was obtained from the velocity field on the GUV generated by applying a point force, based on the hydrodynamic model proposed by Henle and Levine. The velocity field was visualized by a motion of the circular liquid ordered (Lo) domains formed by a phase separation. With increasing PEG density, the membrane viscosity of PEG-grafted GUVs increased gradually in the mushroom region and significantly in the brush region. We propose a hydrodynamic model that includes the excluded volume effect of PEG chains to explain the increase in membrane viscosity in the mushroom region. This work provides a basic understanding of how grafted polymers affect the membrane fluidity.

Morphology of vesicle triplets: shape transformation at weak and strong adhesion limits.

We investigate the morphologies of adhering vesicle triplets as a function of volume-to-area ratio encoded by the reduced volume in strong and weak adhesion regimes. In the strong adhesion regime, the morphology change of the vesicle triplet depends on the arrangement of vesicles. By decreasing the reduced volume, a triangular triplet composed of three spherical caps with a trifurcated flat contact zone deformed to a compact spherical shape with a sigmoidal contact zone, whereas a linear vesicle triplet composed of pancake-shaped vesicles sandwiched between two spherical-cap vesicles with a flat contact zone deformed into a compact spherical shape with biconvex interfaces. The morphologies of vesicle triplets with flat contact zones are reproduced by the so-called two-tension model based on the total energy consisting of bending energy, adhesion energy and surface energy, where the surface tension in the noncontact zone is different from that in the contact zone. When the flat interface deforms, the two-tension model is no longer applicable. The compact spherical triplets with curved interfaces can be reproduced by introducing geometrical constraints requiring that the total area of the non-contact zones is minimal, thereby confining the aggregate to a spherical cavity; this is referred to as the cavity model. In the weak adhesion regime, vesicle triplets with either a triangular or linear topology deform into prolate-based triplets by decreasing the reduced volume.

Synthesising a minimal cell with artificial metabolic pathways.

A "synthetic minimal cell" is considered here as a cell-like artificial vesicle reproduction system in which a chemical and physico-chemical transformation network is regulated by information polymers. Here we synthesise such a minimal cell consisting of three units: energy production, information polymer synthesis, and vesicle reproduction. Supplied ingredients are converted to energy currencies which trigger the synthesis of an information polymer, where the vesicle membrane plays the role of a template. The information polymer promotes membrane growth. By tuning the membrane composition and permeability to osmolytes, the growing vesicles show recursive reproduction over several generations. Our "synthetic minimal cell" greatly simplifies the scheme of contemporary living cells while keeping their essence. The chemical pathways and the vesicle reproduction pathways are well described by kinetic equations and by applying the membrane elasticity model, respectively. This study provides new insights to better understand the differences and similarities between non-living forms of matter and life.

From vesicles toward protocells and minimal cells.

In contrast to ordinary condensed matter systems, "living systems" are unique. They are based on molecular compartments that reproduce themselves through (i) an uptake of ingredients and energy from the environment, and (ii) spatially and timely coordinated internal chemical transformations. These occur on the basis of instructions encoded in information molecules (DNAs). Life originated on Earth about 4 billion years ago as self-organised systems of inorganic compounds and organic molecules including macromolecules (e.g. nucleic acids and proteins) and low molar mass amphiphiles (lipids). Before the first living systems emerged from non-living forms of matter, functional molecules and dynamic molecular assemblies must have been formed as prebiotic soft matter systems. These hypothetical cell-like compartment systems often are called "protocells". Other systems that are considered as bridging units between non-living and living systems are called "minimal cells". They are synthetic, autonomous and sustainable reproducing compartment systems, but their constituents are not limited to prebiotic substances. In this review, we focus on both membrane-bounded (vesicular) protocells and minimal cells, and provide a membrane physics background which helps to understand how morphological transformations of vesicle systems might have happened and how vesicle reproduction might be coupled with metabolic reactions and information molecules. This research, which bridges matter and life, is a great challenge in which soft matter physics, systems chemistry, and synthetic biology must take joined efforts to better understand how the transformation of protocells into living systems might have occurred at the origin of life.

Vesicle deformation and division induced by flip-flops of lipid molecules.

We investigated the deformation of small unilamellar vesicles (SUVs) induced by flip-flops of lipids using coarse-grained molecular dynamics simulations. In the case of single-component SUVs composed of zero spontaneous curvature lipids (ZLs), the flip-flop of ZLs deformed stomatocyte-shaped SUVs into an oblate shape, whereas pear-shaped SUVs were deformed into a prolate shape. These two equilibrium shapes comply with the local minima of elastic energy. In the case of binary vesicles composed of ZLs and negative spontaneous curvature lipids (NLs), the vesicle deformation pathway depended on the initial NL distribution in the bilayer. If the initial difference in the NL concentration between the outer and inner leaflets was small, the flip-flop of ZLs and NLs rapidly deformed pear-shaped SUVs into an equilibrium prolate shape. On the other hand, when NLs were localised in the inner leaflet, the flip-flop of ZLs and NLs deformed pear-shaped SUVs into a limiting shape and then induced vesicle division. Because the flip-flop rate of NLs is much faster than that of ZLs, the total free energy was first relaxed by the flip-flop of NLs and then by that of ZLs. This kinetic effect is responsible for the observed vesicle division induced by flip-flops.

Morphologies of Vesicle Doublets: Competition among Bending Elasticity, Surface Tension, and Adhesion.

To study the mechanical laws governing the form of multicellular organisms, we examine the morphology of adhering vesicle doublets as the simplest model system. We monitor the morphological transformations of doublets induced by changes of adhesion strength and volume/area ratio, which are controlled by intermembrane interactions and thermal area expansion, respectively. When we increase the temperature in the weak adhesion regime, a dumbbell flat-contact doublet is transformed to a parallel-prolate doublet, whereas in the strong adhesion regime, heating transforms the dumbbell flat-contact doublet into a spherical sigmoid-contact doublet. We reproduce the observed doublet morphologies by numerically minimizing the total energy, including the contact-potential adhesion term as well as the surface and bending terms, using the Surface Evolver package. From the reproduced morphologies, we extract the adhesion strength, the surface tension, and the volume/area ratio of the vesicles, which reveals the detailed mechanisms of the morphological transitions in doublets.

Viscosity Landscape of Phase-Separated Lipid Membrane Estimated from Fluid Velocity Field.

In cell membranes, the functional constituents such as peptides, proteins, and polysaccharides diffuse in a sea of lipids as single molecules and molecular aggregates. Thus, the fluidity of the heterogeneous multicomponent membrane is important for understanding the roles of the membrane in cell functionality. Recently, Henle and Levine described the hydrodynamics of molecular diffusion in a spherical membrane. A tangential point force at the north pole induces a pair of vortices whose centers lie on a line perpendicular to the point force and are symmetrical with respect to the point force. The position of the vortex center depends on ηm/Rηw, where R is the radius of the spherical membrane, and ηm and ηw are the viscosities of the membrane and the surrounding medium, respectively. Based on this theoretical prediction, we applied a point force to a phase-separated spherical vesicle composed of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine/1,2-dioleoyl-sn-glycero-3-phosphocholine/cholesterol by means of a microinjection technique. The pathlines were visualized by trajectories of microdomains. We determined the position of the vortex center and estimated the membrane viscosity using the dependence of the position of the vortex center on ηm/Rηw. The obtained apparent membrane viscosities for various compositions are mapped on the phase diagram. The membrane viscosity is almost constant in the range of 0 <ϕLo ≤ 0.5 (ϕLo: area fraction of the liquid ordered phase), whereas that in the range of 0.5 ≤ ϕLo < 1.0 exponentially increases with increase of ϕLo. The obtained viscosity landscape provides a basic understanding of the fluidity of heterogeneous multicomponent membranes.

Migration of Deformable Vesicles Induced by Ionic Stimuli.

We have investigated the dynamics of phospholipid vesicles composed of 1,2-dioleoyl- sn-glycero-3-phosphocholine triggered by ionic stimuli using electrolytes such as CaCl2, NaCl, and NaOH. The ionic stimuli induce two characteristic vesicle dynamics, deformation due to the ion binding to the lipids in the outer leaflet of the vesicle and migration due to the concentration gradient of ions, that is, diffusiophoresis or the interfacial energy gradient mechanism. We examined the deformation pathway for each electrolyte as a function of time and analyzed it based on the surface dissociation model and the area difference elasticity model, which reveals the change of the cross-sectional area of the phospholipid by the ion binding. The metal ions such as Ca2+ and Na+ encourage inward budding deformation by decreasing the cross-sectional area of a lipid, whereas the hydroxide ion (OH-) encourages outward budding deformation by increasing the cross-sectional area of a lipid. When we microinjected these electrolytes toward the vesicles, a strong coupling between the deformation and the migration of the vesicle was observed for CaCl2 and NaOH, whereas for NaCl, the coupling was very weak. This difference probably originates from the binding constants of the ions.

Molecular mechanism of vesicle division induced by coupling between lipid geometry and membrane curvatures.

We investigated the effects of lipid geometry on vesicle division using coarse grained molecular dynamics simulations. When the vesicle is composed of zero and negative spontaneous curvature lipids (ZSLs and NSLs), the difference in their molecular spontaneous curvatures destabilizes the neck of the limiting shape vesicle. In the vesicle division pathway, the neck developed into the stalk intermediates. The stalk was broken when the NSLs were expelled from the stalk. Free energy analysis shows that the coupling between the lipid geometry and the Gaussian rigidity is responsible for the observed vesicle division.

pH sensing by lipids in membranes: The fundamentals of pH-driven migration, polarization and deformations of lipid bilayer assemblies.

Most biological molecules contain acido-basic groups that modulate their structure and interactions. A consequence is that pH gradients, local heterogeneities and dynamic variations are used by cells and organisms to drive or regulate specific biological functions including energetic metabolism, vesicular traffic, migration and spatial patterning of tissues in development. While the direct or regulatory role of pH in protein function is well documented, the role of hydrogen and hydroxyl ions in modulating the properties of lipid assemblies such as bilayer membranes is only beginning to be understood. Here, we review approaches using artificial lipid vesicles that have been instrumental in providing an understanding of the influence of pH gradients and local variations on membrane vectorial motional processes: migration, membrane curvature effects promoting global or local deformations, crowding generation by segregative polarization processes. In the case of pH induced local deformations, an extensive theoretical framework is given and an application to a specific biological issue, namely the structure and stability of mitochondrial cristae, is described. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Migration of Phospholipid Vesicles Can Be Selectively Driven by Concentration Gradients of Metal Chloride Solutions.

We have investigated the migrations of phospholipid vesicles under the concentration gradients of metal ions. We microinjected metal chloride solutions, monovalent (NaCl and KCl), divalent (CaCl2 and MgCl2), and trivalent (LaCl3) salts, toward phospholipid giant vesicles (GVs) composed of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). For NaCl, CaCl2, and MgCl2 solutions, the GVs migrated straight toward the tip of the micropipette in response to the concentration gradients, whereas for KCl and LaCl3, GVs moved to the opposite direction. Our motion tracking of lipid domains in a vesicle membrane showed no unidirectional flow in the membrane during the vesicle migration, indicating that the Marangoni mechanism is not responsible for the observed vesicle migration. We calculated the diffusiophoretic velocities for symmetric and asymmetrical electrolytes by solving the Stokes' equation numerically. The theoretical diffusiophoretic velocities described the observed migration velocities well. Thus, we can control the migration of vesicle in response to the concentration gradient by adapting the electrolytes and the lipids.

Role of Inverse-Cone-Shape Lipids in Temperature-Controlled Self-Reproduction of Binary Vesicles.

We investigate a temperature-driven recursive division of binary giant unilamellar vesicles (GUVs). During the heating step of the heating-cooling cycle, the spherical mother vesicle deforms to a budded limiting shape using up the excess area produced by the chain melting of the lipids and then splits off into two daughter vesicles. Upon cooling, the daughter vesicle opens a pore and recovers the spherical shape of the mother vesicle. Our GUVs are composed of DLPE (1,2-dilauroyl-sn-glycero-3-phosphoethanolamine) and DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphocholine). During each cycle, vesicle deformation is monitored by a fast confocal microscope and the images are analyzed to obtain the time evolution of reduced volume and reduced monolayer area difference as the key geometric parameters that quantify vesicle shape. By interpreting the deformation pathway using the area-difference elasticity theory, we conclude that vesicle division relies on (1) a tiny asymmetric distribution of DLPE within the bilayer, which controls the observed deformation from the sphere to the budded shape; and (2) redistribution of DLPE during the deformation-division stage, which ensures that the process is recursive. The spontaneous coupling between membrane curvature and PE lipid distribution is responsible for the observed recursive division of GUVs. These results shed light on the mechanisms of vesicle self-reproduction.

Migration of phospholipid vesicles in response to OH(-) stimuli.

We demonstrate migration of phospholipid vesicles in response to a pH gradient. Upon simple micro-injection of a NaOH solution, the vesicles linearly moved to the tip of the micro-pipette and the migration velocity was proportional to the gradient of OH(-) concentration. Vesicle migration was characteristic of OH(-) ions and no migration was observed for monovalent salts or nonionic sucrose solutions. The migration of vesicles is quantitatively described by the surface tension gradient model where the hydrolysis of the phospholipids by NaOH solution decreases the surface tension of the vesicle. The vesicles move toward a direction where the surface energy decreases. Thus the chemical modification of lipids produces a mechanical force to drive vesicles.

Dynamics of fatty acid vesicles in response to pH stimuli.

We investigate the dynamics of decanoic acid/decanoate (DA) vesicles in response to pH stimuli. Two types of dynamic processes induced by the micro-injection of NaOH solutions are sequentially observed: deformations and topological transitions. In the deformation stage, DA vesicles show a series of shape deformations, i.e., prolate-oblate-stomatocyte-sphere. In the topological transition stage, spherical DA vesicles follow either of the two pathways, pore formation and vesicle fusion. The pH stimuli modify a critical aggregation concentration of DA molecules, which causes the solubilization of DA molecules in the outer leaflet of the vesicle bilayers. This solubilization decreases the outer surface area of the vesicle, thereby increasing surface tension. A kinetic model based on area difference elasticity theory can accurately describe the dynamics of DA vesicles triggered by pH stimuli.

From vesicles to protocells: the roles of amphiphilic molecules.

It is very challenging to construct protocells from molecular assemblies. An important step in this challenge is the achievement of vesicle dynamics that are relevant to cellular functions, such as membrane trafficking and self-reproduction, using amphiphilic molecules. Soft matter physics will play an important role in the development of vesicles that have these functions. Here, we show that simple binary phospholipid vesicles have the potential to reproduce the relevant functions of adhesion, pore formation and self-reproduction of vesicles, by coupling the lipid geometries (spontaneous curvatures) and the phase separation. This achievement will elucidate the pathway from molecular assembly to cellular life.

Tubular membrane formation of binary giant unilamellar vesicles composed of cylinder and inverse-cone-shaped lipids.

We have succeeded in controlling tubular membrane formations in binary giant unilamellar vesicles (GUVs) using a simple temperature changing between the homogeneous one-phase region and the two-phase coexistence region. The binary GUV is composed of inverse-cone (bulky hydrocarbon chains and a small headgroup) and cylinder-shaped lipids. When the temperature was set in the two-phase coexistence region, the binary GUV had a spherical shape with solidlike domains. By increasing the temperature to the homogeneous one-phase region, the excess area created by the chain melting of the lipid produced tubes inside the GUV. The tubes had a radius on the micrometer scale and were stable in the one-phase region. When we again decreased the temperature to the two-phase coexisting region, the tubes regressed and the GUVs recovered their phase-separated spherical shape. We infer that the tubular formation was based on the mechanical balance of the vesicle membrane (spontaneous tension) coupled with the asymmetric distribution of the inverse-cone-shaped lipids between the inner and outer leaflets of the vesicle (lipid sorting).

Model system of self-reproducing vesicles.

Development of self-reproducing vesicle systems is the first step for autopoietic cycles. We established a model self-reproducing vesicle system without the membrane molecule synthesis route. The model vesicle composed of cylinder- and inverse-cone-shaped lipids formed inclusion vesicles inside the mother vesicle, and the inclusion vesicles were then expelled by a temperature cycling. By changing the vesicle composition, the mother vesicles showed a budding-type self-reproduction pathway. A key concept of this system is the coupling of the main-chain transition and the shape of lipids.

Pore formation in a binary giant vesicle induced by cone-shaped lipids.

We have investigated shape deformations of binary giant unilamellar vesicles (GUVs) composed of cone- and cylinder-shaped lipids. By coupling the spontaneous curvature of lipids with the phase separation, we demonstrated pore opening and closing in GUVs. When the temperature was set below the chain melting transition temperature of the cylinder-shaped lipid, the GUVs burst and then formed a single large pore, where the cone shape lipids form a cap at the edge of the bilayer to stabilize the pore. The pore closed when we increased the temperature above the transition temperature. The pore showed three types of shapes depending on the cone-shaped lipid concentration: simple circular, rolled-rim, and wrinkled-rim pores. These pore shape changes indicate that the distribution of the cone- and cylinder-shaped lipids is asymmetric between the inner and outer leaflets in the bilayer. We have proposed a theoretical model for a two-component membrane with an edge of bilayer where lipids can transfer between two leaflets. Using this model, we have reproduced numerically the observed shape deformations at the rim of pore.