Delineating the shape of COat Protein complex-II coated membrane bud.

Curvature-generating proteins that direct membrane trafficking assemble on the surface of lipid bilayers to bud transport intermediates, which move protein and lipid cargoes from one cellular compartment to another. However, it remains unclear what controls the overall shape of the membrane bud once curvature induction has begun. In vitro experiments showed that excessive concentrations of the COPII protein Sar1 promoted the formation of membrane tubules from synthetic vesicles, while COPII-coated transport intermediates in cells are generally more spherical or lobed in shape. To understand the origin of these morphological differences, we employ atomistic, coarse-grained (CG), and continuum mesoscopic simulations of membranes in the presence of multiple curvature-generating proteins. We first characterize the membrane-bending ability of amphipathic peptides derived from the amino terminus of Sar1, as a function of interpeptide angle and concentration using an atomistic bicelle simulation protocol. Then, we employ CG simulations to reveal that Sec23 and Sec24 control the relative spacing between Sar1 protomers and form the inner-coat unit through an attachment with Sar1. Finally, using dynamical triangulated surface simulations based on the Helfrich Hamiltonian, we demonstrate that the uniform distribution of spacer molecules among curvature-generating proteins is crucial to the spherical budding of the membrane. Overall, our analyses suggest a new role for Sec23, Sec24, and cargo proteins in COPII-mediated membrane budding process in which they act as spacers to preserve a dispersed arrangement of Sar1 protomers and help determine the overall shape of the membrane bud.

Modeling the mechanochemical feedback for membrane-protein interactions using a continuum mesh model.

The Helfrich free energy is widely used to model the generation of membrane curvature due to different physical and chemical components. The governing equations resulting from the energy minimization procedure are a system of coupled higher order partial differential equations. Simulations of membrane deformation for obtaining quantitative comparisons against experimental observations require computational schemes that will allow us to solve these equations without restrictions to axisymmetric coordinates. Here, we describe one such tool that we developed in our group based on discrete differential geometry to solve these equations along with examples.

Structure, interfacial film properties, and thermal fluctuations of microemulsions as seen by scattering experiments.

The physics of microemulsions and in particular Dominique Langevin's contributions to the understanding of microemulsion structure and bending properties using scattering techniques are reviewed. Among the many methods used by her and her co-workers, we particularly emphasize optical techniques and small angle neutron scattering (SANS), but also neutron spin echo spectroscopy (NSE). The review is then extended to more recent studies of properties of microemulsions close to surfaces, using reflectometry and grazing-incidence small angle neutron scattering (GISANS).

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.

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.

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.