Abstract
It is imperative to explore the gigantic available chemical space to identify new scaffolds for drug lead discovery. Identifying potent hits from virtual screening of large chemical databases is challenging and computationally demanding. Rather than the traditional two-dimensional (2D)/three-dimensional (3D) approaches on smaller chemical libraries of a few hundred thousand compounds, we screened a ZINC library of 15 million compounds using multiple computational methods. Here, we present the successful application of a virtual screening methodology that identifies several chemotypes as starting hits against lactate dehydrogenase-A (LDHA). From 29 compounds identified from virtual screening, 17 (58%) showed IC50 values < 63 μM, two showed single-digit micromolar inhibition, and the most potent hit compound had IC50 down to 117 nM. We enriched the database and employed an ensemble approach by combining 2D fingerprint similarity searches, pharmacophore modeling, molecular docking, and molecular dynamics. WaterMap calculations were carried out to explore the thermodynamics of surface water molecules and gain insights into the LDHA binding pocket. The present work has led to the discovery of two new chemical classes, including compounds with a succinic acid monoamide moiety or a hydroxy pyrimidinone ring system. Selected hits block lactate production in cells and inhibit pancreatic cancer cell lines with cytotoxicity IC50 down to 12.26 μM against MIAPaCa-2 cells and 14.64 μM against PANC-1, which, under normoxic conditions, is already comparable or more potent than most currently available known LDHA inhibitors.