Membrane proteins diffuse as dynamic complexes with lipids

Biophysical Journal, 2008

As charged macromolecules adsorb and diffuse on cell membranes in a large variety of cell signaling processes, they can attract or repel oppositely charged lipids. This results in lateral membrane rearrangement and affects the dynamics of protein function. To address such processes quantitatively we introduce a dynamic mean-field scheme that allows selfconsistent calculations of the equilibrium state of membrane-protein complexes after such lateral reorganization of the membrane components, and serves to probe kinetic details of the process. Applicable to membranes with heterogeneous compositions containing several types of lipids, this comprehensive method accounts for mobile salt ions and charged macromolecules in three dimensions, as well as for lateral demixing of charged and net-neutral lipids in the membrane plane. In our model, the mobility of membrane components is governed by the diffusion-like Cahn-Hilliard equation, while the local electrochemical potential is based on nonlinear Poisson-Boltzmann theory. We illustrate the method by applying it to the adsorption of the anionic polypeptide poly-Lysine on negatively charged lipid membranes composed of binary mixtures of neutral and monovalent lipids, or onto ternary mixtures of neutral, monovalent, and multivalent lipids. Consistent with previous calculations and experiments, our results show that at steady-state multivalent lipids (such as PIP 2 ), but not monovalent lipid (such as phosphatidylserine), will segregate near the adsorbing macromolecules. To address the corresponding diffusion of the adsorbing protein in the membrane plane, we couple lipid mobility with the propagation of the adsorbing protein through a dynamic Monte Carlo scheme. We find that due to their higher mobility dictated by the electrochemical potential, multivalent lipids such as PIP 2 more quickly segregate near oppositely charged proteins than do monovalent lipids, even though their diffusion constants may be similar. The segregation, in turn, slows protein diffusion, as lipids introduce an effective drag on the motion of the adsorbate. In contrast, monovalent lipids such as phosphatidylserine only weakly segregate, and the diffusions of protein and lipid remain largely uncorrelated.

Positively charged polybasic domains are essential for recruiting multiple signaling proteins, such as Ras GTPases and Src kinase, to the negatively charged cellular membranes. Much less, however, is known about the influence of electrostatic interactions on the lateral dynamics of these proteins. We developed a dynamic Monte-Carlo automaton that faithfully simulates lateral diffusion of the adsorbed positively charged oligopeptides as well as the dynamics of mono-(phosphatidylserine) and polyvalent (PIP 2 ) anionic lipids within the bilayer. In agreement with earlier results, our simulations reveal lipid demixing that leads to the formation of a lipid shell associated with the peptide. The computed association times and average numbers of bound lipids demonstrate that tetravalent PIP 2 interacts with the peptide much more strongly than monovalent lipid. On the spatially homogeneous membrane, the lipid shell affects the behavior of the peptide only by weakly reducing its lateral mobility. However, spatially heterogeneous distributions of monovalent lipids are found to produce peptide drift, the velocity of which is determined by the total charge of the peptide-lipid complex. We hypothesize that this predicted phenomenon may affect the spatial distribution of proteins with polybasic domains in the context of cell-signaling events that alter the local density of monovalent anionic lipids.

Biological membranes are complex environments whose physico-chemical properties are of utmost importance for the understanding of many crucial biological processes. Much attention has been given in the literature to the description of membranes along the z-axis perpendicular to the membrane. Here, we instead consider the lateral dynamics of lipids and peripheral proteins due to their electrostatic interaction. Previously, we constructed a Monte Carlo automaton capable of simulating mutual diffusive dynamics of charged lipids and associated positively charged peptides. Here, we derive and numerically analyze a system of Poisson-Boltzmann-Nernst-Planck (PBNP) equations that provide a mean-field approximation compatible with our Monte Carlo model. The thorough comparison between the mean-field PBNP equations and Monte Carlo simulations demonstrates that both the approaches are in a good qualitative agreement in all tested scenarios. We find that the two methods quantitatively deviate when the local charge density is high, presumably because the Poisson-Boltzmann formalism is applicable in the so-called weak coupling limit, whose validity is restricted to low charge densities. Nevertheless, we conclude that the mean-field PBNP approach provides a good approximation for the considerably more detailed Monte Carlo model at only a fraction of the associated computational cost and allows simulation of the membrane lateral dynamics on the space and time scales relevant for the realistic biological problems. to 129.215.4.23. Redistribution subject to AIP license or copyright; see http://jcp.aip.org/about/rights_and_permissions 155103-2 Kiselev et al.

The Journal of Chemical Physics, 2006

Molecular dynamics (MD) simulations of dipalmitoylphosphatidylcholine (DPPC) bilayers composed of 72 and 288 lipids are used to examine system size dependence on dynamical properties associated with the particle-mesh Ewald (PME) treatment of electrostatic interactions. The lateral diffusion constant, D ℓ , is 2.92×10 −7 and 0.95×10 −7 cm 2 /s for 72 and 288 lipids, respectively. This dramatic finite size effect originates from the correlation length of lipid diffusion, which extends to next-nearest neighbors in the 288 lipid system. Consequently, diffusional events in smaller systems can propagate across the boundaries of the periodic box. The internal dynamics of lipids calculated from the PME simulations are independent of the system size. Specifically, reorientational correlation functions for the slowly relaxing phosphorous-glycerol hydrogen, phosphorous-nitrogen vectors, and more rapidly relaxing CH vectors in the aliphatic chains are equivalent for the 72 and 288 lipid simulations. A third MD simulation of a bilayer with 72 lipids using spherical force-shift electrostatic cutoffs resulted in interdigitated chains, thereby rendering this cutoff method inappropriate.

Faraday Discussions, 2014

Interactions between lipids and membrane proteins play a key role in determining the nanoscale dynamic and structural properties of biological membranes. Molecular dynamics (MD) simulations provide a valuable tool for studying membrane models, complementing experimental approaches. It is now possible to simulate large membrane systems, such as simplified models of 15 bacterial and viral envelope membranes. Consequently, there is a pressing need to develop tools to visualize and quantify the dynamics of these immense systems, which typically are comprised of millions of particles. To tackle this issue, we have developed visual and quantitative analyses of molecular positions and their velocity field using path line, vector field and streamline techniques. This allows us to highlight large, transient flow-like movements of lipids and to better understand crowding within the lipid bilayer. The current study focuses on visualization and analysis of lipid dynamics. However, the methods are flexible and can be readily applied to e.g. proteins and nanoparticles within large complex membranes. The protocols developed here are readily accessible 25 both as a plugin for the molecular visualization program VMD and as a module for the MDAnalysis library. 40 current.

Molecular Simulation, 2013

The accuracy of the Content should not be relied upon and should be independently verified with primary sources of information. Taylor and Francis shall not be liable for any losses, actions, claims, proceedings, demands, costs, expenses, damages, and other liabilities whatsoever or howsoever caused arising directly or indirectly in connection with, in relation to or arising out of the use of the Content. This article may be used for research, teaching, and private study purposes. Any substantial or systematic reproduction, redistribution, reselling, loan, sub-licensing, systematic supply, or distribution in any form to anyone is expressly forbidden.

Biochimica Et Biophysica Acta-reviews on Biomembranes, 1998

The fundamental physical principles of the lateral organization of trans-membrane proteins and peptides as well as peripheral membrane proteins and enzymes are considered from the point of view of the lipid-bilayer membrane, its structure, dynamics, and cooperative phenomena. Based on a variety of theoretical considerations and model calculations, the nature of lipid-protein interactions is considered both for a single protein and an assembly of proteins that can lead to aggregation and protein crystallization in the plane of the membrane. Phenomena discussed include lipid sorting and selectivity at protein surfaces, protein-lipid phase equilibria, lipid-mediated protein-protein interactions, wetting and capillary condensation as means of protein organization, mechanisms of two-dimensional protein crystallization, as well as non-equilibrium organization of active proteins in membranes. The theoretical findings are compared with a variety of experimental data. ß 1998 Elsevier Science B.V. All rights reserved.

The Journal of Chemical Physics, 2011

This communication presents a molecular dynamics simulation study of a bilayer consisting of 128 dioleoyl-sn-glycero-3-phosphocholine molecules, which focusses on the center-of-mass diffusion of the lipid molecules parallel to the membrane plane. The analysis of the simulation results is performed within the framework of the generalized Langevin equation and leads to a consistent picture of subdiffusion. The mean square displacement of the lipid molecules evolves as ∝ t α , with α between 0.5 and 0.6, and the fractional diffusion coefficient is close to the experimental value for a similar system obtained by fluorescence correlation spectroscopy. We show that the long-time tails of the lateral velocity autocorrelation function and the associated memory function agree well with exact results which have been recently derived by asymptotic analysis [G. Kneller, J. Chem. Phys. 134, 224106 ]. In this context, we define characteristic time scales for these two quantities.

The Journal of chemical physics, 2016

The protein transport inside a cell is a complex phenomenon that goes through several difficult steps. The facilitated transport requires sophisticated machineries involving protein assemblies. In this work, we developed a diffusion-reaction model to simulate co-transport kinetics of proteins and lipids. We assume the following: (a) there is always a small lipid concentration of order of the Critical Micellar Concentration (CMC) in equilibrium with the membrane; (b) the binding of lipids to proteins modulates the hydrophobicity of the complexes and, therefore, their ability to interact and merge with the bilayer; and (c) some lipids leave the bilayer to replenish those bound to proteins. The model leads to a pair of integral equations for the time-evolution of the adsorbed proteins in the lipid bilayer. Relationships between transport kinetics, CMC, and lipid-protein binding constants were found. Under particular conditions, a perturbation analysis suggests the onset of kinks in the...

Biophysical Journal, 2000

This report addresses the following problems associated with the generation of computer models of phospholipid bilayer membranes using molecular dynamics simulations: arbitrary initial structures and short equilibration periods, an Ewald-induced strong coupling of phospholipids, uncertainty regarding which value should be used for surface tension to alleviate the problem of the small size of the membrane, and simultaneous realization of both order parameters and the surface area. We generated a computer model of the liquid-crystalline L-␣-dimyristoylphosphatidylcholine (DMPC) bilayer, starting from a configuration based on a crystal structure (rather than from an arbitrary structure). To break the crystalline structure, a 20-ps high-temperature pulse of 510 K (but not 450 or 480 K) was effective. The system finally obtained is an all-atom model, with Ewald summation to evaluate Coulombic interactions and a constant surface tension of 35 dynes/cm/ water-membrane interface, equilibrated for 12 ns (over 50 ns total calculation time), which reproduces all of the experimentally observed parameters examined in this work. Furthermore, this model shows the presence of significant orientational correlations between neighboring alkyl chains and between shoulder vectors (which show the orientations of the lipids about their long axes) of neighboring DMPCs.