Biomembranes in atomistic and coarse-grained simulations

Influence of phospholipid head and tail molecular structures on cell membrane mechanical response under tension

Biological cell membranes are primarily comprised of a diverse lipid bilayer with multiple phospholipid (lipid) types, each of which is comprised of a hydrophilic headgroup and two hydrophobic hydrocarbon tails. The lipid type determines the molecular structure of head and tail groups, which can affect membrane mechanics at nanoscale and subsequently cell viability under mechanical loading. Hence, using molecular dynamics simulations, the current study investigated seven membrane phospholipids and the effect of their structural differences on physical deformation, mechanoporation damage, and mechanical failure of the membranes under tension. The inspected phospholipids showed similar yield stresses and strains, as well as pore evolution and damage, but significantly different failure strains. In general, failure occurred at a lower strain for lipids with a larger equilibrium area per lipid. The obtained results suggest that larger headgroup structure, greater degree of unsaturation, and tail-length asymmetry influenced the phospholipids' ability to pack against each other, increased the fluidity and equilibrium area per lipid of the membrane, and resulted in lower failure strain. Overall, this study provides insights on how different phospholipid structures affect membrane physical responses at the molecular level and serves as a reference for future studies of more complex membrane systems with intricate biophysical properties.

Effect of Force Field Resolution on Membrane Mechanical Response and Mechanoporation Damage under Deformation Simulations

Damage induced by transient disruption and mechanoporation in an intact cell membrane is a vital nanoscale biomechanical mechanism that critically affects cell viability. To complement experimental studies of mechanical membrane damage and disruption, molecular dynamics (MD) simulations have been performed at different force field resolutions, each of which follows different parameterization strategies and thus may influence the properties and dynamics of membrane systems. Therefore, the current study performed tensile deformation MD simulations of bilayer membranes using all-atom (AA), united-atom (UA), and coarse-grained Martini (CG-M) models to investigate how the damage biomechanics differs across atomistic and coarse-grained (CG) simulations. The mechanical response and mechanoporation damage were qualitatively similar but quantitatively different in the three models, including some progressive changes based on the coarse-graining level. The membranes yielded and reached ultimate strength at similar strains; however, the coarser systems exhibited lower average yield stresses and failure strains. The average failure strain in the UA model was approximately 7% lower than the AA, and the CG-M was 20% lower than UA and 27% lower than AA. The CG systems also nucleated a higher number of pores and larger pores, which resulted in higher damage during the deformation process. Overall, the study provides insight on the impact of force field-a critical factor in modeling biomolecular systems and their interactions-in inspecting membrane mechanosensitive responses and serves as a reference for justifying the appropriate force field for future studies of more complex membranes and more diverse biomolecular assemblies.

Experimental and Theoretical Brownian Dynamics Analysis of Ion Transport During Cellular Electroporation of E. coli Bacteria

Escherichia coli bacterium is a rod-shaped organism composed of a complex double membrane structure. Knowledge of electric field driven ion transport through both membranes and the evolution of their induced permeabilization has important applications in biomedical engineering, delivery of genes and antibacterial agents. However, few studies have been conducted on Gram-negative bacteria in this regard considering the contribution of all ion types. To address this gap in knowledge, we have developed a deterministic and stochastic Brownian dynamics model to simulate in 3D space the motion of ions through pores formed in the plasma membranes of E. coli cells during electroporation. The diffusion coefficient, mobility, and translation time of Ca2+, Mg2+, Na+, K+, and Cl- ions within the pore region are estimated from the numerical model. Calculations of pore's conductance have been validated with experiments conducted at Gustave Roussy. From the simulations, it was found that the main driving force of ionic uptake during the pulse is the one due to the externally applied electric field. The results from this work provide a better understanding of ion transport during electroporation, aiding in the design of electrical pulses for maximizing ion throughput, primarily for application in cancer treatment.

Molecular assessment of drug-phospholipid interactions consequent to cancer treatment: a study of anthracycline-induced cardiotoxicity

Cardiotoxicity limits the use of anthracyclines as potent chemotherapeutics. We employ classical molecular dynamics to explore anthracycline interactions with a realistic myocardial membrane and compare to an ideal membrane widely used in literature. The interaction of these two membranes with four anthracyclines; doxorubicin, epirubicin, daunorubicin, and idarubicin are studied. Careful analysis was conducted on three forms of each drug; pristine, primary metabolite, and cationic salt. By examining the molecular residence time near the membrane's surface, the average number of molecule/membrane hydrogen bonds, the immobilization of the molecules near the membrane, and the location of those molecules relative to the mid-plane of the membrane we found out that salt forms exhibit the highest cardiotoxic probability, followed by the metabolites and pristine forms. Additionally, all forms have more affinity to the upper layer of the realistic myocardial membrane. Meanwhile, an ideal membrane consisting of a single type of phospholipids is not capable of capturing the specific interactions of each drug form. These findings confirm that cardiotoxic mechanisms are membrane-layer and drug-form dependent.

How Do Cyclopropane Fatty Acids Protect the Cell Membrane of Escherichia coli in Cold Shock?

The cyclopropanation of unsaturated lipid acyl chains of some bacterial cell membranes is an important survival strategy to protect the same against drastic cooling. To elucidate the role of cyclopropane ring-containing lipids, we have simulated the lipid membrane of Escherichia coli (E. coli) and two modified membranes by replacing the cyclopropane rings with either single or double bonds at widely different temperatures. It has been observed that the cyclopropane rings provide more rigid kinks in the lipid acyl chain compared to the double bonds and therefore further reduce the packing density of the membrane and subsequently enhance the membrane fluidity at low temperatures. They also inhibit the close packing of other lipids and deleterious phase separation by strongly interacting with them. Therefore, this study has explained why E. coli bacterial strain, susceptible to freezing environments, relies on the cyclopropanation of an unsaturated chain.

Homeoviscous Adaptation of the Lipid Membrane of a Soil Bacterium Surviving under Diurnal Temperature Variation: A Molecular Simulation Perspective

A recent experiment has reported the lipidome remodeling of a soil-based plant-associated bacterium Methylobacterium extorquens due to diurnal temperature variations. The key adaptation strategy is the headgroup-specific remodeling of the acyl chain. To understand the idiosyncratic adaptation at the molecular level, we simulate the model membrane of the same bacterium using the reported lipidome compositions at four different experimental temperatures. We investigate the temperature-dependent packing density and fluidity of the membrane, the constancy of which is key to the homeoviscous adaptation. The results show that complex lipidome remodeling approximately preserves membrane properties under heat and cold stress. The headgroup-specific remodeling of the acyl chain serves to fine-tune the packing density and fluidity of the membrane at different temperatures. While lipids with strongly interacting headgroups are more abundant at higher temperatures, the lipidome is more dominated by lipids with weaker interacting headgroups at lower temperatures. This adaptation alleviates lipid membrane disruption caused by heat and cold stress. This study provides a molecular picture of the homeoviscous adaptation of the realistic lipid membrane of a soil-based bacterium.

Antimicrobial Peptide Mechanism Studied by Scattering-Guided Molecular Dynamics Simulation

In an effort to combat rising antimicrobial resistance, our labs have rationally designed cationic, helical, amphipathic antimicrobial peptides (AMPs) as alternatives to traditional antibiotics since AMPs incur bacterial resistance in weeks, rather than days. One highly positively charged AMP, WLBU2 (+13e), (RRWV RRVR RWVR RVVR VVRR WVRR), has been shown to be effective in killing both Gram-negative (G(-)) and Gram-positive (G(+)) bacteria by directly perturbing the bacterial membrane nonspecifically. Previously, we used two equilibrium experimental methods: synchrotron X-ray diffuse scattering (XDS) providing lipid membrane thickness and neutron reflectometry (NR) providing WLBU2 depth of penetration into three lipid model membranes (LMMs). The purpose of the present study is to use the results from the scattering experiments to guide molecular dynamics (MD) simulations to investigate the detailed biophysics of the interactions of WLBU2 with LMMs of Gram-negative outer and inner membranes, and Gram-positive cell membranes, to elucidate the mechanisms of bacterial killing. Instead of coarse-graining, backmapping, or simulating without bias for several microseconds, all-atom (AA) simulations were guided by the experimental results and then equilibrated for ∼0.5 μs. Multiple replicas of the inserted peptide were run to probe stability and reach a combined time of at least 1.2 μs for G(-) and also 2.0 μs for G(+). The simulations with experimental comparisons help rule out certain structures and orientations and propose the most likely set of structures, orientations, and effects on the membrane. The simulations revealed that water, phosphates, and ions enter the hydrocarbon core when WLBU2 is positioned there. For an inserted peptide, the three types of amino acids, arginine, tryptophan, and valine (R, W, V), are arranged with the 13 Rs extending from the hydrocarbon core to the phosphate group, Ws are located at the interface, and Vs are more centrally located. For a surface state, R, W, and V are positioned relative to the bilayer interface as expected from their hydrophobicities, with Rs closest to the phosphate group, Ws close to the interface, and Vs in between. G(-) and G(+) LMMs are thinned ∼1 Å by the addition of WLBU2. Our results suggest a dual anchoring mechanism for WLBU2 both in the headgroup and in the hydrocarbon region that promotes a defect region where water and ions can flow across the slightly thinned bacterial cell membrane.

Organelle Interaction and Drug Discovery: Towards Correlative Nanoscopy and Molecular Dynamics Simulation

The inter-organelle interactions, including the cytomembrane, endoplasmic reticulum, mitochondrion, lysosome, dictyosome, and nucleus, play the important roles in maintaining the normal function and homeostasis of cells. Organelle dysfunction can lead to a range of diseases (e.g., Alzheimer's disease (AD), Parkinson's disease (PD), and cancer), and provide a new perspective for drug discovery. With the development of imaging techniques and functional fluorescent probes, a variety of algorithms and strategies have been developed for the ever-improving estimation of subcellular structures, organelle interaction, and organelle-related drug discovery with accounting for the dynamic structures of organelles, such as the nanoscopy technology and molecular dynamics (MD) simulations. Accordingly, this work summarizes a series of state-of-the-art examples of the recent progress in this rapidly changing field and uncovering the drug screening based on the structures and interactions of organelles. Finally, we propose the future outlook for exciting applications of organelle-related drug discovery, with the cooperation of nanoscopy and MD simulations.

Anisotropic Electrostatic Interactions in Coarse-Grained Water Models to Enhance the Accuracy and Speed-Up Factor of Mesoscopic Simulations

Water models with realistic physical–chemical properties are essential to study a variety of biomedical processes or engineering technologies involving molecules or nanomaterials. Atomistic models of water are constrained by the feasible computational capacity, but calibrated coarse-grained (CG) ones can go beyond these limits. Here, we compare three popular atomistic water models with their corresponding CG model built using finite-size particles such as ellipsoids. Differently from previous approaches, short-range interactions are accounted for with the generalized Gay–Berne potential, while electrostatic and long-range interactions are computed from virtual charges inside the ellipsoids. Such an approach leads to a quantitative agreement between the original atomistic models and their CG counterparts. Results show that a timestep of up to 10 fs can be achieved to integrate the equations of motion without significant degradation of the physical observables extracted from the computed trajectories, thus unlocking a significant acceleration of water-based mesoscopic simulations at a given accuracy.

The Plasma Membrane at the Cornerstone Between Flexibility and Adaptability: Implications for Saccharomyces cerevisiae as a Cell Factory

In the last decade, microbial-based biotechnological processes are paving the way toward sustainability as they implemented the use of renewable feedstocks. Nonetheless, the viability and competitiveness of these processes are often limited due to harsh conditions such as: the presence of feedstock-derived inhibitors including weak acids, non-uniform nature of the substrates, osmotic pressure, high temperature, extreme pH. These factors are detrimental for microbial cell factories as a whole, but more specifically the impact on the cell’s membrane is often overlooked. The plasma membrane is a complex system involved in major biological processes, including establishing and maintaining transmembrane gradients, controlling uptake and secretion, intercellular and intracellular communication, cell to cell recognition and cell’s physical protection. Therefore, when designing strategies for the development of versatile, robust and efficient cell factories ready to tackle the harshness of industrial processes while delivering high values of yield, titer and productivity, the plasma membrane has to be considered. Plasma membrane composition comprises diverse macromolecules and it is not constant, as cells adapt it according to the surrounding environment. Remarkably, membrane-specific traits are emerging properties of the system and therefore it is not trivial to predict which membrane composition is advantageous under certain conditions. This review includes an overview of membrane engineering strategies applied to Saccharomyces cerevisiae to enhance its fitness under industrially relevant conditions as well as strategies to increase microbial production of the metabolites of interest.

Lipid Dynamics in Membranes Slowed Down by Transmembrane Proteins

Lipids and proteins, as essential components of biological cell membranes, exhibit a significant degree of freedom for different kinds of motions including lateral long-range mobility. Due to their interactions, they not only preserve the cellular membrane but also contribute to many important cellular functions as e.g., signal transport or molecular exchange of the cell with its surrounding. Many of these processes take place on a short time (up to some nanoseconds) and length scale (up to some nanometers) which is perfectly accessible by quasielastic neutron scattering (QENS) experiments and molecular dynamics (MD) simulations. In order to probe the influence of a peptide, a transmembrane sequence of the transferrin receptor (TFRC) protein, on the dynamics of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) large unilamellar vesicles (LUVs) on a nanosecond time scale, high-resolution QENS experiments and complementary MD simulations have been utilized. By using different scattering contrasts in the experiment (chain-deuterated lipids and protonated lipids, respectively), a model could be developed which allows to examine the lipid and peptide dynamics separately. The experimental results revealed a restricted lipid lateral mobility in the presence of the TFRC transmembrane peptides. Also the apparent self-diffusion coefficient of the lateral movement of the peptide molecules could be determined quantitatively for the probed short-time regime. The findings could be confirmed very precisely by MD simulations. Furthermore, the article presents an estimation for the radius of influence of the peptides on the lipid long-range dynamics which could be determined by consistently combining results from experiment and simulation.

Serotonin Alters the Phase Equilibrium of a Ternary Mixture of Phospholipids and Cholesterol

Unsaturated and saturated phospholipids tend to laterally segregate, especially in the presence of cholesterol. Small molecules such as neurotransmitters, toxins, drugs etc. possibly modulate this lateral segregation. The small aromatic neurotransmitter serotonin (5-HT) has been found to bind to membranes. We studied the lipid structure and packing of a ternary membrane mixture consisting of palmitoyl-oleoyl-phosphatidylcholine, palmitoyl-sphingomyelin, and cholesterol at a molar ratio of 4/4/2 in the absence and in the presence of 5-HT, using a combination of solid-state 2H NMR, atomic force microscopy, and atomistic molecular dynamics (MD) simulations. Both NMR and MD report formation of a liquid ordered (L o ) and a liquid disordered (L d ) phase coexistence with small domains. Lipid exchange between the domains was fast such that single component 2H NMR spectra are detected over a wide temperature range. A drastic restructuring of the domains was induced when 5-HT is added to the membranes at a 9 mol% concentration relative to the lipids. 2H NMR spectra of all components of the mixture showed two prominent contributions indicative of molecules of the same kind residing both in the disordered and the ordered phase. Compared to the data in the absence of 5-HT, the lipid chain order in the disordered phase was further decreased in the presence of 5-HT. Likewise, addition of serotonin increased lipid chain order within the ordered phase. These characteristic lipid chain order changes were confirmed by MD simulations. The 5-HT-induced larger difference in lipid chain order results in more pronounced differences in the hydrophobic thickness of the individual membrane domains. The correspondingly enlarged hydrophobic mismatch between ordered and disordered phases is assumed to increase the line tension at the domain boundary, which drives the system into formation of larger size domains. These results not only demonstrate that small membrane binding molecules such as neurotransmitters have a profound impact on essential membrane properties. It also suggests a mechanism by which the interaction of small molecules with membranes can influence the function of membrane proteins and non-cognate receptors. Altered membrane properties may modify lateral sorting of membrane protein, membrane protein conformation, and thus influence their function as suspected for neurotransmitters, local anesthetics, and other small drug molecules.

Molecular Dynamics Simulation of Small Molecules Interacting with Biological Membranes

Cell membranes protect and compartmentalise cells and their organelles. The semi-permeable nature of these membranes controls the exchange of solutes across their structure. Characterising the interaction of small molecules with biological membranes is critical to understanding of physiological processes, drug action and permeation, and many biotechnological applications. This review provides an overview of how molecular simulations are used to study the interaction of small molecules with biological membranes, with a particular focus on the interactions of water, organic compounds, drugs and short peptides with models of plasma cell membrane and stratum corneum lipid bilayers. This review will not delve on other types of membranes which might have different composition and arrangement, such as thylakoid or mitochondrial membranes. The application of unbiased molecular dynamics simulations and enhanced sampling methods such as umbrella sampling, metadynamics and replica exchange are described using key examples. This review demonstrates how state-of-the-art molecular simulations have been used successfully to describe the mechanism of binding and permeation of small molecules with biological membranes, as well as associated changes to the structure and dynamics of these membranes. The review concludes with an outlook on future directions in this field.

The Emerging Role of the Mammalian Glycocalyx in Functional Membrane Organization and Immune System Regulation

All cells in the human body are covered by a dense layer of sugars and the proteins and lipids to which they are attached, collectively termed the "glycocalyx." For decades, the organization of the glycocalyx and its interplay with the cellular state have remained enigmatic. This changed in recent years. Latest research has shown that the glycocalyx is an organelle of vital significance, actively involved in and functionally relevant for various cellular processes, that can be directly targeted in therapeutic contexts. This review gives a brief introduction into glycocalyx biology and describes the specific challenges glycocalyx research faces. Then, the traditional view of the role of the glycocalyx is discussed before several recent breakthroughs in glycocalyx research are surveyed. These results exemplify a currently unfolding bigger picture about the role of the glycocalyx as a fundamental cellular agent.

Molecular mechanism for bidirectional regulation of CD44 for lipid raft affiliation by palmitoylations and PIP2

The co-localization of Cluster-of-Differentiation-44 protein (CD44) and cytoplasmic adaptors in specific membrane environments is crucial for cell adhesion and migration. The process is controlled by two different pathways: On the one hand palmitoylation keeps CD44 in lipid raft domains and disables the linking to the cytoplasmic adaptor, whereas on the other hand, the presence of phosphatidylinositol-4,5-biphosphate (PIP2) lipids accelerates the formation of the CD44-adaptor complex. The molecular mechanism explaining how CD44 is migrating into and out of the lipid raft domains and its dependence on both palmitoylations and the presence of PIP2 remains, however, elusive. In this study, we performed extensive molecular dynamics simulations to study the raft affinity and translocation of CD44 in phase separated model membranes as well as more realistic plasma membrane environments. We observe a delicate balance between the influence of the palmitoylations and the presence of PIP2 lipids: whereas the palmitoylations of CD44 increases the affinity for raft domains, PIP2 lipids have the opposite effect. Additionally, we studied the association between CD44 and the membrane adaptor FERM in dependence of these factors. We find that the presence of PIP2 lipids allows CD44 and FERM to associate in an experimentally observed binding mode whereas the highly palmitoylated species shows no binding affinity. Together, our results shed light on the sophisticated mechanism on how membrane translocation and peripheral protein association can be controlled by both protein modifications and membrane composition.

Cytoskeleton protein complex involving with association of CD44 and ERMs is critical for cancer-related cellular adhesion and migration. The protein interactions are found to be modulated by chemical modification and membrane microenvironments, but the inherent mechanism is unclear. We obtained molecular dynamic details of CD44 localization switching between raft/non-raft subdomains regulated by palmitoylations and PIP2 molecules. Binding of PIP2 on the palmitoylated CD44 enables it to release from lipid raft, revealing an exceptional role of PIP2 in mediating protein translocation. PIP2 is beneficial for CD44 to associate with the active domain of ERM, in a nearly crystal structure mode. The molecular information will enhance our understanding for PIP2 regulation to protein translocation and membrane association.

The orientation of a membrane probe from structural analysis by enhanced Raman scattering

Small fluorescent molecules are widely used as probes of biomembranes. Different probes optically indicate membrane properties such as the lipid phase, thickness, viscosity, and electrical potential. The detailed molecular mechanisms behind probe signals are not well understood, in part due to the lack of tools to determine probe position and orientation in the membrane. Optical measurements on aligned biomembranes and lipid bilayers provide some degree of orientational information based on anisotropy in absorption, fluorescence, or nonlinear optical properties. These methods typically find the polar tilt angle between the membrane normal and the long axis of the molecule. Here we show that solution-phase surface enhanced Raman scattering (SERS) spectra of lipid membranes on gold nanorods can be used to determine molecular orientation of molecules within the membrane. The voltage sensitive dye 4-(2-(6-(dibutylamino)-2-naphthalenyl)ethenyl)-1-(3-sulfopropyl)-hydroxide, known as di-4-ANEPPS, is studied. Through the analysis of several peaks in the SERS spectrum, the polar angle from the membrane normal is found to be 66°, and the roll angle around the long axis of the molecule to be 305° from the original orientation. This structural analysis method could help elucidate the meaning of fluorescent membrane probe signals, and how they are affected by different lipid compositions.

Molecular Dynamics of the Association of L-Selectin and FERM Regulated by PIP2

Phosphatidylinositol 4,5-bisphosphate (PIP2) acts as a signaling lipid, mediating membrane trafficking and recruitment of proteins to membranes. A key example is the PIP2-dependent regulation of the adhesion of L-selectin to the cytoskeleton adaptors of the N-terminal subdomain of ezrin-radixin-moesin (FERM). The molecular details of the mediating behavior of multivalent anionic PIP2 lipids in this process, however, remain unclear. Here, we use coarse-grained molecular dynamics simulation to explore the mechanistic details of PIP2 in the transformation, translocation, and association of the FERM/L-selectin complex. We compare membranes of different compositions and find that anionic phospholipids are necessary for both FERM and the cytoplasmic domain of L-selectin to absorb on the membrane surface. The subsequent formation of the FERM/L-selectin complex is strongly favored by the presence of PIP2, which clusters around both proteins and triggers a conformational transition in the cytoplasmic domain of L-selectin. We are able to quantify the effect of PIP2 on the association free energy of the complex by means of a potential of mean force. We conclude that PIP2 behaves as an adhesive agent to enhance the stability of the FERM/L-selectin complex and identify key residues involved. The molecular information revealed in this study highlights the specific role of membrane lipids such as PIP2 in protein translocation and potential signaling.

Structural Model of the mIgM B-Cell Receptor Transmembrane Domain From Self-Association Molecular Dynamics Simulations

Antigen binding to B-cell antigen receptors (BCRs) followed by signaling initiates the humoral immune response. The signaling is intimately coupled to nanoclustering of BCRs and their sorting to specific membrane domains, a process that is ruled by interactions between the BCR transmembrane domain and lipids. While the structure of the extracellular domains of BCRs has been resolved, little is known about the configuration of the constituting four immunoglobulin domains spanning the membrane. Here, we modeled the structure of the transmembrane (TM) domain of the IgM B-cell receptor using self-assembly coarse-grained molecular dynamics simulations. The obtained quaternary structure was validated against available experimental data and atomistic simulations. The IgM-BCR-TM domain configuration shows a 1:1 stoichiometry between the homodimeric membrane-bound domain of IgM (mIgM) and a Ig-α/Ig-β heterodimer. The mIgM homodimer is based on an asymmetric association of two mIgM domains. We show that a specific site of the Ig-α/Ig-β heterodimer is responsible for the association of IgM-BCRs with lipid rafts. Our results further suggest that this site is blocked in small-sized IgM-BCR clusters. The BCR TM structure provides a molecular basis for the previously suggested dissociation activation model of B-cell receptors. Self-assembly molecular dynamics simulations at the coarse-grained scale here proved as a versatile tool in the study of receptor complexes.

Plant lipids: Key players of plasma membrane organization and function

The plasma membrane (PM) is the biological membrane that separates the interior of all cells from the outside. The PM is constituted of a huge diversity of proteins and lipids. In this review, we will update the diversity of molecular species of lipids found in plant PM. We will further discuss how lipids govern global properties of the plant PM, explaining that plant lipids are unevenly distributed and are able to organize PM in domains. From that observation, it emerges a complex picture showing a spatial and multiscale segregation of PM components. Finally, we will discuss how lipids are key players in the function of PM in plants, with a particular focus on plant-microbe interaction, transport and hormone signaling, abiotic stress responses, plasmodesmata function. The last chapter is dedicated to the methods that the plant membrane biology community needs to develop to get a comprehensive understanding of membrane organization in plants.

Phosphatidylinositol-3,5-bisphosphate lipid-binding-induced activation of the human two-pore channel 2

Mammalian two-pore channels (TPCs) are activated by the low-abundance membrane lipid phosphatidyl-(3,5)-bisphosphate (PI(3,5)P2) present in the endo-lysosomal system. Malfunction of human TPC1 or TPC2 (hTPC) results in severe organellar storage diseases and membrane trafficking defects. Here, we compared the lipid-binding characteristics of hTPC2 and of the PI(3,5)P2-insensitive TPC1 from the model plant Arabidopsis thaliana. Combination of simulations with functional analysis of channel mutants revealed the presence of an hTPC2-specific lipid-binding pocket mutually formed by two channel regions exposed to the cytosolic side of the membrane. We showed that PI(3,5)P2 is simultaneously stabilized by positively charged amino acids (K203, K204, and K207) in the linker between transmembrane helices S4 and S5 and by S322 in the cytosolic extension of S6. We suggest that PI(3,5)P2 cross links two parts of the channel, enabling their coordinated movement during channel gating.

Impact of Electric Fields on the Nanoscale Behavior of Lipid Monolayers at the Surface of Graphite in Solution

The nanoscale organization and dynamics of lipid molecules in self-assembled membranes is central to the biological function of cells and in the technological development of synthetic lipid structures as well as in devices such as biosensors. Here, we explore the nanoscale molecular arrangement and dynamics of lipids assembled in monolayers at the surface of highly ordered pyrolytic graphite (HOPG), in different ionic solutions, and under electrical potentials. Using a combination of atomic force microscopy and fluorescence recovery after photobleaching, we show that HOPG is able to support fully formed and fluid lipid membranes, but mesoscale order and corrugations can be observed depending on the type of the lipid considered (1,2-dioleoyl- sn-glycero-3-phosphocholine, 1,2-dioleoyl- sn-glycero-3-phospho-l-serine (DOPS), and 1,2-dioleoyl-3-trimethylammoniumpropane) and the ion present (Na⁺, Ca²⁺, Cl^-). Interfacial solvation forces and ion-specific effects dominate over the electrostatic changes induced by moderate electric fields (±1.0 V vs Ag/AgCl reference electrode) with particularly marked effects in the presence of calcium, and for DOPS. Our results provide insights into the interplay between the molecular, ionic, and electrostatic interactions and the formation of dynamical ordered structures in fluid lipid membranes.

Molecular Mechanism of Lipid Nanodisk Formation by Styrene-Maleic Acid Copolymers

Experimental characterization of membrane proteins often requires solubilization. A recent approach is to use styrene-maleic acid (SMA) copolymers to isolate membrane proteins in nanometer-sized membrane disks, or so-called SMA lipid particles (SMALPs). The approach has the advantage of allowing direct extraction of proteins, keeping their native lipid environment. Despite the growing popularity of using SMALPs, the molecular mechanism behind the process remains poorly understood. Here, we unravel the molecular details of the nanodisk formation by using coarse-grained molecular dynamics simulations. We show how SMA copolymers bind to the lipid bilayer interface, driven by the hydrophobic effect. Due to the concerted action of multiple adsorbed copolymers, large membrane defects appear, including small, water-filled pores. The copolymers can stabilize the rim of these pores, leading to pore growth and membrane disruption. Although complete solubilization is not seen on the timescale of our simulations, self-assembly experiments show that small nanodisks are the thermodynamically preferred end state. Our findings shed light on the mechanism of SMALP formation and on their molecular structure. This can be an important step toward the design of optimized extraction tools for membrane protein research.

Hydrodynamics of Diffusion in Lipid Membrane Simulations

By performing molecular dynamics simulations with up to 132 million coarse-grained particles in half-micron sized boxes, we show that hydrodynamics quantitatively explains the finite-size effects on diffusion of lipids, proteins, and carbon nanotubes in membranes. The resulting Oseen correction allows us to extract infinite-system diffusion coefficients and membrane surface viscosities from membrane simulations despite the logarithmic divergence of apparent diffusivities with increasing box width. The hydrodynamic theory of diffusion applies also to membranes with asymmetric leaflets and embedded proteins, and to a complex plasma-membrane mimetic.

Lipid-Protein Interactions Are Unique Fingerprints for Membrane Proteins

Cell membranes contain hundreds of different proteins and lipids in an asymmetric arrangement. Our current understanding of the detailed organization of cell membranes remains rather elusive, because of the challenge to study fluctuating nanoscale assemblies of lipids and proteins with the required spatiotemporal resolution. Here, we use molecular dynamics simulations to characterize the lipid environment of 10 different membrane proteins. To provide a realistic lipid environment, the proteins are embedded in a model plasma membrane, where more than 60 lipid species are represented, asymmetrically distributed between the leaflets. The simulations detail how each protein modulates its local lipid environment in a unique way, through enrichment or depletion of specific lipid components, resulting in thickness and curvature gradients. Our results provide a molecular glimpse of the complexity of lipid-protein interactions, with potentially far-reaching implications for our understanding of the overall organization of real cell membranes.

A program for automatically predicting supramolecular aggregates and its application to urea and porphin

Not only the molecular structure but also the presence or absence of aggregates determines many properties of organic materials. Theoretical investigation of such aggregates requires the prediction of a suitable set of diverse structures. Here, we present the open-source program EnergyScan for the unbiased prediction of geometrically diverse sets of small aggregates. Its bottom-up approach is complementary to existing ones by performing a detailed scan of an aggregate's potential energy surface, from which diverse local energy minima are selected. We crossvalidate this approach by predicting both literature-known and heretofore unreported geometries of the urea dimer. We also predict a diverse set of dimers of the less intensely studied case of porphin, which we investigate further using quantum chemistry. For several dimers, we find strong deviations from a reference absorption spectrum, which we explain using computed transition densities. This proof of principle clearly shows that EnergyScan successfully predicts aggregates exhibiting large structural and spectral diversity. © 2018 Wiley Periodicals, Inc.

Statistically derived asymmetric membrane potentials from α-helical and β-barrel membrane proteins

Modeling membrane protein (MP) folding, insertion, association and their interactions with other proteins, lipids, and drugs requires accurate transfer free energies (TFEs). Various TFE scales have been derived to quantify the energy required or released to insert an amino acid or protein into the membrane. Experimental measurement of TFEs is challenging, and only few scales were extended to depth-dependent energetic profiles. Statistical approaches can be used to derive such potentials; however, this requires a sufficient number of MP structures. Furthermore, MPs are tightly coupled to bilayers that are heterogeneous in terms of lipid composition, asymmetry, and protein content between organisms and organelles. Here we derived asymmetric implicit membrane potentials from β-barrel and α-helical MPs and use them to predict topology, depth and orientation of proteins in the membrane. Our data confirm the ‘charge-outside’ and ‘positive-inside’ rules for β-barrels and α-helical proteins, respectively. We find that the β-barrel profiles have greater asymmetry than the ones from α-helical proteins, as a result of the different membrane architecture of gram-negative bacterial outer membranes and the existence of lipopolysaccharide in the outer leaflet. Our data further suggest that pore-facing residues in β-barrels have a larger contribution to membrane insertion and stability than previously suggested.

Closely related, yet unique: Distinct homo- and heterodimerization patterns of G protein coupled chemokine receptors and their fine-tuning by cholesterol

Chemokine receptors, a subclass of G protein coupled receptors (GPCRs), play essential roles in the human immune system, they are involved in cancer metastasis as well as in HIV-infection. A plethora of studies show that homo- and heterodimers or even higher order oligomers of the chemokine receptors CXCR4, CCR5, and CCR2 modulate receptor function. In addition, membrane cholesterol affects chemokine receptor activity. However, structural information about homo- and heterodimers formed by chemokine receptors and their interplay with cholesterol is limited. Here, we report homo- and heterodimer configurations of the chemokine receptors CXCR4, CCR5, and CCR2 at atomistic detail, as obtained from thousands of molecular dynamics simulations. The observed homodimerization patterns were similar for the closely related CC chemokine receptors, yet they differed significantly between the CC receptors and CXCR4. Despite their high sequence identity, cholesterol modulated the CC homodimer interfaces in a subtype-specific manner. Chemokine receptor heterodimers display distinct dimerization patterns for CXCR4/CCR5 and CXCR4/CCR2. Furthermore, associations between CXCR4 and CCR5 reveal an increased cholesterol-sensitivity as compared to CXCR4/CCR2 heterodimerization patterns. This work provides a first comprehensive structural overview over the complex interaction network between chemokine receptors and indicates how heterodimerization and the interaction with the membrane environment diversifies the function of closely related GPCRs.

Exciton Localization on Ru-Based Photosensitizers Induced by Binding to Lipid Membranes

The characterization of electronic properties of metal complexes embedded in membrane environments is of paramount importance to develop efficient photosensitizers in optogenetic applications. Molecular dynamics and QM/MM simulations together with quantitative wave function analysis reveal a directional electronic redistribution of the exciton formed upon excitation of [Ru(bpy)₂(bpy-C17)]²⁺ when going from water to a lipid bilayer, despite the fact that the media influence neither the metal-to-ligand charge-transfer character nor the excitation energy of the absorption spectra. When the photosensitizer is embedded into the DOPC lipid membrane, exciton population is mainly located in the bypyridyl sites proximal to the positively charged surface of the bilayer due to electrostatic interactions. This behavior shows that the electronic structure of metal complexes can be controlled through the binding to external species, underscoring the crucial role of the environment in directing the electronic flow upon excitation and thus helping rational tuning of optogenetic agents.

Prediction of Thylakoid Lipid Binding Sites on Photosystem II

The thylakoid membrane has a unique lipid composition, consisting mostly of galactolipids. These thylakoid lipids have important roles in photosynthesis. Here, we investigate to what extent these lipids bind specifically to the Photosystem II complex. To this end, we performed coarse-grain MD simulations of the Photosystem II complex embedded in a thylakoid membrane with realistic composition. Based on >85 μs simulation time, we find that monogalactosyldiacylglycerol and sulfoquinovosyldiacylglycerol lipids are enriched in the annular shell around the protein, and form distinct binding sites. From the analysis of residue contacts, we conclude that electrostatic interactions play an important role in stabilizing these binding sites. Furthermore, we find that chlorophyll a has a prevalent role in the coordination of the lipids. In addition, we observe lipids to diffuse in and out of the plastoquinone exchange cavities, allowing exchange of cocrystallized lipids with the bulk membrane and suggesting a more open nature of the plastoquinone exchange cavity. Together, our data provide a wealth of information on protein-lipid interactions for a key protein in photosynthesis.

Scalable and fast heterogeneous molecular simulation with predictive parallelization schemes

Multiscale and inhomogeneous molecular systems are challenging topics in the field of molecular simulation. In particular, modeling biological systems in the context of multiscale simulations and exploring material properties are driving a permanent development of new simulation methods and optimization algorithms. In computational terms, those methods require parallelization schemes that make a productive use of computational resources for each simulation and from its genesis. Here, we introduce the heterogeneous domain decomposition approach, which is a combination of an heterogeneity-sensitive spatial domain decomposition with an a priori rearrangement of subdomain walls. Within this approach, the theoretical modeling and scaling laws for the force computation time are proposed and studied as a function of the number of particles and the spatial resolution ratio. We also show the new approach capabilities, by comparing it to both static domain decomposition algorithms and dynamic load-balancing schemes. Specifically, two representative molecular systems have been simulated and compared to the heterogeneous domain decomposition proposed in this work. These two systems comprise an adaptive resolution simulation of a biomolecule solvated in water and a phase-separated binary Lennard-Jones fluid.

Molecular dynamics simulation strategies for designing carbon-nanotube-based targeted drug delivery

The carbon nanotube (CNT)-based target-specific delivery of drugs, or other molecular cargo, has emerged as one of the most promising biomedical applications of nanotechnology. To achieve efficient CNT-based drug delivery, the interactions between the drug, CNT and biomolecular target need to be properly optimized. Recent advances in the computer-aided molecular design tools, in particular molecular dynamics (MD) simulation studies, offer an appropriate low-cost approach for such optimization. This review highlights the various potential MD approaches for the simulation of CNT interactions with cell membranes while emphasizing various methods of cellular internalization and toxicities of CNTs to build new strategies for designing rational CNT-based targeted drug delivery to circumvent the limitations associated with the various clinically available nonspecific therapeutic agents.

Efficient potential of mean force calculation from multiscale simulations: Solute insertion in a lipid membrane

The determination of potentials of mean force for solute insertion in a lipid membrane by means of all-atom molecular dynamics simulations is often hampered by sampling issues. Recently, a multiscale method has been proposed to leverage the conformational ensemble of a lower-resolution model as starting point for higher resolution simulations. In this work, we analyze the efficiency of this method by comparing its predictions for propanol insertion into a lipid membrane against conventional atomistic umbrella sampling simulation results. The multiscale approach is confirmed to provide accurate results with a gain of one order of magnitude in computational time. We then investigate the role of the coarse-grained representation. We find that the accuracy of the results is tightly connected to the presence of a good configurational overlap between the coarse-grained and atomistic models-a general requirement when developing multiscale simulation methods.

Membrane phase transition during heating and cooling: molecular insight into reversible melting

With increasing temperature, lipid bilayers undergo a gel-fluid phase transition, which plays an essential role in many physiological phenomena. In the present work, this first-order phase transition was investigated for variable heating and cooling rates for a dipalmitoylphosphatidylcholine (DPPC) lipid bilayer by means of atomistic molecular dynamics simulations. Alternative methods to track the melting temperature [Formula: see text] are compared. The resulting [Formula: see text] is shown to be independent of the scan rate for small heating rates (0.05-0.3 K/ns) implying reversible melting, and increases for larger heating (0.3-4 K/ns) or cooling rates (2-0.1 K/ns). The reported dependency of the melting temperature on the heating rate is in perfect agreement with a two-state kinetic rate model as suggested previously. Expansion and shrinkage, as well as the dynamics of melting seeds is described. The simulations further exhibit a relative shift between melting seeds in opposing membrane leaflets as predicted from continuum elastic theory.

Unperturbed hydrocarbon chains and liquid phase bilayer lipid chains: a computer simulation study

In this work, the properties of saturated and unsaturated fatty acid acyl chains 16:0, 18:0, 18:1(n-9)cis, 18:2(n-6)cis, 18:3(n-3)cis, 18:4(n-3)cis, 18:5(n-3)cis, 20:4(n-6)cis, 20:5(n-3)cis and 22:6(n-3)cis in a bilayer liquid crystalline state and similar hydrocarbon chains (with CH[Formula: see text] terminal groups instead of C=O groups) in the unperturbed state characterised by a lack of long-range interaction were investigated. The unperturbed hydrocarbon chains were modelled by Monte Carlo simulations at temperature [Formula: see text] K; sixteen fully hydrated homogeneous liquid crystalline phosphatidylcholine bilayers containing these chains were studied by molecular dynamics simulations at the same temperature. To eliminate effects of the simulation parameters, the molecular dynamics and Monte Carlo simulations were carried out using the same structural data and force field coefficients. From these computer simulations, the average distances between terminal carbon atoms of the chains (end-to-end distances) were calculated and compared. The trends in the end-to-end distances obtained for the unperturbed chains were found to be qualitatively similar to those obtained for the same lipid chains in the bilayers. So, for understanding of a number of processes in biological membranes (e.g., changes in fatty acid composition caused by environmental changes such as temperature and pressure), it is possible to use, at least as a first approximation, the relationships between the structure and properties for unperturbed or isolated hydrocarbon chains.

Investigating the Structure of Multicomponent Gel-Phase Lipid Bilayers

Single- and multicomponent lipid bilayers of 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC), 1,2-distearoyl-sn-glycero-3-phosphatidylcholine (DSPC), isostearyl isostearate, and heptadecanoyl heptadecanoate in the gel phase are studied via molecular dynamics simulations. It is shown that the structural properties of multicomponent bilayers can deviate strongly from the structures of their single-component counterparts. Specifically, the lipid mixtures are shown to adopt a compact packing by offsetting the positioning depths at which different lipid species are located in the bilayer. This packing mechanism affects the area per lipid, the bilayer height, and the chain tilt angles and has important consequences for other bilayer properties, such as interfacial hydrogen bonding and bilayer permeability. In particular, the simulations suggest that bilayers containing isostearyl isostearate or heptadecanoyl heptadecanoate are less permeable than pure 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine or DSPC bilayers. Furthermore, hydrogen-bond analysis shows that the residence times of lipid-water hydrogen bonds depend strongly on the bilayer composition, with longer residence times for bilayers that have a higher DSPC content. The findings illustrate and explain the fundamental differences between the properties of single- and multicomponent bilayers.

The Multifaceted Role of SNARE Proteins in Membrane Fusion

Membrane fusion is a key process in all living organisms that contributes to a variety of biological processes including viral infection, cell fertilization, as well as intracellular transport, and neurotransmitter release. In particular, the various membrane-enclosed compartments in eukaryotic cells need to exchange their contents and communicate across membranes. Efficient and controllable fusion of biological membranes is known to be driven by cooperative action of SNARE proteins, which constitute the central components of the eukaryotic fusion machinery responsible for fusion of synaptic vesicles with the plasma membrane. During exocytosis, vesicle-associated v-SNARE (synaptobrevin) and target cell-associated t-SNAREs (syntaxin and SNAP-25) assemble into a core trans-SNARE complex. This complex plays a versatile role at various stages of exocytosis ranging from the priming to fusion pore formation and expansion, finally resulting in the release or exchange of the vesicle content. This review summarizes current knowledge on the intricate molecular mechanisms underlying exocytosis triggered and catalyzed by SNARE proteins. Particular attention is given to the function of the peptidic SNARE membrane anchors and the role of SNARE-lipid interactions in fusion. Moreover, the regulatory mechanisms by synaptic auxiliary proteins in SNARE-driven membrane fusion are briefly outlined.

Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are central to many fundamental cellular signaling pathways. They transduce signals from the outside to the inside of cells in physiological processes ranging from vision to immune response. It is extremely challenging to look at them individually using conventional experimental techniques. Recently, a pseudo atomistic molecular model has emerged as a valuable tool to access information on GPCRs, more specifically on their interactions with their environment in their native cell membrane and the consequences on their supramolecular organization. This approach uses the Martini coarse grain (CG) model to describe the receptors, lipids, and solvent in molecular dynamics (MD) simulations and in enough detail to allow conserving the chemical specificity of the different molecules. The elimination of unnecessary degrees of freedom has opened up large-scale simulations of the lipid-mediated supramolecular organization of GPCRs. Here, after introducing the Martini CGMD method, we review these studies carried out on various members of the GPCR family, including rhodopsin (visual receptor), opioid receptors, adrenergic receptors, adenosine receptors, dopamine receptor, and sphingosine 1-phosphate receptor. These studies have brought to light an interesting set of novel biophysical principles. The insights range from revealing localized and heterogeneous deformations of the membrane bilayer at the surface of the protein, specific interactions of lipid molecules with individual GPCRs, to the effect of the membrane matrix on global GPCR self-assembly. The review ends with an overview of the lessons learned from the use of the CGMD method, the biophysical-chemical findings on lipid-protein interplay.

A Coiled-Coil Peptide Shaping Lipid Bilayers upon Fusion

A system based on two designed peptides, namely the cationic peptide K, (KIAALKE)3, and its complementary anionic counterpart called peptide E, (EIAALEK)3, has been used as a minimal model for membrane fusion, inspired by SNARE proteins. Although the fact that docking of separate vesicle populations via the formation of a dimeric E/K coiled-coil complex can be rationalized, the reasons for the peptides promoting fusion of vesicles cannot be fully explained. Therefore it is of significant interest to determine how the peptides aid in overcoming energetic barriers during lipid rearrangements leading to fusion. In this study, investigations of the peptides' interactions with neutral PC/PE/cholesterol membranes by fluorescence spectroscopy show that tryptophan-labeled K∗ binds to the membrane (KK∗ ∼6.2 103 M-1), whereas E∗ remains fully water-solvated. 15N-NMR spectroscopy, depth-dependent fluorescence quenching, CD-spectroscopy experiments, and MD simulations indicate a helix orientation of K∗ parallel to the membrane surface. Solid-state 31P-NMR of oriented lipid membranes was used to study the impact of peptide incorporation on lipid headgroup alignment. The membrane-immersed K∗ is found to locally alter the bilayer curvature, accompanied by a change of headgroup orientation relative to the membrane normal and of the lipid composition in the vicinity of the bound peptide. The NMR results were supported by molecular dynamics simulations, which showed that K reorganizes the membrane composition in its vicinity, induces positive membrane curvature, and enhances the lipid tail protrusion probability. These effects are known to be fusion relevant. The combined results support the hypothesis for a twofold role of K in the mechanism of membrane fusion: 1) to bring opposing membranes into close proximity via coiled-coil formation and 2) to destabilize both membranes thereby promoting fusion.

Membrane-Mediated Oligomerization of G Protein Coupled Receptors and Its Implications for GPCR Function

The dimerization or even oligomerization of G protein coupled receptors (GPCRs) causes ongoing, controversial debates about its functional role and the coupled biophysical, biochemical or biomedical implications. A continously growing number of studies hints to a relation between oligomerization and function of GPCRs and strengthens the assumption that receptor assembly plays a key role in the regulation of protein function. Additionally, progress in the structural analysis of GPCR-G protein and GPCR-ligand interactions allows to distinguish between actively functional and non-signaling complexes. Recent findings further suggest that the surrounding membrane, i.e., its lipid composition may modulate the preferred dimerization interface and as a result the abundance of distinct dimeric conformations. In this review, the association of GPCRs and the role of the membrane in oligomerization will be discussed. An overview of the different reported oligomeric interfaces is provided and their capability for signaling discussed. The currently available data is summarized with regard to the formation of GPCR oligomers, their structures and dependency on the membrane microenvironment as well as the coupling of oligomerization to receptor function.

Efficient preparation and analysis of membrane and membrane protein systems

Molecular dynamics (MD) simulations have become a highly important technique to consider lipid membrane systems, and quite often they provide considerable added value to laboratory experiments. Rapid development of both software and hardware has enabled the increase of time and size scales reachable by MD simulations to match those attainable by several accurate experimental techniques. However, until recently, the quality and maturity of software tools available for building membrane models for simulations as well as analyzing the results of these simulations have seriously lagged behind. Here, we discuss the recent developments of such tools from the end-users' point of view. In particular, we review the software that can be employed to build lipid bilayers and other related structures with or without embedded membrane proteins to be employed in MD simulations. Additionally, we provide a brief critical insight into force fields and MD packages commonly used for membrane and membrane protein simulations. Finally, we list analysis tools that can be used to study the properties of membrane and membrane protein systems. In all these points we comment on the respective compatibility of the covered tools. We also share our opinion on the current state of the available software. We briefly discuss the most commonly employed tools and platforms on which new software can be built. We conclude the review by providing a few ideas and guidelines on how the development of tools can be further boosted to catch up with the rapid pace at which the field of membrane simulation progresses. This includes improving the compatibility between software tools and promoting the openness of the codes on which these applications rely. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Encapsulated membrane proteins: A simplified system for molecular simulation

Over the past 50years there has been considerable progress in our understanding of biomolecular interactions at an atomic level. This in turn has allowed molecular simulation methods employing full atomistic modelling at ever larger scales to develop. However, some challenging areas still remain where there is either a lack of atomic resolution structures or where the simulation system is inherently complex. An area where both challenges are present is that of membranes containing membrane proteins. In this review we analyse a new practical approach to membrane protein study that offers a potential new route to high resolution structures and the possibility to simplify simulations. These new approaches collectively recognise that preservation of the interaction between the membrane protein and the lipid bilayer is often essential to maintain structure and function. The new methods preserve these interactions by producing nano-scale disc shaped particles that include bilayer and the chosen protein. Currently two approaches lead in this area: the MSP system that relies on peptides to stabilise the discs, and SMALPs where an amphipathic styrene maleic acid copolymer is used. Both methods greatly enable protein production and hence have the potential to accelerate atomic resolution structure determination as well as providing a simplified format for simulations of membrane protein dynamics. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Force Field Development for Lipid Membrane Simulations

With the rapid development of computer power and wide availability of modelling software computer simulations of realistic models of lipid membranes, including their interactions with various molecular species, polypeptides and membrane proteins have become feasible for many research groups. The crucial issue of the reliability of such simulations is the quality of the force field, and many efforts, especially in the latest several years, have been devoted to parametrization and optimization of the force fields for biomembrane modelling. In this review, we give account of the recent development in this area, covering different classes of force fields, principles of the force field parametrization, comparison of the force fields, and their experimental validation. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Membrane pore formation in atomistic and coarse-grained simulations

Biological cells and their organelles are protected by ultra thin membranes. These membranes accomplish a broad variety of important tasks like separating the cell content from the outer environment, they are the site for cell-cell interactions and many enzymatic reactions, and control the in- and efflux of metabolites. For certain physiological functions e.g. in the fusion of membranes and also in a number of biotechnological applications like gene transfection the membrane integrity needs to be compromised to allow for instance for the exchange of polar molecules across the membrane barrier. Mechanisms enabling the transport of molecules across the membrane involve membrane proteins that form specific pores or act as transporters, but also so-called lipid pores induced by external fields, stress, or peptides. Recent progress in the simulation field enabled to closely mimic pore formation as supposed to occur in vivo or in vitro. Here, we review different simulation-based approaches in the study of membrane pores with a focus on lipid pore properties such as their size and energetics, poration mechanisms based on the application of external fields, charge imbalances, or surface tension, and on pores that are induced by small molecules, peptides, and lipids. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.