Partitioning into phosphatidylcholine-cholesterol membranes: liposome measurements, coarse-grained simulations, and implications for bioaccumulation.

Membrane-water partitioning is an important physical property for the assessment of bioaccumulation and environmental impact. Here, we advance simulation methodology for predicting the partitioning of small molecules into lipid membranes and compare the computational predictions to experimental measurements in liposomes. As a step towards high-throughput screening, we present an automated mapping and parametrization procedure to produce coarse-grained models compatible with the Martini 3 force field. The methodology is general and can also be used for other applications where coarse-grained simulations are appropriate. This article addresses the effect on membrane-water partitioning of adding cholesterol to POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) membranes. Nine contrasting neutral, zwitterionic and charged solutes are tested. Agreement between experiment and simulation is generally good, with the most challenging cases being permanently charged solutes. For all solutes, partitioning is found to be insensitive to membrane cholesterol concentration up to 25% mole fraction. Hence, for assessment of bioaccumulation into a range of membranes (such as those found in fish), partitioning data measured in pure lipid membranes are still informative.

Automated Coarse-Grained Mapping Algorithm for the Martini Force Field and Benchmarks for Membrane-Water Partitioning.

With a view to high-throughput simulations, we present an automated system for mapping and parameterizing organic molecules for use with the coarse-grained Martini force field. The method scales to larger molecules and a broader chemical space than existing schemes. The core of the mapping process is a graph-based analysis of the molecule's bonding network, which has the advantages of being fast, general, and preserving symmetry. The parameterization process pays special attention to coarse-grained beads in aromatic rings. It also includes a method for building efficient and stable frameworks of constraints for molecules with structural rigidity. The performance of the method is tested on a diverse set of 87 neutral organic molecules and the ability of the resulting models to capture octanol-water and membrane-water partition coefficients. In the latter case, we introduce an adaptive method for extracting partition coefficients from free-energy profiles to take into account the interfacial region of the membrane. We also use the models to probe the response of membrane-water partitioning to the cholesterol content of the membrane.

Self-assembly and mesophase formation in a non-ionic chromonic liquid crystal: insights from bottom-up and top-down coarse-grained simulation models.

New coarse-grained models are introduced for a non-ionic chromonic molecule, TP6EO2M, in aqueous solution. The multiscale coarse-graining (MS-CG) approach is used, in the form of hybrid force matching (HFM), to produce a bottom-up CG model that demonstrates self-assembly in water and the formation of a chromonic stack. However, the high strength of binding in stacks is found to limit the transferability of the HFM model at higher concentrations. The MARTINI 3 framework is also tested. Here, a top-down CG model is produced which shows self-assembly in solution in good agreement with atomistic studies and transfers well to higher concentrations, allowing the full phase diagram of TP6EO2M to be studied. At high concentration, both self-assembly of molecules into chromonic stacks and self-organisation of stacks into mesophases occurs, with the formation of nematic (N) and hexagonal (M) chromonic phases. This CG-framework is suggested as a suitable way of studying a range of chromonic-type drug and dye molecules that exhibit complex self-assembly and solubility behaviour in solution.

Assessing the transferability of common top-down and bottom-up coarse-grained molecular models for molecular mixtures.

The performance of three methods for developing new coarse-grained models for molecular simulation is critically assessed. Two bottom-up approaches are employed: iterative Boltzmann inversion (IBI) and the multiscale coarse-graining method (MS-CG), using an atomistic n-octane-benzene reference system. Results are compared to a top-down coarse-graining approach employing the SAFT-γ Mie equation of state. The performance of each methodology is assessed against the twin criteria of local structure prediction and accurate free energy representation. In addition, the transferability of the generated potentials is compared across state points. We examine the extent to which the IBI methodology can be improved by using a multi-reference approach (MS-IBI), and demonstrate how a pressure correction can be employed to improve the results for the MS-CG approach. Additionally, we look at the effect of including angle-terms in the SAFT-γ Mie model. Finally, we discuss in detail the strengths and weaknesses of each method and suggest possible ways forward for coarse-graining, which may eventually address the problems of structure prediction, thermodynamic consistency and improved transferability within a single model.