Protein-induced membrane curvature changes membrane tension

Adsorption of proteins onto membranes can alter the local membrane curvature. This phenomenon has been observed in biological processes such as endocytosis, tubulation and vesiculation. However, it is not clear how the local surface properties of the membrane, such as membrane tension, change in response to protein adsorption. In this paper, we show that the classical elastic model of lipid membranes cannot account for simultaneous changes in shape and membrane tension due to protein adsorption in a local region, and a viscous-elastic formulation is necessary to fully describe the system. Therefore, we develop a viscous-elastic model for inhomogeneous membranes of the Helfrich type. Using the new viscous-elastic model, we find that the lipids flow to accommodate changes in membrane curvature during protein adsorption. We show that, at the end of protein adsorption process, the system sustains a residual local tension to balance the difference between the actual mean curvature and the imposed spontaneous curvatures. This change in membrane tension can have a functional impact in many biological phenomena where proteins interact with membranes.Comment: 15 pages, 5 figure

Multiscale modeling for fluid transport in nanosystems.

Atomistic-scale behavior drives performance in many micro- and nano-fluidic systems, such as mircrofludic mixers and electrical energy storage devices. Bringing this information into the traditionally continuum models used for engineering analysis has proved challenging. This work describes one such approach to address this issue by developing atomistic-to-continuum multi scale and multi physics methods to enable molecular dynamics (MD) representations of atoms to incorporated into continuum simulations. Coupling is achieved by imposing constraints based on fluxes of conserved quantities between the two regions described by one of these models. The impact of electric fields and surface charges are also critical, hence, methodologies to extend finite-element (FE) MD electric field solvers have been derived to account for these effects. Finally, the continuum description can have inconsistencies with the coarse-grained MD dynamics, so FE equations based on MD statistics were derived to facilitate the multi scale coupling. Examples are shown relevant to nanofluidic systems, such as pore flow, Couette flow, and electric double layer

Theoretical and experimental studies of electrified interfaces relevant to energy storage

Advances in technology for electrochemical energy storage require increased understanding of electrolyte/electrode interfaces, including the electric double layer structure, and processes involved in charging of the interface, and the incorporation of this understanding into quantitative models. Simplified models such as Helmholtz's electric double-layer (EDL) concept don't account for the molecular nature of ion distributions, solvents, and electrode surfaces and therefore cannot be used in predictive, high-fidelity simulations for device design. This report presents theoretical results from models that explicitly include the molecular nature of the electrical double layer and predict critical electrochemical quantities such as interfacial capacitance. It also describes development of experimental tools for probing molecular properties of electrochemical interfaces through optical spectroscopy. These optical experimental methods are designed to test our new theoretical models that provide descriptions of the electric double layer in unprecedented detail

Homogeneous Non-Equilibrium Molecular Dynamics Methods for Calculating the Heat Transport Coefficient of Solids and Mixtures

This thesis presents a class of homogeneous non-equilibrium molecular dynamics (HNEMD) methods for obtaining the heat transport coefficient that relates the heat flux and temperature gradient in the linear irreversible regime. These methods are based on the linear response theory of statistical mechanics. The proposed HNEMD methods are parallelizable, and yield better statistical averages at lower overall computational cost than the existing direct and Green-Kubo methods. The HNEMD method, as it was initially proposed, is applicable only to single-species systems with two-body interactions. In this thesis, the HNEMD method is extended to single species systems with many-body interactions, and is applied to silicon systems where three-body interactions are taken into account. The HNEMD method developed for single-species systems is inadequate for obtaining the heat transport coefficient of multi-species systems. A further development of the HNEMD method, the Mixture-HNEMD (M-HNEMD) method, is presented for multi-species systems with many-body interactions. This M-HNEMD method satisfies all the requirements of linear response theory and is compatible with periodic boundary conditions. Applications of the M-HNEMD method to liquid argon-krypton systems with two-body interactions and to perfectly crystalline gallium-nitride systems with three-body interactions are presented, and the results are consistent with the results from the Green-Kubo method. This is the first HNEMD method which can be used for calculating the heat-transport coefficient of multi-species systems.The expressions for stress tensor and heat-flux vector needed for the development of HNEMD method for single-species systems and of the M-HNEMD method for multi-species systems with many-body interactions require an extension of the statistical mechanical theory of transport processes proposed by Irving and Kirkwood. This extension forms an integral part of the thesis

Mechanisms for bacterial gliding motility on soft substrates.

The motility mechanism of certain prokaryotes has long been a mystery, since their motion, known as gliding, involves no external appendages. The physical principles behind gliding still remain poorly understood. Using myxobacteria as an example of such organisms, we identify here the physical principles behind gliding motility and develop a theoretical model that predicts a 2-regime behavior of the gliding speed as a function of the substrate stiffness. Our theory describes the elasto-capillary-hydrodynamic interactions between the membrane of the bacteria, the slime it secretes, and the soft substrate underneath. Defining gliding as the horizontal translation under zero net force, we find the 2-regime behavior is due to 2 distinct mechanisms of motility thrust. On mildly soft substrates, the thrust arises from bacterial shape deformations creating a flow of slime that exerts a pressure along the bacterial length. This pressure in conjunction with the bacterial shape provides the necessary thrust for propulsion. On very soft substrates, however, we show that capillary effects must be considered that lead to the formation of a ridge at the slime-substrate-air interface, thereby creating a thrust in the form of a localized pressure gradient at the bacterial leading edge. To test our theory, we perform experiments with isolated cells on agar substrates of varying stiffness and find the measured gliding speeds in good agreement with the predictions from our elasto-capillary-hydrodynamic model. The mechanisms reported here serve as an important step toward an accurate theory of friction and substrate-mediated interactions between bacteria proliferating in soft media

Multiscale modeling for fluid transport in nanosystems.

Atomistic-scale behavior drives performance in many micro- and nano-fluidic systems, such as mircrofludic mixers and electrical energy storage devices. Bringing this information into the traditionally continuum models used for engineering analysis has proved challenging. This work describes one such approach to address this issue by developing atomistic-to-continuum multi scale and multi physics methods to enable molecular dynamics (MD) representations of atoms to incorporated into continuum simulations. Coupling is achieved by imposing constraints based on fluxes of conserved quantities between the two regions described by one of these models. The impact of electric fields and surface charges are also critical, hence, methodologies to extend finite-element (FE) MD electric field solvers have been derived to account for these effects. Finally, the continuum description can have inconsistencies with the coarse-grained MD dynamics, so FE equations based on MD statistics were derived to facilitate the multi scale coupling. Examples are shown relevant to nanofluidic systems, such as pore flow, Couette flow, and electric double layer

A comprehensive model of Drosophila epithelium reveals the role of embryo geometry and cell topology in mechanical responses.

In order to understand morphogenesis, it is necessary to know the material properties or forces shaping the living tissue. In spite of this need, very few in vivo measurements are currently available. Here, using the early Drosophila embryo as a model, we describe a novel cantilever-based technique which allows for the simultaneous quantification of applied force and tissue displacement in a living embryo. By analyzing data from a series of experiments in which embryonic epithelium is subjected to developmentally relevant perturbations, we conclude that the response to applied force is adiabatic and is dominated by elastic forces and geometric constraints, or system size effects. Crucially, computational modeling of the experimental data indicated that the apical surface of the epithelium must be softer than the basal surface, a result which we confirmed experimentally. Further, we used the combination of experimental data and comprehensive computational model to estimate the elastic modulus of the apical surface and set a lower bound on the elastic modulus of the basal surface. More generally, our investigations revealed important general features that we believe should be more widely addressed when quantitatively modeling tissue mechanics in any system. Specifically, different compartments of the same cell can have very different mechanical properties; when they do, they can contribute differently to different mechanical stimuli and cannot be merely averaged together. Additionally, tissue geometry can play a substantial role in mechanical response, and cannot be neglected