Experimental validation of a phase-field model to predict coarsening dynamics of lipid domains in multicomponent membranes

Making the cut: Multiscale simulation of membrane remodeling

Biological membranes are dynamic heterogeneous materials, and their shape and organization are tightly coupled to the properties of the proteins in and around them. However, the length scales of lipid and protein dynamics are far below the size of membrane-bound organelles, much less an entire cell. Therefore, multiscale modeling approaches are often necessary to build a comprehensive picture of the interplay of these factors, and have provided critical insights into our understanding of membrane dynamics. Here, we review computational methods for studying membrane remodeling, as well as passive and active examples of protein-driven membrane remodeling. As the field advances towards the modeling of key aspects of organelles and whole cells - an increasingly accessible regime of study - we summarize here recent successes and offer comments on future trends.

Modelling lipid rafts formation through chemo-mechanical interplay triggered by receptor-ligand binding

Cell membranes, mediator of many biological mechanisms from adhesion and metabolism up to mutation and infection, are highly dynamic and heterogeneous environments exhibiting a strong coupling between biochemical events and structural re-organisation. This involves conformational changes induced, at lower scales, by lipid order transitions and by the micro-mechanical interplay of lipids with transmembrane proteins and molecular diffusion. Particular attention is focused on lipid rafts, ordered lipid microdomains rich of signalling proteins, that co-localise to enhance substance trafficking and activate different intracellular biochemical pathways. In this framework, the theoretical modelling of the dynamic clustering of lipid rafts implies a full multiphysics coupling between the kinetics of phase changes and the mechanical work performed by transmembrane proteins on lipids, involving the bilayer elasticity. This mechanism produces complex interspecific dynamics in which membrane stresses and chemical potentials do compete by determining different morphological arrangements, alteration in diffusive walkways and coalescence phenomena, with a consequent influence on both signalling potential and intracellular processes. Therefore, after identifying the leading chemo-mechanical interactions, the present work investigates from a modelling perspective the spatio-temporal evolution of raft domains to theoretically explain co-localisation and synergy between proteins' activation and raft formation, by coupling diffusive and mechanical phenomena to observe different morphological patterns and clustering of ordered lipids. This could help to gain new insights into the remodelling of cell membranes and could potentially suggest mechanically based strategies to control their selectivity, by orienting intracellular functions and mechanotransduction.

On Fusogenicity of Positively Charged Phased-Separated Lipid Vesicles: Experiments and Computational Simulations

This paper studies the fusogenicity of cationic liposomes in relation to their surface distribution of cationic lipids and utilizes membrane phase separation to control this surface distribution. It is found that concentrating the cationic lipids into small surface patches on liposomes, through phase-separation, can enhance liposome's fusogenicity. Further concentrating these lipids into smaller patches on the surface of liposomes led to an increased level of fusogenicity. These experimental findings are supported by numerical simulations using a mathematical model for phase-separated charged liposomes. Findings of this study may be used for design and development of highly fusogenic liposomes with minimal level of toxicity.

Charged Lipids Modulate the Phase Separation in Multicomponent Membranes

Phase separation in lipid membranes controls the organization of membrane components and thus regulates membrane-mediated processes. Membrane phase behavior is influenced by the molecular properties of its components and their relative concentrations. Charged lipid species are among the most essential components of lipid membranes, and their impact on the membrane phase behavior is yet to be fully understood. Aiming to provide insight into this impact, this paper investigates how the presence and amount of anionic and cationic lipids affect the phase behavior of multicomponent membranes. Membranes of ternary composition DOPC:DPPC:Chol with two distinct molar ratios were used to test the hypothesis that inclusion of charged lipids with saturated tails, beyond a certain concentration, would impede phase separation in an otherwise phase-separating membrane. Fluorescence microscopy examination of electroformed giant liposomes revealed that when more than half of DOPC in the examined mixtures was replaced with DOPA or DOTAP, phase separation in liposomes was somewhat suppressed, and this effect increased with increasing charged lipid content. This effect depended on the membrane surface charge density as the half-maximal effect was observed at around 0.0072 C Å-2 in all examined cases. The phase-separation suppressing effect of DOPA was neutralized when oppositely charged lipid DOTAP was included in the mixture. Likewise, presence of divalent cation Ca2+ in the solution neutralized the impact of negatively charged DOPA. These results underline the detrimental influence of surface charge density on membrane phase behavior. More importantly, these findings suggest that the charged lipid content in membranes may be a regulator of their phase behavior and open new opportunities for the design of synthetic lipid membranes.

Systematic analysis of curvature-dependent lipid dynamics in a stochastic 3D membrane model

To minimize the free energy of the system, lipid membranes display curvature-dependent rearrangements at the local and global scale. The optimal membrane shape is generally approximated by averaging the curvature preference of individual lipids across the whole surface. Potential stress due to imperfections in lipid packing caused by local lipid inhomogeneities, however, is frequently neglected. Here, we developed a stochastic 3D membrane model to investigate the relevance of this parameter for shape-dependent lipid and membrane dynamics. A systematic analysis of the discretized Helfrich type Hamiltonian indicates that stress-energy arising from imperfections in packing is analogous to van der Waals interactions, jointly determining membrane shape and localization of curvature-sensitive lipids based on their relative strengths. Insights from this work can be used to characterize natural and design synthetic agents for membrane-shape changes.

Domain dynamics of phase-separated lipid membranes under shear flow

The dynamical behaviour of lateral domains on phase-separated lipid vesicles under external flow is reported. A microfluidic chamber was used for the immobilization of vesicles and the application of shear. Microscopic observation revealed that domains tended to be localized at the vortex center and to exhibit a stripe morphology as the flow speed increased. We clarified the dependency of domain behaviors on the flow speed and lipid mixing fraction. The cholesterol ratio in the membrane affected these domain behaviors. Next, we investigated the growth of domains under flow. We discuss the mechanism of these trends by considering the free energy of phase separation, and reproduce the experimental results by numerical simulations. These findings may lead to a better understanding of the dynamical properties of the membrane under nonequilibrium situations and the biophysical mechanism of cellular mechanotransduction.

Lipid domain coarsening and fluidity in multicomponent lipid vesicles: A continuum based model and its experimental validation

Liposomes that achieve a heterogeneous and spatially organized surface through phase separation have been recognized to be a promising platform for delivery purposes. However, their design and optimization through experimentation can be expensive and time-consuming. To assist with the design and reduce the associated cost, we propose a computational platform for modeling membrane coarsening dynamics based on the principles of continuum mechanics and thermodynamics. This model couples phase separation to lateral flow and accounts for different membrane fluidity within the different phases, which is known to affect the coarsening dynamics on lipid membranes. The simulation results are in agreement with the experimental data in terms of liquid ordered domains area fraction, total domains perimeter over time, and total number of domains over time for two different membrane compositions (DOPC:DPPC with a 1:1 M ratio with 15% Chol and DOPC:DPPC with a 1:2 M ratio with 25% Chol) that yield opposite and nearly inverse phase behavior. This quantitative validation shows that the developed platform can be a valuable tool in complementing experimental practice.

Conformational state diagram of DOPC/DPPCd62/cholesterol mixtures

Raman spectra of aqueous suspensions of vesicles composed of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), deuterated 1,2-dipalmitoyl-d62-sn-glycero-3-phosphocholine (DPPC_d62), and cholesterol (Chol) were studied at room temperature to determine the conformational states of the phospholipid hydrocarbon chains. Deuteration of DPPC_d62 allowed us to characterize the conformational states of DOPC and DPPC_d62 independently. The parameters of Raman peaks, which are sensitive to the conformational order, were studied in a wide range of compositions. It was found that the DOPC molecules are conformationally disordered for all compositions. The conformational state of the DPPC_d62 molecules changes with composition. Their conformational state is influenced by cholesterol-induced partial disordering and DOPC solvation, transforming the DPPC molecules into the disordered state. The conformational state diagram from the Raman experiment was compared with outcomes from the differential scanning calorimetry (DSC) experiment. The Raman spectra also revealed that the DPPC molecules coexist in the disordered and all-trans ordered states for the DOPC/DPPC_d62/Chol mixtures except for the pure liquid-disordered phase.