Membrane-Protein Interactions in a Generic Coarse-Grained Model for Lipid Bilayers

We study membrane-protein interactions and membrane-mediated protein-protein interactions by Monte Carlo simulations of a generic coarse-grained model for lipid bilayers with cylindrical hydrophobic inclusions. The strength of the hydrophobic force and the hydrophobic thickness of the proteins are systematically varied. The results are compared with analytical predictions of two popular analytical theories: The Landau-de Gennes theory and the elastic theory. The elastic theory provides an excellent description of the fluctuation spectra of pure membranes and successfully reproduces the deformation profiles of membranes around single proteins. However, its prediction for the potential of mean force between proteins is not compatible with the simulation data for large distances. The simulations show that the lipid-mediated interactions are governed by five competing factors: Direct interactions, lipid-induced depletion interactions, lipid bridging, lipid packing, and a smooth long-range contribution. The mechanisms leading to "hydrophobic mismatch" interactions are critically analyzed.Comment: 16 pages, 8 figures, accepted for publication in Biophysical Journa

Lipid Bilayer Deformation and the Free Energy of Interaction of a Kv Channel Gating-Modifier Toxin

A number of membrane proteins act via binding at the water/lipid bilayer interface. An important example of such proteins is provided by the gating-modifier toxins that act on voltage-gated potassium (Kv) channels. They are thought to partition to the headgroup region of lipid bilayers, and so provide a good system for probing the nature of interactions of a protein with the water/bilayer interface. We used coarse-grained molecular dynamics simulations to compute the one-dimensional potential of mean force (i.e., free energy) profile that governs the interaction between a Kv channel gating-modifier toxin (VSTx1) and model phospholipid bilayers. The reaction coordinate sampled corresponds to the position of the toxin along the bilayer normal. The course-grained representation of the protein and lipids enabled us to explore extended time periods, revealing aspects of toxin/bilayer dynamics and energetics that would be difficult to observe on the timescales currently afforded by atomistic molecular dynamics simulations. In particular, we show for this model system that the bilayer deforms as it interacts with the toxin, and that such deformations perturb the free energy profile. Bilayer deformation therefore adds an additional layer of complexity to be addressed in investigations of membrane/protein systems. In particular, one should allow for local deformations that may arise due to the spatial array of charged and hydrophobic elements of an interfacially located membrane protein

Continuum descriptions of membranes and their interaction with proteins: Towards chemically accurate models

Biological membranes deform in response to resident proteins leading to a coupling between membrane shape and protein localization. Additionally, the membrane influences the function of membrane proteins. Here we review contributions to this field from continuum elastic membrane models focusing on the class of models that couple the protein to the membrane. While it has been argued that continuum models cannot reproduce the distortions observed in fully-atomistic molecular dynamics simulations, we suggest that this failure can be overcome by using chemically accurate representations of the protein. We outline our recent advances along these lines with our hybrid continuum-atomistic model, and we show the model is in excellent agreement with fully-atomistic simulations of the nhTMEM16 lipid scramblase. We believe that the speed and accuracy of continuum-atomistic methodologies will make it possible to simulate large scale, slow biological processes, such as membrane morphological changes, that are currently beyond the scope of other computational approaches

Multiscale (re)modeling of lipid bilayer membranes

Many phenomena of biological membranes are inherently multi-scalar: their observation and description requires very different length and time-scales. Especially membrane remodeling processes that are essential for many important cellular activities, such as endo- and exocytosis, cell division, infection, immune response or cell-cell communication, involve large scale morphological changes of the membrane, which are initiated and controlled by molecular scale interactions. The large-scale behavior of the membrane is intimately coupled to the atomic detail of the system, so that for a successful model, the contributions at all scales have to be included. This chapter gives an overview of the computational methods used to model lipid bilayer membranes and their interactions with proteins and other molecules at different levels of resolution. Then strategies to connect the models at different scales in order to provide a multi-resolution picture are described. Methods to quantify free energy changes associated with complex collective rearrangement are outlined. The last sections summarize examples for the application of these methods to pore formation, reshaping membranes into buds and membrane tubes as well as membrane fusion