PackMem: A Versatile Tool to Compute and Visualize Interfacial Packing Defects in Lipid Bilayers

From complex data to clear insights: visualizing molecular dynamics trajectories

Advances in simulations, combined with technological developments in high-performance computing, have made it possible to produce a physically accurate dynamic representation of complex biological systems involving millions to billions of atoms over increasingly long simulation times. The analysis of these computed simulations is crucial, involving the interpretation of structural and dynamic data to gain insights into the underlying biological processes. However, this analysis becomes increasingly challenging due to the complexity of the generated systems with a large number of individual runs, ranging from hundreds to thousands of trajectories. This massive increase in raw simulation data creates additional processing and visualization challenges. Effective visualization techniques play a vital role in facilitating the analysis and interpretation of molecular dynamics simulations. In this paper, we focus mainly on the techniques and tools that can be used for visualization of molecular dynamics simulations, among which we highlight the few approaches used specifically for this purpose, discussing their advantages and limitations, and addressing the future challenges of molecular dynamics visualization.

Analyzing Lipid Membrane Defects via a Coarse-Grained to Triangulated Surface Map: The Role of Lipid Order and Local Curvature in Molecular Binding

Lipid packing defects are known to serve as quantitative indicators for protein binding to lipid membranes. This paper presents a protocol for mapping molecular lipid detail onto a triangulated continuum leaflet representation. Besides establishing the desired forward counterpart to the existing inverse TS2CG map, this coarse-grained to triangulated surface (CG2TS) map enables straightforward extraction of the defect characteristics for any membrane geometry found in nature. We have applied our protocol to investigate the role of local curvature and varying lipid packing on the defect constant π. We find that the defect size is greatly influenced by both factors, arguing strongly against the usual assignment of a single defect constant in the case of more realistic membrane conditions. An important discovery is that lipids in the gel phase produce larger defects, or a higher π, in domains of high (local) curvature than the same lipid in a liquid phase of any curvature. This finding suggests that membranes featuring very ordered lipid packing can bind proteins via large defects in curved regions. Finally, we propose a route for estimating defect constants directly from the standard membrane properties. Identifying the precise role of composition, lipid (tail) order, and (local) curvature in defects for the irregular lipid structures that are (temporally) present in many biological processes is instrumental for obtaining fundamental insight as well as for a rational design of membrane binding targets.

MemPrep, a new technology for isolating organellar membranes provides fingerprints of lipid bilayer stress

Biological membranes have a stunning ability to adapt their composition in response to physiological stress and metabolic challenges. Little is known how such perturbations affect individual organelles in eukaryotic cells. Pioneering work has provided insights into the subcellular distribution of lipids in the yeast Saccharomyces cerevisiae, but the composition of the endoplasmic reticulum (ER) membrane, which also crucially regulates lipid metabolism and the unfolded protein response, remains insufficiently characterized. Here, we describe a method for purifying organelle membranes from yeast, MemPrep. We demonstrate the purity of our ER membrane preparations by proteomics, and document the general utility of MemPrep by isolating vacuolar membranes. Quantitative lipidomics establishes the lipid composition of the ER and the vacuolar membrane. Our findings provide a baseline for studying membrane protein biogenesis and have important implications for understanding the role of lipids in regulating the unfolded protein response (UPR). The combined preparative and analytical MemPrep approach uncovers dynamic remodeling of ER membranes in stressed cells and establishes distinct molecular fingerprints of lipid bilayer stress.

Structural diversity in the membrane-bound hIAPP dimer correlated with distinct membrane disruption mechanisms

Amyloid deposits of the human islet amyloid polypeptide (hIAPP) have been identified in 90% of patients with type II diabetes. Cellular membranes accelerate the hIAPP fibrillation, and the integrity of membranes is also disrupted at the same time, leading to the apoptosis of β cells in pancreas. The molecular mechanism of hIAPP-induced membrane disruption, especially during the initial membrane disruption stage, has not been well understood yet. Herein, we carried out extensive all-atom molecular dynamics simulations investigating the hIAPP dimerization process in the anionic POPG membrane, to provide the detailed molecular mechanisms during the initial hIAPP aggregation stage in the membrane environment. Compared to the hIAPP monomer on the membrane, we observed not only an increase of α-helical structures, but also a substantial increase of β-sheet structures upon spontaneous dimerization. Moreover, the random coiled and α-helical dimer structures insert deep into the membrane interior with a few inter-chain contacts at the C-terminal region, while the β-sheet-rich structures reside on the membrane surface accompanied by strong inter-chain hydrophobic interactions. The coexistence of α and β structures constitutes a diverse structural ensemble of the membrane-bound hIAPP dimer. From α-helical to β-sheet structures, the degree of membrane disruption decreases gradually, and thus the membrane damage induced by random coiled and α-helical structures precedes that induced by β-sheet structures. We speculate that insertion of random coiled and α-helical structures contributes to the initial stage of membrane damage, while β-sheet structures on the membrane surface are more involved in the later stage of fibril-induced membrane disruption.

Phase-Separated Lipid-Based Nanoparticles: Selective Behavior at the Nano-Bio Interface

The membrane-protein interface on lipid-based nanoparticles influences their in vivo behavior. Better understanding may evolve current drug delivery methods toward effective targeted nanomedicine. Previously, the cell-selective accumulation of a liposome formulation in vivo is demonstrated, through the recognition of lipid phase-separation by triglyceride lipases. This exemplified how liposome morphology and composition can determine nanoparticle-protein interactions. Here, the lipase-induced compositional and morphological changes of phase-separated liposomes-which bear a lipid droplet in their bilayer- are investigated, and the mechanism upon which lipases recognize and bind to the particles is unravelled. The selective lipolytic degradation of the phase-separated lipid droplet is observed, while nanoparticle integrity remains intact. Next, the Tryptophan-rich loop of the lipase is identified as the region with which the enzymes bind to the particles. This preferential binding is due to lipid packing defects induced on the liposome surface by phase separation. In parallel, the existing knowledge that phase separation leads to in vivo selectivity, is utilized to generate phase-separated mRNA-LNPs that target cell-subsets in zebrafish embryos, with subsequent mRNA delivery and protein expression. Together, these findings can expand the current knowledge on selective nanoparticle-protein communications and in vivo behavior, aspects that will assist to gain control of lipid-based nanoparticles.

Role of Lipid Packing Defects in Determining Membrane Interactions of Antimicrobial Polymers

Understanding the emergence and role of lipid packing defects in the detection and subsequent partitioning of antimicrobial agents into bacterial membranes is essential for gaining insights into general antimicrobial mechanisms. Herein, using methacrylate polymers as a model platform, we investigate the effects of inclusion of various functional groups in the biomimetic antimicrobial polymer design on the aspects of lipid packing defects in model bacterial membranes. Two antimicrobial polymers are considered: ternary polymers composed of cationic, hydrophobic, and polar moieties and binary polymers with only cationic and hydrophobic moieties. We find that differing modes of insertion of these two polymers lead to different packing defects in the bacterial membrane. While insertion of both binary and ternary polymers leads to an enhanced number of deep defects in the upper leaflet, shallow defects are moderately enhanced upon interaction with ternary polymers only. We provide conclusive evidence that insertion of antimicrobial polymers in bacterial membrane is preceded by sensing of interfacial lipid packing defects. Our simulation results show that the hydrophobic groups are inserted at a single colocalized deep defect site for both binary and ternary polymers. However, the presence of polar groups in the ternary polymers use the shallow defects close to the lipid-water interface, in addition, to insert into the membrane, which leads to a more folded conformation of the ternary polymer in the membrane environment, and hence a different membrane partitioning mechanism compared to the binary polymer, which acquires an amphiphilic conformation.

Influence of Cholesterol on the Membrane Binding and Conformation of α-Synuclein

The α-Synuclein (α-Syn) plays an important role in the pathology of Parkinson's disease (PD), and its oligomers and fibrils are toxic to the nervous system. As organisms age, the cholesterol content in biological membranes increases, which is a potential cause of PD. Cholesterol may affect the membrane binding of α-Syn and its abnormal aggregation, but the mechanism remains unclear. Here, we present our molecular dynamics simulation studies on the interaction between α-Syn and lipid membranes, with or without cholesterol. It is demonstrated that cholesterol provides additional hydrogen bond interaction with α-Syn; however, the coulomb interaction and hydrophobic interaction between α-Syn and lipid membranes could be weakened by cholesterol. In addition, cholesterol leads to the shrinking of lipid packing defects and the decrease of lipid fluidity, thereby shortening the membrane binding region of α-Syn. Under these multifaceted effects of cholesterol, membrane-bound α-Syn shows signs of forming a β-sheet structure, which may further induce the formation of abnormal α-Syn fibrils. These results provide important information for the understanding of membrane binding of α-Syn, and they are expected to promote the bridging between cholesterol and the pathological aggregation of α-Syn.

Confinement effect on hydrolysis in small lipid vesicles

In living organisms most chemical reactions take place within the confines of lipid-membrane bound compartments, while confinement within the bounds of a lipid membrane is thought to be a key step in abiogenesis. In previous work we demonstrated that confinement in the aqueous cavity of a lipid vesicle affords protection against hydrolysis, a phenomenon that we term here confinement effect (C _e) and that we attributed to the interaction with the lipid membrane. Here, we show that both the size and the shape of the cavity of the vesicle modulate the C _e. We link this observation to the packing of the lipid following changes in membrane curvature, and formulate a mathematical model that relates the C _e to the radius of a spherical vesicle and the packing parameter of the lipids. These results suggest that the shape of the compartment where a molecule is located plays a major role in controlling the chemical reactivity of non-enzymatic reactions. Moreover, the mathematical treatment we propose offers a useful tool for the design of vesicles with predictable reaction rates of the confined molecules, e.g., drug delivery vesicles with confined prodrugs. The results also show that a crude form of signal transduction, devoid of complex biological machinery, can be achieved by any external stimuli that drastically changes the structure of the membrane, like the osmotic shocks used in the present work.

“Interplay of lipid-head group and packing defects in driving Amyloid-beta mediated myelin-like model membrane deformation”

Accumulating evidence suggests that amyloid plaque associated myelin lipid loss as a result of elevated amyloid burden might also contribute to Alzheimer's disease. The amyloid fibrils though closely associated with lipids under physiological conditions, however, the progression of membrane remodeling events leading to lipid-fibril assembly remains unknown. Here we first reconstitute the interaction of Aβ-40 with myelin-like model membrane and show that the binding of Aβ-40 induces extensive tubulation. To look into the mechanism of membrane tubulation we chose a set of membrane conditions varying in lipid packing density and net charge that allows us to dissect the contribution of lipid specificity of Aβ-40 binding, aggregation kinetics, and subsequent changes in membrane parameters such as fluidity, diffusion, and compressibility modulus. We show that the binding of Aβ-40 depends predominantly on the lipid packing defect densities and electrostatic interactions and results in rigidification of the myelin-like model membrane during the early phase of amyloid aggregation. Furthermore, elongation of Aβ-40 into higher oligomeric and fibrillar species leads to eventual fluidization of the model membrane followed by extensive lipid membrane tubulation observed in the late phase. Taken together, our results capture mechanistic insights into snapshots of temporal dynamics of Aβ-40 - myelin-like model membrane interaction and demonstrate how short timescale, local phenomena of binding, and fibril-mediated load generation results in the consequent association of lipids with growing amyloid fibrils.

Lessons for Oral Bioavailability: How Conformationally Flexible Cyclic Peptides Enter and Cross Lipid Membranes

Cyclic peptides extend the druggable target space due to their size, flexibility, and hydrogen-bonding capacity. However, these properties impact also their passive membrane permeability. As the "journey" through membranes cannot be monitored experimentally, little is known about the underlying process, which hinders rational design. Here, we use molecular simulations to uncover how cyclic peptides permeate a membrane. We show that side chains can act as "molecular anchors", establishing the first contact with the membrane and enabling insertion. Once inside, the peptides are positioned between headgroups and lipid tails─a unique polar/apolar interface. Only one of two distinct orientations at this interface allows for the formation of the permeable "closed" conformation. In the closed conformation, the peptide crosses to the lower leaflet via another "anchoring" and flipping mechanism. Our findings provide atomistic insights into the permeation process of flexible cyclic peptides and reveal design considerations for each step of the process.

PEO-b-PBD Diblock Copolymers Induce Packing Defects in Lipid/Hybrid Membranes and Improve Insertion Rates of Natively Folded Peptides

Hybrid membranes assembled from biological lipids and synthetic polymers are a promising scaffold for the reconstitution and utilization of membrane proteins. Recent observations indicate that inclusion of small fractions of polymer in lipid membranes can improve protein folding and function, but the exact structural and physical changes a given polymer sequence imparts on a membrane often remain unclear. Here, we use all-atom molecular dynamics simulations to study the structure of hybrid membranes assembled from DOPC phospholipids and PEO-b-PBD diblock copolymers. We verified our computational model using new and existing experimental data and obtained a detailed picture of the polymer conformations in the lipid membrane that we can relate to changes in membrane elastic properties. We find that inclusion of low polymer fractions induces transient packing defects into the membrane. These packing defects act as insertion sites for two model peptides, and in this way, small amounts of polymer content in lipid membranes can lead to large increases in peptide insertion rates. Additionally, we report the peptide conformational space in both pure lipid and hybrid membranes. Both membranes support similar alpha helical peptide structures, exemplifying the biocompatibility of hybrid membranes.

Triglyceride lipolysis triggers liquid crystalline phases in lipid droplets and alters the LD proteome

Lipid droplets (LDs) are reservoirs for triglycerides (TGs) and sterol-esters (SEs), but how these lipids are organized within LDs and influence their proteome remain unclear. Using in situ cryo-electron tomography, we show that glucose restriction triggers lipid phase transitions within LDs generating liquid crystalline lattices inside them. Mechanistically this requires TG lipolysis, which decreases the LD's TG:SE ratio, promoting SE transition to a liquid crystalline phase. Molecular dynamics simulations reveal TG depletion promotes spontaneous TG and SE demixing in LDs, additionally altering the lipid packing of the PL monolayer surface. Fluorescence imaging and proteomics further reveal that liquid crystalline phases are associated with selective remodeling of the LD proteome. Some canonical LD proteins, including Erg6, relocalize to the ER network, whereas others remain LD-associated. Model peptide LiveDrop also redistributes from LDs to the ER, suggesting liquid crystalline phases influence ER-LD interorganelle transport. Our data suggests glucose restriction drives TG mobilization, which alters the phase properties of LD lipids and selectively remodels the LD proteome.

Capturing the Liquid-Crystalline Phase Transformation: Implications for Protein Targeting to Sterol Ester-Rich Lipid Droplets

Lipid droplets are essential organelles that store and traffic neutral lipids. The phospholipid monolayer surrounding their neutral lipid core engages with a highly dynamic proteome that changes according to cellular and metabolic conditions. Recent work has demonstrated that when the abundance of sterol esters increases above a critical concentration, such as under conditions of starvation or high LDL exposure, the lipid droplet core can undergo an amorphous to liquid-crystalline phase transformation. Herein, we study the consequences of this transformation on the physical properties of lipid droplets that are thought to regulate protein association. Using simulations of different sterol-ester concentrations, we have captured the liquid-crystalline phase transformation at the molecular level, highlighting the alignment of sterol esters in alternating orientations to form concentric layers. We demonstrate how ordering in the core permeates into the neutral lipid/phospholipid interface, changing the magnitude and nature of neutral lipid intercalation and inducing ordering in the phospholipid monolayer. Increased phospholipid packing is concomitant with altered surface properties, including smaller area per phospholipid and substantially reduced packing defects. Additionally, the ordering of sterol esters in the core causes less hydration in more ordered regions. We discuss these findings in the context of their expected consequences for preferential protein recruitment to lipid droplets under different metabolic conditions.

Efficient Quantification of Lipid Packing Defect Sensing by Amphipathic Peptides: Comparing Martini 2 and 3 with CHARMM36

In biological systems, proteins can be attracted to curved or stretched regions of lipid bilayers by sensing hydrophobic defects in the lipid packing on the membrane surface. Here, we present an efficient end-state free energy calculation method to quantify such sensing in molecular dynamics simulations. We illustrate that lipid packing defect sensing can be defined as the difference in mechanical work required to stretch a membrane with and without a peptide bound to the surface. We also demonstrate that a peptide's ability to concurrently induce excess leaflet area (tension) and elastic softening─a property we call the "characteristic area of sensing" (CHAOS)─and lipid packing sensing behavior are in fact two sides of the same coin. In essence, defect sensing displays a peptide's propensity to generate tension. The here-proposed mechanical pathway is equally accurate yet, computationally, about 40 times less costly than the commonly used alchemical pathway (thermodynamic integration), allowing for more feasible free energy calculations in atomistic simulations. This enabled us to directly compare the Martini 2 and 3 coarse-grained and the CHARMM36 atomistic force fields in terms of relative binding free energies for six representative peptides including the curvature sensor ALPS and two antiviral amphipathic helices (AH). We observed that Martini 3 qualitatively reproduces experimental trends while producing substantially lower (relative) binding free energies and shallower membrane insertion depths compared to atomistic simulations. In contrast, Martini 2 tends to overestimate (relative) binding free energies. Finally, we offer a glimpse into how our end-state-based free energy method can enable the inverse design of optimal lipid packing defect sensing peptides when used in conjunction with our recently developed evolutionary molecular dynamics (Evo-MD) method. We argue that these optimized defect sensors─aside from their biomedical and biophysical relevance─can provide valuable targets for the development of lipid force fields.

Computational Studies of Lipid Droplets

Lipid droplets (LDs) are intracellular organelles whose primary function is energy storage. Known to emerge from the endoplasmic reticulum (ER) bilayer, LDs have a unique structure with a core consisting of neutral lipids, triacylglycerol (TG) or sterol esters (SE), surrounded by a phospholipid (PL) monolayer and decorated by proteins that come and go throughout their complex lifecycle. In this Feature Article, we review recent developments in computational studies of LDs, a rapidly growing area of research. We highlight how molecular dynamics (MD) simulations have provided valuable molecular-level insight into LD targeting and LD biogenesis. Additionally, we review the physical properties of TG from different force fields compared with experimental data. Possible future directions and challenges are discussed.

Role of Disulphide Bonds in Membrane Partitioning of a Viral Peptide

The importance of disulphide bond in mediating viral peptide entry into host cells is well known. In the present work, we elucidate the role of disulphide (SS) bond in partitioning mechanism of membrane-active Hepatitis A Virus-2B (HAV-2B) peptide, which harbours three cysteine residues promoting formation of multiple SS-bonded states. The inclusion of SS-bond not only results in a compact conformation but also induces distorted α-helical hairpin geometry in comparison to SS-free state. Owing to these, the hydrophobic residues get buried, restricting the insertion of SS-bonded HAV-2B peptide into lipid packing defects and thus the partitioning of the peptide is completely or partly abolished. In this way, the disulphide bond can potentially regulate the partitioning of HAV-2B peptide such that the membrane remodelling effects of this viral peptide are significantly reduced. The current findings may have potential implications in drug designing, targeting the HAV-2B protein by promoting disulphide bond formation within its membrane-active region.

The role of size and nature in nanoparticle binding to a model lung membrane: an atomistic study

Understanding the uptake of nanoparticles (NPs) by different types of cellular membranes plays a pivotal role in the design of NPs for medical applications and in avoiding adverse effects that result in nanotoxicity. Yet, the role of key design parameters, such as the bare NP material, NP size and surface reactivity, and the nature of NP coatings, in membrane remodelling and uptake mechanisms is still very poorly understood, particularly towards the lower range of NP dimensions that are beyond the experimental imaging resolution. The same can be said about the role of a particular membrane composition. Here, we systematically employ biased and unbiased molecular dynamics simulations to calculate the binding energy for three bare materials (Ag/SiO₂/TiO₂) and three NP sizes (1/3/5 nm diameter) with a representative lung surfactant membrane, and to study their binding kinetics. The calculated binding energies show that irrespective of size, Ag nanoparticles bind very strongly to the bilayer, while the NPs made of SiO₂ or TiO₂ experience very low to no binding. The unbiased simulations provide insight into how the NPs and membrane affect each other in terms of the solvent-accessible surface area (SASA) of the NPs and the defect types and fluidity of the membrane. Using these systematic fine-grained results in coarsening procedures will pave the way for simulations considering NP sizes that are well beyond the membrane thickness, i.e. closer to experimental dimensions, for which different binding characteristics and more significant membrane remodelling are expected.

Cooperation of Conical and Polyunsaturated Lipids to Regulate Initiation and Processing of Membrane Fusion

The shape of lipids has long been suspected to be a critical determinant for the control of membrane fusion. To experimentally test this assertion, we used conical and malleable lipids and measured their influence on the fusion kinetics. We found that, as previously suspected, both types of lipids accelerate fusion. However, the implicated molecular mechanisms are strikingly different. Malleable lipids, with their ability to change shape with low energy cost, favor fusion by decreasing the overall activation energy. On the other hand, conical lipids, with their small polar head relative to the area occupied by the hydrophobic chains, tend to make fusion less energetically advantageous because they tend to migrate towards the most favorable lipid leaflet, hindering fusion pore opening. They could however facilitate fusion by generating hydrophobic defects on the membranes; this is suggested by the similar trend observed between the experimental rate of fusion nucleation and the surface occupied by hydrophobic defects obtained by molecular simulations. The synergy of dual-process, activation energy and nucleation kinetics, could facilitate membrane fusion regulation in vivo.

Effect of cholesterol on the membrane partitioning dynamics of hepatitis A virus-2B peptide

Understanding viral peptide detection and partitioning and the subsequent host membrane composition-based response is essential for gaining insights into the viral mechanism. Here, we probe the crucial role of the presence of membrane lipid packing defects, depending on the membrane composition, in allowing the viral peptide belonging to C-terminal Hepatitis A Virus-2B (HAV-2B) to detect, attach and subsequently partition into host cell membrane mimics. Using molecular dynamics simulations, we conclusively show that the hydrophobic residues in the viral peptide detect transiently present lipid packing defects, insert themselves into such defects, form anchor points and facilitate the partitioning of the peptide, thereby inducing membrane disruption. We also show that the presence of cholesterol significantly alters such lipid packing defects, both in size and in number, thus mitigating the partitioning of the membrane active viral peptide into cholesterol-rich membranes. Our results are in excellent agreement with previously published experimental data and further explain the role of lipid defects in understanding such data. These results show differential ways in which the presence and absence of cholesterol can alter the permeability of the host membranes to the membrane active peptide component of HAV-2B virus, via lipid packing defects, and can possibly be a part of the general membrane detection mechanism for viroporins.

Membrane hydrophobicity determines the activation free energy of passive lipid transport

The collective behavior of lipids with diverse chemical and physical features determines a membrane's thermodynamic properties. Yet, the influence of lipid physicochemical properties on lipid dynamics, in particular interbilayer transport, remains underexplored. Here, we systematically investigate how the activation free energy of passive lipid transport depends on lipid chemistry and membrane phase. Through all-atom molecular dynamics simulations of 11 chemically distinct glycerophospholipids, we determine how lipid acyl chain length, unsaturation, and headgroup influence the free energy barriers for two elementary steps of lipid transport: lipid desorption, which is rate limiting, and lipid insertion into a membrane. Consistent with previous experimental measurements, we find that lipids with longer, saturated acyl chains have increased activation free energies compared to lipids with shorter, unsaturated chains. Lipids with different headgroups exhibit a range of activation free energies; however, no clear trend based solely on chemical structure can be identified, mirroring difficulties in the interpretation of previous experimental results. Compared to liquid-crystalline phase membranes, gel phase membranes exhibit substantially increased free energy barriers. Overall, we find that the activation free energy depends on a lipid's local hydrophobic environment in a membrane and that the free energy barrier for lipid insertion depends on a membrane's interfacial hydrophobicity. Both of these properties can be altered through changes in lipid acyl chain length, lipid headgroup, and membrane phase. Thus, the rate of lipid transport can be tuned through subtle changes in local membrane composition and order, suggesting an unappreciated role for nanoscale membrane domains in regulating cellular lipid dynamics.

NMR identification of a conserved Drp1 cardiolipin-binding motif essential for stress-induced mitochondrial fission

Mitochondria form tubular networks that undergo coordinated cycles of fission and fusion. Emerging evidence suggests that a direct yet unresolved interaction of the mechanoenzymatic GTPase dynamin-related protein 1 (Drp1) with mitochondrial outer membrane-localized cardiolipin (CL), externalized under stress conditions including mitophagy, catalyzes essential mitochondrial hyperfragmentation. Here, using a comprehensive set of structural, biophysical, and cell biological tools, we have uncovered a CL-binding motif (CBM) conserved between the Drp1 variable domain (VD) and the unrelated ADP/ATP carrier (AAC/ANT) that intercalates into the membrane core to effect specific CL interactions. CBM mutations that weaken VD-CL interactions manifestly impair Drp1-dependent fission under stress conditions and induce "donut" mitochondria formation. Importantly, VD membrane insertion and GTP-dependent conformational rearrangements mediate only transient CL nonbilayer topological forays and high local membrane constriction, indicating that Drp1-CL interactions alone are insufficient for fission. Our studies establish the structural and mechanistic bases of Drp1-CL interactions in stress-induced mitochondrial fission.

Stressed Lipid Droplets: How Neutral Lipids Relieve Surface Tension and Membrane Expansion Drives Protein Association

Lipid droplets (LDs) are intracellular storage organelles composed of neutral lipids, such as triacylglycerol (TG), surrounded by a phospholipid (PL) monolayer decorated with specific proteins. Herein, we investigate the mechanism of protein association during LD and bilayer membrane expansion. We find that the neutral lipids play a dynamic role in LD expansion by further intercalating with the PL monolayer to create more surface-oriented TG molecules (SURF-TG). This interplay both reduces high surface tension incurred during LD budding or growth and also creates expansion-specific surface features for protein recognition. We then show that the autoinhibitory (AI) helix of CTP:phosphocholine cytidylyltransferase, a protein known to target expanding monolayers and bilayers, preferentially associates with large packing defects in a sequence-specific manner. Despite the presence of three phenylalanines, the initial binding with bilayers is predominantly mediated by the sole tryptophan due to its preference for membrane interfaces. Subsequent association is dependent on the availability of large, neighboring defects that can accommodate the phenylalanines, which are more probable in the stressed systems. Tryptophan, once fully associated, preferentially interacts with the glycerol moiety of SURF-TG in LDs. The calculation of AI binding free energy, hydrogen bonding and depth analysis, and in silico mutation experiments support the findings. Hence, SURF-TG can both reduce surface tension and mediate protein association, facilitating class II protein recruitment during LD expansion.

Investigating the structural properties of hydrophobic solvent-rich lipid bilayers

In vitro reconstitutions of lipid membranes have proven to be an indispensable tool to rationalize their molecular complexity and to understand their role in countless cellular processes. However, amongst the various techniques used to reconstitute lipid bilayers in vitro, several approaches are not solvent-free, but rather contain residual hydrophobic solvents in between the two bilayer leaflets, generally as a consequence of the procedure used to generate the bilayer. To what extent the presence of these hydrophobic solvents modifies bilayer properties with respect to native, solvent-free, conditions remains an open question that has important implications for the appropriate interpretation of numerous experimental observations. Here, we thorouhgly characterize hydrophobic solvent-rich lipid bilayers using atomistic molecular dynamics simulations. Our data indicate that while the presence of hydrophobic solvents at high concentrations, such as hexadecane, has a significant effect on membrane thickness, their effects on surface properties, membrane order and lateral stress are quite moderate. Our results corroborate the validity of in vitro approaches as model systems for the investigations of biological membranes but raise a few cautionary aspects that must be considered when investigating specific membrane properties.

Polyunsaturated Fatty Acid Modulates Membrane-Bound Monomeric α-Synuclein by Modulating Membrane Microenvironment through Preferential Interactions

There is ample evidence that both native functions and pathogenic aggregation of α-synuclein are intimately dependent on lipid interactions and fatty acid type; the regulatory mechanism however remains unclear. In the present work, using extensive atomistic molecular dynamics simulations and enhanced-sampling, we have focused on exploring the mechanism of fatty acid dependent regulation of monomeric α-Syn₁₀₀ in a native synaptic vesicle-like membrane. Our results show that α-Syn₁₀₀ spontaneously binds to the membrane through its N-terminal region (residues 1-34), where the depth of membrane insertion, the structure, and orientation of the membrane-bound α-Syn₁₀₀ and its impact on membrane structure are modulated by docosahexaenoic acid (DHA). DHA is a polyunsaturated fatty acid abundantly found in the brain and known to promote the oligomerization of α-synuclein. We found that DHA exhibits marked propensity to interact with monomeric α-Syn₁₀₀ and modulates the microenvironment of the protein by preferentially sorting DHA-containing phospholipids, depleting other phospholipids and cholesterol as well as increasing the proportion of anionic to neutral lipids in the immediate vicinity of the protein. Owing to the unique conformational flexibility, DHA chains form more lipid-packing defects in the membrane and efficiently coat the membrane-embedded surface of the protein, compared to the saturated and monounsaturated fatty acids. DHA thus makes the bilayer more amiable to protein adsorption and less prone to α-synuclein-induced perturbation associated with cytotoxicity. Indeed, in the absence of DHA, we observed significant thinning of the local bilayer membrane induced by α-Syn₁₀₀. Though α-Syn₁₀₀ is predominantly α-helical in membranes studied here, in the presence of DHA we observe formation of β-sheet/β-strands in the C-terminal region (residues 35-100) of α-Syn₁₀₀, which is extended out from the membrane surface. Notably, DHA induces β structure in the NAC domain of α-Syn₁₀₀ and promotes extended conformations as well as large solvent exposure of this hydrophobic domain, properties that are known to facilitate self-assembly of α-synuclein. To the best of our knowledge, this study for the first time provides the atomistic insights into DHA-induced regulatory mechanism of monomeric α-synuclein, having implications in protein structure and its physiological/pathological functions.

In Silico Prediction of the Binding, Folding, Insertion, and Overall Stability of Membrane-Active Peptides

Membrane-active peptides (MAPs) are short-length peptides used for potential biomedical applications in diagnostic imaging of tissues, targeted drug delivery, gene delivery, and antimicrobials and antibiotics. The broad appeal of MAPs is that they are infinitely variable, relatively low cost, and biocompatible. However, experimentally characterizing the specific properties of a MAP or its many variants is a low-resolution and potentially time-consuming endeavor; molecular dynamics (MD) simulations have emerged as an invaluable tool in identifying the biophysical interactions that are fundamental to the function of MAPs. In this chapter, a step-by-step approach to discreetly model the binding, folding, and insertion of a membrane-active peptide to a model lipid bilayer using MD simulations is described. Detailed discussion is devoted to the critical aspects of running these types of simulations: prior knowledge of the system, understanding the strengths and weaknesses of molecular mechanics force fields, proper construction and equilibration of the system, realistically estimating both experimental and computational timescales, and leveraging analysis to make direct comparisons to experimental results as often as possible.

Generation of antimicrobial peptides Leg1 and Leg2 from chickpea storage protein, active against food spoilage bacteria and foodborne pathogens

Contamination with bacteria leads to food waste and foodborne diseases with severe consequences for the environment and human health. Aiming to reduce food spoilage and infection, the present study developed novel highly active food-grade antimicrobial peptides affecting a wide range of bacteria. After extraction from chickpea, the storage protein legumin was hydrolyzed by the digestive protease chymotrypsin. Subsequent analysis by ultrahigh-performance micro-liquid chromatography-triple quadrupole time-of-flight tandem mass spectrometry determined the resulting peptide profiles. Virtual screening identified 21 potential antimicrobial peptides in the hydrolysates. Among those, the peptides Leg1 (RIKTVTSFDLPALRFLKL) and Leg2 (RIKTVTSFDLPALRWLKL) exhibited antimicrobial activity against 16 different bacteria, including pathogens, spoilage-causing bacteria and two antibiotic-resistant strains. Leg1/Leg2 showed minimum inhibitory concentrations (MIC) down to 15.6 µmol/L and were thus 10-1,000-fold more active compared to conventional food preservatives. Moreover, Leg1 and Leg2 showed bactericidal activity in contrast to the bacteriostatic activity of conventional preservatives.

Three-Dimensional Packing Defects in Lipid Membrane as a Function of Membrane Order

Lipid membrane packing defects are considered to be an essential parameter that regulates specific membrane binding of several peripheral proteins. In the absence of direct experimental characterization, lipid packing defects and their role in the binding of peripheral proteins are generally investigated through computational studies, which have been immensely successful in unraveling the key steps of the membrane-binding process. However, packing defects are calculated using two-dimensional (2D) projections and the crucial information on their depths is generally overlooked. Here, we present a simple yet computationally efficient algorithm, which identifies these defects in three dimensions. We validate the algorithm on a number of model membrane systems that are previously studied using 2D defect calculations and find that the defect size and the defect depth may not always be directly correlated. We employ the algorithm to understand the nature of packing defects in flat bilayer membranes exhibiting liquid-ordered (L_o), liquid-disordered (L_d), and co-existing (L_o/L_d) phases. Our results indicate the presence of shallower, smaller, and spatially localized defects in the L_o phase membranes as compared to the defects in L_d and mixed L_o/L_d phase membranes. Such analyses can elucidate the molecular-scale mechanisms that drive the preferential localization of certain proteins to either of the liquid phases or their interface. We also analyze the membrane sensing and anchoring process of a peptide in terms of the three-dimensional defects, which provides additional insights into the process with respect to depth distributions across the bilayer leaflets.

The Surface and Hydration Properties of Lipid Droplets

Lipid droplets (LDs) are energy storage organelles composed of neutral lipids, such as triacylglycerol (TG) and sterol esters, surrounded by a phospholipid (PL) monolayer. Their central role in metabolism, complex life cycle, and unique lipid monolayer surface have garnered great attention over the last decade. In this article, results from the largest and longest all-atom simulations to date suggest that 5-8% of the LD surface is occupied by TG molecules, a number that exceeds the maximal solubility reported for TGs in PL bilayers (2.8%). Two distinct classes of TG molecules that interact with the LD monolayer are found. Those at the monolayer surface (SURF-TG) are ordered like PLs with the glycerol moiety exposed to water, creating a significant amount of chemically unique packing defects, and the acyl chains extended toward the LD center. In contrast, the TGs that intercalate just into the PL tail region (CORE-TG) are disordered and increase the amount of PL packing defects and the PL tail order. The degree of interdigitation caused by CORE-TG is stable and determines the width of the TG-PL overlap, whereas that caused by SURF-TG fluctuates and is highly correlated with the area per PL or the expansion of the monolayer. Thus, when the supply of PLs is limited, SURF-TG may reduce surface tension by behaving as a secondary membrane component. The hydration properties of the simulated LD systems demonstrate ∼10 times more water in the LD core than previously reported. Collectively, the reported surface and hydration properties are expected to play a direct role in the mechanisms by which proteins target and interact with LDs.

Membrane Thinning Induces Sorting of Lipids and the Amphipathic Lipid Packing Sensor (ALPS) Protein Motif

Heterogeneities (e.g., membrane proteins and lipid domains) and deformations (e.g., highly curved membrane regions) in biological lipid membranes cause lipid packing defects that may trigger functional sorting of lipids and membrane-associated proteins. To study these phenomena in a controlled and efficient way within molecular simulations, we developed an external field protocol that artificially enhances packing defects in lipid membranes by enforcing local thinning of a flat membrane region. For varying lipid compositions, we observed strong thinning-induced depletion or enrichment, depending on the lipid's intrinsic shape and its effect on a membrane's elastic modulus. In particular, polyunsaturated and lysolipids are strongly attracted to regions high in packing defects, whereas phosphatidylethanolamine (PE) lipids and cholesterol are strongly repelled from it. Our results indicate that externally imposed changes in membrane thickness, area, and curvature are underpinned by shared membrane elastic principles. The observed sorting toward the thinner membrane region is in line with the sorting expected for a positively curved membrane region. Furthermore, we have demonstrated that the amphipathic lipid packing sensor (ALPS) protein motif, a known curvature and packing defect sensor, is effectively attracted to thinner membrane regions. By extracting the force that drives amphipathic molecules toward the thinner region, our thinning protocol can directly quantify and score the lipid packing sensing of different amphipathic molecules. In this way, our protocol paves the way toward high-throughput exploration of potential defect- and curvature-sensing motifs, making it a valuable addition to the molecular simulation toolbox.

Structure and Dynamics of the Acyl Chains in the Membrane Trafficking and Enzymatic Processing of Lipids

The regulatory chemical mechanisms of lipid trafficking and degradation are involved in many pathophysiological processes, being implicated in severe pain, inflammation, and cancer. In addition, the processing of lipids is also relevant for industrial and environmental applications. However, there is poor understanding of the chemical features that control lipid membrane trafficking and allow lipid-degrading enzymes to efficiently select and hydrolyze specific fatty acids from a complex cellular milieu of bioactive lipids. This is particularly true for lipid acyl chains, which have diverse structures that can critically affect the many complex reactions needed to elongate, desaturate, or transport fatty acids. Building upon our own contributions in this field, we will discuss how molecular simulations, integrated with experimental evidence, have revealed that the structure and dynamics of the lipid tail are actively involved in modulating membrane trafficking at cellular organelles, and enzymatic reactions at cell membranes. Further evidence comes from recent crystal structures of lipid receptors and remodeling enzymes. Taken together, these recent works have identified those structural features of the lipid acyl chain that are crucial for the regioselectivity and stereospecificity of essential desaturation reactions. In this context, we will first illustrate how atomistic and coarse-grained simulations have elucidated the structure-function relationships between the chemical composition of the lipid's acyl chains and the molecular properties of lipid bilayers. Particular emphasis will be given to the prominent chemical role of the number of double carbon-carbon bonds along the lipid acyl chain, that is, discriminating between saturated, monounsaturated, and polyunsaturated lipids. Different levels of saturation in fatty acid molecules dramatically influence the biophysical properties of lipid assemblies and their interaction with proteins. We will then discuss the processing of lipids by membrane-bound enzymes. Our focus will be on lipids such as anandamide and 2-arachidonoylglycerol. These are the main molecules that act as neurotransmitters in the endocannabinoid system. Specifically, recent findings indicate a crucial interplay between the level of saturation of the lipid tail, its energetically and sterically favored conformations, and the hydrophobic accessory cavities in lipid-degrading enzymes, which help form catalytically active conformations of the selected substrate. This Account will emphasize how the specific chemical structure of acyl chains affects the molecular mechanisms for modulating membrane trafficking and selective hydrolysis. The results examined here show that, by using molecular simulations to investigate lipid plasticity and substrate flexibility, researchers can enrich their interpretation of experimental results about the structure-function relationships of lipids. This could positively impact chemical and biological studies in the field and ultimately support protein engineering studies and structure-based drug discovery to target lipid-processing enzymes.

Lipids or Proteins: Who Is Leading the Dance at Membrane Contact Sites?

Understanding the mode of action of membrane contact sites (MCSs) across eukaryotic organisms at the near-atomic level to infer function at the cellular and tissue levels is a challenge scientists are currently facing. These peculiar systems dedicated to inter-organellar communication are perfect examples of cellular processes where the interplay between lipids and proteins is critical. In this mini review, we underline the link between membrane lipid environment, the recruitment of proteins at specialized membrane domains and the function of MCSs. More precisely, we want to give insights on the crucial role of lipids in defining the specificity of plant endoplasmic reticulum (ER)-plasma membrane (PM) MCSs and we further propose approaches to study them at multiple scales. Our goal is not so much to go into detailed description of MCSs, as there are numerous focused reviews on the subject, but rather try to pinpoint the critical elements defining those structures and give an original point of view by considering the subject from a near-atomic angle with a focus on lipids. We review current knowledge as to how lipids can define MCS territories, play a role in the recruitment and function of the MCS-associated proteins and in turn, how the lipid environment can be modified by proteins.

Formation of α-helical and β-sheet structures in membrane-bound human IAPP monomer and the resulting membrane deformation

Upon adsorption on membrane, human IAPP monomer takes conformational changes from coils to α-helices and β-sheets. The helices inserted and β on surface cause different types of membrane deformation, implying two distinct aggregation mechanisms.

GROmaρs: A GROMACS-Based Toolset to Analyze Density Maps Derived from Molecular Dynamics Simulations

We introduce a computational toolset, named GROmaρs, to obtain and compare time-averaged density maps from molecular dynamics simulations. GROmaρs efficiently computes density maps by fast multi-Gaussian spreading of atomic densities onto a three-dimensional grid. It complements existing map-based tools by enabling spatial inspection of atomic average localization during the simulations. Most importantly, it allows the comparison between computed and reference maps (e.g., experimental) through calculation of difference maps and local and time-resolved global correlation. These comparison operations proved useful to quantitatively contrast perturbed and control simulation data sets and to examine how much biomolecular systems resemble both synthetic and experimental density maps. This was especially advantageous for multimolecule systems in which standard comparisons like RMSDs are difficult to compute. In addition, GROmaρs incorporates absolute and relative spatial free-energy estimates to provide an energetic picture of atomistic localization. This is an open-source GROMACS-based toolset, thus allowing for static or dynamic selection of atoms or even coarse-grained beads for the density calculation. Furthermore, masking of regions was implemented to speed up calculations and to facilitate the comparison with experimental maps. Beyond map comparison, GROmaρs provides a straightforward method to detect solvent cavities and average charge distribution in biomolecular systems. We employed all these functionalities to inspect the localization of lipid and water molecules in aquaporin systems, the binding of cholesterol to the G protein coupled chemokine receptor type 4, and the identification of permeation pathways through the dermicidin antimicrobial channel. Based on these examples, we anticipate a high applicability of GROmaρs for the analysis of molecular dynamics simulations and their comparison with experimentally determined densities.