Mixing small proteins with lipids and cholesterol

Many ternary mixtures composed of saturated and unsaturated lipids with cholesterol (Chol) exhibit a region of coexistence between liquid-disordered (Ld) and liquid-ordered (Lo) domains, bearing some similarities to lipid rafts in biological membranes. However, biological rafts also contain many proteins that interact with the lipids and modify the distribution of lipids. Here, we extend a previously published lattice model of ternary DPPC/DOPC/Chol mixtures by introducing a small amount of small proteins (peptides). We use Monte Carlo simulations to explore the mixing phase behavior of the components as a function of the interaction parameter representing the affinity between the proteins and the saturated DPPC chains and for different mixture compositions. At moderate fractions of DPPC, the system is in a two-phase Ld + Lo coexistence, and the proteins exhibit a simple partition behavior between the phases that depends on the protein-lipid affinity parameter. At low DPPC compositions, the mixture is in Ld phase with local nanoscopic ordered domains. The addition of proteins with sufficiently strong attraction to the saturated lipids can induce the separation of a distinct Lo large domain with tightly packed gel-like clusters of proteins and saturated lipids. Consistent with the theory of phase transitions, we observe that the domain sizes grow when the mixture composition is in the vicinity of the critical point. Our simulations show that the addition of a small amount of proteins to such mixtures can cause their size to grow even further and lead to the formation of metastable dynamic Lo domains with sizes comparable to biological rafts.

Molecular dynamics simulations of lipid composition and its impact on structural and dynamic properties of skin membrane

The stratum corneum (SC) plays the most important role in the absorption of topical and transdermal drugs. In this study, we developed a multi-layered SC model using coarse-grained molecular dynamics (CGMD) simulations of ceramides, cholesterol, and fatty acids in equimolar proportions, starting from two different initial configurations. In the first approach, all ceramide molecules were initially in the hairpin conformation, and the membrane bilayers were pre-formed. In the second approach, ceramide molecules were introduced in either the hairpin or splayed conformation, with the lipid molecules randomly oriented at the start of the simulation. The aim was to evaluate the effects of lipid chain length on the structural and dynamic properties of SC. By incorporating ceramides and fatty acids of different chain lengths, we simulated the SC membrane in healthy and diseased states. We calculated key structural properties including the thickness, normalized lipid area, lipid tail order parameters, and spatial ordering of the lipids from each system. The results showed that systems with higher ordering and structural integrity contained an equimolar ratio of ceramides (chain length of 24 carbon atoms), fatty acids with chain lengths ≥ of 20 carbon atoms, and cholesterol. In these systems, strong apolar interactions between the ceramide and fatty acid long acyl chains restricted the mobility of the lipid molecules, thereby maintaining a compact lipid headgroup region and high order in the lipid tail region. The simulations also revealed distinct flip-flop mechanisms for cholesterol and fatty acid within the multi-layered membrane. Cholesterol is mostly diffused through the tail-tail interface region of the membrane and could flip-flop in the same bilayer. In contrast, fatty acids flip-flopped between adjacent leaflets of two bilayers in which the tails crossed the thinner headgroup region of the membrane. To conclude, our SC model provides mechanistic insights into lipid mobility and is flexible in its design and composition of different lipids, enabling studies of varying skin conditions.

The density-threshold affinity: Calculating lipid binding affinities from unbiased coarse-grained molecular dynamics simulations

Many membrane proteins are sensitive to their local lipid environment. As structural methods for membrane proteins have improved, there is growing evidence of direct, specific binding of lipids to protein surfaces. Unfortunately the workhorse of understanding protein-small molecule interactions, the binding affinity for a given site, is experimentally inaccessible for these systems. Coarse-grained molecular dynamics simulations can be used to bridge this gap, and are relatively straightforward to learn. Such simulations allow users to observe spontaneous binding of lipids to membrane proteins and quantify localized densities of individual lipids or lipid fragments. In this chapter we outline a protocol for extracting binding affinities from these localized distributions, known as the "density threshold affinity." The density threshold affinity uses an adaptive and flexible definition of site occupancy that alleviates the need to distinguish between "bound'' lipids and bulk lipids that are simply diffusing through the site. Furthermore, the method allows "bead-level" resolution that is suitable for the case where lipids share binding sites, and circumvents ambiguities about a relevant reference state. This approach provides a convenient and straightforward method for comparing affinities of a single lipid species for multiple sites, multiple lipids for a single site, and/or a single lipid species modeled using multiple forcefields.

Multimillion Atom Simulations of Di-8-ANEPPS Chromophores Embedded in a Model Plasma Membrane: Toward the Investigation of Realistic Dyed Cell Membranes

A multistep computational approach has been employed to study a multimillion all-atom dyed plasma membrane, with no less than 42 different lipid species spanning the major head groups and a variety of fatty acids, as well as cholesterol, with the objective of investigating its structure and dynamics, as well as its impact on the embedded di-8-ANEPPS dyes. The latter are commonly used as bioimaging probes and serve as local microscopes. So, they provide information on membrane morphology via their second harmonic nonlinear optical (NLO) responses, which have the advantage of being specific to interface regions and sensitive to the chromophore environment. In previous studies, this chromophore has only been studied in simpler membrane models, far from the complexity of real lipid bilayers, while, owing to the ever-increasing computational resources, multimillion lipid bilayers have been studied, giving access to the effects of its heterogeneity. First, using molecular dynamics (MD) simulations, it is found that the combination of lipids produces a more ordered and denser membrane compared to its homogeneous model counterparts, while the local environment of the embedded dyes becomes enriched in phosphatidylcholine. Subsequently, the second harmonic first hyperpolarizability of the probes was calculated at the TDDFT level on selected frames of MD, highlighting the influence of the lipid environment. Due to the complexity of the system, machine learning (ML) tools have been employed to establish relationships between the membrane structural parameters, the orientation of the probes, and their NLO responses. These ML approaches have revealed influential features, including the presence of diacylglycerol lipids close to the dye. On the whole, this work provides a first step toward understanding the cooperation, synergy, and interactions that occur in such complex guest-host environments, which have emerged as new targets for drug design and membrane lipid therapy.

Martini 3 Force Field Parameters for Protein Lipidation Post-Translational Modifications

Protein lipidations are vital co/post-translational modifications that tether lipid tails to specific protein amino acids, allowing them to anchor to biological membranes, switch their subcellular localization, and modulate association with other proteins. Such lipidations are thus crucial for multiple biological processes including signal transduction, protein trafficking, and membrane localization and are implicated in various diseases as well. Examples of lipid-anchored proteins include the Ras family of proteins that undergo farnesylation; actin and gelsolin that are myristoylated; phospholipase D that is palmitoylated; glycosylphosphatidylinositol-anchored proteins; and others. Here, we develop parameters for cysteine-targeting farnesylation, geranylgeranylation, and palmitoylation, as well as glycine-targeting myristoylation for the latest version of the Martini 3 coarse-grained force field. The parameters are developed using the CHARMM36m all-atom force field parameters as reference. The behavior of the coarse-grained models is consistent with that of the all-atom force field for all lipidations and reproduces key dynamical and structural features of lipid-anchored peptides, such as the solvent-accessible surface area, bilayer penetration depth, and representative conformations of the anchors. The parameters are also validated in simulations of the lipid-anchored peripheral membrane proteins Rheb and Arf1, after comparison with independent all-atom simulations. The parameters, along with mapping schemes for the popular martinize2 tool, are available for download at 10.5281/zenodo.7849262 and also as supporting information.

Understanding the free-energy landscape of phase separation in lipid bilayers using molecular dynamics

Liquid-liquid phase separation inside the cell often results in biological condensates that can critically affect cell homeostasis. Such phase separation events occur in multiple parts of cells, including the cell membranes, where the "lipid raft" hypothesis posits the formation of ordered domains floating in a sea of disordered lipids. The resulting lipid domains often have functional roles. However, the thermodynamics of lipid phase separation and their resulting mechanistic effects on cell function and dysfunction are poorly understood. Understanding such complex phenomena in cell membranes, with their diverse lipid compositions, is exceptionally difficult. For these reasons, simple model systems that can recapitulate similar behavior are widely used to study this phenomenon. Despite these simplifications, the timescale and length scales of domain formation pose a challenge for molecular dynamics (MD) simulations. Thus, most MD studies focus on spontaneous lipid phase separation-essentially measuring the sign (but not the amplitude) of the free-energy change upon separation-rather than directly interrogating the thermodynamics. Here, we propose a proof-of-concept pipeline that can directly measure this free energy by combining coarse-grained MD with enhanced sampling protocols using a novel collective variable. This approach will be a useful tool to help connect the thermodynamics of phase separation with the mechanistic insights already available from MD simulations.

Understanding the Free Energy Landscape of Phase Separation in Lipid Bilayers using Molecular Dynamics

Liquid-liquid phase separation (LLPS) inside the cell often results in biological condensates that can critically impact cell homeostasis. Such phase separation events occur in multiple parts of cells, including the cell membranes, where the so-called "lipid raft" hypothesis posits the formation of ordered domains floating in a sea of disordered lipids. The resulting lipid domains often have functional roles. However, the thermodynamics of lipid phase separation and their resulting mechanistic effects on cell function and dysfunction are poorly understood. Understanding such complex phenomena in cell membranes, with their diverse lipid compositions, is exceptionally difficult. For this reasons, simple model systems that can recapitulate similar behavior are widely used to study this phenomenon. Despite these simplifications, the timescale and and length scales of domain formation pose a challenge for molecular dynamics (MD) simulations. Thus, most MD studies focus on spontaneous lipid phase separation - essentially measuring the sign (but not the amplitude) of the free energy change upon separation - rather than directly interrogating the thermodynamics. Here, we propose a proof-of-concept pipeline that can directly measure this free energy by combining coarse-grained MD with enhanced sampling protocols using a novel collective variable. This approach will be a useful tool to help connect the thermodynamics of phase separation with the mechanistic insights already available from molecular dynamics simulations.

SIGNIFICANCE

Standard molecular dynamics simulations can determine the sign the free energy change upon phase separation, but not the amplitude. We present a new method to determine the phase separation free energy for lipid membranes, based on a enhanced sampling using the weighted ensemble method combined with a novel collective variable, validated using coarse-grained simulations applied to several simple systems. The new method will be valuable as a way to develop models that connect molecular-level structural features to the thermodynamics of phase separation.

Lipid packing in biological membranes governs protein localization and membrane permeability

Plasma membrane (PM) heterogeneity has long been implicated in various cellular functions. However, mechanistic principles governing functional regulations of lipid environment are not well understood due to the inherent complexities associated with the relevant length and timescales that limit both direct experimental measurements and their interpretation. In this context, computer simulations hold immense potential to investigate molecular-level interactions and mechanisms that lead to PM heterogeneity and its functions. Herein, we investigate spatial and dynamic heterogeneity in model membranes with coexisting liquid ordered and liquid disordered phases and characterize the membrane order in terms of the local topological changes in lipid environment using the nonaffine deformation framework. Furthermore, we probe the packing defects in these membranes, which can be considered as the conjugate of membrane order assessed in terms of the nonaffine parameter. In doing so, we formalize the connection between membrane packing and local membrane order and use that to explore the mechanistic principles behind their functions. Our observations suggest that heterogeneity in mixed phase membranes is a consequence of local lipid topology and its temporal evolution, which give rise to disparate lipid packing in ordered and disordered domains. This in turn governs the distinct nature of packing defects in these domains that can play a crucial role in preferential localization of proteins in mixed phase membranes. Furthermore, we observe that lipid packing also leads to contrasting distribution of free volume in the membrane core region in ordered and disordered membranes, which can lead to distinctive membrane permeability of small molecules. Our results, thus, indicate that heterogeneity in mixed phase membranes closely governs the membrane functions that may emerge from packing-related basic design principles.

Identification of possible binding modes of SARS-CoV-2 spike N-terminal domain for ganglioside GM1

Coarse-grained molecular dynamics simulations of the lipid bilayer mixture of POPC and cholesterol were carried out in the presence and absence of ganglioside monosialo 1 (GM1) with N - terminal domain (NTD) of SARS-CoV-2 spike glycoprotein. The interactions of GM1 with two different NTD orientations were compared. NTD orientation I compactly bind GM1 predominantly through the sialic acid and the external galactose moieties providing more restriction to GM1 mobility whereas orientation II is more distributed on the lipid surface and due to the relaxed mobility of GM1 there, presumably, the NTD receptor penetrates more through the membrane.

Cofilin-Membrane Interactions: Electrostatic Effects in Phosphoinositide Lipid Binding

The actin cytoskeleton interacts with the cell membrane primarily through the indirect interactions of actin-binding proteins such as cofilin-1. The molecular mechanisms underlying the specific interactions of cofilin-1 with membrane lipids are still unclear. Here, we performed coarse-grain molecular dynamics simulations of cofilin-1 with complex lipid bilayers to analyze the specificity of protein-lipid interactions. We observed the maximal interactions with phosphoinositide (PIP) lipids, especially PIP₂ and PIP₃ lipids. A good match was observed between the residues predicted to interact and previous experimental studies. The clustering of PIP lipids around the membrane bound protein leads to an overall lipid demixing and gives rise to persistent membrane curvature. Further, through a series of control simulations, we observe that both electrostatics and geometry are critical for specificity of lipid binding. Our current study is a step towards understanding the physico-chemical basis of cofilin-PIP lipid interactions.

Martini 3 Coarse-Grained Force Field for Carbohydrates

The Martini 3 force field is a full reparametrization of the Martini coarse-grained model for biomolecular simulations. Due to the improved interaction balance, it allows for a more accurate description of condensed phase systems. In the present work, we develop a consistent strategy to parametrize carbohydrate molecules accurately within the framework of Martini 3. In particular, we develop a canonical mapping scheme which decomposes arbitrarily large carbohydrates into a limited number of fragments. Bead types for these fragments have been assigned by matching physicochemical properties of mono- and disaccharides. In addition, guidelines for assigning bonds, angles, and dihedrals were developed. These guidelines enable a more accurate description of carbohydrate conformations than in the Martini 2 force field. We show that models obtained with this approach are able to accurately reproduce osmotic pressures of carbohydrate water solutions. Furthermore, we provide evidence that the model differentiates correctly the solubility of the polyglucoses dextran (water-soluble) and cellulose (water insoluble but soluble in ionic liquids). Finally, we demonstrate that the new building blocks can be applied to glycolipids. We show they are able to reproduce membrane properties and induce binding of peripheral membrane proteins. These test cases demonstrate the validity and transferability of our approach.

Molecular Mechanism of CD44 Homodimerization Modulated by Palmitoylation and Membrane Environments

The homodimerization of CD44 plays a key role in an intercellular-to-extracellular signal transduction and tumor progression. Acylated modification and specific membrane environments have been reported to mediate translocation and oligomerization of CD44, however, the underlying molecular mechanism remains elusive. In this study, extensive molecular dynamics simulations are performed to characterize the dimerization of palmitoylated CD44 variants in different bilayer environments. CD44 forms homodimer depending on the cysteines on the juxta-membrane domains, and the dimerization efficiency and packing configurations are defected by their palmitoylated modifications. In the phase-segregated (raft included) membrane, homodimerization of the palmitoylated CD44 is hardly observed, whereas PIP2 addition compensates to realize dimerization. However, the structure of CD44 homodimer formed in the phase-segregated bilayer turns susceptive and PIP2 addition allows for an extensive conformation of the cytoplasmic domain, a proposal prerequisite to access the cytoskeleton linker proteins. The results unravel a delicate competitive relationship between PIP2, lipid raft and palmitoylation in mediating protein homodimerization, which helps to clarify the dynamic dimer conformations and involved cellular signaling of the CD44 likewise proteins.

Molecular Mechanisms Underlying Caveolin-1 Mediated Membrane Curvature

Caveolin-1 is one of the main protein components of caveolae that acts as a mechanosensor at the cell membrane. The interactions of caveolin-1 with membranes have been shown to lead to complex effects such as curvature and the clustering of specific lipids. Here, we review the emerging concepts on the molecular interactions of caveolin-1, with a focus on insights from coarse-grain molecular dynamics simulations. Consensus structural models of caveolin-1 report a helix-turn-helix core motif with flanking domains of higher disorder that could be membrane composition dependent. Caveolin-1 appears to be mainly surface-bound and does not embed very deep in the membrane to which it is bound. The most interesting aspect of caveolin-1 membrane binding is the interplay of cholesterol clustering and membrane curvature. Although cholesterol has been reported to cluster in the vicinity of caveolin-1 by several approaches, simulations show that the clustering is maximal in membrane leaflet opposing the surface-bound caveolin-1. The intrinsic negative curvature of cholesterol appears to stabilize the negative curvature in the opposing leaflet. In fact, the simulations show that blocking cholesterol clustering (through artificial position restraints) blocks membrane curvature, and vice versa. Concomitant with cholesterol clustering is sphingomyelin clustering, again in the opposing leaflet, but in a concentration-dependent manner. The differential stress due to caveolin-1 binding and the inherent asymmetry of the membrane leaflets could be the determinant for membrane curvature and needs to be further probed. The review is an important step to reconcile the molecular level details emerging from simulations with the mesoscopic details provided by state of the art experimental approaches.

In Silico-Based Experiments on Mechanistic Interactions between Several Intestinal Permeation Enhancers with a Lipid Bilayer Model

Oral administration of drugs is generally considered convenient and patient-friendly. However, oral administration of biological drugs exhibits low oral bioavailability (BA) due to enzymatic degradation and low intestinal absorption. A possible approach to circumvent the low BA of oral peptide drugs is to coformulate the drugs with permeation enhancers (PEs). PEs have been studied since the 1960s and are molecules that enhance the absorption of hydrophilic molecules with low permeability over the gastrointestinal epithelium. In this study, we investigated the impact of six PEs on the structural properties of a model membrane using molecular dynamics (MD) simulations. The PEs included were the sodium salts of the medium chain fatty acids laurate, caprate, and caprylate and the caprylate derivative SNAC─all with a negative charge─and neutral caprate and neutral sucrose monolaurate. Our results indicated that the PEs, once incorporated into the membrane, could induce membrane leakiness in a concentration-dependent manner. Our simulations suggest that a PE concentration of at least 70-100 mM is needed to strongly affect transcellular permeability. The increased aggregation propensity seen for neutral PEs might provide a molecular-level mechanism for the membrane disruptions seen at higher concentrations in vivo. The ability for neutral PEs to flip-flop across the lipid bilayer is also suggestive of possible intracellular modes of action other than increasing membrane fluidity. Taken together, our results indicate that MD simulations are useful for gaining insights relevant to the design of oral dosage forms based around permeability enhancer molecules.

Characteristic lengths of transmembrane peptides controlling their tilt and lateral distribution between membrane domains

Lipids and proteins of plasma membranes of eukaryotic cells are supposed to form protein-lipid domains, characterized by a different molecular order, bilayer thickness, and elastic parameters. Several mechanisms of preferable distribution of transmembrane proteins to the ordered or disordered membrane domains have been revealed. The mismatch between the length of the protein transmembrane domain and hydrophobic thickness of the lipid bilayer is considered to be an important driving force of protein lateral sorting. Utilizing the continuum theory of elasticity, we analyzed optimal configurations and preferable membrane domains for single-pass transmembrane peptides of various hydrophobic lengths and effective molecular shapes. We obtained that short transmembrane peptides stand perpendicularly to the membrane plane. The exceedance of a certain characteristic length leads to the tilt of the peptide. This length depends on the bilayer thickness. Thus, in the membrane with coexisting ordered (thicker) and disordered (thinner) phases tilting of the peptide in each phase is governed by its individual characteristic length. The lateral distribution of the peptides between ordered and disordered membrane domains is shown to be described by two additional characteristic lengths. The exceedance of the smaller one drives the peptide towards a more ordered and thicker membrane, while the exceedance of the larger characteristic length switches the preferable membrane domain from ordered and thicker to disordered and thinner. Thus, membrane proteins with long enough transmembrane domains are predicted to accumulate in the thinner disordered membrane as compared to the thicker ordered bilayer. For hourglass-like and barrel-like shaped transmembrane peptides the specific regime of sorting was obtained: the peptides distributed almost equally between the phases in a wide range of peptide lengths. This finding allowed explaining the experimental data on lateral distribution of transmembrane peptide tLAT.

Unveiling the Dynamics of KRAS4b on Lipid Model Membranes

Small GTPase proteins are ubiquitous and responsible for regulating several processes related to cell growth and differentiation. Mutations that stabilize their active state can lead to uncontrolled cell proliferation and cancer. Although these proteins are well characterized at the cellular scale, the molecular mechanisms governing their functions are still poorly understood. In addition, there is limited information about the regulatory function of the cell membrane which supports their activity. Thus, we have studied the dynamics and conformations of the farnesylated KRAS4b in various membrane model systems, ranging from binary fluid mixtures to heterogeneous raft mimics. Our approach combines long time-scale coarse-grained (CG) simulations and Markov state models to dissect the membrane-supported dynamics of KRAS4b. Our simulations reveal that protein dynamics is mainly modulated by the presence of anionic lipids and to some extent by the nucleotide state (activation) of the protein. In addition, our results suggest that both the farnesyl and the polybasic hypervariable region (HVR) are responsible for its preferential partitioning within the liquid-disordered (Ld) domains in membranes, potentially enhancing the formation of membrane-driven signaling platforms.

Advances in understanding and in multi-disciplinary methodology used to assess lipid regulation of signalling cascades from the cancer cell plasma membrane

The lipid bilayer is a functional component of cells, forming a stable platform for the initiation of key biological processes, including cell signalling. There are distinct changes in the lipid composition of cell membranes during oncogenic transformation resulting in aberrant activation and inactivation of signalling transduction pathways. Studying the role of the cell membrane in cell signalling is challenging, since techniques are often limited to by timescale, resolution, sensitivity, and averaging. To overcome these limitations, combining 'computational', 'wet-lab' and 'semi-dry' approaches offers the best opportunity to resolving complex biological processes involved in membrane organisation. In this review, we highlight analytical tools that have been applied for the study of cell signalling initiation from the cancer cell membranes through computational microscopy, biological assays, and membrane biophysics. The cancer therapeutic potential of extracellular membrane-modulating agents, such as cholesterol-reducing agents is also discussed, as is the need for future collaborative inter-disciplinary research for studying the role of the cell membrane and its components in cancer therapy.

A Trifunctional Linker for Palmitoylation and Peptide and Protein Localization in Biological Membranes

Attachment of lipophilic groups is an important post-translational modification of proteins, which involves the coupling of one or more anchors such as fatty acids, isoprenoids, phospholipids, or glycosylphosphatidyl inositols. To study its impact on the membrane partitioning of hydrophobic peptides or proteins, we designed a tyrosine-based trifunctional linker. The linker allows the facile incorporation of two different functionalities at a cysteine residue in a single step. We determined the effect of the lipid modification on the membrane partitioning of the synthetic α-helical model peptide WALP with or without here and in all cases below; palmitoyl groups in giant unilamellar vesicles that contain a liquid-ordered (Lo ) and liquid-disordered (Ld ) phase. Introduction of two palmitoyl groups did not alter the localization of the membrane peptides, nor did the membrane thickness or lipid composition. In all cases, the peptide was retained in the Ld phase. These data demonstrate that the Lo domain in model membranes is highly unfavorable for a single membrane-spanning peptide.

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.

Elastic deformations mediate interaction of the raft boundary with membrane inclusions leading to their effective lateral sorting

Liquid-ordered lipid domains represent a lateral inhomogeneity in cellular membranes. These domains have elastic and physicochemical properties different from those of the surrounding membrane. In particular, their thickness exceeds that of the disordered membrane. Thus, elastic deformations arise at the domain boundary in order to compensate for the thickness mismatch. In equilibrium, the deformations lead to an incomplete register of monolayer ordered domains: the elastic energy is minimal if domains in opposing monolayers lie on the top of each other, and their boundaries are laterally shifted by about 3 nm. This configuration introduces a region, composed of one ordered and one disordered monolayers, with an intermediate bilayer thickness. Besides, a jump in a local monolayer curvature takes place in this intermediate region, concentrating here most of the elastic stress. This region can participate in a lateral sorting of membrane inclusions by offering them an optimal bilayer thickness and local curvature conditions. In the present study, we consider the sorting of deformable lipid inclusions, undeformable peripheral and deeply incorporated peptide inclusions, and undeformable transmembrane inclusions of different molecular geometry. With rare exceptions, all types of inclusions have an affinity to the ordered domain boundary as compared to the bulk phases. The optimal lateral distribution of inclusions allows relaxing the elastic stress at the boundary of domains.

Post-Translational Modifications at the Coarse-Grained Level with the SIRAH Force Field

Post-translational modifications (PTMs) on proteins significantly enlarge the physicochemical diversity present in biological macromolecules, altering function, localization, and interactions. Despite their critical role in regulating cellular processes, theoretical methods are not yet fully capable of coping with this diversity. These limitations are particularly more marked for coarse-grained (CG) models, in which comprehensive and self-consistent parametrizations are less frequent. Here we present a set of topologies and interaction parameters for the most common PTMs, fully compatible with the SIRAH force field. The PTMs introduced here reach the same level of structural description of the existing SIRAH force field, expanding the chemical spectrum with promising applications in dynamical protein-protein interactions in large and complex cellular environments.

Interaction of Antimicrobial Lipopeptides with Bacterial Lipid Bilayers

The resistance of pathogens to traditional antibiotics is currently a global issue of enormous concern. As the discovery and development of new antibiotics become increasingly challenging, synthetic antimicrobial lipopeptides (AMLPs) are now receiving renewed attention as a new class of antimicrobial agents. In contrast to traditional antibiotics, AMLPs act by physically disrupting the cell membrane (rather than targeting specific proteins), thus reducing the risk of inducing bacterial resistance. In this study, we use microsecond-timescale atomistic molecular dynamics simulations to quantify the interaction of a short AMLP (C16-KKK) with model bacterial lipid bilayers. In particular, we investigate how fundamental transmembrane properties change in relation to a range of lipopeptide concentrations. A number of structural, mechanical, and dynamical features are found to be significantly altered in a non-linear fashion. At 10 mol% concentration, lipopeptides have a condensing effect on bacterial bilayers, characterized by a decrease in the area per lipid and an increase in the bilayer order. Higher AMLP concentrations of 25 and 40 mol% destabilize the membrane by disrupting the bilayer core structure, inducing membrane thinning and water leakage. Important transmembrane properties such as the lateral pressure and dipole potential profiles are also affected. Potential implications on membrane function and associated proteins are discussed.

Gangliosides Destabilize Lipid Phase Separation in Multicomponent Membranes

Gangliosides (GMs) form an important class of lipids found in the outer leaflet of the plasma membrane. Typically, they colocalize with cholesterol and sphingomyelin in ordered membrane domains. However, detailed understanding of the lateral organization of GM-rich membranes is still lacking. To gain molecular insight, we performed molecular dynamics simulations of GMs in model membranes composed of coexisting liquid-ordered and liquid-disordered domains. We found that GMs indeed have a preference to partition into the ordered domains. At higher concentrations (>10 mol %), we observed a destabilizing effect of GMs on the phase coexistence. Further simulations with modified GMs show that the structure of the GM headgroup affects the phase separation, whereas the nature of the tail determines the preferential location. Together, our findings provide a molecular basis to understand the lateral organization of GM-rich membranes.

Antimicrobial peptide ROAD-1 triggers phase change in local membrane environment to execute its activity

Emergence of antibiotic-resistant pathogens has paved way for development of newer class of drugs that would not be susceptible to resistance. Antimicrobial peptides such as defensins that target the microbial membrane are promising candidates. ROAD-1 is an alpha-defensin present in the oral cavity of rhesus macaque and shares very high sequence similarity to human enteric defensin 5. In this study we have performed microsecond long all atom molecular dynamic simulations to understand the mechanism of action of ROAD-1. We find that ROAD-1 is able to adopt an energetically stable conformation predominantly stabilized by electrostatic interactions only in presence of bacterial membranes. In mammalian membrane even though it gets absorbed onto the bilayer, it is unable to adopt an equilibrium conformation. Binding of ROAD-1 to bilayer induces clustering of POPG molecules up to 15 Å around the peptide. POPG molecules show higher order parameters than the neighboring POPE implying coexistence of different phases. Analysis of binding free energy of ROAD-1-membrane complex indicates Arg1, Arg2, Arg7, and Arg25 to play key role in its antimicrobial activity. Unlike its homolog HD5, ROAD-1 is not observed to form a dimer. Our study gives insight into the membrane-bound conformation of ROAD-1 and its mechanism of action that can aid in designing defensin-based therapeutics. Graphical abstract Antimicrobial peptide ROAD-1 adopts a different membrane-bound conformation as compared with HD5 even though they belong to the same family implying a different mechanism of action.

Protein Partitioning into Ordered Membrane Domains: Insights from Simulations

Cellular membranes are laterally organized into domains of distinct structures and compositions by the differential interaction affinities between various membrane lipids and proteins. A prominent example of such structures are lipid rafts, which are ordered, tightly packed domains that have been widely implicated in cellular processes. The functionality of raft domains is driven by their selective recruitment of specific membrane proteins to regulate their interactions and functions; however, there have been few general insights into the factors that determine the partitioning of membrane proteins between coexisting liquid domains. In this work, we used extensive coarse-grained and atomistic molecular dynamics simulations, potential of mean force calculations, and conceptual models to describe the partitioning dynamics and energetics of a model transmembrane domain from the linker of activation of T cells. We find that partitioning between domains is determined by an interplay between protein-lipid interactions and differential lipid packing between raft and nonraft domains. Specifically, we show that partitioning into ordered domains is promoted by preferential interactions between peptides and ordered lipids, mediated in large part by modification of the peptides by saturated fatty acids (i.e., palmitoylation). Ordered phase affinity is also promoted by elastic effects, specifically hydrophobic matching between the membrane and the peptide. Conversely, ordered domain partitioning is disfavored by the tight molecular packing of the lipids therein. The balance of these dominant drivers determines partitioning. In the case of the wild-type linker of activation of T cells transmembrane domain, these factors combine to yield enrichment of the peptide at L_o/L_d interfaces. These results define some of the general principles governing protein partitioning between coexisting membrane domains and potentially explain previous disparities among experiments and simulations across model systems.

Exploring the impact of proteins on the line tension of a phase-separating ternary lipid mixture

The separation of lipid mixtures into thermodynamically stable phase-separated domains is dependent on lipid composition, temperature, and system size. Using molecular dynamics simulations, the line tension between thermodynamically stable lipid domains formed from ternary mixtures of di-C16:0 PC:di-C18:2 PC:cholesterol at 40:40:20 mol. % ratio was investigated via two theoretical approaches. The line tension was found to be 3.1 ± 0.2 pN by capillary wave theory and 4.7 ± 3.7 pN by pressure tensor anisotropy approaches for coarse-grained models based on the Martini force field. Using an all-atom model of the lipid membrane based on the CHARMM36 force field, the line tension was found to be 3.6 ± 0.9 pN using capillary wave theory and 1.8 ± 2.2 pN using pressure anisotropy approaches. The discrepancy between estimates of the line tension based on capillary wave theory and pressure tensor anisotropy methods is discussed. Inclusion of protein in Martini membrane lipid mixtures was found to reduce the line tension by 25%-35% as calculated by the capillary wave theory approach. To further understand and predict the behavior of proteins in phase-separated membranes, we have formulated an analytical Flory-Huggins model and parameterized it against the simulation results. Taken together these results suggest a general role for proteins in reducing the thermodynamic cost associated with domain formation in lipid mixtures and quantifies the thermodynamic driving force promoting the association of proteins to domain interfaces.

Computational Modeling of Realistic Cell Membranes

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

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.

Understanding Membrane Domain-Partitioning Thermodynamics of Transmembrane Domains with Potential of Mean Force Calculations

The transmembrane domain (TMD) of membrane proteins plays an essential role in their dynamics and functions. Certain properties of TMDs, such as raft affinity and orientation, have been studied extensively both experimentally and computationally. However, the extent to which specific physicochemical properties of TMDs determine their membrane domain-partitioning thermodynamics is still far from clear. In this work, we propose an approach based on umbrella sampling molecular dynamics simulations of model membranes and idealized TMDs to quantify the effect of TMD physicochemical properties, namely, length, degree of hydrophobicity, and size of TMDs, on their membrane domain-partitioning thermodynamics. The results, which are fully consistent with previous experimental and simulation data, indicate that the concept of "hydrophobic mismatch" should go beyond differences in hydrophobic thickness to include mismatch in the degree of hydrophobicity between the TMD and the surrounding hydrocarbon lipid chains. Our method provides quantitative insights into the role of specific physicochemical features of TMDs in membrane localization and orientation, which will be broadly useful for predicting the raft affinity and membrane partitioning of any transmembrane protein.

Interplay between Membrane Curvature and Cholesterol: Role of Palmitoylated Caveolin-1

Caveolin-1 (cav-1) is an important player in cell signaling and endocytosis that has been shown to colocalize with cholesterol-rich membrane domains. Experimental studies with varying cav-1 constructs have suggested that it can induce both cholesterol clustering and membrane curvature. Here, we probe the molecular origin of membrane curvature and cholesterol clustering by cav-1 by using coarse-grain molecular dynamics simulations. We have performed a series of simulations of a functionally important cav-1 construct, comprising the membrane-interacting domains and a C-terminal palmitoyl tail. Our results suggest that cav-1 is able to induce cholesterol clustering in the membrane leaflet to which it is bound as well as the opposing leaflet. A positive membrane curvature is observed upon cav-1 binding in cholesterol-containing bilayers. Interestingly, we observe an interplay between cholesterol clustering and membrane curvature such that cav-1 is able to induce higher membrane curvature in cholesterol-rich membranes. The role of the cav-1 palmitoyl tail is less clear and appears to increase the membrane contacts. Further, we address the importance of the secondary structure of cav-1 domains and show that it could play an important role in membrane curvature and cholesterol clustering. Our work is an important step toward a molecular picture of caveolae and vesicular endocytosis.

Theoretical Insights into the Interactions between Star-Shaped Antimicrobial Polypeptides and Bacterial Membranes

A structurally nanoengineered antimicrobial polypeptide consisting of lysine and valine residues is a new class of antimicrobial agent with superior antibacterial activity against multidrug-resistant bacteria and low toxicity toward mammalian cells. Utilizing coarse-grained models, we studied the interactions of microbial cytoplasmic membranes with polypeptides of either (K₂V₁)₅ (star-KV) or CM15 (star-CM15). Our computational results verify the low toxicity of polypeptides of (K₂V₁)₅ toward the dipalmitoyl phosphatidylcholine bilayer. This low toxicity is demonstrated to originate from weakened hydrophobicity combined with its random coil conformation for (K₂V₁)₅ because of the highly abundant valine residues, compared with the typical antimicrobial peptides, such as CM15. In the interactions with a palmitoyl-oleoyl-phosphatidylethanolamine/palmitoyl-oleoyl-phosphatidylglycerol bilayer, star-KV has greater ability in phase separation and generation of phase boundary defects not only in lipid redistribution but also in lateral dynamic movements, although both star-KV and star-CM15 can extract the phosphatidylglycerol lipids and purify the phosphatidylethanolamine lipids into continuum domains. We suggest that the polypeptide of (K₂V₁)₅ can nondisruptively kill bacteria by hampering bacterial metabolism through reorganizing lipid domain distribution and simultaneously "freezing" lipid movement.

Phase Behavior of GM1-Containing DMPC-Cholesterol Monolayer: Experimental and Theoretical Study

Organization and distribution of lipids in cellular membranes play an important role in a diverse range of biological processes, such as membrane trafficking and signaling. Here, we present the combined experimental and simulated results to elucidate the phase behavioral features of ganglioside monosialo 1 (GM1)-containing mixed monolayer of the lipids 1,2-dimyristoyl- sn-glycero-3-phosphocholine (DMPC) and cholesterol (CHOL). Two monolayers having compositions DMPC-CHOL and GM1-DMPC-CHOL are investigated at air-water and air-solid interfaces using Langmuir-Blodgett experiments and scanning electron microscopy (SEM), respectively, to ascertain the phase behavior change of the monolayers. Surface pressure isotherms and SEM imaging of domain formation indicate that addition of GM1 to the monolayer at low surface pressure causes a fluidization of the system but once the system attains the surface pressure corresponding to its liquid-condensed phase, the monolayer becomes more ordered than the system devoid of GM1 and interacts among each other more cooperatively. Besides, the condensing effect of cholesterol on the DMPC monolayer was also verified by our experiments. Apart from these, the effects induced by GM1 on the phase behavior of the binary mixture of DMPC-CHOL were studied with and without applying liquid-expanded (L_E)-liquid-condensed (L_C) equilibrium surface pressure using molecular dynamics (MD) simulation. Our molecular dynamics (MD) simulation results give an atomistic-level explanation of our experimental findings and furnish a similar conclusion.

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.

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.

Parameterization of Palmitoylated Cysteine, Farnesylated Cysteine, Geranylgeranylated Cysteine, and Myristoylated Glycine for the Martini Force Field

Peripheral membrane proteins go through various post-translational modifications that covalently bind fatty acid tails to specific amino acids. These post-translational modifications significantly alter the lipophilicity of the modified proteins and allow them to anchor to biological membranes. Over 1000 different proteins have been identified to date that require such membrane-protein interactions to carry out their biological functions, including members of the Src and Ras superfamilies that play key roles in cell signaling and carcinogenesis. We have used all-atom simulations with the CHARMM36 force field to parameterize four of the most common post-translational modifications for the Martini 2.2 force field: palmitoylated cysteine, farnesylated cysteine, geranylgeranylated cysteine, and myristoylated glycine. The parameters reproduce the key features of clusters of configurations of the different anchors in lipid membranes as well as the water-octanol partitioning free energies of the anchors, which are crucial for the correct reproduction of the expected biophysical behavior of peripheral membrane proteins at the membrane-water interface. Implementation in existing Martini setup tools facilitates the use of the new parameters.

Characterizing the interactions of two lipid modifications with lipid rafts: farnesyl anchors vs. palmitoyl anchors

Farnesyl (Far) and palmitoyl (Pal) anchors play important roles in the traffic of many lipidated proteins. Herein, we show the distinctive interactions and influences of the two lipid modifications on lipid rafts (LRs) and non-raft-like membranes using molecular dynamics simulations. Palmitoyl anchors behave in a more ordered fashion, pack tighter with the lipids of LRs and diffuse at a slower rate than farnesyl anchors in LRs. When interacting with non-raft-like membranes these two types of anchors become less ordered, pack more loosely with lipids, and diffuse at a higher rate. By calculating both the number of contacts per chain and the number of contact atoms per carbon of the two anchors with the lipid components, we found that the palmitoyl chains preferred to associate with the saturated chains of lipids and cholesterol molecules in LRs, while farnesyl chains favored association with saturated chains and unsaturated chains. For non-raft-like membranes, these two lipid anchors had roughly the same preference for the three types of contact lipid chains. Additionally, palmitoyl anchors caused cholesterol to orient more perpendicular to the membrane surface, surrounding lipids to become more ordered, and lipid lateral fluidity to reduce significantly, compared to farnesyl anchors in LRs. By contrast, the POPE and DSPC became much less ordered, cholesterol became more tilted, and lipids became more fluid, when the two types lipid anchors were inserted in non-raft-like membranes. These findings are useful for understanding the traffic mechanisms of lipidated proteins with farnesyl and palmitoyl modifications in cell membranes.

Protein post-translational modifications: In silico prediction tools and molecular modeling

Post-translational modifications (PTMs) occur in almost all proteins and play an important role in numerous biological processes by significantly affecting proteins' structure and dynamics. Several computational approaches have been developed to study PTMs (e.g., phosphorylation, sumoylation or palmitoylation) showing the importance of these techniques in predicting modified sites that can be further investigated with experimental approaches. In this review, we summarize some of the available online platforms and their contribution in the study of PTMs. Moreover, we discuss the emerging capabilities of molecular modeling and simulation that are able to complement these bioinformatics methods, providing deeper molecular insights into the biological function of post-translational modified proteins.

Effect of the Southeast Asian Ovalocytosis Deletion on the Conformational Dynamics of Signal-Anchor Transmembrane Segment 1 of Red Cell Anion Exchanger 1 (AE1, Band 3, or SLC4A1)

The first transmembrane (TM1) helix in the red cell anion exchanger (AE1, Band 3, or SLC4A1) acts as an internal signal anchor that binds the signal recognition particle and directs the nascent polypeptide chain to the endoplasmic reticulum (ER) membrane where it moves from the translocon laterally into the lipid bilayer. The sequence N-terminal to TM1 forms an amphipathic helix that lies at the membrane interface and is connected to TM1 by a bend at Pro403. Southeast Asian ovalocytosis (SAO) is a red cell abnormality caused by a nine-amino acid deletion (Ala400-Ala408) at the N-terminus of TM1. Here we demonstrate, by extensive (∼4.5 μs) molecular dynamics simulations of TM1 in a model 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine membrane, that the isolated TM1 peptide is highly dynamic and samples the structure of TM1 seen in the crystal structure of the membrane domain of AE1. The SAO deletion not only removes the proline-induced bend but also causes a "pulling in" of the part of the amphipathic helix into the hydrophobic phase of the bilayer, as well as the C-terminal of the peptide. The dynamics of the SAO peptide very infrequently resembles the structure of TM1 in AE1, demonstrating the disruptive effect the SAO deletion has on AE1 folding. These results provide a precise molecular view of the disposition and dynamics of wild-type and SAO TM1 in a lipid bilayer, an important early biosynthetic intermediate in the insertion of AE1 into the ER membrane, and extend earlier results of cell-free translation experiments.

The Ebola virus protein VP40 hexamer enhances the clustering of PI(4,5)P2 lipids in the plasma membrane

The Ebola virus is a lipid-enveloped virus that obtains its lipid coat from the plasma membrane of the host cell it infects during the budding process. The Ebola virus protein VP40 localizes to the inner leaflet of the plasma membrane and forms the viral matrix, which provides the major structure for the Ebola virus particles. VP40 is initially a dimer that rearranges to a hexameric structure that mediates budding. VP40 hexamers and larger filaments have been shown to be stabilized by PI(4,5)P2 in the plasma membrane inner leaflet. Reduction in the plasma membrane levels of PI(4,5)P2 significantly reduce formation of VP40 oligomers and virus-like particles. We investigated the lipid-protein interactions in VP40 hexamers at the plasma membrane. We quantified lipid-lipid self-clustering by calculating the fractional interaction matrix and found that the VP40 hexamer significantly enhances the PI(4,5)P2 clustering. The radial pair distribution functions suggest a strong interaction between PI(4,5)P2 and the VP40 hexamer. The cationic Lys side chains are found to mediate the PIP2 clustering around the protein, with cholesterol filling the space between the interacting PIP2 molecules. These computational studies support recent experimental data and provide new insights into the mechanisms by which VP40 assembles at the plasma membrane inner leaflet, alters membrane curvature, and forms new virus-like particles.

Ganglioside-Lipid and Ganglioside-Protein Interactions Revealed by Coarse-Grained and Atomistic Molecular Dynamics Simulations

Gangliosides are glycolipids in which an oligosaccharide headgroup containing one or more sialic acids is connected to a ceramide. Gangliosides reside in the outer leaflet of the plasma membrane and play a crucial role in various physiological processes such as cell signal transduction and neuronal differentiation by modulating structures and functions of membrane proteins. Because the detailed behavior of gangliosides and protein-ganglioside interactions are poorly known, we investigated the interactions between the gangliosides GM1 and GM3 and the proteins aquaporin (AQP1) and WALP23 using equilibrium molecular dynamics simulations and potential of mean force calculations at both coarse-grained (CG) and atomistic levels. In atomistic simulations, on the basis of the GROMOS force field, ganglioside aggregation appears to be a result of the balance between hydrogen bond interactions and steric hindrance of the headgroups. GM3 clusters are slightly larger and more ordered than GM1 clusters due to the smaller headgroup of GM3. The different structures of GM1 and GM3 clusters from atomistic simulations are not observed at the CG level based on the Martini model, implying a difference in driving forces for ganglioside interactions in atomistic and CG simulations. For protein-ganglioside interactions, in the atomistic simulations, GM1 lipids bind to specific sites on the AQP1 surface, whereas they are depleted from WALP23. In the CG simulations, the ganglioside binding sites on the AQP1 surface are similar, but ganglioside aggregation and protein-ganglioside interactions are more prevalent than in the atomistic simulations. Using the polarizable Martini water model, results were closer to the atomistic simulations. Although experimental data for validation is lacking, we proposed modified Martini parameters for gangliosides to more closely mimic the sizes and structures of ganglioside clusters observed at the atomistic level.

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.

Effects of Transmembrane α-Helix Length and Concentration on Phase Behavior in Four-Component Lipid Mixtures: A Molecular Dynamics Study

We used coarse-grained molecular dynamics simulations to examine the effects of transmembrane α-helical WALP peptides on the behavior of four-component lipid mixtures. These mixtures contain a high-melting temperature (high-Tm) lipid, a nanodomain-inducing low-Tm lipid, a macrodomain-inducing low-Tm lipid and cholesterol to model the outer leaflet of cell plasma membranes. In a series of simulations, we incrementally replace the nanodomain-inducing low-Tm lipid by the macrodomain-inducing low-Tm lipid and measure how lipid and phase properties are altered by the addition of WALPs of different length. Regardless of the ratio of the two low-Tm lipids, shorter WALPs increase domain size and all WALPs increase domain alignment between the two leaflets. These effects are smallest for the longest WALP tested, and increase with increasing WALP concentration. Thus, our simulations explain the experimental observation that WALPs induce macroscopic domains in otherwise nanodomain-forming lipid-only mixtures (unpublished). Since the cell plasma membrane contains a large fraction of transmembrane proteins, these findings link the behavior of lipid-only model membranes in vitro to phase behavior in vivo.

Lateral pressure-mediated protein partitioning into liquid-ordered/liquid-disordered domains

We have studied the contributions of stored elastic energies in liquid-ordered (Lo) and liquid-disordered (Ld) domains to transmembrane proteins using the lateral pressure concept. In particular we applied previously reported experimental data for the membrane thickness, intrinsic curvature and bending elasticities of coexisting Lo/Ld domains to calculate whether proteins of simple geometric shapes would preferentially diffuse into Lo or Ld domains and form oligomers of a certain size. For the studied lipid mixture we generally found that proteins with convex shapes prefer sorting to Ld phases and the formation of large clusters. Lo domains in turn would be enriched in monomers of concave shaped proteins. We further observed that proteins which are symmetric with respect to the bilayer center prefer symmetric Lo or Ld domains, while asymmetric proteins favor a location in domains with Lo/Ld asymmetry. In the latter case we additionally retrieved a strong dependence on protein directionality, thus providing a mechanism for transmembrane protein orientation.

Molecular simulations of glycolipids: Towards mammalian cell membrane models

Glycolipids are key components of mammalian cell membranes, influencing a diverse range of cellular functions. For example, a number of receptor tyrosine kinases, including the epidermal growth factor receptor (EGFR), are allosterically regulated by the glycolipid monosialodihexosylganglioside (GM3). Recent advances in molecular dynamics methods, especially the development of coarse-grained models, have enabled simulations of increasingly complex models of cell membranes. We demonstrate these methodological developments via a case study of a coarse-grained model for the ganglioside GM3. This glycolipid is included in simulations of a mixed lipid bilayer model reflecting the compositional complexity of a mammalian cell membrane. The resultant membrane model is used to simulate the interactions of GM3 with the transmembrane domain of the EGFR.

•

We review recent progress in molecular dynamics simulations of glycolipids.

•

A coarse-grained model of the ganglioside GM3 is described.

•

The GM3 model is used in coarse-grained simulations of a mammalian cell membrane.

•

We describe the interactions of GM3 with the transmembrane domain of the EGFR.

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

Cross-linking of specific lipid components by proteins mediates transmembrane signaling and material transport. In this work, we conducted coarse-grained simulation to investigate the interactions of binding units of chorela toxin (CTB) with mixed ganglioside GM1 and dipalmitoylphosphatidylcholine (DPPC) lipid bilayer membrane. We determine that the binding of CTB pentamers cross-links GM1 molecules into protein-sized nanodomains that have distinct lipid order compared with the bulk. The toxin in the nanodomain partially penetrates into the membrane. The local disordering can also transmit across the membrane via lipid coupling. Comparison simulations on CTB binding to a membrane that is composed of various lipid components demonstrate that several factors are responsible for the nanodomain formation: (a) the negatively charged headgroup of a GM1 receptor is responsible for the multivalent binding; (b) the head groups being full of hydrogen-bonding donors and receptors stabilize the GM1 cluster itself and ensure the toxin binding with high affinity; and

Biomembranes in atomistic and coarse-grained simulations

The architecture of biological membranes is tightly coupled to the localization, organization, and function of membrane proteins. The organelle-specific distribution of lipids allows for the formation of functional microdomains (also called rafts) that facilitate the segregation and aggregation of membrane proteins and thus shape their function. Molecular dynamics simulations enable to directly access the formation, structure, and dynamics of membrane microdomains at the molecular scale and the specific interactions among lipids and proteins on timescales from picoseconds to microseconds. This review focuses on the latest developments of biomembrane force fields for both atomistic and coarse-grained molecular dynamics (MD) simulations, and the different levels of coarsening of biomolecular structures. It also briefly introduces scale-bridging methods applicable to biomembrane studies, and highlights selected recent applications.

Dynamics of the Glycophorin A Dimer in Membranes of Native-Like Composition Uncovered by Coarse-Grained Molecular Dynamics Simulations

Membranes are central for cells as borders to the environment or intracellular organelle definition. They are composed of and harbor different molecules like various lipid species and sterols, and they are generally crowded with proteins. The membrane system is very dynamic and components show lateral, rotational and translational diffusion. The consequence of the latter is that phase separation can occur in membranes in vivo and in vitro. It was documented that molecular dynamics simulations of an idealized plasma membrane model result in formation of membrane areas where either saturated lipids and cholesterol (liquid-ordered character, L_o) or unsaturated lipids (liquid-disordered character, L_d) were enriched. Furthermore, current discussions favor the idea that proteins are sorted into the liquid-disordered phase of model membranes, but experimental support for the behavior of isolated proteins in native membranes is sparse. To gain insight into the protein behavior we built a model of the red blood cell membrane with integrated glycophorin A dimer. The sorting and the dynamics of the dimer were subsequently explored by coarse-grained molecular dynamics simulations. In addition, we inspected the impact of lipid head groups and the presence of cholesterol within the membrane on the dynamics of the dimer within the membrane. We observed that cholesterol is important for the formation of membrane areas with L_o and L_d character. Moreover, it is an important factor for the reproduction of the dynamic behavior of the protein found in its native environment. The protein dimer was exclusively sorted into the domain of L_d character in the model red blood cell plasma membrane. Therefore, we present structural information on the glycophorin A dimer distribution in the plasma membrane in the absence of other factors like e.g. lipid anchors in a coarse grain resolution.