Two novel approaches to construct empirical schemes for partial atomic charge calculation were proposed. The charge schemes possess important benefits. First, they produce both topologically symmetrical and environment dependent charges. Second, they can be parameterised to reasonably reproduce ab initio molecular electrostatic potential (MEP), which guarantees their successful use in molecular modelling. To validate the approaches, the parameters of the proposed charge schemes were fitted to best reproduce MEP simultaneously on grids around a set of 227 diverse organic compounds. The residual errors in MEP reproduction due to calculated atomic charges were compared to those due to charges from known charge schemes.

The notion of a contribution of a specific group in an organic molecule’s property and/or activity is both common in our thinking and is still not strictly correct due to the inherent non-additivity of free energy with respect to molecular fragments composing a molecule. The fragment- based drug discovery (FBDD) approach has proven to be fruitful in addressing the above notions. The main difficulty of the FBDD, however, is in its reliance on the low throughput and expensive experimental means of determining the fragment-sized molecules binding. In this article we propose a way to enhance the throughput and availability of the FBDD methods by judiciously using an in silico means of assessing the contribution to ligand-receptor binding energy of fragments of a molecule under question using a previously developed in silico Reverse Fragment Based Drug Discovery (R-FBDD) approach. It has been shown that the proposed structure-based drug discovery (SBDD) type of approach fills in the vacant niche among the existing in silico approaches, which mainly stem from the ligand-based drug discovery (LBDD) counterparts. In order to illustrate the applicability of the approach, our work retrospectively repeats the findings of the use case of an FBDD hit-to-lead project devoted to the experimentally based determination of additive group efficiency (GE)—an analog of ligand efficiency (LE) for a group in the molecule—using the Free-Wilson (FW) decomposition. It is shown that in using our in silico approach to evaluate fragment contributions of a ligand and to estimate GE one can arrive at similar decisions as those made using the experimentally determined activity-based FW decomposition. It is also shown that the approach is rather robust to the choice of the scoring function, provided the latter demonstrates a decent scoring power. We argue that the proposed approach of in silico assessment of GE has a wider applicability domain and expect that it will be widely applicable to enhance the net throughput of drug discovery based on the FBDD paradigm.

In the drug design process, a frequent task is the decomposition of small molecules into fragments. There exist a number of approaches and methods to break molecules into fragments. However, a method that allows the decomposition of molecules into non-overlapping fragments that is meaningful in terms of medicinal chemistry is absent, and in this work, we present a new simple approach for the decomposition of molecules—MedChemFrag. It aims to break drug-like molecules into a set of rings and linkers, which are close to the perception of “fragments” by medicinal chemists. In contrast to most previous efforts aimed at breaking molecules using retrosynthetic feasible rules, our approach strives to preserve the functional groups, which may reveal the specific interaction pattern, e.g., the amide groups.