Identifying the Onset of Phase Separation in Quaternary Lipid Bilayer Systems from Coarse-Grained Simulations

Dynamic Simulations of Interaction of the PEG-DPPE Micelle-Encapsulated Short-Chain Ceramides with the Raft-Included Membrane

The lipid raft subdomains in cancer cell membranes play a key role in signal transduction, biomolecule recruitment, and drug transmembrane transport. Augmented membrane rigidity due to the formation of a lipid raft is unfavorable for the entry of drugs, a limiting factor in clinical oncology. The short-chain ceramide (CER) has been reported to promote drug entry into membranes and disrupt lipid raft formation, but the underlying mechanism is not well understood. We recently explored the carrier-membrane fusion dynamics of PEG-DPPE micelles in delivering doxorubicin (DOX). Based on the phase-segregated membrane model composed of DPPC/DIPC/CHOL/GM1/PIP2, we aim to explore the dynamic mechanism of the PEG-DPPE micelle-encapsulating DOXs in association with the raft-included cell membrane modulated by C8 acyl tail CERs. The results show that the lipid raft remains integrated and DOX-resistant subjected to free DOXs and the micelle-encapsulating ones. Addition of CERs disorganizes the lipid raft by pushing CHOL aside from DPPC. It subsequently allows for a good permeability for PEG-DPPE micelle-encapsulated DOXs, which penetrate deeper as CER concentration increases. GM1 is significant in guiding drugs' redistributing between bilayer phases, and the anionic PIP2 further helps DOXs attain the inner bilayer surface. These results elaborate on the perturbing effect of CERs on lipid raft stability, which provides a new comprehensive approach for further design of drug delivery systems.

Atomistic Insight into the Lipid Nanodomains of Synaptic Vesicles

Membrane curvature, once regarded as a passive consequence of membrane composition and cellular architecture, has been shown to actively modulate various properties of the cellular membrane. These changes could also lead to segregation of the constituents of the membrane, generating nanodomains with precise biological properties. Proteins often linked with neurodegeneration (e.g., tau, alpha-synuclein) exhibit an unintuitive affinity for synaptic vesicles in neurons, which are reported to lack distinct, ordered nanodomains based on their composition. In this study, all-atom molecular dynamics simulations are used to study a full-scale synaptic vesicle of realistic Gaussian curvature and its effect on the membrane dynamics and lipid nanodomain organization. Compelling indicators of nanodomain formation, from the perspective of composition, surface areas per lipid, order parameter, and domain lifetime, are identified in the vesicle membrane, which are absent in a flat bilayer of the same lipid composition. Therefore, our study supports the idea that curvature may induce phase separation in an otherwise fluid, disordered membrane.

Martini 3 Coarse-Grained Force Field for Cholesterol

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. The proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

Simulation study of domain formation in a model bacterial membrane

Recent experimental studies revealed that functional membrane microdomains (FMMs) are formed in prokaryotic cells which are structurally and functionally similar to the lipid rafts formed in eukaryotic cells. In this study, we employ coarse-grained molecular dynamics simulations to investigate the mechanism of domain formation and its physiochemical properties in a model methicillin-resistant staphylococcus aureus (MRSA) cell membrane. We find that domains are formed through lateral segregation of staphyloxanthin (STX), a carotenoid which shields the bacteria from the host's immune because of its antioxidant nature. Simulation results suggest that membrane integrity increases with the size of the domain, which is assessed by computing bond order parameter of the lipid tails, membrane expansion modulus and water permeability across the membrane. Various membrane domain proteins such as flotillin-like protein floA and penicillin binding protein (PBP2a) preferentially bind with the STX and accumulate in the membrane domain which is consistent with the recent experimental results.

The effects of mutation on the drug binding affinity of Neuraminidase: case study of Capsaicin using steered molecular dynamics simulation

The influenza virus is an important respiratory pathogen that causes many incidences of diseases and even death each year. One of the primary factors of this virus is the Neuraminidase surface protein, which causes the virus to leave the host cell and spread to new target cells. The main antiviral medication for influenza is designed as a protein inhibitor ligand that prevents further spread of the disease, and eventually relieves the emerged symptoms. The effectiveness of such inhibitory drugs is highly associated with their binding affinity. In this paper, the binding affinity of an herbal ligand of Capsaicin bound to Neuraminidase of the influenza virus is investigated using steered molecular dynamics (SMD) simulation. Since mutations of the virus directly impact the binding affinity of the inhibitory drugs, different mutations were generated by using Mutagenesis module. The rapid spread of infection during the avian influenza A/H5N1 epidemic has raised concerns about far more dangerous consequences if the virus becomes resistant to current drugs. Currently, oseltamivir (Tamiflu), zanamivir (Relenza), pramivir (Rapivab), and laninamivir (Inavir) are increasingly used to treat the flu. However, with the rapid evolution of the virus, some drug-resistant strains are emerging. Therefore, it is very important to seek alternative therapies and identify the roots of drug resistance. Obtained results demonstrated a reduced binding affinity for the applied mutations. This reduction in binding affinity will cause the virus mutation to become resistant to the drug, which will spread the disease and make it more difficult to treat. From a molecular prospect, this decrease in binding affinity is due to the loss of a number of effective bonds between the ligand and the receptor, which occurs with mutations of the wild-type (WT) species. The results of the present study can be used in the rational design of novel drugs that are compatible with specific mutations.

Study of lipid heterogeneity on bilayer membranes using molecular dynamics simulations

Human cell membranes consist of various lipids that are essential for their structure and function. It typically comprises phosphatidylcholine (POPC), phosphatidylethanolamine (POPE), phosphatidylserine (POPS), sphingomyelin (PSM), and cholesterol (CHL). Several experimental and computational techniques have been employed to characterize the composition of human cell membranes, however, CHL enriched membrane is still not clearly understood through these techniques. Molecular dynamics simulation results illustrated the biophysical properties of heterogeneous membranes based on the lipid composition as well as the concentration of lipids, exclusively for CHL and PSM. Herein, we have investigated the structure-function relationships of lipids comparatively to delineate the effect of heterogeneity on the biophysical properties of different membranes. It has been observed that the significant fraction of CHL (i.e., ~33% in ternary, ~25% in quaternary, and ~16% in senary type bilayers) in combination with other lipids introduced compactness, and increased the thickness of the membrane. The analysis of lipid mass density stated that the density of lipid head group, phosphate_, and glycerol-ester in presence of CHL with or without PSM is an underlying reason for membrane ordering. Results also revealed that the presence of POPI and POPS are the reasons for an adequate drop in the ordering of lipid chain, particularly on POPE chain. The self-interaction of CHL, PSM, POPE and the interaction of CHL and POPC with POPE seem to determine the structure and function of the heterogeneous membrane. Our findings provide a qualitative understanding of the effect of membrane heterogeneity on the physiological properties of membranes. The structures inspected in this study would help to select the heterogeneous bilayer model to mimic the human cell membranes to analyse or characterize the membrane-associated phenomena.

Effect of alcohol on the phase separation in model membranes

The discovery of coexisting liquid-ordered and liquid-disordered phases in multicomponent lipid bilayers has received widespread attention due to its potential relevance for biological systems. One of the many open questions is how the presence of additional components affects the nature of the coexisting phases. Of particular interest is the addition of alcohols because their anesthetic properties may arise from modulating bilayer behavior. We use coarse-grained Molecular Dynamics simulations to gain insight into the partitioning preferences of linear n-alcohols into ordered and disordered bilayers alongside their effects on local membrane structure. We find that alcohols cause only small changes to membrane composition alongside a lack of significant effects on membrane thickness and lipid tail order. Cholesterol and n-alcohol trans-bilayer motion is measured and found to be near or within the range of previous atomistic results. The cholesterol flip-flop rates increase with both n-alcohol length and concentration for octanol, dodecanol, and hexadecanol, indicating a decrease in lipid order. Umbrella sampling simulations of removing cholesterol from tertiary membranes find no significant difference with or without n-alcohols at various concentrations. Simulations of a phase-separated bilayer show that octanol preferentially partitions into the liquid-disordered phase in a ratio of approximately 3:1 over the liquid-ordered phase. Furthermore, partition coefficients of alcohol in single-phase membranes show a preference for longer alcohols (dodecanol and hexadecanol) to partition preferentially into the liquid-ordered phase, while decreasing the length of the alcohol reverses this trend. Our work tests experimental results while also investigating the ability for coarse-grained MARTINI simulations to capture minute differences in model membrane spatial arrangements on the nanoscale level.

Modulated and spiral surface patterns on deformable lipid vesicles

We investigate the behavior of two-dimensional systems that exhibit a transition between homogeneous and spatially inhomogeneous phases, which have spherical topology, and whose mechanical properties depend on the local value of the order parameter. One example of such a system is multicomponent lipid bilayer vesicles, which serve as a model to study cellular membranes. Under certain conditions, such bilayers separate into coexisting liquid-ordered and liquid-disordered regions. When arranged into the shape of small vesicles, this phase coexistence can result in spatial patterns that are more complex than the basic two-domain configuration encountered in typical bulk systems. The difference in bending rigidity between the liquid-ordered and liquid-disordered regions couples the shape of the vesicle to the local composition. We show that this interplay gives rise to a rich phase diagram that includes homogeneous, separated, and axisymmetric modulated phases that are divided by regions of spiral patterns in the surface morphology.

Laterally Resolved Small-Angle Scattering Intensity from Lipid Bilayer Simulations: An Exact and a Limited-Range Treatment

When combined, molecular simulations and small-angle scattering experiments are able to provide molecular-scale resolution of structure. Separately, scattering experiments provide only intermingled pair correlations between atoms, while molecular simulations are limited by model quality and the relatively short time scales that they can access. Their combined strength relies on agreement between the experimental spectra and those computed by simulation. To date, computing the neutron spectra from a molecular simulation of a lipid bilayer is straightforward only if the structure is approximated by laterally averaging the in-plane bilayer structure. However, this neglects all information about lateral heterogeneity, e.g., clustering of components in a lipid mixture. This paper presents two methods for computing the scattering intensity of simulated bilayers with in-plane heterogeneity, enabling a full treatment of both the transverse and lateral bilayer structure for the first time. The first method, termed the Dirac Brush, computes the exact spectra including spurious artifacts resulting from using information from neighboring periodic cells to account for the long-range structure of the bilayer. The second method, termed PFFT, applies a mean-field treatment in the field far from a scattering element, resulting in a correlation range that can be tuned (eliminating correlations with neighboring periodic images), but with computational cost that prohibits obtaining the exact (Dirac Brush) spectra. Following their derivation, the two methods are applied to a coarse-grained molecular simulation of a bilayer inhomogeneity, demonstrating the contributions of lateral correlations to the resulting spectra.

Inserting Small Molecules across Membrane Mixtures: Insight from the Potential of Mean Force

Small solutes have been shown to alter the lateral organization of cell membranes and reconstituted phospholipid bilayers; however, the mechanisms by which these changes happen are still largely unknown. Traditionally, both experiment and simulation studies have been restricted to testing only a few compounds at a time, failing to identify general molecular descriptors or chemical properties that would allow extrapolating beyond the subset of considered solutes. In this work, we probe the competing energetics of inserting a solute in different membrane environments by means of the potential of mean force. We show that these calculations can be used as a computationally efficient proxy to establish whether a solute will stabilize or destabilize domain phase separation. Combined with umbrella-sampling simulations and coarse-grained molecular dynamics simulations, we are able to screen solutes across a wide range of chemistries and polarities. Our results indicate that for the system under consideration, preferential partitioning and therefore effectiveness in altering membrane phase separation are strictly linked to the location of insertion in the bilayer (i.e., midplane or interface). Our approach represents a fast and simple tool for obtaining structural and thermodynamic insight into the partitioning of small molecules between lipid domains and its relation to phase separation, ultimately providing a platform for identifying the key determinants of this process.

Deciphering Melatonin-Stabilized Phase Separation in Phospholipid Bilayers

Lipid bilayers are fundamental building blocks of cell membranes, which contain the machinery needed to perform a range of biological functions, including cell-cell recognition, signal transduction, receptor trafficking, viral budding, and cell fusion. Importantly, many of these functions are thought to take place in the laterally phase-separated regions of the membrane, commonly known as lipid rafts. Here, we provide experimental evidence for the "stabilizing" effect of melatonin, a naturally occurring hormone produced by the brain's pineal gland, on phase-separated model membranes mimicking the outer leaflet of plasma membranes. Specifically, we show that melatonin stabilizes the liquid-ordered/liquid-disordered phase coexistence over an extended range of temperatures. The melatonin-mediated stabilization effect is observed in both nanometer- and micrometer-sized liposomes using small angle neutron scattering (SANS), confocal fluorescence microscopy, and differential scanning calorimetry. To experimentally detect nanoscopic domains in 50 nm diameter phospholipid vesicles, we developed a model using the Landau-Brazovskii approach that may serve as a platform for detecting the existence of nanoscopic lateral heterogeneities in soft matter and biological materials with spherical and planar geometries.

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.

Regimes of Complex Lipid Bilayer Phases Induced by Cholesterol Concentration in MD Simulation

Cholesterol is essential to the formation of phase-separated lipid domains in membranes. Lipid domains can exist in different thermodynamic phases depending on the molecular composition and play significant roles in determining structure and function of membrane proteins. We investigate the role of cholesterol in the structure and dynamics of ternary lipid mixtures displaying phase separation using molecular dynamics simulations, employing a physiologically relevant span of cholesterol concentration. We find that cholesterol can induce formation of three regimes of phase behavior: 1) miscible liquid-disordered bulk, 2) phase-separated, domain-registered coexistence of liquid-disordered and liquid-ordered domains, and 3) phase-separated, domain-antiregistered coexistence of liquid-disordered and newly identified nanoscopic gel domains composed of cholesterol threads we name "cholesterolic gel" domains. These findings are validated and discussed in the context of current experimental knowledge, models of cholesterol spatial distributions, and models of ternary lipid-mixture phase separation.

Tuning Length Scales of Small Domains in Cell-Derived Membranes and Synthetic Model Membranes

Micron-scale, coexisting liquid-ordered (L_o) and liquid-disordered (L_d) phases are straightforward to observe in giant unilamellar vesicles (GUVs) composed of ternary lipid mixtures. Experimentally, uniform membranes undergo demixing when temperature is decreased: domains subsequently nucleate, diffuse, collide, and coalesce until only one domain of each phase remains. The sizes of these two domains are limited only by the size of the system. Under different conditions, vesicles exhibit smaller-scale domains of fixed sizes, leading to the question of what sets the length scale. In membranes with excess area, small domains are expected when coarsening is hindered or when a microemulsion or modulated phase arises. Here, we test predictions of how the size, morphology, and fluorescence levels of small domains vary with the membrane's temperature, tension, and composition. Using GUVs and cell-derived giant plasma membrane vesicles, we find that 1) the characteristic size of domains decreases when temperature is increased or membrane tension is decreased, 2) stripes are favored over circular domains for lipid compositions with low energy per unit interface, 3) fluorescence levels are consistent with domain registration across both monolayer leaflets of the bilayer, and 4) small domains form in GUVs composed of lipids both with and without ester-linked lipids. Our experimental results are consistent with several elements of current theories for microemulsions and modulated phases and inconsistent with others, suggesting a motivation to modify or enhance current theories.

Relating the structure factors of two-dimensional materials in planar and spherical geometries

Scattering structure factors provide essential insight into material properties and are routinely obtained in experiments, computer simulations, and theoretical analyses. Different approaches favor different geometries of the material. In case of lipid bilayers, scattering experiments can be performed on spherical vesicles, while simulations and theory often consider planar membrane patches. We derive an approximate relationship between the structure functions of such a material in planar and spherical geometries. We illustrate the usefulness of this relationship in a case study of a Gaussian material that supports both homogeneous and microemulsion phases. Within its range of applicability, this relationship enables a model-free comparison of structure factors of the same material in different geometries.