We present a method for simulating fluid vesicles with in-plane orientational ordering. The method involves computation of local curvature tensor and parallel transport of the orientational field on a randomly triangulated surface. It is shown that the model reproduces the known equilibrium conformation of fluid membranes and work well for a large range of bending rigidities. Introduction of nematic ordering leads to stiffening of the membrane. Nematic ordering can also result in anisotropic rigidity on the surface leading to formation of membrane tubes.

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 set...

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 towards 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.

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.

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-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.

Inflamed organs display marked spatial heterogeneity of inflammation, with patches of inflamed tissue adjacent to healthy tissue. To investigate how nanocarriers (NCs) distribute between such patches, we created a mouse model that recapitulates the spatial heterogeneity of the inflammatory lung disease ARDS. NCs targeting the epitope PECAM strongly accumulated in the lungs, but were shunted away from inflamed lung regions due to hypoxic vasoconstriction (HVC). In contrast, ICAM-targeted NCs, which had lower whole-lung uptake than PECAM/NCs in inflamed lungs, displayed markedly higher NC levels in inflamed regions than PECAM/NCs, due to increased regional ICAM. Regional HVC, epitope expression, and capillary leak were sufficient to predict intra-organ of distribution of NCs, antibodies, and drugs. Importantly, these effects were not observable with traditional spatially-uniform models of ARDS, nor when examining only whole-organ uptake. This study underscores how examining NCs’ intra-organ distribution in spatially heterogeneous animal models can guide rational NC design.

We study the equilibrium shapes of vesicles, with an in-plane nematic order, using a Monte-Carlo scheme and show that highly curved shapes, like tubes and discs, with a striking similarity to the structures engendered by certain curvature sensing peripheral membrane proteins, can be spontaneously generated by anisotropic directional curvature with nematic disclinations playing and important role. We show that the coupling between nematic order and local curvature could lead to like defects moving towards each other and unlike defects moving away, in turn leading to tube formation. Thermally induced defect pair production lead to branched tubular structures. It is also shown that helical arrangement of the membrane tubes, with nematic field spiraling around it, is a dominant soft mode of the system.