AFMoC enhances predictivity of 3D QSAR: a case study with DOXP-reductoisomerase

The advancement of multidimensional QSAR for novel drug discovery - where are we headed?

INTRODUCTION

The Multidimensional quantitative structure-activity relationship (multidimensional-QSAR) method is one of the most popular computational methods employed to predict interesting biochemical properties of existing or hypothetical molecules. With continuous progress, the QSAR method has made remarkable success in various fields, such as medicinal chemistry, material science and predictive toxicology. Areas covered: In this review, the authors cover the basic elements of multidimensional -QSAR including model construction, validation and application. It includes and emphasizes the very recent developments of multidimensional -QSAR such as: HQSAR, G-QSAR, MIA-QSAR, multi-target QSAR. The advantages and disadvantages of each method are also discussed and typical examples of their application are detailed. Expert opinion: Although there are defects in multidimensional-QSAR modeling, it is still of enormous help to chemists, biologists and other researchers in various fields. In the authors' opinion, the latest more precise and feasible QSAR models should be further developed by integrating new descriptors, algorithms and other relevant computational techniques. Apart from being applied in traditional fields (e.g. lead optimization and predictive risk assessment), QSAR should be used more widely as a routine method in other emerging research fields including the modeling of nanoparticles(NPs), mixture toxicity and peptides.

Molecular docking as a popular tool in drug design, an in silico travel

New molecular modeling approaches, driven by rapidly improving computational platforms, have allowed many success stories for the use of computer-assisted drug design in the discovery of new mechanism-or structure-based drugs. In this overview, we highlight three aspects of the use of molecular docking. First, we discuss the combination of molecular and quantum mechanics to investigate an unusual enzymatic mechanism of a flavoprotein. Second, we present recent advances in anti-infectious agents' synthesis driven by structural insights. At the end, we focus on larger biological complexes made by protein-protein interactions and discuss their relevance in drug design. This review provides information on how these large systems, even in the presence of the solvent, can be investigated with the outlook of drug discovery.

Chalcone-based Selective Inhibitors of a C4 Plant Key Enzyme as Novel Potential Herbicides

Weeds are a challenge for global food production due to their rapidly evolving resistance against herbicides. We have identified chalcones as selective inhibitors of phosphoenolpyruvate carboxylase (PEPC), a key enzyme for carbon fixation and biomass increase in the C4 photosynthetic pathway of many of the world's most damaging weeds. In contrast, many of the most important crop plants use C3 photosynthesis. Here, we show that 2',3',4',3,4-Pentahydroxychalcone (IC50 = 600 nM) and 2',3',4'-Trihydroxychalcone (IC50 = 4.2 μM) are potent inhibitors of C4 PEPC but do not affect C3 PEPC at a same concentration range (selectivity factor: 15-45). Binding and modeling studies indicate that the active compounds bind at the same site as malate/aspartate, the natural feedback inhibitors of the C4 pathway. At the whole plant level, both substances showed pronounced growth-inhibitory effects on the C4 weed Amaranthus retroflexus, while there were no measurable effects on oilseed rape, a C3 plant. Growth of selected soil bacteria was not affected by these substances. Our chalcone compounds are the most potent and selective C4 PEPC inhibitors known to date. They offer a novel approach to combat C4 weeds based on a hitherto unexplored mode of allosteric inhibition of a C4 plant key enzyme.

Fosmidomycin as an antimalarial drug: a meta-analysis of clinical trials

With first indications of resistance against artemisinin compounds, the development of novel alternative antimalarials remains an urgent need. One candidate is fosmidomycin (Fos), a phosphonic acid derivative. This PRISMA guideline-adhering and PROSPERO-registered systematic review and meta-analysis provides an overview of the state-of-the-art of the clinical development of Fos as an antimalarial. Pooling six clinical trials of Fos against uncomplicated malaria in African children yielded an overall day 28 cure rate of 85% (95% CI: 71-98%); a parasite clearance time of 39 h; and a fever clearance time of 30 h. In four adult cohorts, the corresponding values were 70% (95% CI: 40-100%), 49 and 42 h, respectively. Data suggest that besides the partner drug, formulation determines efficacy. We advocate further clinical development of Fos-combinations. PROSPERO registration number: CRD42014013688.

Exploring Drug Targets in Isoprenoid Biosynthetic Pathway for Plasmodium falciparum

Emergence of rapid drug resistance to existing antimalarial drugs in Plasmodium falciparum has created the need for prediction of novel targets as well as leads derived from original molecules with improved activity against a validated drug target. The malaria parasite has a plant plastid-like apicoplast. To overcome the problem of falciparum malaria, the metabolic pathways in parasite apicoplast have been used as antimalarial drug targets. Among several pathways in apicoplast, isoprenoid biosynthesis is one of the important pathways for parasite as its multiplication in human erythrocytes requires isoprenoids. Therefore targeting this pathway and exploring leads with improved activity is a highly attractive approach. This report has explored progress towards the study of proteins and inhibitors of isoprenoid biosynthesis pathway. For more comprehensive analysis, antimalarial drug-protein interaction has been covered.

Expression, characterization and inhibition of Toxoplasma gondii 1-deoxy-D-xylulose-5-phosphate reductoisomerase

The apicomplexan parasite Toxoplasma gondii, the causative agent of toxoplasmosis, is an important human pathogen. 1-Deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) in the non-mevalonate isoprene biosynthesis pathway is essential to the organism and therefore a target for developing anti-toxoplasmosis drugs. In order to find potent inhibitors, we expressed and purified recombinant T. gondii DXR (TgDXR). Biochemical properties of this enzyme were characterized and an enzyme activity/inhibition assay was developed. A collection of 11 compounds with a broad structural diversity were tested against TgDXR and several potent inhibitors were identified with Ki values as low as 48 nM. Analysis of the results as well as those of Escherichia coli and Plasmodium falciparum DXR enzymes revealed a different structure-activity relationship profile for the inhibition of TgDXR.

Antimalarial and Structural Studies of Pyridine-containing Inhibitors of 1-Deoxyxylulose-5-phosphate Reductoisomerase

1-Deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) in the non-mevalonate isoprene biosynthesis pathway is a target for developing antimalarial drugs. Fosmidomycin, a potent DXR inhibitor, showed safety as well as efficacy against P. falciparum malaria in clinical trials. Based on our previous quantitative structure activity relationship (QSAR) and crystallographic studies, several novel pyridine-containing fosmidomycin derivatives were designed, synthesized and found to be highly potent inhibitors of P. falciparum DXR (PfDXR) having Ki values of 1.9 - 13 nM, with the best one being ~11× more active than fosmidomycin. These compounds also potently block the proliferation of multi-drug resistant P. falciparum with EC50 values as low as 170 nM. A 2.3 Å crystal structure of PfDXR in complex with one of the inhibitors is reported, showing the flexible loop of the protein undergoes conformational changes upon ligand binding and a hydrogen bond and favorable hydrophobic interactions between the pyridine group and PfDXR account for the enhanced activity.

Design of Potential Bisubstrate Inhibitors against Mycobacterium tuberculosis (Mtb) 1-Deoxy-D-Xylulose 5-Phosphate Reductoisomerase (Dxr)-Evidence of a Novel Binding Mode

In most bacteria, the nonmevalonate pathway is used to synthesize isoprene units. Dxr, the second step in the pathway, catalyzes the NADPH-dependent reductive isomerization of 1-deoxy-D-xylulose-5-phosphate (DXP) to 2-C-methyl-D-erythritol-4-phosphate (MEP). Dxr is inhibited by natural products fosmidomycin and FR900098, which bind in the DXP binding site. These compounds, while potent inhibitors of Dxr, lack whole cell activity against Mycobacterium tuberculosis (Mtb) due to their polarity. Our goal was to use the Mtb Dxr-fosmidomycin co-crystal structure to design bisubstrate ligands to bind to both the DXP and NADPH sites. Such compounds would be expected to demonstrate improved whole cell activity due to increased lipophilicity. Two series of compounds were designed and synthesized. Compounds from both series inhibited Mtb Dxr. The most potent compound (8) has an IC₅₀ of 17.8 µM. Analysis shows 8 binds to Mtb Dxr via a novel, non-bisubstrate mechanism. Further, the diethyl ester of 8 inhibits Mtb growth making this class of compounds interesting lead molecules in the search for new antitubercular agents.

From Hansch-Fujita Analysis to AFMoC: A Road to Structure-Based QSAR

Since the pioneering effort of Hansch and Fujita, quantitative structure-activity relationships (QSAR) have proved valuable in optimizing lead structures. Enriching classical 3D-QSAR analysis, which exploits the three-dimensional structure of ligands, with structural information of the target has helped to improve the interpretability of the derived models and to increase their predictive power. One such method is the Adaption of Fields for Molecular Comparison (AFMoC) approach where protein-specifically adapted knowledge-based pair-potentials are tailored to one particular protein by considering additional structural and energetic information about ligands. Here, we summarize applications of AFMoC, describe recent developments, and provide an outlook on how to improve the method.

Thermodynamic Investigation of Inhibitor Binding to 1-Deoxy-D-Xylulose-5-Phosphate Reductoisomerase

Isothermal titration calorimetry (ITC) was used to investigate the binding of six inhibitors to 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR), a target for developing novel anti-infectives. The binding of hydroxamate inhibitors to E. coli DXR is Mg(2+)-dependent, highly endothermic (ΔH: 22.7-24.3 kJ/mol) and entropy-driven, while that of non-hydroxamate compounds is metal ion independent and exothermic (ΔH: -19.4- -13.8 kJ/mol), showing hydration/dehydration of the enzyme metal ion binding pocket account for the drastic ΔH change. However, for DXRs from Plasmodium falciparum and Mycobacterium tuberculosis, the binding of all inhibitors is exothermic (ΔH: -24.9 - -9.2 kJ/mol), suggesting the metal ion binding sites of these two enzymes are considerably less hydrated. The dissociation constants measured by ITC are well correlated with those obtained by enzyme inhibition assays (R(2) = 0.75). Given the rapid rise of antibiotic resistance, this work is of interest since it provides novel structural implications for rational development of potent DXR inhibitors.

Improvement of the Prediction Power of the CoMFA and CoMSIA Models on Histamine H3 Antagonists by Different Variable Selection Methods

The aim of this study is to enhance the predictivity power of CoMFA and CoMSIA models by means of different variable selection algorithms. The genetic algorithm (GA), successive projection algorithm (SPA), stepwise multiple linear regression (SW-MLR), and the enhanced replacement method (ERM) were used and tested as variable selection algorithms. Then, the selected variables were used to generate a simple and predictive model by the multilinear regression algorithm. A set of 74 histamine H₃ antagonists were split into 40 compounds as a training set, and 17 compounds as a test set, by the Kennard-Stone algorithm. Before splitting the data, 17 compounds were randomly selected from the pool of the whole data set as an evaluation set without any supervision, pretreatment, or visual inspection. Among applied variable selection algorithms, ERM had noticeable improvement on the statistical parameters. The r² values of training, test, and evaluation sets for the ERM-MLR model using CoMFA fields were 0.9560, 0.8630, and 0.8460 and using the CoMSIA fields were 0.9800, 0.8521, and 0.9080, respectively. In this study, the principles of organization for economic cooperation and development (OECD) for regulatory acceptability of QSARs are considered.

Inhibition of 1-deoxy-D-xylulose-5-phosphate reductoisomerase by lipophilic phosphonates: SAR, QSAR, and crystallographic studies

1-Deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) is a novel target for developing new antibacterial (including antituberculosis) and antimalaria drugs. Forty-one lipophilic phosphonates, representing a new class of DXR inhibitors, were synthesized, among which 5-phenylpyridin-2-ylmethylphosphonic acid possesses the most activity against E. coli DXR (EcDXR) with a K(i) of 420 nM. Structure-activity relationships (SAR) are discussed, which can be rationalized using our EcDXR:inhibitor structures, and a predictive quantitative SAR (QSAR) model is also developed. Since inhibition studies of DXR from Mycobacterium tuberculosis (MtDXR) have not been performed well, 48 EcDXR inhibitors with a broad chemical diversity were found, however, to generally exhibit considerably reduced activity against MtDXR. The crystal structure of a MtDXR:inhibitor complex reveals the flexible loop containing the residues 198-208 has no strong interactions with the 3,4-dichlorophenyl group of the inhibitor, representing a structural basis for the reduced activity. Overall, these results provide implications in the future design and development of potent DXR inhibitors.

Structures of 1-Deoxy-D-Xylulose-5-Phosphate Reductoisomerase/Lipophilic Phosphonate Complexes

Fosmidomycin, a potent inhibitor of 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR), has antibacterial and antimalaria activity. Due to its poor pharmacokinetics, more lipophilic DXR inhibitors are needed. However, the hydrophobic binding site(s) in DXR remains elusive. Here, pyridine/quinoline containing phosphonates are identified to be DXR inhibitors with IC(50) values as low as 840 nM. We also report three DXR:inhibitor structures, revealing a novel binding mode. The indole group of Trp211 is found to move ~4.6 Å to open up a mainly hydrophobic pocket, where the pyridine/quinoline rings of the inhibitors are located and have strong π-π stacking/charge-transfer interactions with the indole. Docking studies demonstrate our structures could be used to predict the binding modes of other lipophilic DXR inhibitors. Overall, this work shows an important role of Trp211 in inhibitor recognition and provides a structural basis for future drug design and development.

Free-energy-based methods for binding profile determination in a congeneric series of CDK2 inhibitors

Free-energy pathway methods show great promise in computing the mode of action and the free energy profile associated with the binding of small molecules with proteins, but are generally very computationally demanding. Here we apply a novel approach based on metadynamics and path collective variables. We show that this combination is able to find an optimal reaction coordinate and the free energy profile of binding with explicit solvent and full flexibility, while minimizing human intervention and computational costs. We apply it to predict the binding affinity of a congeneric series of 5 CDK2 inhibitors. The predicted binding free energy profiles are in accordance with experiment.

Addressing the malaria drug resistance challenge using flow cytometry to discover new antimalarials

A new flow cytometry method that uses an optimized DNA and RNA staining strategy to monitor the growth and development of the Plasmodium falciparum strain W2mef has been used in a pilot study and has identified Bay 43-9006 1, SU 11274 2, and TMC 125 5 as compounds that exhibit potent (<1 microM) overall and ring stage in vitro antimalarial activity.

Towards new antimalarial drugs: synthesis of non-hydrolyzable phosphate mimics as feed for a predictive QSAR study on 1-deoxy-D-xylulose-5-phosphate reductoisomerase inhibitors

The conversion of 1-deoxy-D-xylulose-5-phosphate (DOXP) to 2-C-methyl-D-erythritol-4-phosphate (MEP) is effectively blocked by 1-deoxy-D-xylulose-5-phosphate reductoisomerase (Dxr) inhibitors such as the natural antibiotic fosmidomycin. Prediction of binding affinities for closely related Dxr ligands as well as estimation of the affinities of structurally more distinct inhibitors within this class of non-hydrolyzable phosphate mimics relies on the synthesis of fosmidomycin derivatives with a broad range of target affinity. Maintaining the phosphonic acid moiety, linear modifications of the lead structure were carried out in an effort to expand the SAR of this physicochemically challenging class of compounds. Synthetic access to a set of phosphonic acids with inhibitory activity (IC(50)) in the range from 1 to >30 microM vs. E. coli Dxr and 0.4 to 20 microM against P. falciparum Dxr is reported.

Novel Deoxyxylulosephosphate-Reductoisomerase Inhibitors: Fosmidomycin Derivatives with Spacious Acyl Residues

1-deoxy-D-xylulose-5-phosphate reductoisomerase (Dxr) represents an essential enzyme of the mevalonate-independent pathway of the isoprenoid biosynthesis. Using fosmidomycin as a specific inhibitor of Dxr, this enzyme was previously validated as target for the treatment of malaria and bacterial infections. The replacement of the formyl residue of fosmidomycin by spacious acyl residues yielded inhibitors active in the micromolar range. As predicted by flexible docking, evidence was obtained for the formation of a hydrogen bond between an appropriately placed carbonyl group in the acyl residue and the main-chain NH of Met214 located in the flexible catalytic loop of the enzyme.

Integrated approach using protein and ligand information to analyze selectivity- and affinity-determining features of carbonic anhydrase isozymes

The application and comparison of selected protein- and ligand-based approaches to elucidate factors important for affinity and selectivity towards the carbonic anhydrase isozymes I, II, and IV are described. Carbonic anhydrases are abundant in pro- and eukaryotes. These enzymes catalyze the reversible hydration of carbon dioxide to bicarbonate and H(+) ions and are thus involved in many important physiological and pathophysiological processes. Due to the fact that the human carbonic anhydrase family consists of 16 closely related isozymes, the design of selective inhibitors is a special challenge for medicinal chemists. In order to extract selectivity-determining features, we applied purely ligand-based 3D QSAR techniques as well as qualitative comparative molecular field analyses of the targets' binding sites using consensus principal component analysis (CPCA). The dataset for the QSAR studies was deliberately compiled from 1,748 inhibitors and comprises about 140 ligands, mainly of the sulfonamide type. Additionally, we employed the novel AFMoC approach, which intrinsically combines protein and ligand information. The simultaneous use of these different techniques gives deeper insight into selectivity and affinity-determining features and provides quantitative models for prediction.

Comparative protein modeling of 1-deoxy-D-xylulose-5-phosphate reductoisomerase enzyme from Plasmodium falciparum: a potential target for antimalarial drug discovery

Plasmodium falciparum 1-deoxy-D-xylulose-5-phosphate reductoisomerase (Pf-DXR) is a potential target for antimalarial chemotherapy. The three-dimensional model (3D) of this enzyme was determined by means of comparative modeling through multiple alignment followed by intensive optimization, minimization, and validation. The resulting model demonstrates a reasonable topology as gauged from the Ramachandran plot and acceptable three-dimensional structure compatibility as assessed by the Profiles-3D score. The modeled monomeric subunit consists of three domains: (1) N-terminal NADPH binding domain, (2) connective or linker domain (with most of the active site residues located in this domain), and (3) a C-terminal domain. This structure proved to be consistent with known DXR crystal structures from other species. The predicted active site compared favorably with those of the templates and appears to have an active site with a highly conserved architecture. Additionally, the model explains several site-directed mutagenesis data. Besides using several protein structure-checking programs to validate the model, a set of known inhibitors of DXR were also docked into the active site of the modeled Pf-DXR. The docked scores correlated reasonably well with experimental pIC50 values with a regression coefficient (R2) equal to 0.84. Results of the current study should prove useful in the early design and development of inhibitors by either de novo drug design or virtual screening of large small-molecule databases leading to development of new antimalarial agents.

Evaluation of fosmidomycin analogs as inhibitors of the Synechocystis sp. PCC6803 1-deoxy-d-xylulose 5-phosphate reductoisomerase

Analogs of the antibiotic fosmidomycin, an inhibitor of the methylerythritol phosphate pathway to isoprenoids, were synthesized and evaluated against the recombinant Synechocystis sp. PCC6803 1-deoxy-d-xylulose 5-phosphate reductoisomerase (DXR). Fosfoxacin, the phosphate analog of fosmidomycin, and its acetyl congener were found to be more potent inhibitors of DXR than fosmidomycin.

Synthesis and Biological Evaluation of Cyclopropyl Analogues of Fosmidomycin as Potent Plasmodium falciparum Growth Inhibitors

A series of fosmidomycin analogues featuring restricted conformational mobility has been synthesized and evaluated as inhibitors of 1-deoxy-D-xylulose 5-phosphate (DOXP) reductoisomerase and as growth inhibitors of P. falciparum. The enantiomerically pure trans-cyclopropyl N-acetyl analogue 3b showed comparable inhibitory activity as fosmidomycin toward E. coli DOXP reductoisomerase and proved equally active when tested in vitro for P. falciparum growth inhibition. Conversely, the alpha-phenyl cis-cyclopropyl analogue 4 showed virtually no inhibition of the enzyme.

Synthesis and Evaluation of 1-Deoxy-d-xylulose 5-Phosphate Analogues as Chelation-Based Inhibitors of Methylerythritol Phosphate Synthase

[structures: see text] A series of 1-deoxy-D-xylulose 5-phosphate (DXP) analogues were synthesized and evaluated as inhibitors of E. coli methylerythritol phosphate (MEP) synthase. In analogues 1-4, the methyl group in DXP was replaced by hydroxyl, hydroxylamino, methoxy, and amino moieties, respectively. In analogues 5 and 6, the acetyl moiety in DXP was replaced by hydroxymethyl and aminomethyl groups. These compounds were designed to coordinate to the active site divalent metal in MEP synthase. The carboxylate (1), methyl ester (3), amide (4), and alcohol (5) analogues were inhibitors with IC50's ranging from 0.25 to 1.0 mM. The hydroxamic acid (2) and amino (6) analogues did not inhibit the enzyme.