Nanodomain Formation of Ganglioside GM1 in Lipid Membrane: Effects of Cholera Toxin-Mediated Cross-Linking

Specific interactions of peripheral membrane proteins with lipids: what can molecular simulations show us?

Peripheral membrane proteins (PMPs) can reversibly and specifically bind to biological membranes to carry out functions such as cell signalling, enzymatic activity, or membrane remodelling. Structures of these proteins and of their lipid-binding domains are typically solved in a soluble form, sometimes with a lipid or lipid headgroup at the binding site. To provide a detailed molecular view of PMP interactions with the membrane, computational methods such as molecular dynamics (MD) simulations can be applied. Here, we outline recent attempts to characterise these binding interactions, focusing on both intracellular proteins, such as phosphatidylinositol phosphate (PIP)-binding domains, and extracellular proteins such as glycolipid-binding bacterial exotoxins. We compare methods used to identify and analyse lipid-binding sites from simulation data and highlight recent work characterising the energetics of these interactions using free energy calculations. We describe how improvements in methodologies and computing power will help MD simulations to continue to contribute to this field in the future.

Structure and Formation Mechanism of Antimicrobial Peptides Temporin B- and L-Induced Tubular Membrane Protrusion

Temporins are a family of antimicrobial peptides (AMPs) isolated from frog skin, which are very short, weakly charged, and highly hydrophobic. They execute bactericidal activities in different ways from many other AMPs. This work investigated morphological changes of planar bilayer membranes composed of mixed zwitterionic and anionic phospholipids induced by temporin B and L (TB and TL) using all-atom and coarse-grained molecular dynamics simulations. We found that TB and TL fold to α-helices at the membrane surface and penetrate shallowly into the bilayer. These short AMPs have low propensity to induce membrane pore formation but possess high ability to extract lipids out. At relatively high peptide concentrations, the strong hydrophobicity of TB and TL promotes them to aggregate into clusters on the membrane surface. These aggregates attract a large amount of lipids out of the membrane to release compression induced by other dispersed peptides binding to the membrane. The extruded lipids mix evenly with the peptides in the cluster and form tubule-like protrusions. Certain water molecules follow the movement of lipids, which not only fill the cavities of the protrusion but also assist in maintaining the tubular structures. In contrast, the peptide-free leaflet remains intact. The present results unravel distinctive antimicrobial mechanisms of temporins disturbing membranes.

Pore Forming Protein Induced Biomembrane Reorganization and Dynamics: A Focused Review

Pore forming proteins are a broad class of pathogenic proteins secreted by organisms as virulence factors due to their ability to form pores on the target cell membrane. Bacterial pore forming toxins (PFTs) belong to a subclass of pore forming proteins widely implicated in bacterial infections. Although the action of PFTs on target cells have been widely investigated, the underlying membrane response of lipids during membrane binding and pore formation has received less attention. With the advent of superresolution microscopy as well as the ability to carry out molecular dynamics (MD) simulations of the large protein membrane assemblies, novel microscopic insights on the pore forming mechanism have emerged over the last decade. In this review, we focus primarily on results collated in our laboratory which probe dynamic lipid reorganization induced in the plasma membrane during various stages of pore formation by two archetypal bacterial PFTs, cytolysin A (ClyA), an α-toxin and listeriolysin O (LLO), a β-toxin. The extent of lipid perturbation is dependent on both the secondary structure of the membrane inserted motifs of pore complex as well as the topological variations of the pore complex. Using confocal and superresolution stimulated emission depletion (STED) fluorescence correlation spectroscopy (FCS) and MD simulations, lipid diffusion, cholesterol reorganization and deviations from Brownian diffusion are correlated with the oligomeric state of the membrane bound protein as well as the underlying membrane composition. Deviations from free diffusion are typically observed at length scales below ∼130 nm to reveal the presence of local dynamical heterogeneities that emerge at the nanoscale-driven in part by preferential protein binding to cholesterol and domains present in the lipid membrane. Interrogating the lipid dynamics at the nanoscale allows us further differentiate between binding and pore formation of β- and α-PFTs to specific domains in the membrane. The molecular insights gained from the intricate coupling that occurs between proteins and membrane lipids and receptors during pore formation are expected to improve our understanding of the virulent action of PFTs.

Multiscale Molecular Dynamics Studies Reveal Different Modes of Receptor Clustering by Gb3-Binding Lectins

The recognition of carbohydrate receptors on host cell membranes by pathogenic lectins is a crucial step in the microbial invasion. Two bacterial lectins, the B-subunit of Shiga toxin from Shigella dysenteria (StxB) and lectin I from Pseudomonas aeruginosa (LecA), are specific to the same galactolipid-globotriaosylceramide (Gb3). In this study we present a coarse-grained (cg) model of Gb3, which we further apply to unravel the molecular details of glycolipid binding by two lectins on the surface of a DOPC/cholesterol/Gb3 bilayer. In cg molecular dynamics simulations with time scales of dozens of microseconds, Gb3 was randomly distributed. The binding of both StxB or LecA is accompanied by Gb3 clustering in a cholesterol environment and with exclusion of DOPC in protein vicinity. StxB being bound by all 15 binding sites induced membrane bending, while LecA interacted with two out of four binding sites for most of the time causing a smaller inward curvature of the model membrane. Stable interactions occurred preferably when LecA was normal to the membrane surface. Furthermore, all-atom simulations revealed that LecA bound Gb3's headgroup at only one out of two possible conformations of the carbohydrate moiety observed at protein-free conditions. The results shed light on the mechanism of interactions between two lectins and Gb3 on the membrane surface and offer a coarse-grained model to study more complex systems at large spatiotemporal scales.

Single-lipid dynamics in phase-separated supported lipid bilayers

Phase separation is a fundamental organizing mechanism on cellular membranes. Lipid phases have complex dependencies on the membrane composition, curvature, tension, and temperature. Lipid diffusion rates vary by up to ten-fold between liquid-disordered (L_d) and liquid-ordered (L_o) phases depending on the membrane composition, measurement technique, and the surrounding environment. This manuscript reports the lipid diffusion on phase-separated supported lipid bilayers (SLBs) with varying temperature, composition, and lipid phase. Lipid diffusion is measured by single-particle tracking (SPT) and fluorescence correlation spectroscopy (FCS) via custom data acquisition and analysis protocols that apply to diverse membranes systems. Traditionally, SPT is sensitive to diffuser aggregation, whereas the diffusion rates reported by FCS are unaffected by the presence of immobile aggregates. Within this manuscript, we report (1) improved single-particle tracking analysis of lipid diffusion, (2) comparison and consistency between diffusion measurement methods for non-Brownian diffusers, and (3) the application of these methods to measure the phase, temperature, and composition dependencies in lipid diffusion. We demonstrate improved SPT analysis methods that yield consistent FCS and SPT diffusion results even when most fluorescent lipids are frequently confined within aggregates within the membrane. With varying membrane composition and temperature, we demonstrate differences in diffusion between the L_d and L_o phases of SLBs.

Bringing rafts to life: Lessons learned from lipid organization across diverse biological membranes

The ability of lipids to drive lateral organization is a remarkable feature of membranes and has been hypothesized to underlie the architecture of cells. Models for lipid rafts and related domains were originally based on the mammalian plasma membrane, but the nature of heterogeneity in this system is still not fully resolved. However, the concept of lipid-driven organization has been highly influential across biology, and has led to discoveries in organisms that feature a diversity of lipid chemistries and physiological needs. Here we review several emerging and instructive cases of membrane organization in non-mammalian systems. In bacteria, several types of membrane domains that act in metabolism and signaling have been elucidated. These widen our view of what constitutes a raft, but also introduce new questions about the relationship between organization and function. In yeast, observable membrane organization is found in both the plasma membrane and the vacuole. The latter serves as the best example of classic membrane phase partitioning in a living system to date, suggesting that internal organelles are important membranes to investigate across eukaryotes. Finally, we highlight plants as powerful model systems for complex membrane interactions in multicellular organisms. Plant membranes are organized by unique glycosphingolipids, supporting the importance of carbohydrate interactions in organizing lateral domains. These examples demonstrate that membrane organization is a potentially universal phenonenon in biology and argue for the continued broadening of lipid physical chemistry research into a wide range of systems.

Structured clustering of the glycosphingolipid GM1 is required for membrane curvature induced by cholera toxin

AB₅ bacterial toxins and polyomaviruses induce membrane curvature as a mechanism to facilitate their entry into host cells. How membrane bending is accomplished is not yet fully understood but has been linked to the simultaneous binding of the pentameric B subunit to multiple copies of glycosphingolipid receptors. Here, we probe the toxin membrane binding and internalization mechanisms by using a combination of superresolution and polarized localization microscopy. We show that cholera toxin subunit B (CTxB) can induce membrane curvature only when bound to multiple copies of its glycosphingolipid receptor, GM1, and the ceramide structure of GM1 is likely not a determinant of this activity as assessed in model membranes. A mutant CTxB capable of binding only a single GM1 fails to generate curvature either in model membranes or in cells, and clustering the mutant CTxB-single-GM1 complexes by antibody cross-linking does not rescue the membrane curvature phenotype. We conclude that both the multiplicity and specific geometry of GM1 binding sites are necessary for the induction of membrane curvature. We expect this to be a general rule of membrane behavior for all AB₅ toxins and polyomaviruses that bind glycosphingolipids to invade host cells.

Interleaflet Coupling of Lipid Nanodomains - Insights From in vitro Systems

The plasma membrane is a complex system, consisting of two layers of lipids and proteins compartmentalized into small structures called nanodomains. Despite the asymmetric composition of both leaflets, coupling between the layers is surprisingly strong. This can be evidenced, for example, by recent experimental studies performed on phospholipid giant unilamellar vesicles showing that nanodomains formed in the outer layer are perfectly registered with those in the inner leaflet. Similarly, microscopic phase separation in one leaflet can induce phase separation in the opposing leaflet that would otherwise be homogeneous. In this review, we summarize the current theoretical and experimental knowledge that led to the current view that domains are – irrespective of their size – commonly registered across the bilayer. Mechanisms inducing registration of nanodomains suggested by theory and calculations are discussed. Furthermore, domain coupling is evidenced by experimental studies based on the sparse number of methods that can resolve registered from independent nanodomains. Finally, implications that those findings using model membrane studies might have for cellular membranes are discussed.

Peripheral Antimicrobial Peptide Gomesin Induces Membrane Protrusion, Folding, and Laceration

Optical microscopy shows that the peripheral antimicrobial peptide (AMP) gomesin does not disrupt the bacterial membrane by forming stable transmembrane pores but induces lipid accumulation domains, which is followed by a sudden burst near the domains. The molecular action mechanisms of gomesin on vesicle and planar bilayer membranes are investigated in this work using coarse-grained molecular dynamics simulations. By comparing the membrane morphology and property changes induced by gomesin and the pore-forming AMP melittin, we determined that the amphiphilic shape of the AMPs is a key factor affecting the mechanism of cell death. The binding of wedge-shaped gomesin, with a small hydrophobic surface, onto the membrane induces protrusion and folding of the outer monolayer followed by sudden membrane lacerations at the axillae of the protuberances. Alternatively, cylinder-shaped melittins with comparable hydrophilic and hydrophobic surfaces destroy membranes by forming stable pores coexisting with exocytosis-like buddings and endocytosis-like invaginations. The multiple actions of AMPs on the bacterial membrane suggest diverse paradigms for designing molecular carriers for delivering drugs to the cell.

Role of Transmembrane Proteins for Phase Separation and Domain Registration in Asymmetric Lipid Bilayers

It is well known that the formation and spatial correlation of lipid domains in the two apposed leaflets of a bilayer are influenced by weak lipid–lipid interactions across the bilayer’s midplane. Transmembrane proteins span through both leaflets and thus offer an alternative domain coupling mechanism. Using a mean-field approximation of a simple bilayer-type lattice model, with two two-dimensional lattices stacked one on top of the other, we explore the role of this “structural” inter-leaflet coupling for the ability of a lipid membrane to phase separate and form spatially correlated domains. We present calculated phase diagrams for various effective lipid–lipid and lipid–protein interaction strengths in membranes that contain a binary lipid mixture in each leaflet plus a small amount of added transmembrane proteins. The influence of the transmembrane nature of the proteins is assessed by a comparison with “peripheral” proteins, which result from the separation of one single integral protein into two independent units that are no longer structurally connected across the bilayer. We demonstrate that the ability of membrane-spanning proteins to facilitate domain formation requires sufficiently strong lipid–protein interactions. Weak lipid–protein interactions generally tend to inhibit phase separation in a similar manner for transmembrane as for peripheral proteins.

Dissipative Particle Dynamics Simulations of a Protein-Directed Self-Assembly of Nanoparticles

Design and fabrication of multifunctional porous structures play key roles in the development of high-performance energy storage devices. Our experiments demonstrated that nanostructured porous components, such as electrodes and interlayers, generated from the protein-directed self-assembly of nanoparticles can significantly improve the battery performances. The protein-directed assembly of nanoparticles in solution is a complex process involving the complicated interactions among proteins, particles, and solvent molecules. In this paper, we investigate the effects of coating proteins and specific solvent environments on the assembled porous structures. Comprehensive dissipative particle dynamics (DPD) simulations have been implemented to explore the molecular interactions and uncover the fundamental mechanisms in a gelatin-directed self-assembly of carbon black particles under different solvent conditions. Our simulations show that compact triple-strand "rod-like" structures are formed in water while loose curved "sheet-like" structures are formed in an acetic acid/water mixture. The structural difference is mainly due to the redistribution of the charges on the gelatin side chains under specific acid-solvent conditions. The strong and flexible "sheet-like" structures lead to a homogenous porous structure with high porosity and with large functionalized surfaces. Our simulations results can reasonably explain the experimental observations; this work demonstrates the great potential of DPD as a powerful tool in guiding future experimental design and optimization.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Binding of SV40's Viral Capsid Protein VP1 to Its Glycosphingolipid Receptor GM1 Induces Negative Membrane Curvature: A Molecular Dynamics Study

The binding of the pentameric capsid protein VP1 of simian virus 40 to its glycosphingolipid receptor GM1 is a key step for the entry of the virus into the host cell. Recent experimental studies have shown that the interaction of variants of soluble VP1 pentamers with giant unilamellar vesicles composed of GM1, DOPC, and cholesterol leads to the formation of tubular membrane invaginations to the inside of the vesicles, mimicking the initial steps of endocytosis. We have used coarse-grained and atomistic molecular dynamics (MD) simulations to study the interaction of VP1 with GM1/DOPC/cholesterol bilayers. In the presence of one VP1 protein, we monitor the formation of small local negative curvature and membrane thinning at the protein binding site as well as reduction of area per lipid. These membrane deformations are also observed under cholesterol-free conditions. However, here, the number of GM1 molecules attached to the VP1 binding pockets increases. The membrane curvature is slightly increased for asymmetric GM1 distribution that mimics conditions in vivo, compared to symmetric GM1 distributions which are often applied in experiments. Slightly smaller inward curvature was observed in atomistic control simulations. Binding of four VP1 proteins leads to an increase of the average intrinsic area per lipid in the protein binding leaflet. Membrane fluctuations appear to be the driving force of VP1 aggregation, as was previously shown for membrane-adhering particles because no VP1 aggregation is observed in the absence of a lipid membrane.

Implicit-solvent dissipative particle dynamics force field based on a four-to-one coarse-grained mapping scheme

A new set of efficient solvent-free dissipative particle dynamics (DPD) force fields was developed for phospholipids and peptides. To enhance transferability, this model maps around four heavy atoms and their connected hydrogen atoms into a coarse-grained elementary bead based on functional group. The effective hybrid potential between any pair of beads is composed of a short-range repulsive soft-core potential that directly adopts the form of an explicit-solvent DPD model and a long-range attractive hydrophobic potential. The parameters of the attractive potentials for lipid molecules were obtained by fitting the explicit-solvent DPD simulation of one bead of any type in a water box, then finely tuning it until the bilayer membrane properties obtained in the explicit-solvent model were matched. These parameters were further extended to amino acids according to bead type. The structural and elastic properties of bilayer membranes, free energy profiles for a lipid flip-flop and amino acid analogues translocating across the membrane, and membrane pore formation induced by antimicrobial peptides obtained from this solvent-free DPD force field considerably agreed with the explicit-solvent DPD results. Importantly, the efficiency of this method is guaranteed to accelerate the assembly of vesicles composed of several thousand lipids by up to 50-fold, rendering the experimental liposome dynamics as well as membrane-peptide interactions feasible at accessible computational expense.

Surface Chemistry and Spectroscopic Study of a Cholera Toxin B Langmuir Monolayer

In this article, we explored the surface chemistry properties of a cholera toxin B (CTB) monolayer at the air-subphase interface and investigated the change in interfacial properties through in situ spectroscopy. The study showed that the impact of the blue shift was negligible, suggesting that the CTB molecules were minimally affected by the subphase molecules. The stability of the CTB monolayer was studied by maintaining the constant surface pressure for a long time and also by using the compression-decompression cycle experiments. The high stability of the Langmuir monolayer of CTB clearly showed that the driving force of CTB going to the amphiphilic membrane was its amphiphilic nature. In addition, no major change was detected in the various in situ spectroscopy results (such as UV-vis, fluorescence, and IR ER) of the CTB Langmuir monolayer with the increase in surface pressure. This indicates that no aggregation occurs in the Langmuir monolayer of CTB.

Nanoscale Membrane Budding Induced by CTxB and Detected via Polarized Localization Microscopy

For endocytosis and exocytosis, membranes transition among planar, budding, and vesicular topographies through nanoscale reorganization of lipids, proteins, and carbohydrates. However, prior attempts to understand the initial stages of nanoscale bending have been limited by experimental resolution. Through the implementation of polarized localization microscopy, this article reports the inherent membrane bending capability of cholera toxin subunit B (CTxB) in quasi-one-component-supported lipid bilayers. Membrane buds were first detected with <50 nm radius, grew to >200 nm radius, and extended into longer tubules with dependence on the membrane tension and CTxB concentration. Compared to the concentration of the planar-supported lipid bilayers, CTxB was (12 ± 4)× more concentrated on the positive curvature top and (26 ± 11)× more concentrated on the negative Gaussian curvature neck of the nanoscale membrane buds. CTxB is frequently used as a marker for liquid-ordered lipid phases; however, the coupling between CTxB and membrane bending provides an alternate understanding of CTxB-induced membrane reorganization. These findings allow for the reinterpretation of prior observations by correlating CTxB clustering and diffusion to CTxB-induced membrane bending. Single-particle tracking was performed on single lipids and CTxB to reveal the correlations among single-molecule diffusion, CTxB accumulation, and membrane topography. Slowed lipid and CTxB diffusion was observed at the nanoscale bud locations, suggesting a local increase in the effective membrane viscosity or molecular crowding upon membrane bending. These results suggest inherent CTxB-induced membrane bending as a mechanism for initiating CTxB internalization in cells that could be independent of clathrin, caveolin, actin, and lipid phase separation.

Cholera toxin B subunit induces local curvature on lipid bilayers

The B subunit of the bacterial cholera toxin (CTxB) is responsible for the toxin binding to the cell membrane and its intracellular trafficking. CTxB binds to the monosialotetrahexosyl ganglioside at the plasma membrane of the target cell and mediates toxin internalization by endocytosis. CTxB induces a local membrane curvature that is essential for its clathrin-independent uptake. Using all-atom molecular dynamics, we show that CTxB induces local curvature, with the radius of curvature around 36 nm. The main feature of the CTxB molecular structure that causes membrane bending is the protruding alpha helices in the middle of the protein. Our study points to a generic protein design principle for generating local membrane curvature through specific binding to their lipid anchors.

GM1 Softens POPC Membranes and Induces the Formation of Micron-Sized Domains

The influence of the glycolipid GM1 on the physical properties of POPC membranes was studied systematically by using different methods applied to giant and large unilamellar vesicles. The charge per GM1 molecule in the membrane was estimated from electrophoretic mobility measurements. Optical microscopy and differential scanning calorimetry were employed to construct a partial phase diagram of the GM1/POPC system. At room temperature, phase separation in the membrane was detected for GM1 fractions at and above ∼5 mol %, whereby GM1-rich gel-like domains were observed by fluorescent microscopy. Fluctuation analysis, vesicle electrodeformation, and micropipette aspiration were used to assess the bending rigidity of the membrane as a function of GM1 content. In the fluid phase, GM1 was shown to strongly soften the bilayer. In the region of coexistence of fluid and gel-like domains, the micropipette aspiration technique allowed measurements of the bending rigidity of the fluid phase only, whereas electrodeformation and fluctuation analysis were affected by the presence of the gel-phase domains. The observation that GM1 decreased the bilayer bending rigidity is important for understanding the role of this ganglioside in the flexibility of neuronal membranes.

Interleaflet Coupling, Pinning, and Leaflet Asymmetry-Major Players in Plasma Membrane Nanodomain Formation

The plasma membrane has a highly asymmetric distribution of lipids and contains dynamic nanodomains many of which are liquid entities surrounded by a second, slightly different, liquid environment. Contributing to the dynamics is a continuous repartitioning of components between the two types of liquids and transient links between lipids and proteins, both to extracellular matrix and cytoplasmic components, that temporarily pin membrane constituents. This make plasma membrane nanodomains exceptionally challenging to study and much of what is known about membrane domains has been deduced from studies on model membranes at equilibrium. However, living cells are by definition not at equilibrium and lipids are distributed asymmetrically with inositol phospholipids, phosphatidylethanolamines and phosphatidylserines confined mostly to the inner leaflet and glyco- and sphingolipids to the outer leaflet. Moreover, each phospholipid group encompasses a wealth of species with different acyl chain combinations whose lateral distribution is heterogeneous. It is becoming increasingly clear that asymmetry and pinning play important roles in plasma membrane nanodomain formation and coupling between the two lipid monolayers. How asymmetry, pinning, and interdigitation contribute to the plasma membrane organization is only beginning to be unraveled and here we discuss their roles and interdependence.

Glycolipid Crosslinking Is Required for Cholera Toxin to Partition Into and Stabilize Ordered Domains

Current models of lipid rafts propose that lipid domains exist as nanoscale compositional fluctuations and these fluctuations can potentially be stabilized into larger domains, consequently better compartmentalizing cellular functions. However, the mechanisms governing stabilized raft assembly and function remain unclear. Here, we test the role of glycolipid crosslinking as a raft targeting and ordering mechanism using the well-studied raft marker cholera toxin B pentamer (CTxB) that binds up to five GM1 glycosphingolipids to enter host cells. We show that when applied to cell-derived giant plasma membrane vesicles, a variant of CTxB containing only a single functional GM1 binding site exhibits significantly reduced partitioning to the ordered phase compared to wild-type CTxB with five binding sites. Moreover, monovalent CTxB does not stabilize membrane domains, unlike wild-type CTxB. These results support the long-held hypothesis that CTxB stabilizes raft domains via a lipid crosslinking mechanism and establish a role for crosslinking in the partitioning of CTxB to ordered domains.

Multi-scale molecular dynamics study of cholera pentamer binding to a GM1-phospholipid membrane

The AB5 type toxin produced by the Vibrio cholerae bacterium is the causative agent of the cholera disease. The cholera toxin (CT) has been shown to bind specifically to GM1 glycolipids on the membrane surface. This binding of CT to the membrane is the initial step in its endocytosis and has been postulated to cause significant disruption to the membrane structure. In this work, we have carried out a combination of coarse-grain and atomistic simulations to study the binding of CT to a membrane modelled as an asymmetrical GM1-DPPC bilayer. Simulation results indicate that the toxin binds to the membrane through only three of its five B subunits, in effect resulting in a tilted bound configuration. Additionally, the binding of the CT can increase the area per lipid of GM1 leaflet, which in turn can cause the membrane regions interacting with the bound subunits to experience significant bilayer thinning and lipid tail disorder across both the leaflets.

Dissipative Particle Dynamics Simulations for Phospholipid Membranes Based on a Four-To-One Coarse-Grained Mapping Scheme

In this article, a new set of parameters compatible with the dissipative particle dynamics (DPD) force field is developed for phospholipids. The coarse-grained (CG) models of these molecules are constructed by mapping four heavy atoms and their attached hydrogen atoms to one bead. The beads are divided into types distinguished by charge type, polarizability, and hydrogen-bonding capacity. First, we derive the relationship between the DPD repulsive force and Flory-Huggins χ-parameters based on this four-to-one CG mapping scheme. Then, we optimize the DPD force parameters for phospholipids. The feasibility of this model is demonstrated by simulating the structural and thermodynamic properties of lipid bilayer membranes, including the membrane thickness, the area per lipid, the lipid tail orientation, the bending rigidity, the rupture behavior, and the potential of mean force for lipid flip-flop.

Molecular Dynamics Simulation of Cholera Toxin A-1 Polypeptide

A molecular dynamics (MD) simulation study of the enzymatic portion of cholera toxin; cholera toxin A-1 polypeptide (CTA1) was performed at 283, 310 and 323 K. From total energy analysis it was observed that this toxin is stable thermodynamically and these outcomes were likewise confirmed by root mean square deviations (RMSD) investigations. The Cα root mean square fluctuation (RMSF) examinations revealed that there are a number of residues inside CTA1, which can be used as target for designing and synthesizing inhibitory drugs, in order to inactivate cholera toxin inside the human body. The fluctuations in the radius of gyration and hydrogen bonding in CTA1 proved that protein unfolding and refolding were normal routine phenomena in its structure at all temperatures. Solvent accessible surface area study identified the hydrophilic nature of the CTA1, and due to this property it can be a potential biological weapon. The structural identification (STRIDE) algorithm for proteins was successfully used to determine the partially disordered secondary structure of CTA1. On account of this partially disordered secondary structure, it can easily deceive the proteolytic enzymes of the endoplasmic reticulum of host cells.

Interleaflet Coupling, Pinning, and Leaflet Asymmetry-Major Players in Plasma Membrane Nanodomain Formation

The plasma membrane has a highly asymmetric distribution of lipids and contains dynamic nanodomains many of which are liquid entities surrounded by a second, slightly different, liquid environment. Contributing to the dynamics is a continuous repartitioning of components between the two types of liquids and transient links between lipids and proteins, both to extracellular matrix and cytoplasmic components, that temporarily pin membrane constituents. This make plasma membrane nanodomains exceptionally challenging to study and much of what is known about membrane domains has been deduced from studies on model membranes at equilibrium. However, living cells are by definition not at equilibrium and lipids are distributed asymmetrically with inositol phospholipids, phosphatidylethanolamines and phosphatidylserines confined mostly to the inner leaflet and glyco- and sphingolipids to the outer leaflet. Moreover, each phospholipid group encompasses a wealth of species with different acyl chain combinations whose lateral distribution is heterogeneous. It is becoming increasingly clear that asymmetry and pinning play important roles in plasma membrane nanodomain formation and coupling between the two lipid monolayers. How asymmetry, pinning, and interdigitation contribute to the plasma membrane organization is only beginning to be unraveled and here we discuss their roles and interdependence.