Computer-Aided Drug Design of Falcipain Inhibitors: Virtual Screening, Structure–Activity Relationships, Hydration Site Thermodynamics, and Reactivity Analysis

Falcipains (FPs) are hemoglobinases of Plasmodium falciparum that are validated targets for the development of antimalarial chemotherapy. A combined ligand- and structure-based virtual screening of commercial databases was performed to identify structural analogs of virtual screening hits previously discovered in our laboratory. A total of 28 low micromolar inhibitors of FP-2 and FP-3 were identified and the structure–activity relationship (SAR) in each series was elaborated. The SAR of the compounds was unusually steep in some cases and could not be explained by a traditional analysis of the ligand–protein interactions (van der Waals, electrostatics, and hydrogen bonds). To gain further insights, a statistical thermodynamic analysis of explicit solvent in the ligand binding domains of FP-2 and FP-3 was carried out to understand the roles played by water molecules in binding of these inhibitors. Indeed, the energetics associated with the displacement of water molecules upon ligand binding explained some of the complex trends in the SAR. Furthermore, low potency of a subset of FP-2 inhibitors that could not be understood by the water energetics was explained in the context of poor chemical reactivity of the reactive centers of these compounds. The present study highlights the importance of considering energetic contributors to binding beyond traditional ligand–protein interactions.