Computer-Aided Design of Orally Bioavailable Pyrrolidine Carboxamide Inhibitors of Enoyl-Acyl Carrier Protein Reductase of Mycobacterium tuberculosis with Favorable Pharmacokinetic Profiles

Structure-Based Design and Pharmacophore-Based Virtual Screening of Combinatorial Library of Triclosan Analogues Active against Enoyl-Acyl Carrier Protein Reductase of Plasmodium falciparum with Favourable ADME Profiles

Cost-effective therapy of neglected and tropical diseases such as malaria requires everlasting drug discovery efforts due to the rapidly emerging drug resistance of the plasmodium parasite. We have carried out computational design of new inhibitors of the enoyl-acyl carrier protein reductase (ENR) of Plasmodium falciparum (PfENR) using computer-aided combinatorial and pharmacophore-based molecular design. The Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) complexation QSAR model was developed for triclosan-based inhibitors (TCL) and a significant correlation was established between the calculated relative Gibbs free energies of complex formation (∆∆Gcom) between PfENR and TCL and the observed inhibitory potencies of the enzyme (IC50exp) for a training set of 20 known TCL analogues. Validation of the predictive power of the MM-PBSA QSAR model was carried out with the generation of 3D QSAR pharmacophore (PH4). We obtained a reasonable correlation between the relative Gibbs free energy of complex formation ∆∆Gcom and IC50exp values, which explained approximately 95% of the PfENR inhibition data: pIC50exp=-0.0544×∆∆Gcom+6.9336,R2=0.95. A similar agreement was established for the PH4 pharmacophore model of the PfENR inhibition (pIC50exp=0.9754×pIC50pre+0.1596, R2=0.98). Analysis of enzyme-inhibitor binding site interactions suggested suitable building blocks to be used in a virtual combinatorial library of 33,480 TCL analogues. Structural information derived from the complexation model and the PH4 pharmacophore guided us through in silico screening of the virtual combinatorial library of TCL analogues to finally identify potential new TCL inhibitors effective at low nanomolar concentrations. Virtual screening of the library by PfENR-PH4 led to a predicted IC50pre value for the best inhibitor candidate as low as 1.9 nM. Finally, the stability of PfENR-TCLx complexes and the flexibility of the active conformation of the inhibitor for selected top-ranking TCL analogues were checked with the help of molecular dynamics. This computational study resulted in a set of proposed new potent inhibitors with predicted antimalarial effects and favourable pharmacokinetic profiles that act on a novel pharmacological target, PfENR.

Anti-tuberculosis drug development via targeting the cell envelope of Mycobacterium tuberculosis

Mycobacterium tuberculosis possesses a dynamic cell envelope, which consists of a peptidoglycan layer, a mycolic acid layer, and an arabinogalactan polysaccharide. This envelope possesses a highly complex and unique structure representing a barrier that protects and assists the growth of M. tuberculosis and allows its adaptation to the host. It regulates the immune response of the host cells, causing their damage. Therefore, the cell envelope of M. tuberculosis is an attractive target for vaccine and drug development. The emergence of multidrug-resistant as well as extensively drug resistant tuberculosis and co-infection with HIV prevented an effective control of this disease. Thus, the discovery and development of new drugs is a major keystone for TB treatment and control. This review mainly summarizes the development of drug enzymes involved in the biosynthesis of the cell wall in M. tuberculosis, and other potential drug targets in this pathway, to provide more effective strategies for the development of new drugs.

Computer-Aided Design of Peptidomimetic Inhibitors of Falcipain-3: QSAR and Pharmacophore Models

In this work, antiparasitic peptidomimetics inhibitors (PEP) of falcipain-3 (FP3) of Plasmodium falciparum (Pf) are proposed using structure-based and computer-aided molecular design. Beginning with the crystal structure of PfFP3-K11017 complex (PDB ID: 3BWK), three-dimensional (3D) models of FP3-PEPx complexes with known activities ( IC50exp) were prepared by in situ modification, based on molecular mechanics and implicit solvation to compute Gibbs free energies (GFE) of inhibitor-FP3 complex formation. This resulted in a quantitative structure–activity relationships (QSAR) model based on a linear correlation between computed GFE (ΔΔGcom) and the experimentally measured  IC50exp. Apart from the structure-based relationship, a ligand-based quantitative pharmacophore model (PH4) of novel PEP analogues where substitutions were directed by comparative analysis of the active site interactions was derived using the proposed bound conformations of the PEPx. This provided structural information useful for the design of virtual combinatorial libraries (VL), which was virtually screened based on the 3D-QSAR PH4. The end results were predictive inhibitory activities falling within the low nanomolar concentration range.

Structure-Based Design and in Silico Screening of Virtual Combinatorial Library of Benzamides Inhibiting 2-trans Enoyl-Acyl Carrier Protein Reductase of Mycobacterium tuberculosis with Favorable Predicted Pharmacokinetic Profiles

Background: During the previous decade a new class of benzamide-based inhibitors of 2-trans enoyl-acyl carrier protein reductase (InhA) of Mycobacterium tuberculosis (Mt) with unusual binding mode have emerged. Here we report in silico design and evaluation of novel benzamide InhA-Mt inhibitors with favorable predicted pharmacokinetic profiles. Methods: By using in situ modifications of the crystal structure of N-benzyl-4-((heteroaryl)methyl) benzamide (BHMB)-InhA complex (PDB entry 4QXM), 3D models of InhA-BHMBx complexes were prepared for a training set of 19 BHMBs with experimentally determined inhibitory potencies (half-maximal inhibitory concentrations IC50exp). In the search for active conformation of the BHMB1-19, linear QSAR model was prepared, which correlated computed gas phase enthalpies of formation (∆∆HMM) of InhA-BHMBx complexes with the IC50exp. Further, taking into account the solvent effect and entropy changes upon ligand, binding resulted in a superior QSAR model correlating computed complexation Gibbs free energies (∆∆Gcom). The successive pharmacophore model (PH4) generated from the active conformations of BHMBs served as a virtual screening tool of novel analogs included in a virtual combinatorial library (VCL) of compounds containing benzamide scaffolds. The VCL filtered by Lipinski’s rule-of-five was screened by the PH4 model to identify new BHMB analogs. Results: Gas phase QSAR model: −log₁₀(IC₅₀^exp) = pIC₅₀^exp = −0.2465 × ∆∆H_MM + 7.95503, R² = 0.94; superior aqueous phase QSAR model: pIC₅₀^exp = −0.2370 × ∆∆G_com + 7.8783, R² = 0.97 and PH4 pharmacophore model: pIC50exp = 1.0013 × pIC50exp − 0.0085, R² = 0.95. The VCL of more than 114 thousand BHMBs was filtered down to 73,565 analogs Lipinski’s rule. The five-point PH4 screening retained 90 new and potent BHMBs with predicted inhibitory potencies IC₅₀^pre up to 65 times lower than that of BHMB1 (IC₅₀^exp = 20 nM). Predicted pharmacokinetic profile of the new analogs showed enhanced cell membrane permeability and high human oral absorption compared to current anti-tuberculotics. Conclusions: Combined use of QSAR models that considered binding of the BHMBs to InhA, pharmacophore model, and ADME properties helped to recognize bound active conformation of the benzamide inhibitors, permitted in silico screening of VCL of compounds sharing benzamide scaffold and identification of new analogs with predicted high inhibitory potencies and favorable pharmacokinetic profiles.

Virtual design of novel Plasmodium falciparum cysteine protease falcipain-2 hybrid lactone-chalcone and isatin-chalcone inhibitors probing the S2 active site pocket

We report computer-aided design of new lactone–chalcone and isatin–chalcone (HLCIC) inhibitors of the falcipain-2 (PfFP-2). 3D models of 15 FP-2:HLCIC1-15 complexes with known observed activity (IC₅₀^exp) were prepared to establish a quantitative structure–activity (QSAR) model and linear correlation between relative Gibbs free energy of enzyme:inhibitor complex formation (ΔΔG_com) and IC₅₀^exp: pIC₅₀^exp = −0.0236 × ΔΔG_com+5.082(#); R² = 0.93. A 3D pharmacophore model (PH4) derived from the QSAR directed our effort to design novel HLCIC analogues. During the design, an initial virtual library of 2621440 HLCIC was focused down to 18288 drug-like compounds and finally, PH4 screened to identify 81 promising compounds. Thirty-three others were added from an intuitive substitution approach intended to fill better the enzyme S2 pocket. One hundred and fourteen theoretical IC₅₀ (IC₅₀^pre) values were predicted by means of (#) and their pharmacokinetics (ADME) profiles. More than 30 putative HLCICs display IC₅₀^pre 100 times superior to that of the published most active training set inhibitor HLCIC1.

Target discovery focused approaches to overcome bottlenecks in the exploitation of antimycobacterial natural products

Tuberculosis is a major global health hazard. The search for new antimycobacterials has focused on such as screening combinational chemistry libraries or designing chemicals to target predefined pockets of essential bacterial proteins. The relative ineffectiveness of these has led to a reappraisal of natural products for new antimycobacterial drug leads. However, progress has been limited, we suggest through a failure in many cases to define the drug target and optimize the hits using this information. We highlight methods of target discovery needed to develop a drug into a candidate for clinical trials. We incorporate these into suggested analysis pipelines which could inform the research strategies to accelerate the development of new drug leads from natural products.