Unraveling the relation between the chemical structure of small druglike compounds and their rate of passive permeation across lipid membranes is of fundamental importance for pharmaceutical applications. The elucidation of a comprehensive structure–permeability relationship expressed in terms of a few molecular descriptors is unfortunately hampered by the overwhelming number of possible compounds. In this work, we reduce a priori the size and diversity of chemical space to solve an analogous—but smoothed out—structure–property relationship problem. This is achieved by relying on a physics-based coarse-grained model that reduces the size of chemical space, enabling a comprehensive exploration of this space with greatly reduced computational cost. We perform high-throughput coarse-grained (HTCG) simulations to derive a permeability surface in terms of two simple molecular descriptors—bulk partitioning free energy and p K a . The surface is constructed by exhaustively simulating all coarse-grained compounds that are representative of small organic molecules (ranging from 30 to 160 Da) in a high-throughput scheme. We provide results for acidic, basic, and zwitterionic compounds. Connecting back to the atomic resolution, the HTCG predictions for more than 500 000 compounds allow us to establish a clear connection between specific chemical groups and the resulting permeability coefficient, enabling for the first time an inverse design procedure. Our results have profound implications for drug synthesis: the predominance of commonly employed chemical moieties narrows down the range of permeabilities.

The partitioning of small molecules in cell membranes-a key parameter for pharmaceutical applicationstypically relies on experimentally-available bulk partitioning coefficients. Computer simulations provide a structural resolution of the insertion thermodynamics via the potential of mean force, but require significant sampling at the atomistic level. Here, we introduce high-throughput coarse-grained molecular dynamics simulations to screen thermodynamic properties. This application of physics-based models in a large-scale study of small molecules establishes linear relationships between partitioning coefficients and key features of the potential of mean force. This allows us to predict the structure of the insertion from bulk experimental measurements for more than 400,000 compounds. The potential of mean force hereby becomes an easily accessible quantity-already recognized for its high predictability of certain properties, e.g., passive permeation. Further, we demonstrate how coarse graining helps reduce the size of chemical space, enabling a hierarchical approach to screening small molecules.

Increasing the efficiency of materials design and discovery remains a significant challenge, especially given the prohibitively large size of chemical compound space. The use of a chemically transferable coarse-grained model enables different molecular fragments to map to the same bead type, while also reducing computational expense. These properties further increase screening efficiency, as many compounds are screened through the use of a single coarse-grained simulation, effectively reducing the size of chemical compound space. Here, we propose new criteria for the rational design of coarse-grained models that allows for the optimization of their chemical transferability and evaluate the Martini model within this framework. We further investigate the scope of this chemical transferability by parameterizing three Martini-like models, in which the number of bead types ranges from five to sixteen for the different force fields. We then implement a Bayesian approach to determining which chemical groups are more likely to be present on fragments corresponding to specific bead types for each model. We demonstrate that a level of performance and accuracy comparable to Martini can be obtained by using a force field with fewer bead types. However, the advantage of including more bead types is a reduction of uncertainty with respect to back-mapping these bead types to specific chemistries. Just as reducing the size of the coarse-grained particles leads to a finer mapping of conformational space, increasing the number of bead types yields a finer mapping of chemical compound space. Finally, we note that, due to the relatively large size of the chemical fragments that map to a single martini bead, a clear resolution limit arises when using ∆G W→Ol as the only descriptor when coarse-graining chemical compound space.

An honest discussion on the strengths and limitations of molecular dynamics force fields for P3HT through neutron scattering.