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 molar ratio with 15% Chol and DOPC:DPPC with a 1:2 molar 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.

[1]  T. Hughes,et al.  Isogeometric analysis of the Cahn–Hilliard phase-field model , 2008 .

[2]  Charles M. Elliott,et al.  Error estimates with smooth and nonsmooth data for a finite element method for the Cahn-Hilliard equation , 1992 .

[3]  H. Amenitsch,et al.  In Situ Determination of Structure and Fluctuations of Coexisting Fluid Membrane Domains , 2015, Biophysical journal.

[4]  Maxim A. Olshanskii,et al.  Trace Finite Element Methods for PDEs on Surfaces , 2016, 1612.00054.

[5]  P. Schwille,et al.  Partitioning, diffusion, and ligand binding of raft lipid analogs in model and cellular plasma membranes. , 2012, Biochimica et biophysica acta.

[6]  K. Gawrisch,et al.  Critical fluctuations in domain-forming lipid mixtures , 2007, Proceedings of the National Academy of Sciences.

[7]  Annalisa Quaini,et al.  Experimental validation of a phase-field model to predict coarsening dynamics of lipid domains in multicomponent membranes , 2020, Biochimica et biophysica acta. Biomembranes.

[8]  Stavroula Sofou,et al.  Masking and triggered unmasking of targeting ligands on liposomal chemotherapy selectively suppress tumor growth in vivo. , 2013, Molecular pharmaceutics.

[9]  Q. Du,et al.  Modelling and simulations of multi-component lipid membranes and open membranes via diffuse interface approaches , 2006, Journal of Mathematical Biology.

[10]  V. Santen.  The Ostwald step rule , 1984 .

[11]  H. Abels,et al.  Thermodynamically Consistent, Frame Indifferent Diffuse Interface Models for Incompressible Two-Phase Flows with Different Densities , 2011, 1104.1336.

[12]  A. Quaini,et al.  A computational study of lateral phase separation in biological membranes , 2018, International journal for numerical methods in biomedical engineering.

[13]  Stavroula Sofou,et al.  Floret-shaped solid domains on giant fluid lipid vesicles induced by pH. , 2012, Langmuir : the ACS journal of surfaces and colloids.

[14]  Kai Simons,et al.  Lipid Rafts As a Membrane-Organizing Principle , 2010, Science.

[15]  Maxim A. Olshanskii,et al.  A Finite Element Method for Elliptic Equations on Surfaces , 2009, SIAM J. Numer. Anal..

[16]  Frederick A. Heberle,et al.  Lipid Rafts: Controversies Resolved, Mysteries Remain. , 2020, Trends in cell biology.

[17]  Axel Voigt,et al.  Dynamics of multicomponent vesicles in a viscous fluid , 2010, J. Comput. Phys..

[18]  G. Pabst,et al.  Monolayer spontaneous curvature of raft-forming membrane lipids , 2013, Soft matter.

[19]  I. Kevrekidis,et al.  Sticky Patches on Lipid Nanoparticles Enable the Selective Targeting and Killing of Untargetable Cancer Cells. , 2016, Langmuir : the ACS journal of surfaces and colloids.

[20]  D. Haar,et al.  Statistical Physics , 1971, Nature.

[21]  Sheereen Majd,et al.  A Simple and Versatile Method for the Formation of Arrays of Giant Vesicles with Controlled Size and Composition , 2013, Advanced materials.

[22]  M. Mayer,et al.  Hydrogel stamping of arrays of supported lipid bilayers with various lipid compositions for the screening of drug-membrane and protein-membrane interactions. , 2005, Angewandte Chemie.

[23]  Sarah L Veatch,et al.  Separation of liquid phases in giant vesicles of ternary mixtures of phospholipids and cholesterol. , 2003, Biophysical journal.

[24]  Soohyung Park,et al.  Reconstitution and functional studies of hamster P-glycoprotein in giant liposomes , 2018, PloS one.

[25]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[26]  Annalisa Quaini,et al.  A decoupled, stable, and linear FEM for a phase-field model of variable density two-phase incompressible surface flow , 2021, Computer Methods in Applied Mechanics and Engineering.

[27]  K. Michalak,et al.  Giant unilamellar vesicles - a perfect tool to visualize phase separation and lipid rafts in model systems. , 2009, Acta biochimica Polonica.

[28]  Jie Shen,et al.  Coarsening kinetics from a variable-mobility Cahn-Hilliard equation: application of a semi-implicit Fourier spectral method. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[29]  Michael Loran Dustin,et al.  Localizing order to boost signaling , 2017, eLife.

[30]  J. Lowengrub,et al.  Locomotion, wrinkling, and budding of a multicomponent vesicle in viscous fluids , 2012 .

[31]  H. Lekkerkerker,et al.  Insights into phase transition kinetics from colloid science , 2002, Nature.

[32]  Petra Schwille,et al.  Probing Lipid Mobility of Raft-exhibiting Model Membranes by Fluorescence Correlation Spectroscopy* , 2003, Journal of Biological Chemistry.

[33]  Deborah A. Brown,et al.  Structure and Function of Sphingolipid- and Cholesterol-rich Membrane Rafts* , 2000, The Journal of Biological Chemistry.

[34]  F. Solis,et al.  Dynamics of coarsening in multicomponent lipid vesicles with non-uniform mechanical properties. , 2014, The Journal of chemical physics.

[35]  T. Taniguchi,et al.  Viscosity Landscape of Phase-Separated Lipid Membrane Estimated from Fluid Velocity Field. , 2020, Biophysical journal.

[36]  L. Bagatolli,et al.  Lipid domains in giant unilamellar vesicles and their correspondence with equilibrium thermodynamic phases: a quantitative fluorescence microscopy imaging approach. , 2009, Biochimica et biophysica acta.

[37]  G. Putzel,et al.  Coarsening dynamics of domains in lipid membranes. , 2012, Biophysical journal.

[38]  James H. Davis,et al.  Fluorescent probe partitioning in giant unilamellar vesicles of 'lipid raft' mixtures. , 2010, The Biochemical journal.

[39]  C. Hayden,et al.  Liposome-Mediated Chemotherapeutic Delivery is Synergistically Enhanced by Ternary Lipid Compositions and Cationic Lipids. , 2019, Langmuir : the ACS journal of surfaces and colloids.

[40]  J. Zimmerberg,et al.  Line tension and interaction energies of membrane rafts calculated from lipid splay and tilt. , 2005, Biophysical journal.

[41]  S. Semrau,et al.  Membrane heterogeneity – from lipid domains to curvature effects , 2009 .

[42]  Siewert J. Marrink,et al.  The molecular face of lipid rafts in model membranes , 2008, Proceedings of the National Academy of Sciences.

[43]  R. Brown,et al.  Sphingolipid organization in biomembranes: what physical studies of model membranes reveal. , 1998, Journal of cell science.

[44]  A. Klymchenko,et al.  Fluorescent probes for lipid rafts: from model membranes to living cells. , 2014, Chemistry & biology.

[45]  D. Tieleman,et al.  Computer simulations of lipid membrane domains. , 2013, Biochimica et biophysica acta.

[46]  Axel Voigt,et al.  Phase-field modeling of the dynamics of multicomponent vesicles: Spinodal decomposition, coarsening, budding, and fission. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.