ABSTRACTSmall solutes have been shown to alter the lateral organization of cell membranes and reconstituted phospholipid bilayers; however, the mechanisms by which these changes happen are still largely unknown. Traditionally, both experiment and simulation studies have been restricted to testing only a few compounds at a time, failing to identify general molecular descriptors or chemical properties that would allow extrapolating beyond the subset of considered solutes. In this work, we probe the competing energetics of inserting a solute in different membrane environments by means of the potential of mean force. We show that these calculations can be used as a computationally-efficient proxy to establish whether a solute will stabilize or destabilize domain phase separation. Combined with umbrella sampling simulations and coarse-grained molecular dynamics simulations, we are able to screen solutes across a wide range of chemistries and polarities. Our results indicate that, for the system under consideration, preferential partitioning and therefore effectiveness in altering membrane phase separation are strictly linked to the location of insertion in the bilayer (i.e., midplane or interface). Our approach represents a fast and simple tool for obtaining structural and thermodynamic insight into the partitioning of small molecules between lipid domains and its relation to phase separation, ultimately providing a platform for identifying the key determinants of this process.SIGNIFICANCEIn this work we explore the relationship between solute chemistry and the thermodynamics of insertion in a mixed lipid membrane. By combining a coarse-grained resolution and umbrella-sampling simulations we efficiently sample conformational space to study the thermodynamics of phase separation. We demonstrate that measures of the potential of mean force—a computationally-efficient quantity—between different lipid environments can serve as a proxy to predict a compound’s ability to alter the thermodynamics of the lipid membrane. This efficiency allows us to set up a computational screening across many compound chemistries, thereby gaining insight beyond the study of a single or a handful of compounds.
Potential of mean force and umbrella sampling simulation for the transport of
5‐oxazolidinone
in heterotetrameric sarcosine oxidase
