An Ab Initio Method for Designing Multi-Target Specific Pharmacophores using Complementary Interaction Field of Aspartic Proteases

A probable means to an end: exploring P131 pharmacophoric scaffold to identify potential inhibitors of Cryptosporidium parvum inosine monophosphate dehydrogenase

Compound P131 has been established to inhibit Cryptosporidium parvum's inosine monophosphate dehydrogenase (CpIMPDH). Its inhibitory activity supersedes that of paromomycin, which is extensively used in treating cryptosporidiosis. Through the per-residue energy decomposition approach, crucial moieties of P131 were identified and subsequently adopted to create a pharmacophore model for virtual screening in the ZINC database. This search generated eight ADMET-compliant hits that were examined thoroughly to fit into the active site of CpIMPDH via molecular docking. Three compounds ZINC46542062, ZINC58646829, and ZINC89780094, with favorable docking scores of - 8.3 kcal/mol, - 8.2 kcal/mol, and - 7.5 kcal/mol, were selected. The potential inhibitory mechanism of these compounds was probed using molecular dynamics simulation and Molecular Mechanics Generalized Poisson Boltzmann Surface Area (MM/PBSA) analyses. Results revealed that one of the hits (ZINC46542062) exhibited a lower binding free energy of - 39.52 kcal/mol than P131, which had - 34.6 kcal/mol. Conformational perturbation induced by the binding of the identified hits to CpIMPDH was similar to P131, suggesting a similarity in inhibitory mechanisms. Also, in silico investigation of the properties of the hit compounds implied superior physicochemical properties with regards to their synthetic accessibility, lipophilicity, and number of hydrogen bond donors and acceptors in comparison with P131. ZINC46542062 was identified as a promising hit compound with the highest binding affinity to the target protein and favorable physicochemical and pharmacokinetic properties relative to P131. The identified compounds can serve as a basis for conducting further experimental investigations toward the development of anticryptosporidials, which can overcome the challenges of existing therapeutic options. Graphical abstract P131 and the identified compounds docked in the NAD⁺ binding site of Cryptosporidium parvum IMPDH.

AutoPH4: An Automated Method for Generating Pharmacophore Models from Protein Binding Pockets

Pharmacophore models are widely used in computational drug discovery (e.g., in the virtual screening of drug molecules) to capture essential information about interactions between ligands and a target protein. Generating pharmacophore models from protein structures is typically a manual process, but there has been growing interest in automated pharmacophore generation methods. Automation makes feasible the processing of large numbers of protein conformations, such as those generated by molecular dynamics (MD) simulations, and thus may help achieve the longstanding goal of incorporating protein flexibility into virtual screening workflows. Here, we present AutoPH4, a new automated method for generating pharmacophore models based on protein structures; we show that a virtual screening workflow incorporating AutoPH4 ranks compounds more accurately than any other pharmacophore-based virtual screening workflow for which results on a public benchmark have been reported. The strong performance of the virtual screening workflow indicates that the AutoPH4 component of the workflow generates high-quality pharmacophores, making AutoPH4 promising for use in future virtual screening workflows as well, such as ones that use conformations generated by MD simulations.

Multiple target-based pharmacophore design from active site structures

Health care systems have benefitted from rational drug discovery processes like vHTS, virtual high throughput screening pharmacophores and quantitative structure-activity relationships, and many challenges have been explored using such techniques: decisions on specificity and selectivity are made after screening millions of molecules for multiple targets. Recent challenges in drug research emphasize the design of drugs that bind with more than one target of interest (multi-target) and do not bind with undesirable targets. This work attempts to use a three-dimensional interaction profile of the active site of a class of proteins, identify selective positions for the binding of functional groups, called features, and develop ensembles of multi-targeted pharmacophores that retain specificity and selectivity. The goal of this study is to develop multi-target pharmacophores by computational methods using protein structures alone to guide the discovery of novel inhibitors of plasmepsins, displaying selectivity over their human homologs, cathepsin D and pepsin. The development of such novel tools is attempted using a combination of different approaches such as the molecular interaction field, clique graph and inductive logic programming to identify and compare specific and selective complementary features. The identification of selective combinations of features has resulted in the design of multi-featured specific and selective pharmacophores that are validated using antimalarial compounds in ChEMBL that are known for their anti-plasmepsin II activity. This novel method is computationally less intensive and is applicable to any known class of target structures for finding specific and selective binders simultaneously.