Landscape of finite-temperature equilibrium behaviour of curvature-inducing proteins on a bilayer membrane explored using a linearized elastic free energy model

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.

Multiscale computational models in physical systems biology of intracellular trafficking

In intracellular trafficking, a definitive understanding of the interplay between protein binding and membrane morphology remains incomplete. The authors describe a computational approach by integrating coarse-grained molecular dynamics (CGMD) simulations with continuum Monte Carlo (CM) simulations of the membrane to study protein-membrane interactions and the ensuing membrane curvature. They relate the curvature field strength discerned from the molecular level to its effect at the cellular length-scale. They perform thermodynamic integration on the CM model to describe the free energy landscape of vesiculation in clathrin-mediated endocytosis. The method presented here delineates membrane morphologies and maps out the free energy changes associated with membrane remodeling due to varying coat sizes, coat curvature strengths, membrane bending rigidities, and tensions; furthermore several constraints on mechanisms underlying clathrin-mediated endocytosis have also been identified, Their CGMD simulations have revealed the importance of PIP2 for stable binding of proteins essential for curvature induction in the bilayer and have provided a molecular basis for the positive curvature induction by the epsin N-terminal homology (EIMTH) domain. Calculation of the free energy landscape for vesicle budding has identified the critical size and curvature strength of a clathrin coat required for nucleation and stabilisation of a mature vesicle.

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.

Nanocarrier Hydrodynamics and Binding in Targeted Drug Delivery: Challenges in Numerical Modeling and Experimental Validation

This review discusses current progress and future challenges in the numerical modeling of targeted drug delivery using functionalized nanocarriers (NC). Antibody coated nanocarriers of various size and shapes, also called functionalized nanocarriers, are designed to be injected in the vasculature, whereby they undergo translational and rotational motion governed by hydrodynamic interaction with blood particulates as well as adhesive interactions mediated by the surface antibody binding to target antigens/receptors on cell surfaces. We review current multiscale modeling approaches rooted in computational fluid dynamics and nonequilibrium statistical mechanics to accurately resolve fluid, thermal, as well as adhesive interactions governing nanocarrier motion and their binding to endothelial cells lining the vasculature. We also outline current challenges and unresolved issues surrounding the modeling methods. Experimental approaches in pharmacology and bioengineering are discussed briefly from the perspective of model validation.

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.

Systems biology and physical biology of clathrin-mediated endocytosis

In this review, we describe the application of experimental data and modeling of intracellular endocytic trafficking mechanisms with a focus on the process of clathrin-mediated endocytosis. A detailed parts-list for the protein-protein interactions in clathrin-mediated endocytosis has been available for some time. However, recent experimental, theoretical, and computational tools have proved to be critical in establishing a sequence of events, cooperative dynamics, and energetics of the intracellular process. On the experimental front, total internal reflection fluorescence microscopy, photo-activated localization microscopy, and spinning-disk confocal microscopy have focused on assembly and patterning of endocytic proteins at the membrane, while on the theory front, minimal theoretical models for clathrin nucleation, biophysical models for membrane curvature and bending elasticity, as well as methods from computational structural and systems biology, have proved insightful in describing membrane topologies, curvature mechanisms, and energetics.

Cancer Cell: Linking Oncogenic Signaling to Molecular Structure

A multiscale strategy is presented for constructing models of intracellular signaling networks in which the oncogenic behavior of the network is encoded through alternate parameterization of the kinetic and structural properties of mutant oncoproteins. The approach uses molecular dynamics and docking simulations to quantify altered topologies of interactions as well as to provide the missing parameters for network models of both wild-type and oncogenic signaling. Through simulation of the resulting signaling networks, the global behavior of these networks may then be compared and functional roles may be assigned to the mutant oncoproteins. An example of this approach is presented in which structural alterations found in a mutant form of the epidermal growth factor receptor are represented as kinetic perturbations in a model of growth factor signaling. Based on network parameters estimated from molecular-level simulations, simulations at the network level show that small perturbations in molecular structure can lead to profoundly altered cellular phenotype.

Minimal mesoscale model for protein-mediated vesiculation in clathrin-dependent endocytosis

In eukaryotic cells, the internalization of extracellular cargo via the endocytic machinery is an important regulatory process required for many essential cellular functions. The role of cooperative protein-protein and protein-membrane interactions in the ubiquitous endocytic pathway in mammalian cells, namely the clathrin-dependent endocytosis, remains unresolved. We employ the Helfrich membrane Hamiltonian together with surface evolution methodology to address how the shapes and energetics of vesicular-bud formation in a planar membrane are stabilized by presence of the clathrin-coat assembly. Our results identify a unique dual role for the tubulating protein epsin: multiple epsins localized spatially and orientationally collectively play the role of a curvature inducing capsid; in addition, epsin serves the role of an adapter in binding the clathrin coat to the membrane. Our results also suggest an important role for the clathrin lattice, namely in the spatial- and orientational-templating of epsins. We suggest that there exists a critical size of the coat above which a vesicular bud with a constricted neck resembling a mature vesicle is stabilized. Based on the observed strong dependence of the vesicle diameter on the bending rigidity, we suggest that the variability in bending stiffness due to variations in membrane composition with cell type can explain the experimentally observed variability on the size of clathrin-coated vesicles, which typically range 50-100 nm. Our model also provides estimates for the number of epsins involved in stabilizing a coated vesicle, and without any direct fitting reproduces the experimentally observed shapes of vesicular intermediates as well as their probability distributions quantitatively, in wildtype as well as CLAP IgG injected neuronal cell experiments. We have presented a minimal mesoscale model which quantitatively explains several experimental observations on the process of vesicle nucleation induced by the clathrin-coated assembly prior to vesicle scission in clathrin dependent endocytosis.

Minimal mesoscale model for protein-mediated vesiculation in clathrin-dependent endocytosis

In eukaryotic cells, the internalization of extracellular cargo via the endocytic machinery is an important regulatory process required for many essential cellular functions. The role of cooperative protein-protein and protein-membrane interactions in the ubiquitous endocytic pathway in mammalian cells, namely the clathrin-dependent endocytosis, remains unresolved. We employ the Helfrich membrane Hamiltonian together with surface evolution methodology to address how the shapes and energetics of vesicular-bud formation in a planar membrane are stabilized by presence of the clathrin-coat assembly. Our results identify a unique dual role for the tubulating protein epsin: multiple epsins localized spatially and orientationally collectively play the role of a curvature inducing capsid; in addition, epsin serves the role of an adapter in binding the clathrin coat to the membrane. Our results also suggest an important role for the clathrin lattice, namely in the spatial- and orientational-templating of epsins. We suggest that there exists a critical size of the coat above which a vesicular bud with a constricted neck resembling a mature vesicle is stabilized. Based on the observed strong dependence of the vesicle diameter on the bending rigidity, we suggest that the variability in bending stiffness due to variations in membrane composition with cell type can explain the experimentally observed variability on the size of clathrin-coated vesicles, which typically range 50-100 nm. Our model also provides estimates for the number of epsins involved in stabilizing a coated vesicle, and without any direct fitting reproduces the experimentally observed shapes of vesicular intermediates as well as their probability distributions quantitatively, in wildtype as well as CLAP IgG injected neuronal cell experiments. We have presented a minimal mesoscale model which quantitatively explains several experimental observations on the process of vesicle nucleation induced by the clathrin-coated assembly prior to vesicle scission in clathrin dependent endocytosis.

Calculation of free energies in fluid membranes subject to heterogeneous curvature fields

We present a computational methodology for incorporating thermal effects and calculating relative free energies for elastic fluid membranes subject to spatially dependent intrinsic curvature fields using the method of thermodynamic integration. Based on a simple model for the intrinsic curvature imposed only in a localized region of the membrane, we employ thermodynamic integration to calculate the free-energy change as a function of increasing strength of the intrinsic curvature field and a thermodynamic cycle to compute free-energy changes for different sizes of the localized region. By explicitly computing the free-energy changes and by quantifying the loss of entropy accompanied with increasing membrane deformation, we show that the membrane stiffness increases with increasing intrinsic field, thereby, renormalizing the membrane bending rigidity. The second main conclusion of this work is that the entropy of the membrane decreases with increasing size of the localized region subject to the curvature field. Our results help to quantify the free-energy change when a planar membrane deforms under the influence of curvature-inducing proteins at a finite temperature.