Broad chemical transferability in structure-based coarse-graining.

Compared to top-down coarse-grained (CG) models, bottom-up approaches are capable of offering higher structural fidelity. This fidelity results from the tight link to a higher resolution reference, making the CG model chemically specific. Unfortunately, chemical specificity can be at odds with compound-screening strategies, which call for transferable parameterizations. Here, we present an approach to reconcile bottom-up, structure-preserving CG models with chemical transferability. We consider the bottom-up CG parameterization of 3441 C7O2 small-molecule isomers. Our approach combines atomic representations, unsupervised learning, and a large-scale extended-ensemble force-matching parameterization. We first identify a subset of 19 representative molecules, which maximally encode the local environment of all gas-phase conformers. Reference interactions between the 19 representative molecules were obtained from both homogeneous bulk liquids and various binary mixtures. An extended-ensemble parameterization over all 703 state points leads to a CG model that is both structure-based and chemically transferable. Remarkably, the resulting force field is on average more structurally accurate than single-state-point equivalents. Averaging over the extended ensemble acts as a mean-force regularizer, smoothing out both force and structural correlations that are overly specific to a single-state point. Our approach aims at transferability through a set of CG bead types that can be used to easily construct new molecules while retaining the benefits of a structure-based parameterization.

Resolution limit of data-driven coarse-grained models spanning chemical space.

Increasing the efficiency of materials design remains a significant challenge given the large size of chemical compound space (CCS). The use of a chemically transferable coarse-grained model enables different molecular fragments to map to the same bead type, significantly increasing screening efficiency. Here, we propose new criteria for the design of coarse-grained models allowing for the optimization of their chemical transferability and evaluate the Martini model within this framework. We further investigate the scope of this transferability by parameterizing three Martini-like models in which the number of bead types ranges from 5 to 16. These force fields are fully compatible with existing Martini environments because they are parameterized by interpolating the Martini interaction matrix. We then implement a Bayesian approach to determining which chemical groups are likely to be present on fragments corresponding to specific bead types for each model. We demonstrate that a level of accuracy comparable to Martini is obtained with a force field with fewer bead types, using the water/octanol partitioning free energy (ΔGW→Ol) as our metric for comparison. However, the advantage of including more bead types is a reduction of uncertainty when 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 CCS. Finally, we note that, due to the large size of fragments mapping to a single Martini bead, a resolution limit arises when using ΔGW→Ol as the only descriptor when coarse-graining CCS.

Controlled exploration of chemical space by machine learning of coarse-grained representations.

The size of chemical compound space is too large to be probed exhaustively. This leads high-throughput protocols to drastically subsample and results in sparse and nonuniform datasets. Rather than arbitrarily selecting compounds, we systematically explore chemical space according to the target property of interest. We first perform importance sampling by introducing a Markov chain Monte Carlo scheme across compounds. We then train a machine learning (ML) model on the sampled data to expand the region of chemical space probed. Our boosting procedure enhances the number of compounds by a factor 2 to 10, enabled by the ML model's coarse-grained representation, which both simplifies the structure-property relationship and reduces the size of chemical space. The ML model correctly recovers linear relationships between transfer free energies. These linear relationships correspond to features that are global to the dataset, marking the region of chemical space up to which predictions are reliable; this is a more robust alternative to the predictive variance. Bridging coarse-grained simulations with ML gives rise to an unprecedented database of drug-membrane insertion free energies for 1.3 million compounds.

Assessment of molecular dynamics simulations for amorphous poly(3-hexylthiophene) using neutron and X-ray scattering experiments.

The molecular morphology and dynamics of conjugated polymers in the bulk solid state play a significant role in determining macroscopic charge transport properties. To understand this relationship, molecular dynamics (MD) simulations and quantum mechanical calculations are used to evaluate local electronic properties. In this work, we investigate the importance of system and simulation parameters, such as force fields and equilibration methods, when simulating amorphous poly(3-hexylthiophene) (P3HT), a model semiconducting polymer. An assessment of MD simulations for five different published P3HT force fields is made by comparing results to experimental wide-angle X-ray scattering (WAXS) and to a broad range of quasi-elastic neutron scattering (QENS) data. Moreover, an in silico analysis of force field parameters reveals that atomic partial charges and torsion potentials along the backbone and side chains have the greatest impact on structure and dynamics related to charge transport mechanisms in P3HT.

Drug-Membrane Permeability across Chemical Space.

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

In silico screening of drug-membrane thermodynamics reveals linear relations between bulk partitioning and the potential of mean force.

The partitioning of small molecules in cell membranes-a key parameter for pharmaceutical applications-typically 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.