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.

Two top-down coarse-grained molecular simulation models for a chromonic liquid crystal, 3,6,7,10,11-hexa-(1,4,7-trioxa-octyl)-triphenylene, are tested. We use an extension of the wellknown MARTINI model and develop a new coarse-grained model based on statistical associating fluid theory (SAFT)-γ perturbation theory. For both models, we demonstrate self-assembly in the isotropic phase of the chromonic and we test the effectiveness of both models in terms of the structures of the chromonic aggregates that are produced in solution and the thermodynamics of association. The latter is tested by calculations of the potential of mean force for dimers in solution, which measures the strength of molecular association. SAFT-γ provides valuable insights into the thermodynamics of assembly. Exploration of a range of interactions between unlike sites demonstrates that chromonic self-assembly only occurs in a small parameter space where the hydrophilic-lipophilic balance between aromatic core and ethylene oxide chains is optimal. Outside of this balance, chromonic self-assembly is replaced by the formation of conglomerates of molecules or short stacks. ARTICLE HISTORY

The risk assessment of thousands of chemicals used in our society benefits from adequate grouping of chemicals based on the mode and mechanism of toxic action (MoA).

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...