Monte Carlo simulations of fluid vesicles with in-plane orientational ordering

Mesoscale simulation of biomembranes with FreeDTS

We present FreeDTS software for performing computational research on biomembranes at the mesoscale. In this software, a membrane is represented by a dynamically triangulated surface equipped with vertex-based inclusions to integrate the effects of integral and peripheral membrane proteins. Several algorithms are included in the software to simulate complex membranes at different conditions such as framed membranes with constant tension, vesicles and high-genus membranes with various fixed volumes or constant pressure differences and applying external forces to membrane regions. Furthermore, the software allows the user to turn off the shape evolution of the membrane and focus solely on the organization of proteins. As a result, we can take realistic membrane shapes obtained from, for example, cryo-electron tomography and backmap them into a finer simulation model. In addition to many biomembrane applications, this software brings us a step closer to simulating realistic biomembranes with molecular resolution. Here we provide several interesting showcases of the power of the software but leave a wide range of potential applications for interested users.

Effects of receptor properties on particle internalization through receptor-mediated endocytosis

Receptor-mediated endocytosis (RME) is a highly complex process carried out by bioparticles, such as viruses and drug carriers, to enter cells. The discovery of both clathrin-dependent and clathrin-free pathways makes the RME process even more intriguing. Numerical models have been developed to facilitate the exploration of the process. However, the impacts of the receptor properties on RME have been less studied partially due to the oversimplifications of the receptor models. In this paper, we implement a stochastic model to systematically investigate the effects of mechanical (receptor flexure), geometrical (receptor length) and biochemical (ligand-receptor cutoff) properties of receptors, on RME with and without the existence of clathrin. Our simulation results show that the receptor's flexural rigidity plays an important role in RME with clathrin. There is a threshold beyond which particle internalization will not occur. Without clathrin, it is very difficult to achieve complete endocytosis with ligand-receptor interactions alone. A shorter receptor length and longer ligand-receptor reaction cutoff promote the formation of ligand-receptor bonds and facilitate particle internalization. Complete internalization can only be obtained with an extremely short receptor length and long reaction cutoff. Therefore, there are most likely some additional mechanisms to drive the membrane deformation in clathrin-free RME. Our results yield important fundamental insights into RME and provide crucial guidance when correlating the simulation results with experimental observations.

Synthetic Membrane Shaper for Controlled Liposome Deformation

Shape defines the structure and function of cellular membranes. In cell division, the cell membrane deforms into a "dumbbell" shape, while organelles such as the autophagosome exhibit "stomatocyte" shapes. Bottom-up in vitro reconstitution of protein machineries that stabilize or resolve the membrane necks in such deformed liposome structures is of considerable interest to characterize their function. Here we develop a DNA-nanotechnology-based approach that we call the synthetic membrane shaper (SMS), where cholesterol-linked DNA structures attach to the liposome membrane to reproducibly generate high yields of stomatocytes and dumbbells. In silico simulations confirm the shape-stabilizing role of the SMS. We show that the SMS is fully compatible with protein reconstitution by assembling bacterial divisome proteins (DynaminA, FtsZ:ZipA) at the catenoidal neck of these membrane structures. The SMS approach provides a general tool for studying protein binding to complex membrane geometries that will greatly benefit synthetic cell research.

Crowding-induced membrane remodeling: Interplay of membrane tension, polymer density, architecture

The plasma membrane hosts a wide range of biomolecules, mainly proteins and carbohydrates, that mediate cellular interactions with its environment. The crowding of such biomolecules regulates cellular morphologies and cellular trafficking. Recent discoveries have shown that the structure and density of cell surface polymers and hence the signaling machinery change with the state of the cell, especially in cancer progression. The alterations in membrane-attached glycocalyx and glycosylation of proteins and lipids are common features of cancer cells. The overexpression of glycocalyx polymers, such as mucin and hyaluronan, strongly correlates with cancer metastasis. Here, we present a mesoscale biophysics-based model that accounts for the shape regulation of membranes by crowding of membrane-attached biopolymer-glycocalyx and actin networks. Our computational model is based on the dynamically triangulated Monte Carlo model for membranes and coarse-grained representations of polymer chains. The model allows us to investigate the crowding-induced shape transformations in cell membranes in a tension- and graft polymer density-dependent manner. Our results show that the number of membrane protrusions and their shape depend on membrane tension, with higher membrane tension inducing more tubular protrusions than the vesicular shapes formed at low tension at high surface coverage of polymers. The shape transformations occur above the threshold density predicted by the polymer brush theory, but this threshold also depends on the membrane tension. Increasing the size of the polymer, either by changing the length or by adding side chains, is shown to increase the crowding-induced curvature. The effect of crowding is more prominent for flexible polymers than for semiflexible rigid polymers. We also present an extension of the model that incorporates properties of the actin-like filament networks and demonstrate how tubular structures can be generated by biopolymer crowding on the cytosolic side of cell membranes.

Membrane Remodeling Due to a Mixture of Multiple Types of Curvature Proteins

We present an extension of the Monte Carlo based mesoscopic membrane model, where the membrane is represented as a dynamically triangulated surface and the proteins are modeled as anisotropic inclusions formulated as in-plane nematic field variables adhering to the deformable elastic sheet. In the extended model, we have augmented the Hamiltonian to study membrane deformation due to a mixture of multiple types of curvature generating proteins. This feature opens the door for understanding how multiple kinds of curvature-generating proteins may be working in a coordinated manner to induce desired membrane morphologies. For example, among other things, we study membrane deformations and tubulation due to a mixture of positive and negative curvature proteins as mimics of various proteins from BAR domain family. We also study the effect of membrane anisotropy that manifests as differential binding affinity and organization of curvature proteins, leading to insights into the tightly regulated cargo sorting and transport processes. Our simulation results show different morphology of deformed vesicles that depend on membrane tension, the curvatures and number of the participating proteins as well as on protein-protein and membrane-protein interactions.

Data driven and biophysical insights into the regulation of trafficking vesicles by extracellular matrix stiffness

Biomechanical signals from remodeled extracellular matrix (ECM) promote tumor progression. Here, we show that cell-matrix and cell-cell communication may be inherently linked and tuned through mechanisms of mechanosensitive biogenesis of trafficking vesicles. Pan-cancer analysis of cancer cells' mechanical properties (focusing primarily on cell stiffness) on substrates of varied stiffness and composition elucidated a heterogeneous cellular response to mechanical stimuli. Through machine learning, we identified a fingerprint of cytoskeleton-related proteins that accurately characterize cell stiffness in different ECM conditions. Expression of their respective genes correlates with patient prognosis across different tumor types. The levels of selected cytoskeleton proteins indicated that cortical tension mirrors the increase (or decrease) in cell stiffness with a change in ECM stiffness. A mechanistic biophysical model shows that the tendency for curvature generation by curvature-inducing proteins has an ultrasensitive dependence on cortical tension. This study thus highlights the effect of ECM stiffness, mediated by cortical tension, in modulating vesicle biogenesis.

Quantification of Curvature Sensing Behavior of Curvature-Inducing Proteins on Model Wavy Substrates

Curvature-inducing proteins are involved in a variety of membrane remodeling processes in the cell. Several in vitro experiments have quantified the curvature sensing behavior of these proteins in model lipid systems. One such system consists of a membrane bilayer laid atop a wavy substrate (Hsieh in Langmuir 28:12838-12843, 2012). In these experiments, the bilayer conforms to the wavy substrate, and curvature-inducing proteins show preferential segregation on the wavy membrane. Using a mesoscale computational membrane model based on the Helfrich Hamiltonian, here we present a study which analyzes the curvature sensing characteristics of this membrane-protein system, and elucidates key physical principles governing protein segregation on the wavy substrate and other in vitro systems. In this article we compute the local protein densities from the free energy landscape associated with membrane remodeling by curvature-inducing proteins. In specific, we use the Widom insertion technique to compute the free energy landscape for an inhomogeneous system with spatially varying density and the results obtained with this minimal model show excellent agreement with experimental studies that demonstrate the association between membrane curvature and local protein density. The free energy-based framework employed in this study can be used for different membrane morphologies and varied protein characteristics to gain mechanistic insights into protein sorting on membranes.

Probing lipid membrane bending mechanics using gold nanorod tracking

Lipid bilayer membranes undergo rapid bending undulations with wavelengths from tens of nanometers to tens of microns due to thermal fluctuations. Here, we probe such undulations and the membranes' mechanics by measuring the time-varying orientation of single gold nanorods (GNRs) adhered to the membrane, using high-speed dark field microscopy. In a lipid vesicle, such measurements allow the determination of the membrane's viscosity, bending rigidity, and tension as well as the friction coefficient for sliding of the monolayers over one another. The in-plane rotation of the GNR is hindered by undulations in a tension dependent manner, consistent with simulations. The motion of single GNRs adhered to the plasma membrane of living cultured cells similarly reveals the membrane's complex physics and coupling to the cell's actomyosin cortex.

Sorting of proteins with shape and curvature anisotropy on a lipid bilayer tube

Curvature induced sorting of lipid membrane bound proteins has been widely studied through experiments that induce curvature variation in a giant unilamellar lipid-bilayer vesicle with adsorbed proteins by pulling thin cylindrical tethers. In the theoretical space, this has been supplemented with models that capture curvature dependent interaction between membrane and idealized protein particles, through free energy contributions. Many membrane proteins such as the BAR domain proteins are known to have extremely anisotropic shapes and soft interacting potentials, whereas the idealizations of protein particles explored in models have only assumed them as hard disk-like particles with curvature anisotropy. Here, we present a model of sorting of the proteins while including the effects of softness in their interaction potentials, shape anisotropy in the protein structure, and curvature anisotropy in the interactions with the membrane. This is based on a clean separation of free energy contributions from non-ideal fluid behavior of soft anisotropic particles and curvature interactions between proteins and membranes. We probe the behavior of the sorting function under limiting conditions and show that it converges to the previously derived models. In addition to this, we present a comparison of the variation in sorting ratio due to the observed variation in the shape parameter values in known membrane proteins. Finally, using published experimental data for membrane proteins, we perform fitting and derive model parameters. We observe that shape anisotropy adversely affects the sorting of proteins to a high curvature region, whereas curvature anisotropy and softer interaction between proteins favor sorting.

Deformation of membrane vesicles due to chiral surface proteins

Chiral, rod-like molecules can self-assemble into cylindrical membrane tubules and helical ribbons. They have been successfully modeled using the theory of chiral nematics. Models have also predicted the role of chiral lipids in forming nanometer-sized membrane buds in the cell. However, in most theoretical studies, the membrane shapes are considered fixed (cylinder, sphere, saddle, etc.), and their optimum radii of curvatures are found variationally by minimizing the energy of the composite system consisting of membrane and chiral nematics. Numerical simulations have only recently started to consider membrane deformation and chiral orientation simultaneously. Here we examine how deformable, closed membrane vesicles and chiral nematic rods mutually influence each other's shape and orientation, respectively, using Monte Carlo (MC) simulation on a closed triangulated surface. For this, we adopt a discrete form of chiral interaction between rods, originally proposed by Van der Meer et al., for off-lattice simulations. In our simulation, both conical and short cylindrical tubules emerge, depending on the membrane stiffness and the intrinsic chirality of the molecules. We show that the Helfrich-Prost term, which couples nematic tilt with local membrane curvature in continuum models, can account for most of the observations in the simulation. At higher chirality, our theory also predicts a chiral tweed phase on cones, with varying bandwidths.

Deformation and orientational order of chiral membranes with free edges

Motivated by experiments on colloidal membranes composed of chiral rod-like viruses, we use Monte Carlo methods to simulate these systems and determine the phase diagram for the liquid crystalline order of the rods and the membrane shape. We generalize the Lebwohl-Lasher model for a nematic with a chiral coupling to a curved surface with edge tension and a resistance to bending, and include an energy cost for tilting of the rods relative to the local membrane normal. The membrane is represented by a triangular mesh of hard beads joined by bonds, where each bead is decorated by a director. The beads can move, the bonds can reconnect and the directors can rotate at each Monte Carlo step. When the cost of tilt is small, the membrane tends to be flat, with the rods only twisting near the edge for low chiral coupling, and remaining parallel to the normal in the interior of the membrane. At high chiral coupling, the rods twist everywhere, forming a cholesteric state. When the cost of tilt is large, the emergence of the cholesteric state at high values of the chiral coupling is accompanied by the bending of the membrane into a saddle shape. Increasing the edge tension tends to flatten the membrane. These results illustrate the geometric frustration arising from the inability of a surface normal to have twist.

Creasing of flexible membranes at vanishing tension

The properties of freestanding tensionless interfaces and membranes at low bending rigidity κ are dominated by strong fluctuations and self-avoidance and are thus outside the range of standard perturbative analysis. We analyze this regime by a simple discretized, self-avoiding membrane model on a frame subject to periodic boundary conditions by use of Monte Carlo simulation and dynamically triangulated surface techniques. We find that at low bending rigidities, the membrane properties fall into three regimes: Below the collapse transition κ_{BP} it is subject to branched polymer instability where the framed surface is not defined, in a range below a threshold rigidity κ_{c} the conformational correlation function are characterized by power-law behavior with a continuously varying exponent α, 2<α≤4 and above κ_{c}, α=4 characteristic for linearized bending excitations. Response functions specific heat and area compressibility display pronounced peaks close to κ_{c}. The results may be important for the description of soft interface systems, such as microemulsions and membranes with in-plane cooperative phenomena.

Annexin A4 trimers are recruited by high membrane curvatures in giant plasma membrane vesicles

The plasma membrane (PM) of eukaryotic cells consists of a crowded environment comprised of a high diversity of proteins in a complex lipid matrix. The lateral organization of membrane proteins in the PM is closely correlated with biological functions such as endocytosis, membrane budding and other processes which involve protein mediated shaping of the membrane into highly curved structures. Annexin A4 (ANXA4) is a prominent player in a number of biological functions including PM repair. Its binding to membranes is activated by Ca2+ influx and it is therefore rapidly recruited to the cell surface near rupture sites where Ca2+ influx takes place. However, the free edges near rupture sites can easily bend into complex curvatures and hence may accelerate recruitment of curvature sensing proteins to facilitate rapid membrane repair. To analyze the curvature sensing behavior of curvature inducing proteins in crowded membranes, we quantifify the affinity of ANXA4 monomers and trimers for high membrane curvatures by extracting membrane nanotubes from giant PM vesicles (GPMVs). ANXA4 is found to be a sensor of negative membrane curvatures. Multiscale simulations, in which we extract molecular information from atomistic scale simulations as input to our macroscopic scale simulations, furthermore predicted that ANXA4 trimers generate membrane curvature upon binding and have an affinity for highly curved membrane regions only within a well defined membrane curvature window. Our results indicate that curvature sensing and mobility of ANXA4 depend on the trimer structure of ANXA4 which could provide new biophysical insight into the role of ANXA4 in membrane repair and other biological processes.

Backmapping triangulated surfaces to coarse-grained membrane models

Many biological processes involve large-scale changes in membrane shape. Computer simulations of these processes are challenging since they occur across a wide range of spatiotemporal scales that cannot be investigated in full by any single current simulation technique. A potential solution is to combine different levels of resolution through a multiscale scheme. Here, we present a multiscale algorithm that backmaps a continuum membrane model represented as a dynamically triangulated surface (DTS) to its corresponding molecular model based on the coarse-grained (CG) Martini force field. Thus, we can use DTS simulations to equilibrate slow large-scale membrane conformational changes and then explore the local properties at CG resolution. We demonstrate the power of our method by backmapping a vesicular bud induced by binding of Shiga toxin and by transforming the membranes of an entire mitochondrion to near-atomic resolution. Our approach opens the way to whole cell simulations at molecular detail.

Computer simulations of large-scale changes in membrane shape are challenging since they occur across a wide range of spatiotemporal scales. Here, authors present a multiscale algorithm that backmaps a continuum membrane model represented as a dynamically triangulated surface to its corresponding molecular model based on the coarse-grained Martini force field.

Confined filaments in soft vesicles - the case of sickle red blood cells

Abnormal shapes of red blood cells (RBC) have been associated with various diseases. Diverse RBC shapes have also been intriguing for membrane biophysics. Here we focus on sickle shaped RBC which form due to abnormal growth of semi-rigid hemoglobin (HbS) fibers confined in RBC. Using the area difference elasticity (ADE) model for RBC and worm-like chain model for the confined HbS fibers, we explore shape deformations at equilibrium using Monte-Carlo simulations. We show that while a single HbS fiber is not rigid enough to produce sickle like deformation, a fiber bundle can do so. We also consider multiple disjoint filaments and find that confinement can generate multipolar RBC shapes and can even promote helical filament conformations which have not been discussed before. We show that the same model, when applied to microtubules confined in phospholipid vesicles, predicts vesicle tubulation. In addition we reproduce the tube collapse transition and tennis racket type vesicle shapes, as reported in experiments. We conclude that with a decrease in the surface area to volume ratio, and membrane rigidity, the vesicles prefer tubulation over sickling. The highlight of this work is several important non-axisymmetric RBC and vesicle shapes, which have never been explored in simulations.

Fluctuations and conformational stability of a membrane patch with curvature inducing inclusions

Membranes with curvature inducing inclusions display a range of cooperative phenomena, which can be linked to biomembrane function, e.g. membrane tubulation, vesiculation, softening and spontaneous tension. We investigate how these phenomena are related for a fluctuating, framed membrane through analysis of a descretized membrane model by Monte Carlo simulation techniques. The membrane model is based on a dynamically triangulated surface equipped with non-interacting, up-down symmetry breaking inclusions where only terms coupled linearly to mean-curvature are maintained. We show that the lateral configurational entropy plays a key role for the mechanical properties of the semi-flexible membrane, e.g. a pronounced softening at intermediate inclusion coverages of the membrane and generation of membrane tension. Tensionless framed membranes will remain quasi-flat up to some threshold coverage, where a shape instability occurs with formation of pearling or tubular membranes, which below full coverage is associated with segregation of inclusions between the curved and flat membrane geometries. For inclusions with preference for highly curved membranes the instability appears at dilute inclusion coverages and is accompanied by strong configurational fluctuations.

Thermodynamic analysis of multivalent binding of functionalized nanoparticles to membrane surface reveals the importance of membrane entropy and nanoparticle entropy in adhesion of flexible nanoparticles

We present a quantitative model for multivalent binding of ligand-coated flexible polymeric nanoparticles (NPs) to a flexible membrane expressing receptors. The model is developed using a multiscale computational framework by coupling a continuum field model for the cell membrane with a coarse-grained model for the polymeric NPs. The NP is modeled as a self-avoiding bead-spring polymer chain, and the cell membrane is modeled as a triangulated surface using the dynamically triangulated Monte Carlo method. The nanoparticle binding affinity to a cell surface is mainly determined by the delicate balance between the enthalpic gain due to the multivalent ligand-receptor binding and the entropic penalties of various components including receptor translation, membrane undulation, and NP conformation. We have developed new methods to compute the free energy of binding, which includes these enthalpy and entropy terms. We show that the multivalent interactions between the flexible NP and the cell surface are subject to entropy-enthalpy compensation. Three different entropy contributions, namely, those due to receptor-ligand translation, NP flexibility, and membrane undulations, are all significant, although the first of these terms is the most dominant. However, both NP flexibility and membrane undulations dictate the receptor-ligand translational entropy making the entropy compensation context-specific, i.e., dependent on whether the NP is rigid or flexible, and on the state of the membrane given by the value of membrane tension or its excess area.

A Multi-Scale Approach to Membrane Remodeling Processes

We present a multi-scale simulation procedure to describe membrane-related biological processes that span over a wide range of length scales. At macroscopic length-scale, a membrane is described as a flexible thin film modeled by a dynamic triangulated surface with its spatial conformations governed by an elastic energy containing only a few model parameters. An implicit protein model allows us to include complex effects of membrane-protein interactions in the macroscopic description. The gist of this multi-scale approach is a scheme to calibrate the implicit protein model using finer scale simulation techniques e.g., all atom and coarse grain molecular dynamics. We previously used this approach and properly described the formation of membrane tubular invaginations upon binding of B-subunit of Shiga toxin. Here, we provide a perspective of our multi-scale approach, summarizing its main features and sketching possible routes for future development.

Entry modes of ellipsoidal nanoparticles on a membrane during clathrin-mediated endocytosis

The membrane wrapping and internalization of nanoparticles, such as viruses and drug nanocarriers, through clathrin-mediated endocytosis (CME) are vitally important for intracellular transport. During CME, the shape of the particle plays crucial roles in the determination of particle-membrane interactions, but much of the previous work has been focused on spherical particles. In this work, we develop a stochastic model to study the CME of ellipsoidal particles. In our model, the deformation of the membrane and wrapping of the nanoparticles are driven by the accumulation of clathrin lattices, which is stimulated by the ligand-receptor interactions. Using our model, we systematically investigate the effect of particle shape (ellipsoids with different aspect ratios) on the CME. Our results show three entry modes: tip-first, tilted, and laying-down modes, used by ellipsoidal nanoparticles for internalization depending on the aspect ratio. Certain ellipsoids are able to take multiple entry modes for internalization. Interestingly, the prolate ellipsoid with an aspect ratio of 0.42 can be internalized with a significantly reduced number of ligand-receptor bonds. Particles which can be internalized with fewer bonds are excellent candidates for transcellular drug delivery. Moreover, our results demonstrate that internalization of ellipsoids with intermediate aspect ratios is easier than that of particles with low and high aspect ratios. Our model and simulations provide critical mechanistic insights into CME of ellipsoidal particles, and represent a viable platform for optimal design of nanoparticles for targeted drug delivery applications.

Tubulation pattern of membrane vesicles coated with biofilaments

Narrow membrane tubes are commonly pulled out from the surface of phospholipid vesicles using forces applied either through laser or magnetic tweezers or through the action of processive motor proteins. Recent examples have emerged in which an array of such tubes grows spontaneously from vesicles coated with bioactive cytoskeletal filaments (e.g., FtsZ, microtubule) in the presence GTP or ATP. We show how a soft vesicle deforms as a result of the interplay between its topology, local curvature, and the forces due to filament bundles. We present results from dynamically triangulated Monte Carlo simulations of a closed membrane vesicle coated with a nematic field (the filaments), and we show how the intrinsic curvature of the filaments and their bundling interactions drive membrane tubulation. We predict interesting patterns consisting of a large number of nematic defects that accompany tubulation. A common theme emerges: defect locations on vesicle surfaces are hot spots of membrane deformation activity, which could be useful for vesicle origami. Although our equilibrium model is not applicable to the nonequilibrium shape dynamics exhibited by active microtubule-coated vesicles, we show that some of the features, such as the size-dependent vesicle shape and the number of tubes, can still be understood from our equilibrium model.

Modeling Membrane Curvature Generation due to Membrane⁻Protein Interactions

To alter and adjust the shape of the plasma membrane, cells harness various mechanisms of curvature generation. Many of these curvature generation mechanisms rely on the interactions between peripheral membrane proteins, integral membrane proteins, and lipids in the bilayer membrane. Mathematical and computational modeling of membrane curvature generation has provided great insights into the physics underlying these processes. However, one of the challenges in modeling these processes is identifying the suitable constitutive relationships that describe the membrane free energy including protein distribution and curvature generation capability. Here, we review some of the commonly used continuum elastic membrane models that have been developed for this purpose and discuss their applications. Finally, we address some fundamental challenges that future theoretical methods need to overcome to push the boundaries of current model applications.

Biophysics of membrane curvature remodeling at molecular and mesoscopic lengthscales

At the micron scale, where cell organelles display an amazing complexity in their shape and organization, the physical properties of a biological membrane can be better-understood using continuum models subject to thermal (stochastic) undulations. Yet, the chief orchestrators of these complex and intriguing shapes are a specialized class of membrane associating often peripheral proteins called curvature remodeling proteins (CRPs) that operate at the molecular level through specific protein-lipid interactions. We review multiscale methodologies to model these systems at the molecular as well as at the mesoscopic and cellular scales, and also present a free energy perspective of membrane remodeling through the organization and assembly of CRPs. We discuss the morphological space of nearly planar to highly curved membranes, methods to include thermal fluctuations, and review studies that model such proteins as curvature fields to describe the emergent curved morphologies. We also discuss several mesoscale models applied to a variety of cellular processes, where the phenomenological parameters (such as curvature field strength) are often mapped to models of real systems based on molecular simulations. Much insight can be gained from the calculation of free energies of membranes states with protein fields, which enable accurate mapping of the state and parameter values at which the membrane undergoes morphological transformations such as vesiculation or tubulation. By tuning the strength, anisotropy, and spatial organization of the curvature-field, one can generate a rich array of membrane morphologies that are highly relevant to shapes of several cellular organelles. We review applications of these models to budding of vesicles commonly seen in cellular signaling and trafficking processes such as clathrin mediated endocytosis, sorting by the ESCRT protein complexes, and cellular exocytosis regulated by the exocyst complex. We discuss future prospects where such models can be combined with other models for cytoskeletal assembly, and discuss their role in understanding the effects of cell membrane tension and the mechanics of the extracellular microenvironment on cellular processes.

Stochastic simulations of nanoparticle internalization through transferrin receptor dependent clathrin-mediated endocytosis

BACKGROUND

Receptor dependent clathrin-mediated endocytosis (CME) is one of the most important endocytic pathways for the internalization of bioparticles into cells. During CME, the ligand-receptor interactions, development of clathrin-coated pit (CCP) and membrane evolution all act together to drive the internalization of bioparticles. In this work, we develop a stochastic computational model to investigate the CME based on the Metropolis Monte Carlo simulations.

METHODS

The model is based on the combination of a stochastic particle binding model with a membrane model. The energetic costs of membrane bending, CCP formation and ligand-receptor interactions are systematically linked together.

RESULTS

We implement our model to investigate the effects of particle size, ligand density and membrane stiffness on the overall process of CME from the drug delivery perspectives. Consistent with some experiments, our results show that the intermediate particle size and ligand density favor the particle internalization. Moreover, our results show that it is easier for a particle to enter a cell with softer membrane.

CONCLUSIONS

The model presented here is able to provide mechanistic insights into CME and can be readily modified to include other important factors, such as actins. The predictions from the model will aid in the therapeutic design of intracellular/transcellular drug delivery and antiviral interventions.

Lipid-protein interaction induced domains: Kinetics and conformational changes in multicomponent vesicles

The spatio-temporal organization of proteins and the associated morphological changes in membranes are of importance in cell signaling. Several mechanisms that promote the aggregation of proteins at low cell surface concentrations have been investigated in the past. We show, using Monte Carlo simulations, that the affinity of proteins for specific lipids can hasten their aggregation kinetics. The lipid membrane is modeled as a dynamically triangulated surface with the proteins defined as in-plane fields at the vertices. We show that, even at low protein concentrations, strong lipid-protein interactions can result in large protein clusters indicating a route to lipid mediated signal amplification. At high protein concentrations, the domains form buds similar to that seen in lipid-lipid interaction induced phase separation. Protein interaction induced domain budding is suppressed when proteins act as anisotropic inclusions and exhibit nematic orientational order. The kinetics of protein clustering and resulting conformational changes are shown to be significantly different for the isotropic and anisotropic curvature inducing proteins.

The tension of framed membranes from computer simulations

We have analyzed the behavior of a randomly triangulated, self-avoiding surface model of a flexible, fluid membrane subject to a circular boundary by Wang-Landau Monte Carlo computer simulation techniques. The dependence of the canonical free energy and frame tension on the frame area is obtained for flexible membranes. It is shown that for low bending rigidities the framed membrane is only stable above a threshold tension, suggesting a discontinuous transition from the collapsed (branched polymer) state to a finite tension extended state. In a tension range above this threshold tension the membranes display power-law characteristics for the equation of state, while higher tension levels includes both an extended linear (elastic) as well as a highly non-linear stretching regime. For semi-flexible membranes a transition from extended to buckled conformations takes place at negative frame tensions. Our analysis indicates that at zero frame tension the crumpling transition of fluid membranes show characteristics of both critical behavior and a discontinuous transition at low bending rigidities.

Thermal properties of vortices on curved surfaces

We use Monte Carlo simulations to study the finite temperature behavior of vortices in the XY model for tangent vector order on curved backgrounds. Contrary to naive expectations, we show that the underlying geometry does not affect the proliferation of vortices with temperature respect to what is observed on a flat surface. Long-range order in these systems is analyzed by using two-point correlation functions. As expected, in the case of slightly curved substrates these correlations behave similarly to the plane. However, for high curvatures, the presence of geometry-induced unbounded vortices at low temperatures produces the rapid decay of correlations and an apparent lack of long-range order. Our results shed light on the finite-temperature physics of soft-matter systems and anisotropic magnets deposited on curved substrates.

Computing Curvature Sensitivity of Biomolecules in Membranes by Simulated Buckling

Membrane curvature sensing, where the binding free energies of membrane-associated molecules depend on the local membrane curvature, is a key factor to modulate and maintain the shape and organization of cell membranes. However, the microscopic mechanisms are not well understood, partly due to absence of efficient simulation methods. Here, we describe a method to compute the curvature dependence of the binding free energy of a membrane-associated probe molecule that interacts with a buckled membrane, which has been created by lateral compression of a flat bilayer patch. This buckling approach samples a wide range of curvatures in a single simulation, and anisotropic effects can be extracted from the orientation statistics. We develop an efficient and robust algorithm to extract the motion of the probe along the buckled membrane surface, and evaluate its numerical properties by extensive sampling of three coarse-grained model systems: local lipid density in a curved environment for single-component bilayers, curvature preferences of individual lipids in two-component membranes, and curvature sensing by a homotrimeric transmembrane protein. The method can be used to complement experimental data from curvature partition assays and provides additional insight into mesoscopic theories and molecular mechanisms for curvature sensing.

Excess area dependent scaling behavior of nano-sized membrane tethers

Thermal fluctuations in cell membranes manifest as an excess area ([Formula: see text]) which governs a multitude of physical process at the sub-micron scale. We present a theoretical framework, based on an in silico tether pulling method, which may be used to reliably estimate [Formula: see text] in live cells. We perform our simulations in two different thermodynamic ensembles: (i) the constant projected area and (ii) the constant frame tension ensembles and show the equivalence of our results in the two. The tether forces estimated from our simulations compare well with our experimental measurements for tethers extracted from ruptured GUVs and HeLa cells. We demonstrate the significance and validity of our method by showing that all our calculations performed in the initial tether formation regime (i.e. when the length of the tether is comparable to its radius) along with experiments of tether extraction in 15 different cell types collapse onto two unified scaling relationships mapping tether force, tether radius, bending stiffness κ, and membrane tension σ. We show that [Formula: see text] is an important determinant of the radius of the extracted tether, which is equal to the characteristic length [Formula: see text] for [Formula: see text], and is equal to [Formula: see text] for [Formula: see text]. We also find that the estimated excess area follows a linear scaling behavior that only depends on the true value of [Formula: see text] for the membrane, based on which we propose a self-consistent technique to estimate the range of excess membrane areas in a cell.

On the role of external force of actin filaments in the formation of tubular protrusions of closed membrane shapes with anisotropic membrane components

Biological membranes are composed of different components and there is no a priori reason to assume that all components are isotropic. It was previously shown that the anisotropic properties of membrane components may explain the stability of membrane tubular protrusions even without the application of external force. Our theoretical study focuses on the role of anisotropic membrane components in the stability of membrane tubular structures generated or stabilized by actin filaments. We show that the growth of the actin cytoskeleton inside the vesicle can induce the partial lateral segregation of different membrane components. The entropy of mixing of membrane components hinders the total lateral segregation of the anisotropic and isotropic membrane components. Self-assembled aggregates formed by anisotropic membrane components facilitate the growth of long membrane tubular protrusions. Protrusive force generated by actin filaments favors strong segregation of membrane components by diminishing the opposing effect of mixing entropy.

Numerical insights into the phase diagram of p-atic membranes with spherical topology

The properties of self-avoiding p-atic membranes restricted to spherical topology have been studied by Monte Carlo simulations of a triangulated random surface model. Spherically shaped p-atic membranes undergo a Kosterlitz-Thouless transition as expected with topology induced mutually repelling disclinations of the p-atic ordered phase. For flexible membranes the phase behaviour bears some resemblance to the spherically shaped case with a p-atic disordered crumpled phase and p-atic ordered, conformationally ordered (crinkled) phase separated by a KT-like transition with proliferation of disclinations. We confirm the proposed buckling of disclinations in the p-atic ordered phase, while the expected associated disordering (crumpling) transition at low bending rigidities is absent in the phase diagram.

Anisotropic Membrane Curvature Sensing by Amphipathic Peptides

Many proteins and peptides have an intrinsic capacity to sense and induce membrane curvature, and play crucial roles for organizing and remodeling cell membranes. However, the molecular driving forces behind these processes are not well understood. Here, we describe an approach to study curvature sensing by simulating the interactions of single molecules with a buckled lipid bilayer. We analyze three amphipathic antimicrobial peptides, a class of membrane-associated molecules that specifically target and destabilize bacterial membranes, and find qualitatively different sensing characteristics that would be difficult to resolve with other methods. Our findings provide evidence for direction-dependent curvature sensing mechanisms in amphipathic peptides and challenge existing theories of hydrophobic insertion. The buckling approach is generally applicable to a wide range of curvature-sensing molecules, and our results provide strong motivation to develop new experimental methods to track position and orientation of membrane proteins.

Membrane invagination induced by Shiga toxin B-subunit: from molecular structure to tube formation

The bacterial Shiga toxin is composed of an enzymatically active A-subunit, and a receptor-binding homopentameric B-subunit (STxB) that mediates intracellular toxin trafficking. Upon STxB-mediated binding to the glycolipid globotriaosylceramide (Gb3) at the plasma membrane of target cells, Shiga toxin is internalized by clathrin-dependent and independent endocytosis. The formation of tubular membrane invaginations is an essential step in the clathrin-independent STxB uptake process. However, the mechanism by which STxB induces these invaginations has remained unclear. Using a combination of all-atom molecular dynamics and Monte Carlo simulations we show that the molecular architecture of STxB enables the following sequence of events: the Gb3 binding sites on STxB are arranged such that tight avidity-based binding results in a small increment of local curvature. Membrane-mediated clustering of several toxin molecules then creates a tubular membrane invagination that drives toxin entry into the cell. This mechanism requires: (1) a precise molecular architecture of the STxB binding sites; (2) a fluid bilayer in order for the tubular invagination to form. Although, STxB binding to the membrane requires specific interactions with Gb3 lipids, our study points to a generic molecular design principle for clathrin-independent endocytosis of nanoparticles.

Thermodynamic free energy methods to investigate shape transitions in bilayer membranes

The conformational free energy landscape of a system is a fundamental thermodynamic quantity of importance particularly in the study of soft matter and biological systems, in which the entropic contributions play a dominant role. While computational methods to delineate the free energy landscape are routinely used to analyze the relative stability of conformational states, to determine phase boundaries, and to compute ligand-receptor binding energies its use in problems involving the cell membrane is limited. Here, we present an overview of four different free energy methods to study morphological transitions in bilayer membranes, induced either by the action of curvature remodeling proteins or due to the application of external forces. Using a triangulated surface as a model for the cell membrane and using the framework of dynamical triangulation Monte Carlo, we have focused on the methods of Widom insertion, thermodynamic integration, Bennett acceptance scheme, and umbrella sampling and weighted histogram analysis. We have demonstrated how these methods can be employed in a variety of problems involving the cell membrane. Specifically, we have shown that the chemical potential, computed using Widom insertion, and the relative free energies, computed using thermodynamic integration and Bennett acceptance method, are excellent measures to study the transition from curvature sensing to curvature inducing behavior of membrane associated proteins. The umbrella sampling and WHAM analysis has been used to study the thermodynamics of tether formation in cell membranes and the quantitative predictions of the computational model are in excellent agreement with experimental measurements. Furthermore, we also present a method based on WHAM and thermodynamic integration to handle problems related to end-point-catastrophe that are common in most free energy methods.

Biophysically inspired model for functionalized nanocarrier adhesion to cell surface: roles of protein expression and mechanical factors

In order to achieve selective targeting of affinity-ligand coated nanoparticles to the target tissue, it is essential to understand the key mechanisms that govern their capture by the target cell. Next-generation pharmacokinetic (PK) models that systematically account for proteomic and mechanical factors can accelerate the design, validation and translation of targeted nanocarriers (NCs) in the clinic. Towards this objective, we have developed a computational model to delineate the roles played by target protein expression and mechanical factors of the target cell membrane in determining the avidity of functionalized NCs to live cells. Model results show quantitative agreement with in vivo experiments when specific and non-specific contributions to NC binding are taken into account. The specific contributions are accounted for through extensive simulations of multivalent receptor-ligand interactions, membrane mechanics and entropic factors such as membrane undulations and receptor translation. The computed NC avidity is strongly dependent on ligand density, receptor expression, bending mechanics of the target cell membrane, as well as entropic factors associated with the membrane and the receptor motion. Our computational model can predict the in vivo targeting levels of the intracellular adhesion molecule-1 (ICAM1)-coated NCs targeted to the lung, heart, kidney, liver and spleen of mouse, when the contributions due to endothelial capture are accounted for. The effect of other cells (such as monocytes, etc.) do not improve the model predictions at steady state. We demonstrate the predictive utility of our model by predicting partitioning coefficients of functionalized NCs in mice and human tissues and report the statistical accuracy of our model predictions under different scenarios.

Curvature inducing macroion condensation driven shape changes of fluid vesicles

We study the effect of curvature inducing macroion condensation on the shapes of charged deformable fluid interfaces using dynamically triangulated Monte Carlo simulations. In the weak electrostatic coupling regime, surface charges are weakly screened and the conformations of a vesicle, with fixed spherical topology, depend on the charge-charge interaction on the surface. While in the strong coupling regime, condensation driven curvature induction plays a dominant role in determining the conformations of these surfaces. Condensation itself is observed to be dependent on the induced curvature, with larger induced curvatures favoring increased condensation. We show that both curvature generation and curvature sensing, induced by the interplay of electrostatics and curvature energy, contribute to determination of the vesicle configurations.

Application of a free-energy-landscape approach to study tension-dependent bilayer tubulation mediated by curvature-inducing proteins

We investigate the phenomenon of protein-induced tubulation of lipid bilayer membranes within a continuum framework using Monte Carlo simulations coupled with the Widom insertion technique to compute excess chemical potentials. Tubular morphologies are spontaneously formed when the density and the curvature-field strength of the membrane-bound proteins exceed their respective thresholds and this transition is marked by a sharp drop in the excess chemical potential. We find that the planar to tubular transition can be described by a micellar model and that the corresponding free-energy barrier increases with an increase in the curvature-field strength (i.e., of protein-membrane interactions) and also with an increase in membrane tension.

Monte Carlo simulations of fluid vesicles

Lipid vesicles are closed two dimensional fluid surfaces that are studied extensively as model systems for understanding the physical properties of biological membranes. Here we review the recent developments in the Monte Carlo techniques for simulating fluid vesicles and discuss some of their applications. The technique, which treats the membrane as an elastic sheet, is most suitable for the study of large scale conformations of membranes. The model can be used to study vesicles with fixed and varying topologies. Here we focus on the case of multi-component membranes with the local lipid and protein composition coupled to the membrane curvature leading to a variety of shapes. The phase diagram is more intriguing in the case of fluid vesicles having an in-plane orientational order that induce anisotropic directional curvatures. Methods to explore the steady state morphological structures due to active flux of materials have also been described in the context of Monte Carlo simulations.

Organelle morphogenesis by active membrane remodeling

Intracellular organelles are subject to a steady flux of lipids and proteins through active, energy consuming transport processes. Active fission and fusion are promoted by GTPases, e.g., Arf-Coatamer and the Rab-Snare complexes, which both sense and generate local membrane curvature. Here we investigate, through Dynamical Triangulation Monte Carlo simulations, the role that these active processes play in determining the morphology and composition segregation in closed membranes. We find that the steady state shapes obtained as a result of such active processes, bear a striking resemblance to the ramified morphologies of organelles in vivo, pointing to the relevance of nonequilibrium fission-fusion in organelle morphogenesis.

Phenomenology based multiscale models as tools to understand cell membrane and organelle morphologies

An intriguing question in cell biology is "how do cells regulate their shape?" It is commonly believed that the observed cellular morphologies are a result of the complex interaction among the lipid molecules (constituting the cell membrane), and with a number of other macromolecules, such as proteins. It is also believed that the common biophysical processes essential for the functioning of a cell also play an important role in cellular morphogenesis. At the cellular scale-where typical dimensions are in the order of micrometers-the effects arising from the molecular scale can either be modeled as equilibrium or non-equilibrium processes. In this chapter, we discuss the dynamically triangulated Monte Carlo technique to model and simulate membrane morphologies at the cellular scale, which in turn can be used to investigate several questions related to shape regulation in cells. In particular, we focus on two specific problems within the framework of isotropic and anisotropic elasticity theories: namely, (i) the origin of complex, physiologically relevant, membrane shapes due to the interaction of the membrane with curvature remodeling proteins, and (ii) the genesis of steady state cellular shapes due to the action of non-equilibrium forces that are generated by the fission and fusion of transport vesicles and by the binding and unbinding of proteins from the parent membrane.

Terasaki spiral ramps in the rough endoplasmic reticulum

We present a model describing the morphology as well as the assembly of "Terasaki ramps," the recently discovered helicoidal connections linking adjacent sheets of the rough endoplasmic reticulum (ER). The fundamental unit is a localized symmetric double-ramped "parking garage" formed by two separated gently pitched, approximately helicoidal, ramps of opposite chiralities. This geometry is stabilized by a short-range repulsive interaction between ramps associated with bending energy which opposes the long-range attraction associated with tension. The ramp inner boundaries are themselves stabilized by the condensation of membrane-shaping proteins along their length. A mechanism for parking garage self-assembly is proposed involving the nucleation of dipoles at the center of tubular three-way junctions within the smooth ER. Our predictions are compared with the experimental data.

Mesoscale computational studies of membrane bilayer remodeling by curvature-inducing proteins

Biological membranes constitute boundaries of cells and cell organelles. These membranes are soft fluid interfaces whose thermodynamic states are dictated by bending moduli, induced curvature fields, and thermal fluctuations. Recently, there has been a flood of experimental evidence highlighting active roles for these structures in many cellular processes ranging from trafficking of cargo to cell motility. It is believed that the local membrane curvature, which is continuously altered due to its interactions with myriad proteins and other macromolecules attached to its surface, holds the key to the emergent functionality in these cellular processes. Mechanisms at the atomic scale are dictated by protein-lipid interaction strength, lipid composition, lipid distribution in the vicinity of the protein, shape and amino acid composition of the protein, and its amino acid contents. The specificity of molecular interactions together with the cooperativity of multiple proteins induce and stabilize complex membrane shapes at the mesoscale. These shapes span a wide spectrum ranging from the spherical plasma membrane to the complex cisternae of the Golgi apparatus. Mapping the relation between the protein-induced deformations at the molecular scale and the resulting mesoscale morphologies is key to bridging cellular experiments across the various length scales. In this review, we focus on the theoretical and computational methods used to understand the phenomenology underlying protein-driven membrane remodeling. Interactions at the molecular scale can be computationally probed by all atom and coarse grained molecular dynamics (MD, CGMD), as well as dissipative particle dynamics (DPD) simulations, which we only describe in passing. We choose to focus on several continuum approaches extending the Canham - Helfrich elastic energy model for membranes to include the effect of curvature-inducing proteins and explore the conformational phase space of such systems. In this description, the protein is expressed in the form of a spontaneous curvature field. The approaches include field theoretical methods limited to the small deformation regime, triangulated surfaces and particle-based computational models to investigate the large-deformation regimes observed in the natural state of many biological membranes. Applications of these methods to understand the properties of biological membranes in homogeneous and inhomogeneous environments of proteins, whose underlying curvature fields are either isotropic or anisotropic, are discussed. The diversity in the curvature fields elicits a rich variety of morphological states, including tubes, discs, branched tubes, and caveola. Mapping the thermodynamic stability of these states as a function of tuning parameters such as concentration and strength of curvature induction of the proteins is discussed. The relative stabilities of these self-organized shapes are examined through free-energy calculations. The suite of methods discussed here can be tailored to applications in specific cellular settings such as endocytosis during cargo trafficking and tubulation of filopodial structures in migrating cells, which makes these methods a powerful complement to experimental studies.

Defining the free-energy landscape of curvature-inducing proteins on membrane bilayers

Curvature-sensing and curvature-remodeling proteins, such as Amphiphysin, Epsin, and Exo70, are known to reshape cell membranes, and this remodeling event is essential for key biophysical processes such as tubulation, exocytosis, and endocytosis. Curvature-inducing proteins can act as curvature sensors; they aggregate to membrane regions matching their intrinsic curvature; as well as induce curvature in cell membranes to stabilize emergent high curvature, nonspherical, structures such as tubules, discs, and caveolae. A definitive understanding of the interplay between protein recruitment and migration, the evolution of membrane curvature, and membrane morphological transitions is emerging but remains incomplete. Here, within a continuum framework and using the machinery of Monte Carlo simulations, we introduce and compare three free-energy methods to delineate the free-energy landscape of curvature-inducing proteins on bilayer membranes. We demonstrate the utility of the Widom test particle (or field) insertion methodology in computing the excess chemical potentials associated with curvature-inducing proteins on the membrane-in particular, we use this method to track the onset of morphological transitions in the membrane at elevated protein densities. We validate this approach by comparing the results from the Widom method with those of thermodynamic integration and Bennett acceptance ratio methods. Furthermore, the predictions from the Widom method have been tested against analytical calculations of the excess chemical potential at infinite dilution. Our results are useful in precisely quantifying the free-energy landscape, and also in determining the phase boundaries associated with curvature-induction, curvature-sensing, and morphological transitions. This approach can be extended to studies exploring the role of thermal fluctuations and other external (control) variables, such as membrane excess area, in shaping curvature-mediated interactions on bilayer membranes.

Anisotropic spontaneous curvatures in lipid membranes

Symmetry restrictions due to fluidity require the strain energy in the Helfrich theory of lipid membranes to be locally isotropic in nature. Although this framework is suitable for modeling the interaction of membranes with proteins that generate spherical curvature such as clathrin, there are other important membrane-bending proteins such as BIN-amphiphysin-Rvs proteins that form a cylindrical coat with different curvatures in the longitudinal and the circumferential directions. In this work, we present a detailed mathematical treatment of the theory of lipid membranes incorporating anisotropic spontaneous curvatures. We derive the associated Euler-Lagrange equations and the edge conditions in a generalized setting that allows spatial heterogeneities in the properties of the membrane-protein system. We employ this theory to model the constriction of a membrane tubule by a cylindrical scaffold. In particular, we highlight the role of the equilibrium equation in the tangential plane in regulating the spatial variation of the surface tension field.

Exo70 generates membrane curvature for morphogenesis and cell migration

Dynamic shape changes of the plasma membrane are fundamental to many processes, ranging from morphogenesis and cell migration to phagocytosis and viral propagation. Here, we demonstrate that Exo70, a component of the exocyst complex, induces tubular membrane invaginations toward the lumen of synthetic vesicles in vitro and generates protrusions on the surface of cells. Biochemical analyses using Exo70 mutants and independent molecular dynamics simulations based on Exo70 structure demonstrate that Exo70 generates negative membrane curvature through an oligomerization-based mechanism. In cells, the membrane-deformation function of Exo70 is required for protrusion formation and directional cell migration. Exo70 thus represents a membrane-bending protein that may couple actin dynamics and plasma membrane remodeling for morphogenesis.

Membrane-mediated aggregation of curvature-inducing nematogens and membrane tubulation

The shapes of cell membranes are largely regulated by membrane-associated, curvature-active proteins. Herein, we use a numerical model of the membrane, recently developed by us, with elongated membrane inclusions possessing spontaneous directional curvatures that could be different along, and perpendicular to, the membrane's long axis. We show that, due to membrane-mediated interactions, these curvature-inducing membrane-nematogens can aggregate spontaneously, even at low concentrations, and change the local shape of the membrane. We demonstrate that for a large group of such inclusions, where the two spontaneous curvatures have equal sign, the tubular conformation and sometimes the sheet conformation of the membrane are the common equilibrium shapes. We elucidate the factors necessary for the formation of these protein lattices. Furthermore, the elastic properties of the tubes, such as their compressional stiffness and persistence length, are calculated. Finally, we discuss the possible role of nematic disclination in capping and branching of the tubular membranes.

Multiscale perspectives of virus entry via endocytosis

Most viruses take advantage of endocytic pathways to gain entry into host cells and initiate infections. Understanding of virus entry via endocytosis is critically important for the design of antiviral strategies. Virus entry via endocytosis is a complex process involving hundreds of cellular proteins. The entire process is dictated by events occurring at multiple time and length scales. In this review, we discuss and evaluate the available means to investigate virus endocytic entry, from both experimental and theoretical/numerical modeling fronts, and highlight the importance of multiscale features. The complexity of the process requires investigations at a systems biology level, which involves the combination of different experimental approaches, the collaboration of experimentalists and theorists across different disciplines, and the development of novel multiscale models.

Vesicles and vesicle fusion: coarse-grained simulations

Biological cells are highly dynamic, and continually move material around their own volume and between their interior and exterior. Much of this transport encapsulates the material inside phospholipid vesicles that shuttle to and from, fusing with, and budding from, other membranes. A feature of vesicles that is crucial for this transport is their ability to fuse to target membranes and release their contents to the distal side. In industry, some personal care products contain vesicles to help transport reagents across the skin, and research on drug formulation shows that packaging active compounds inside vesicles delays their clearance from the blood stream. In this chapter, we survey the biological role and physicochemical properties of phospholipids, and describe progress in coarse-grained simulations of vesicles and vesicle fusion. Because coarse-grained simulations retain only those molecular details that are thought to influence the large-scale processes of interest, they act as a model embodying our current understanding. Comparing the predictions of these models with experiments reveals the importance of the retained microscopic details and also the deficiencies that can suggest missing details, thereby furthering our understanding of the complex dynamic world of vesicles.

Mesoscale simulations of curvature-inducing protein partitioning on lipid bilayer membranes in the presence of mean curvature fields

The membrane-surface migration of curvature-inducing proteins in response to membrane curvature gradients has been investigated using Monte Carlo simulations of a curvilinear membrane model based on the Helfrich Hamiltonian. Consistent with theoretical and experimental data, we find the proteins that generate curvature can also sense the background membrane curvature, wherein they preferentially partition to the high curvature regions. The partitioning strength depends linearly on local membrane curvature and the slope (or the coupling constant) of the partitioning probability versus mean curvature depends on the membrane bending rigidity and instantaneous curvature field caused by different proteins. Our simulation study allows us to quantitatively characterize and identify the important factors affecting the coupling constant (slope), which may be difficult to determine in experiments. Furthermore, the membrane model is used to study budding of vesicles where it is found that in order to stabilize a mature vesicle with a stable 'neck-region' (or stable membrane overhangs), the area (extent) of the intrinsic curvature region needs to exceed a threshold-critical value. The migration and partitioning of curvature-inducing proteins in a budding vesicle with a stable neck (with a characteristic negative value of the Gaussian curvature) is investigated.