Search

Liquid-liquid phase separation (LLPS) inside the cell often results in biological condensates that can critically impact cell homeostasis. Such phase separation events occur in multiple parts of cells, including the cell membranes, where the so-called “lipid raft” hypothesis posits the formation of ordered domains floating in a sea of disordered lipids. The resulting lipid domains often have functional roles. However, the thermodynamics of lipid phase separation and their resulting mechanistic effects on cell function and dysfunction are poorly understood. Understanding such complex phenomenon in cell membranes, with their diverse lipid compositions, is exceptionally difficult. For this reasons, simple model systems that can recapitulate similar behavior are widely used to study this phenomenon. Despite these simplifications, the timescale and and length scales of domain formation pose a challenge for molecular dynamics (MD) simulations. Thus, most MD studies focus on spontaneous lipid phase separation — essentially measuring the sign (but not the amplitude) of the free energy change upon separation — rather than directly interrogating the thermodynamics. Here, we propose a proof-of-concept pipeline that can directly measure this free energy by combining coarse-grained MD with enhanced sampling protocols using a novel collective variable. This approach will be a useful tool to help connect the thermodynamics of phase separation with the mechanistic insights already available from molecular dynamics simulations.SIGNIFICANCEStandard molecular dynamics simulations can determine the sign the free energy change upon phase separation, but not the amplitude. We present a new method to determine the phase separation free energy for lipid membranes, based on a enhanced sampling using the weighted ensemble method combined with a novel collective variable, validated using coarse-grained simulations applied to several simple systems. The new method will be valuable as a way to develop models that connect molecular-level structural features to the thermodynamics of phase separation.

Fengycins are a class of antifungal lipopeptides synthesized by the bacteria Bacillus subtilis, commercially available as the primary component of the agricultural fungicide Serenade®. They are toxic to fungi, but far less to mammalian cells. One key difference between mammalian and fungal cell membranes is the presence of cholesterol only in the former; recent experimental work showed that the presence of cholesterol reduces fengycin-induced membrane leakage. (1) Since our previous all-atom and coarse-grained simulations (2, 3) suggested that aggregation of membrane-bound fengycin is central to its ability to disrupt membranes, we hypothesized that cholesterol might reduce fengycin aggregation. Here, we test this hypothesis using coarse-grained molecular dynamics simulations, with sampling enhanced via the weighted ensemble method. The results indicate that cholesterol subtly alters the size distribution for fengycin aggregates, limits the lateral range of their membrane disordering, and reduces the ability of aggregates to bend the membrane. Taken together, these phenomena may account for cholesterol’s affects on fengycin activity.