Ligand Strain Energy in Large Library Docking

Abstract

While small molecule internal strain is crucial to molecular docking, using it in evaluating ligand scores has remained elusive. Here, we investigate a technique that calculates strain using relative torsional populations in the Cambridge Structural Database, enabling fast precalculation of these energies. In retrospective studies of large docking screens of the dopamine D4 receptor and of AmpC β-lactamase, where close to 600 docking hits were tested experimentally, including such strain energies improved hit rates by preferentially reducing the ranks of strained high-scoring decoy molecules. In a 40-target subset of the DUD-E benchmark, we found two thresholds that usefully distinguished between ligands and decoys: one based on the total strain energy of the small molecules and another based on the maximum strain allowed for any given torsion within them. Using these criteria, about 75% of the benchmark targets had improved enrichment after strain filtering. Relying on precalculated population distributions, this approach is rapid, taking less than 0.04 s to evaluate a conformation on a standard core, making it pragmatic for precalculating strain in even ultralarge libraries. Since it is scoring function agnostic, it may be useful to multiple docking approaches; it is openly available at http://tldr.docking.org.

Potential Therapeutic Target Protein Tyrosine Phosphatase-1B for Modulation of Insulin Resistance with Polyphenols and Its Quantitative Structure–Activity Relationship

“…The torsional energy is calculated by multiplying a weight factor (Wtor) with the torsion number of the ligand (Ntor). As the torsional energy (ΔG torsion ) increases, it lowers the binding energy [ 38 , 39 ]. On comparing the torsional energy (ΔG torsion ) of silydianin (+2.49 kcal/mol) with tolbutamide (+1.39 kcal/mol), but also the intermolecular energy as well as the van der Waals, the electrostatic and internal energies were much higher in the silydianin-bound protein complex, which minimized the high torsional energy penalty as compared to tolbutamide.…”

Iterative computational design and crystallographic screening identifies potent inhibitors targeting the Nsp3 Macrodomain of SARS-CoV-2

“…The top 1 million scored compounds from each screen were investigated for intramolecular strain (total strain <7.5 TEU, maximum single torsion strain <2.5 TEU ( 41 )) and hydrogen bonding with Asp22, Ile23, Gly48, Val49, Phe156 and Asp157. Molecules with unsatisfied hydrogen bond donors or more than three unsatisfied acceptors were not considered for experimental evaluation.…”

Structure-based Discovery of Conformationally Selective Inhibitors of the Serotonin Transporter

“…An average of 4358 orientations was sampled, and for each orientation about 187 conformations—over 1.57 x 10 11 ligand configurations in total were sampled in 121,018 core hours (or 5 days over 1000 cores). High ranking molecules were filtered for interactions with Tyr95, Asp98, Tyr176, Ile172, Asn177, Phe335 and Phe341, those adopting strained conformations (Gu et al, 2021) were deprioritized, as were molecules that topologically resembled ∼28,000 annotated aminergic ligands acting at serotonin, dopamine and adrenergic receptors as well as known inhibitors of SERT, DAT or NET with ECFP4-based Tanimoto coefficients (Tcs) < 0.35, based on molecules annotated in ChEMBL20 (Gaulton et al, 2017). Of the remaining molecules, the top ranking 300,000 were clustered for similarity to one-another.…”