Biomembrane simulations of 12 lipid types using the general amber force field in a tensionless ensemble

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid–protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid–protein interactions, and to link lipid–protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid–protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid–protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid–protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid–protein interactions.

What Can and Cannot Be Learned from Molecular Dynamics Simulations of Bacterial Proton-Coupled Oligopeptide Transporter GkPOT?

We have performed an extensive set of all-atom molecular dynamics (MD) simulations of a bacterial proton-coupled oligopeptide transporter (POT) in an explicit membrane environment. We have characterized both the local and global conformational dynamics of the transporter upon the proton and/or substrate binding, within a statistical framework. Our results reveal a clearly distinct behavior for local conformational dynamics in the absence and presence of the proton at the putative proton binding residue E310. Particularly, we find that the substrate binding conformation is drastically different in the two conditions, where the substrate binds to the protein in a lateral/vertical manner, in the presence/absence of the proton. We do not observe any statistically significant distinctive behavior in terms of the global conformational changes in different simulation conditions, within the time scales of our simulations. Our extensive simulations and analyses call into question the implicit assumption of many MD studies that local conformational changes observed in short simulations could provide clues to the global conformational changes that occur on much longer time scales. The linear regression analysis of quantities associated with the global conformational fluctuations, however, provides an indication of a mechanism involving the concerted motion of the transmembrane helices, consistent with the rocker-switch mechanism.

Influence of temperature, anions and size distribution on the zeta potential of DMPC, DPPC and DMPE lipid vesicles

The purpose of the work is to compare the influence of the multilamellarity, phase state, lipid head groups and ionic media on the origin of the surface potential of lipid membranes. With this aim, we present a new analysis of the zeta potential of multilamellar and unilamellar vesicles composed by phosphatidylcholines (PC) and phosphatidylethanolamines (PE) dispersed in water and ionic solutions of polarizable anions, at temperatures below and above the phase transition. In general, the adsorption of anions seems to explain the origin of the zeta potential in vesicles only above the transition temperature (Tc). In this case, the sign of the surface potential is ascribed to a partial orientation of head group moiety toward the aqueous phase. This is noticeable in PC head groups but not in PEs, due to the strong lateral interaction between PO and NH group in PE.

Interfacial water screens the protein-induced transmembrane voltage

Transmembrane proteins are crucial in cellular traffic, signal transduction, and energy storage in a form of transmembrane voltage. These proteins are stabilized by hydrophobic and hydrophilic interactions, in which cytoplasmic and exoplasmic water plays a special role. Water structural ordering generates the dipole potential that typically overcompensates for an intrinsic membrane-protein potential gradient, and thus it modifies and sustains an overall cellular electrostatics. Although the transmembrane voltage has been extensively studied, the dipole potential has attracted very little attention. Here, by using molecular dynamics, we examined water electrostatic response to the transmembrane charge, field, and potential asymmetry introduced by the presence of four integral membrane proteins: typical of inner (α-helix) and outer membrane (β-barrel). In all cases, the protein presence introduces electrostatic directionality in the transmembrane dipole field and voltage. In particular, water generates a deep potential sink if strongly polar residues are densely packed on one side of bilayer, as frequently occurs in a selectivity filter of the K(+) channel. We also found that protein secondary structure is less important than the polar residue distribution along the protein channel. Our findings are relevant for understanding the driving force behind biomembrane conductivity: the ability of biological water to electrostatically screen the transmembrane voltage.