Simulating realistic membrane shapes

Modeling membranes in situ

Molecular dynamics simulations of cellular membranes have come a long way-from simple model lipid bilayers to multicomponent systems capturing the crowded and complex nature of real cell membranes. In this opinionated minireview, we discuss the current challenge to simulate the dynamics of membranes in their native environment, in situ, with the prospect of reaching the level of whole cells and cell organelles using an integrative modeling framework.

Efficient calculation of orientation-dependent lipid dynamics from membrane simulations

Molecular dynamics simulations of lipid membranes have become increasingly impactful in biophysics because they offer atomistic resolution of structural fluctuations in relation to their functional outputs. Yet quantitative characterization of multiscale processes is a formidable challenge due to the distribution of motions that evade analysis of discrete simulation data. Here we investigate the efficient calculation of CH bond relaxation rates from membrane simulations. Widely used computational approaches offer numerical simplicity but fall short of capturing crucial aspects of the orientation dependence of the dynamics. To circumvent this problem, we introduced a robust framework based on liquid crystal theory which considers explicitly the CH bond motions with respect to the director axis (bilayer normal). Analysis of the orientation dependence of the dynamics shows excellent agreement with experiment, illustrating how the ordering potential affects the calculated relaxation rates. Furthermore, a fit-based resampling of the autocorrelation function of the bond fluctuations validates the new approach for low-temporal resolution data. The recovered relaxation rates indicate that at short timescales, both with and without cholesterol, the local motions of CH bonds describe the bilayer microviscosity and resemble liquid hydrocarbons. Our results establish the critical role of the orientational anisotropy in analysis of membrane simulations, explain fundamental aspects of lipid dynamics, and provide guidelines for extracting information that can be compared to experimental data.

Diffusion Analyses along Mean and Gaussian-Curved Membranes with CurD

Curved cellular membranes are both abundant and functionally relevant. While novel tomography approaches reveal the structural details of curved membranes, their dynamics pose an experimental challenge. Curvature especially affects the diffusion of lipids and macromolecules, yet neither experiments nor continuum models distinguish geometric effects from those caused by curvature-induced changes in membrane properties. Molecular simulations could excel here, yet despite community interest toward curved membranes, tools for their analysis are still lacking. Here, we satisfy this demand by introducing CurD, our novel and openly available implementation of the Vertex-oriented Triangle Propagation algorithm to the study of lipid diffusion along membranes with mean and/or Gaussian curvature. This approach, aided by our highly optimized implementation, computes geodetic distances significantly faster than conventional implementations of path-finding algorithms. Our tool, applied to coarse-grained simulations, allows for the first time the analysis of curvature effects on diffusion at size scales relevant to physiological processes such as endocytosis. Our analyses with different membrane geometries reveal that Gaussian curvature plays a surprisingly small role on lipid motion, whereas mean curvature; i.e., the packing of lipid headgroups largely dictates their mobility.

The intricate link between membrane lipid structure and composition and membrane structural properties in bacterial membranes

It is now evident that the cell manipulates lipid composition to regulate different processes such as membrane protein insertion, assembly and function. Moreover, changes in membrane structure and properties, lipid homeostasis during growth and differentiation with associated changes in cell size and shape, and responses to external stress have been related to drug resistance across mammalian species and a range of microorganisms. While it is well known that the biomembrane is a fluid self-assembled nanostructure, the link between the lipid components and the structural properties of the lipid bilayer are not well understood. This perspective aims to address this topic with a view to a more detailed understanding of the factors that regulate bilayer structure and flexibility. We describe a selection of recent studies that address the dynamic nature of bacterial lipid diversity and membrane properties in response to stress conditions. This emerging area has important implications for a broad range of cellular processes and may open new avenues of drug design for selective cell targeting.

Mesoscale simulation of biomembranes with FreeDTS

We present FreeDTS software for performing computational research on biomembranes at the mesoscale. In this software, a membrane is represented by a dynamically triangulated surface equipped with vertex-based inclusions to integrate the effects of integral and peripheral membrane proteins. Several algorithms are included in the software to simulate complex membranes at different conditions such as framed membranes with constant tension, vesicles and high-genus membranes with various fixed volumes or constant pressure differences and applying external forces to membrane regions. Furthermore, the software allows the user to turn off the shape evolution of the membrane and focus solely on the organization of proteins. As a result, we can take realistic membrane shapes obtained from, for example, cryo-electron tomography and backmap them into a finer simulation model. In addition to many biomembrane applications, this software brings us a step closer to simulating realistic biomembranes with molecular resolution. Here we provide several interesting showcases of the power of the software but leave a wide range of potential applications for interested users.

Shocker─A Molecular Dynamics Protocol and Tool for Accelerating and Analyzing the Effects of Osmotic Shocks

The process of osmosis, a fundamental phenomenon in life, drives water through a semipermeable membrane in response to a solute concentration gradient across this membrane. In vitro, osmotic shocks are often used to drive shape changes in lipid vesicles, for instance, to study fission events in the context of artificial cells. While experimental techniques provide a macroscopic picture of large-scale membrane remodeling processes, molecular dynamics (MD) simulations are a powerful tool to study membrane deformations at the molecular level. However, simulating an osmotic shock is a time-consuming process due to slow water diffusion across the membrane, making it practically impossible to examine its effects in classic MD simulations. In this article, we present Shocker, a Python-based MD tool for simulating the effects of an osmotic shock by selecting and relocating water particles across a membrane over the course of several pumping cycles. Although this method is primarily aimed at efficiently simulating volume changes in vesicles, it can also handle membrane tubes and double bilayer systems. Additionally, Shocker is force field-independent and compatible with both coarse-grained and all-atom systems. We demonstrate that our tool is applicable to simulate both hypertonic and hypotonic osmotic shocks for a range of vesicular and bilamellar setups, including complex multicomponent systems containing membrane proteins or crowded internal solutions.

Curvature sensing lipid dynamics in a mitochondrial inner membrane model

Membrane curvature is essential for many cellular structures and processes, and factors such as leaflet asymmetry, lipid composition, and proteins all play important roles. Cardiolipin is the signature lipid of mitochondrial membranes and is essential for maintaining the highly curved shapes of the inner mitochondrial membrane (IMM) and the spatial arrangement of membrane proteins. In this study, we investigate the partitioning behavior of various lipids present in the IMM using coarse-grained molecular dynamics simulations. This study explores curved bilayer systems containing phosphatidylcholine (PC), phosphatidylethanolamine (PE), and cardiolipin (CDL) in binary and ternary component mixtures. Curvature properties such as mean and Gaussian curvatures, as well as the distribution of lipids into the various curved regions of the cristae models, are quantified. Overall, this work represents an advance beyond previous studies on lipid curvature sensing by simulating these systems in a geometry that has the morphological features and scales of curvature consistent with regions of the IMM. We find that CDL has a stronger preference for accumulating in regions of negative curvature than PE lipids, in agreement with previous results. Furthermore, we find lipid partitioning propensity is dominated by sensitivity to mean curvature, while there is a weaker correlation with Gaussian curvature.

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.

Building Asymmetric Lipid Bilayers for Molecular Dynamics Simulations: What Methods Exist and How to Choose One?

The compositional asymmetry of biological membranes has attracted significant attention over the last decade. Harboring more differences from symmetric membranes than previously appreciated, asymmetric bilayers have proven quite challenging to study with familiar concepts and techniques, leaving many unanswered questions about the reach of the asymmetry effects. One particular area of active research is the computational investigation of composition- and number-asymmetric lipid bilayers with molecular dynamics (MD) simulations. Offering a high level of detail into the organization and properties of the simulated systems, MD has emerged as an indispensable tool in the study of membrane asymmetry. However, the realization that results depend heavily on the protocol used for constructing the asymmetric bilayer models has sparked an ongoing debate about how to choose the most appropriate approach. Here we discuss the underlying source of the discrepant results and review the existing methods for creating asymmetric bilayers for MD simulations. Considering the available data, we argue that each method is well suited for specific applications and hence there is no single best approach. Instead, the choice of a construction protocol-and consequently, its perceived accuracy-must be based primarily on the scientific question that the simulations are designed to address.

Mesoscale simulations: An indispensable approach to understand biomembranes

Computer simulation techniques form a versatile tool, a computational microscope, for exploring biological processes. This tool has been particularly effective in exploring different features of biological membranes. In recent years, thanks to elegant multiscale simulation schemes, some fundamental limitations of investigations by distinct simulation techniques have been resolved. As a result, we are now capable of exploring processes spanning multiple scales beyond the capacity of any single technique. In this perspective, we argue that mesoscale simulations require more attention and must be further developed to fill evident gaps in a quest toward simulating and modeling living cell membranes.

Close, but not too close: a mesoscopic description of (a)symmetry and membrane shaping mechanisms

Biomembranes are fundamental to our understanding of the cell, the basic building block of all life. An intriguing aspect of membranes is their ability to assume a variety of shapes, which is crucial for cell function. Here, we review various membrane shaping mechanisms with special focus on the current understanding of how local curvature and local rigidity induced by membrane proteins leads to emerging forces and consequently large-scale membrane deformations. We also argue that describing the interaction of rigid proteins with membranes purely in terms of local membrane curvature is incomplete and that changes in the membrane rigidity moduli must also be considered.

A brief history of visualizing membrane systems in molecular dynamics simulations

Understanding lipid dynamics and function, from the level of single, isolated molecules to large assemblies, is more than ever an intensive area of research. The interactions of lipids with other molecules, particularly membrane proteins, are now extensively studied. With advances in the development of force fields for molecular dynamics simulations (MD) and increases in computational resources, the creation of realistic and complex membrane systems is now common. In this perspective, we will review four decades of the history of molecular dynamics simulations applied to membranes and lipids through the prism of molecular graphics.

Kinetics of antigen cross-presentation assessed experimentally and by a model of the complete endomembrane system

Dendritic cells (DCs) have a specialized endomembrane system capable of presenting exogenous antigens in the context of MHC class I (MHC-I) molecules. This process, named cross-presentation, is crucial to activate CD8+ T lymphocytes and initiate cytotoxic immune responses. In this report, we present an Agent-Based Model in combination with Ordinary Differential Equations with enough complexity to reproduce cross-presentation. The model embraces the secretory and endocytic pathways, in connection with the plasma membrane, the endoplasmic reticulum, and the cytosol. Key molecules required for cross-presentation were included as cargoes. In the simulations, the kinetics of MHC-I uptake and recycling, and cross-presentation accurately reproduced experimental values. The model proved to be a suitable tool to elaborate hypotheses and design experiments. In particular, the model predictions and the experimental results obtained indicate that the rate-limiting step in cross-presentation of soluble ovalbumin is MHC-I loading after proteasomal processing of the antigenic protein.

A Review of Mechanics-Based Mesoscopic Membrane Remodeling Methods: Capturing Both the Physics and the Chemical Diversity

Specialized classes of proteins, working together in a tightly orchestrated manner, induce and maintain highly curved cellular and organelles membrane morphology. Due to the various experimental constraints, including the resolution limits of imaging techniques, it is non-trivial to accurately elucidate interactions among the various components involved in membrane deformation. The spatial and temporal scales of the systems also make it formidable to investigate them using simulations with molecular details. Interestingly, mechanics-based mesoscopic models have been used with great success in recapitulating the membrane deformations observed in experiments. In this review, we collate together and discuss the various mechanics-based mesoscopic models for protein-mediated membrane deformation studies. In particular, we provide an elaborate description of a mesoscopic model where the membrane is modeled as a triangulated sheet and proteins are represented as either nematics or filaments. This representation allows us to explore the various aspects of protein-protein and protein-membrane interactions as well as examine the underlying mechanistic pathways for emergent behavior such as curvature-mediated protein localization and membrane deformation. We also put forward current efforts in the field towards back-mapping these mesoscopic models to finer-grained particle-based models-a framework that could be used to explore how molecular interactions propagate to physical scales and vice-versa. We end the review with an integrative-modeling-based road map where experimental imaging micrograph and biochemical data are combined with mesoscopic and molecular simulations methods in a theoretically consistent manner to faithfully recapitulate the multiple length and time scales in the membrane remodeling processes.

Synthetic Membrane Shaper for Controlled Liposome Deformation

Shape defines the structure and function of cellular membranes. In cell division, the cell membrane deforms into a "dumbbell" shape, while organelles such as the autophagosome exhibit "stomatocyte" shapes. Bottom-up in vitro reconstitution of protein machineries that stabilize or resolve the membrane necks in such deformed liposome structures is of considerable interest to characterize their function. Here we develop a DNA-nanotechnology-based approach that we call the synthetic membrane shaper (SMS), where cholesterol-linked DNA structures attach to the liposome membrane to reproducibly generate high yields of stomatocytes and dumbbells. In silico simulations confirm the shape-stabilizing role of the SMS. We show that the SMS is fully compatible with protein reconstitution by assembling bacterial divisome proteins (DynaminA, FtsZ:ZipA) at the catenoidal neck of these membrane structures. The SMS approach provides a general tool for studying protein binding to complex membrane geometries that will greatly benefit synthetic cell research.

Activation energy and force fields during topological transitions of fluid lipid vesicles

Topological transitions of fluid lipid membranes are fundamental processes for cell life. For example, they are required for endo- and exocytosis or to enable neurotransmitters to cross the neural synapses. Here, inspired by the idea that fusion and fission proteins could have evolved in Nature in order to carry out a minimal work expenditure, we evaluate the minimal free energy pathway for the transition between two spherical large unilamellar vesicles and a dumbbell-shaped one. To address the problem, we propose and successfully use a Ginzburg-Landau type of free energy, which allows us to uniquely describe without interruption the whole, full-scale topological change. We also compute the force fields needed to overcome the involved energy barriers. The obtained forces are in excellent agreement, in terms of intensity, scale, and spatial localization with experimental data on typical fission protein systems, whereas they suggest the presence of additional features in fusion proteins.

Molecular Dynamics Simulations of Curved Lipid Membranes

Eukaryotic cells contain membranes with various curvatures, from the near-plane plasma membrane to the highly curved membranes of organelles, vesicles, and membrane protrusions. These curvatures are generated and sustained by curvature-inducing proteins, peptides, and lipids, and describing these mechanisms is an important scientific challenge. In addition to that, some molecules can sense membrane curvature and thereby be trafficked to specific locations. The description of curvature sensing is another fundamental challenge. Curved lipid membranes and their interplay with membrane-associated proteins can be investigated with molecular dynamics (MD) simulations. Various methods for simulating curved membranes with MD are discussed here, including tools for setting up simulation of vesicles and methods for sustaining membrane curvature. The latter are divided into methods that exploit scaffolding virtual beads, methods that use curvature-inducing molecules, and methods applying virtual forces. The variety of simulation tools allow researcher to closely match the conditions of experimental studies of membrane curvatures.

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

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

Quantitative Analysis of the Interactions of Metal Complexes and Amphiphilic Systems: Calorimetric, Spectroscopic and Theoretical Aspects

Metals and metal-based compounds have many implications in biological systems. They are involved in cellular functions, employed in the formation of metal-based drugs and present as pollutants in aqueous systems, with toxic effects for living organisms. Amphiphilic molecules also play important roles in the above bio-related fields as models of membranes, nanocarriers for drug delivery and bioremediating agents. Despite the interest in complex systems involving both metal species and surfactant aggregates, there is still insufficient knowledge regarding the quantitative aspects at the basis of their binding interactions, which are crucial for extensive comprehension of their behavior in solution. Only a few papers have reported quantitative analyses of the thermodynamic, kinetic, speciation and binding features of metal-based compounds and amphiphilic aggregates, and no literature review has yet addressed the quantitative study of these complexes. Here, we summarize and critically discuss the recent contributions to the quantitative investigation of the interactions of metal-based systems with assemblies made of amphiphilic molecules by calorimetric, spectrophotometric and computational techniques, emphasizing the unique picture and parameters that such an analytical approach may provide, to support a deep understanding and beneficial use of these systems for several applications.

Identification of electroporation sites in the complex lipid organization of the plasma membrane

The plasma membrane of a biological cell is a complex assembly of lipids and membrane proteins, which tightly regulate transmembrane transport. When a cell is exposed to strong electric field, the membrane integrity becomes transiently disrupted by formation of transmembrane pores. This phenomenon termed electroporation is already utilized in many rapidly developing applications in medicine including gene therapy, cancer treatment, and treatment of cardiac arrhythmias. However, the molecular mechanisms of electroporation are not yet sufficiently well understood; in particular, it is unclear where exactly pores form in the complex organization of the plasma membrane. In this study, we combine coarse-grained molecular dynamics simulations, machine learning methods, and Bayesian survival analysis to identify how formation of pores depends on the local lipid organization. We show that pores do not form homogeneously across the membrane, but colocalize with domains that have specific features, the most important being high density of polyunsaturated lipids. We further show that knowing the lipid organization is sufficient to reliably predict poration sites with machine learning. Additionally, by analysing poration kinetics with Bayesian survival analysis we show that poration does not depend solely on local lipid arrangement, but also on membrane mechanical properties and the polarity of the electric field. Finally, we discuss how the combination of atomistic and coarse-grained molecular dynamics simulations, machine learning methods, and Bayesian survival analysis can guide the design of future experiments and help us to develop an accurate description of plasma membrane electroporation on the whole-cell level. Achieving this will allow us to shift the optimization of electroporation applications from blind trial-and-error approaches to mechanistic-driven design.

Coupling Bulk Phase Separation of Disordered Proteins to Membrane Domain Formation in Molecular Simulations on a Bespoke Compute Fabric

Phospholipid membranes surround the cell and its internal organelles, and their multicomponent nature allows the formation of domains that are important in cellular signalling, the immune system, and bacterial infection. Cytoplasmic compartments are also created by the phase separation of intrinsically disordered proteins into biomolecular condensates. The ubiquity of lipid membranes and protein condensates raises the question of how three-dimensional droplets might interact with two-dimensional domains, and whether this coupling has physiological or pathological importance. Here, we explore the equilibrium morphologies of a dilute phase of a model disordered protein interacting with an ideal-mixing, two-component lipid membrane using coarse-grained molecular simulations. We find that the proteins can wet the membrane with and without domain formation, and form phase separated droplets bound to membrane domains. Results from much larger simulations performed on a novel non-von-Neumann compute architecture called POETS, which greatly accelerates their execution compared to conventional hardware, confirm the observations. Reducing the wall clock time for such simulations requires new architectures and computational techniques. We demonstrate here an inter-disciplinary approach that uses real-world biophysical questions to drive the development of new computing hardware and simulation algorithms.

Multiscale simulations of protein and membrane systems

Classical multiscale simulations are perfectly suited to investigate biological soft matter systems. Owing to the bridging between microscopically realistic and lower-resolution models or the integration of a hierarchy of subsystems, one gets access to biologically relevant system sizes and timescales. In recent years, increasingly complex systems and processes have come into focus such as multidomain proteins, phase separation processes in biopolymer solutions, multicomponent biomembranes, or multiprotein complexes up to entire viruses. The review shows factors that have contributed to this progress - from improved models to machine-learning-based analysis and scale-bridging methods.

Interactions between polymyxin B and various bacterial membrane mimics: A molecular dynamics study

Polymyxin B (PMB) is clinically used as a last-line therapy against life-threatening Gram-negative "superbugs". However, thorough understanding of the membrane actions of PMB at a molecular level is still lacking. In this work, a variety of bacterial membrane mimics with varying lipid compositions were built, and their interactions with PMB were systematically investigated using coarse-grained molecular dynamics simulation. PMB demonstrated characteristic preference to specific lipid species during its interaction with different membrane systems, such as the rough mutant lipipolysacchrides (Re LPS) preference in an outer membrane (OM) or the cardiolipin and POPG affinity in an inner membrane (IM). As a result of the lipid-specific actions, complicated membrane interaction states of PMB were observed, including adsorption on the OM surface. In contrast, for the IM or a mutative OM containing "impurity lipids" like POPE, POPG or lipid A, it could insert into the membrane via its acyl chain. Such actions of PMB influence the structure and lipid mobility of the membrane. In particular, the OM-bound PMB breaks the synchronous movement of Re LPS molecules in the outer leaflet and makes them diffuse more randomly, while its insertion into IM blocks the phospholipid diffusion and makes the membrane more homogeneous in the trajectory space. Our results provide insight into the action mechanism of PMB at a membrane level and a foundation for developing novel and safer polymyxin strategies for better clinical use.

Fabrication of Adsorbed Fe(III) and Structurally Doped Fe(III) in Montmorillonite/TiO2 Composite for Photocatalytic Degradation of Phenol

The Fe(III)-doped montmorillonite (Mt)/TiO2 composites were fabricated by adding Fe(III) during or after the aging of TiO2/Ti(OH)4 sol–gel in Mt, named as xFe-Mt/(1 − x)Fe-TiO2 and Fe/Mt/TiO2, respectively. In the xFe-Mt/(1 − x)Fe-TiO2, Fe(III) cations were expected to be located in the structure of TiO2, in the Mt, and in the interface between them, while Fe(III) ions are physically adsorbed on the surfaces of the composites in the Fe/Mt/TiO2. The narrower energy bandgap (Eg) lower photo-luminescence intensity were observed for the composites compared with TiO2. Better photocatalytic performance for phenol degradation was observed in the Fe/Mt/TiO2. The 94.6% phenol degradation was due to greater charge generation and migration capacity, which was confirmed by photocurrent measurements and electrochemical impedance spectroscopy (EIS). The results of the energy-resolved distribution of electron traps (ERDT) suggested that the Fe/Mt/TiO2 possessed a larger amorphous rutile phase content in direct contact with crystal anatase than that of the xFe-Mt/(1 − x)Fe-TiO2. This component is the fraction that is mainly responsible for the photocatalytic phenol degradation by the composites. As for the xFe-Mt/(1 − x)Fe-TiO2, the active rutile phase was followed by isolated amorphous phases which had larger (Eg) and which did not act as a photocatalyst. Thus, the physically adsorbed Fe(III) enhanced light adsorption and avoided charge recombination, resulting in improved photocatalytic performance. The mechanism of the photocatalytic reaction with the Fe(III)-doped Mt/TiO2 composite was proposed.

Langmuir Monolayer Techniques for the Investigation of Model Bacterial Membranes and Antibiotic Biodegradation Mechanisms

The amounts of antibiotics of anthropogenic origin released and accumulated in the environment are known to have a negative impact on local communities of microorganisms, which leads to disturbances in the course of the biodegradation process and to growing antimicrobial resistance. This mini-review covers up-to-date information regarding problems related to the omnipresence of antibiotics and their consequences for the world of bacteria. In order to understand the interaction of antibiotics with bacterial membranes, it is necessary to explain their interaction mechanism at the molecular level. Such molecular-level interactions can be probed with Langmuir monolayers representing the cell membrane. This mini-review describes monolayer experiments undertaken to investigate the impact of selected antibiotics on components of biomembranes, with particular emphasis on the role and content of individual phospholipids and lipopolysaccharides (LPS). It is shown that the Langmuir technique may provide information about the interactions between antibiotics and lipids at the mixed film surface (π–A isotherm) and about the penetration of the active substances into the phospholipid monolayer model membranes (relaxation of the monolayer). Effects induced by antibiotics on the bacterial membrane may be correlated with their bactericidal activity, which may be vital for the selection of appropriate bacterial consortia that would ensure a high degradation efficiency of pharmaceuticals in the environment.

Model Membrane Systems Used to Study Plasma Membrane Lipid Asymmetry

It is well known that the lipid distribution in the bilayer leaflets of mammalian plasma membranes (PMs) is not symmetric. Despite this, model membrane studies have largely relied on chemically symmetric model membranes for the study of lipid–lipid and lipid–protein interactions. This is primarily due to the difficulty in preparing stable, asymmetric model membranes that are amenable to biophysical studies. However, in the last 20 years, efforts have been made in producing more biologically faithful model membranes. Here, we review several recently developed experimental and computational techniques for the robust generation of asymmetric model membranes and highlight a new and particularly promising technique to study membrane asymmetry.

Computational Approaches to Explore Bacterial Toxin Entry into the Host Cell

Many bacteria secrete toxic protein complexes that modify and disrupt essential processes in the infected cell that can lead to cell death. To conduct their action, these toxins often need to cross the cell membrane and reach a specific substrate inside the cell. The investigation of these protein complexes is essential not only for understanding their biological functions but also for the rational design of targeted drug delivery vehicles that must navigate across the cell membrane to deliver their therapeutic payload. Despite the immense advances in experimental techniques, the investigations of the toxin entry mechanism have remained challenging. Computer simulations are robust complementary tools that allow for the exploration of biological processes in exceptional detail. In this review, we first highlight the strength of computational methods, with a special focus on all-atom molecular dynamics, coarse-grained, and mesoscopic models, for exploring different stages of the toxin protein entry mechanism. We then summarize recent developments that are significantly advancing our understanding, notably of the glycolipid–lectin (GL-Lect) endocytosis of bacterial Shiga and cholera toxins. The methods discussed here are also applicable to the design of membrane-penetrating nanoparticles and the study of the phenomenon of protein phase separation at the surface of the membrane. Finally, we discuss other likely routes for future development.