Systematic Parameterization of Flory–Huggins Models from Molecular Dynamics Simulations for Ternary Lipid Mixtures

Ternary lipid mixtures with saturated and unsaturated lipids and cholesterol serve as widely used model systems to study coexisting liquid-ordered (<i>L</i><sub>o</sub>) and liquid-disordered (<i>L</i><sub>d</sub>) domains in biomembranes. Despite advances in both experiments and simulations, a computationally efficient systematic approach to extract effective pairwise lipid-lipid and lipid-cholesterol interactions from molecular dynamics (MD) simulations for mesoscale models remains limited. Moreover, parametrizing Flory-Huggins (FH) or lattice Monte Carlo (MC) models becomes especially challenging for multicomponent systems due to the size and complexity of the pairwise interaction matrix. To address this, we aim to bridge molecular-level interactions and continuum mean-field models by implementing a multiscale framework that integrates coarse-grained (CG-MD) and all-atom (AA-MD) simulations with FH theory and lattice MC models. Pairwise interaction energies (<i>w</i><sub><i>ij</i></sub>) are extracted from radial distribution functions using the reversible work theorem and formulated into FH interaction parameters (χ<sub><i>ij</i></sub>). We apply this multiscale framework to determine phase count and analyze tie-line and bistable phase behavior of (<i>L</i><sub>d</sub>/<i>L</i><sub>o</sub>) domains for ternary lipid mixtures, specifically using the FH model. The extracted χ<sub><i>ij</i></sub> values from MD reproduce domain formation in both MC and FH models. Comparisons across AA-MD, CG-MD, MC, and FH frameworks reveal both quantitative differences and similarities, as well as conserved trends in phase separation behavior. This study establishes a systematic approach to parametrize mesoscale models for multicomponent lipid membranes directly from MD simulations. Furthermore, our approach can serve as a tool for system-specific mesoscale biomembrane modeling studies that embed molecular-level detail and bridge top-down and bottom-up perspectives for future studies.

Authors

Institutions

Publication Details

Journal
Journal of Chemical Theory and Computation
Published
2025-12-26
DOI
https://doi.org/10.1021/acs.jctc.5c01550
Primary Topic
Lipid Membrane Structure and Behavior
Type
article
Field-Weighted Citation Impact
0.00
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

Systematic Parameterization of Flory–Huggins Models from Molecular Dynamics Simulations for Ternary Lipid Mixtures

Sonali Nayak, John E. Straub, David Zwicker
Journal of Chemical Theory and Computation
Lipid Membrane Structure and Behavior
article

Systematic Parameterization of Flory–Huggins Models from Molecular Dynamics Simulations for Ternary Lipid Mixtures

Sonali Nayak, John E. Straub, David Zwicker
article en

Abstract

Ternary lipid mixtures with saturated and unsaturated lipids and cholesterol serve as widely used model systems to study coexisting liquid-ordered (<i>L</i><sub>o</sub>) and liquid-disordered (<i>L</i><sub>d</sub>) domains in biomembranes. Despite advances in both experiments and simulations, a computationally efficient systematic approach to extract effective pairwise lipid-lipid and lipid-cholesterol interactions from molecular dynamics (MD) simulations for mesoscale models remains limited. Moreover, parametrizing Flory-Huggins (FH) or lattice Monte Carlo (MC) models becomes especially challenging for multicomponent systems due to the size and complexity of the pairwise interaction matrix. To address this, we aim to bridge molecular-level interactions and continuum mean-field models by implementing a multiscale framework that integrates coarse-grained (CG-MD) and all-atom (AA-MD) simulations with FH theory and lattice MC models. Pairwise interaction energies (<i>w</i><sub><i>ij</i></sub>) are extracted from radial distribution functions using the reversible work theorem and formulated into FH interaction parameters (χ<sub><i>ij</i></sub>). We apply this multiscale framework to determine phase count and analyze tie-line and bistable phase behavior of (<i>L</i><sub>d</sub>/<i>L</i><sub>o</sub>) domains for ternary lipid mixtures, specifically using the FH model. The extracted χ<sub><i>ij</i></sub> values from MD reproduce domain formation in both MC and FH models. Comparisons across AA-MD, CG-MD, MC, and FH frameworks reveal both quantitative differences and similarities, as well as conserved trends in phase separation behavior. This study establishes a systematic approach to parametrize mesoscale models for multicomponent lipid membranes directly from MD simulations. Furthermore, our approach can serve as a tool for system-specific mesoscale biomembrane modeling studies that embed molecular-level detail and bridge top-down and bottom-up perspectives for future studies.

Journal of Chemical Theory and ComputationVol. 22(1)
Max Planck Institute for Dynamics and Self-Organization (DE), Max Planck School Matter to Life, University of Göttingen (DE)
Openalex Percentile: Top 36%
Lipid Membrane Structure and Behavior
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.