Figures
Abstract
A novel class of phthalimides functionalized with privileged scaffolds was designed, synthesized and evaluated as potential inhibitors of plasmepsin 2 (Ki: 0.99 ± 0.1 μM for 6u) and plasmepsin 4 (Ki: 3.3 ± 0.3 μM for 6t), enzymes found in the digestive vacuole of the plasmodium parasite and considered as crucial drug targets. Three compounds were identified as potential candidates for further development. The listed compounds were also assayed for their antimalarial efficacy against chloroquine (CQ) sensitive strain (3D7) of Plasmodium falciparum. Assay of twenty seven hydroxyethylamine derivatives revealed four (5e, 6j, 6o and 6s) as strongly active, which were further evaluated against CQ resistant strain (7GB) of P. falciparum. Compound 5e possessing the piperidinopiperidine moiety exhibited promising antimalarial activity with an IC50 of 1.16 ± 0.04 μM. Further, compounds 5e, 6j, 6o and 6s exhibited low cytotoxic effect on MCF-7 cell line. Compound 6s possessing C2 symmetry was identified as the least cytotoxic with significant antimalarial activity (IC50: 1.30 ± 0.03 μM). The combined presence of hydroxyethylamine and cyclic amines (piperazines and piperidines) was observed as crucial for the activity. The current studies suggest that hydroxyethylamine based molecules act as potent antimalarial agent and may be helpful in drug development.
Citation: Singh AK, Rathore S, Tang Y, Goldfarb NE, Dunn BM, Rajendran V, et al. (2015) Hydroxyethylamine Based Phthalimides as New Class of Plasmepsin Hits: Design, Synthesis and Antimalarial Evaluation. PLoS ONE 10(10): e0139347. https://doi.org/10.1371/journal.pone.0139347
Editor: Farhat Afrin, Department of Medical Lab Technology, Faculty of Applied Medical Sciences, Taibah University, SAUDI ARABIA
Received: December 2, 2014; Accepted: September 12, 2015; Published: October 26, 2015
Copyright: © 2015 Singh et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: The author Brijesh Rathi thanks the Department of Science and Technology, Ministry of Science and Technology, India (SR/FT/CS-108/2010) for the young scientist project award. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Malaria is known as one of the most dangerous infectious disease influencing almost half of the global population and imposing a large toll on human health [1–6]. In 2010, ~219 million cases of malaria were identified that resulted in 660 thousand deaths, particularly among children [7]. The malarial parasite, Plasmodium falciparum alone owns the majority of lethal cases worldwide. The malarial problems in the endemic regions have reached a critical stage owing to the arrival of extensive resistance to the current antimalarial drugs and insecticide-resistant mosquitoes. The increasing resistance of P. falciparum to the available drugs, e.g., chloroquine (CQ) [8, 9], and even artemisinin derivatives [10, 11] has abridged the drug efficiency that ultimately has affected endemic regions. Only artemisinin-based combination therapies were recently considered as effective cures of malaria [12] but unfortunately the entrance of resistant malaria parasites is a serious matter. Efforts are underway to produce effective alternative drug molecules; however, until now none of the new molecules has emerged as an antimalarial drug after 1996 [13]. These factors mandate the urgent design and progress of new therapeutics with novel modes of action against multiple targets.
In the malaria parasite hemoglobin degradation in the erythrocytic stage is an inevitable process that occurs in an acidic digestive vacuole (DV) [14,15]. The important class of aspartic proteases, plasmepsins (plasmepsins 1, 2 and 4) and histo-aspartic protease (HAP) found in P. falciparum, have been identified as potential key drug targets for novel antimalarial chemotherapy [16–18]. Plasmepsin 2 and 4 have been demonstrated as crucial enzymes involved in hemoglobin degradation [19] with almost similar mode of functions. The inhibitors of plasmepsins possess significant in vitro and in vivo antimalarial effects, which suggest their suitability as potential drug candidates for antimalarial chemotherapy [20,21]. Plasmodium species other than P. falciparum have only plasmepsin 4 as a digestive vacuole aspartic proteinase that further advocates this enzyme as a strong target for the development of new antimalarials [22]. Recent studies revealed that not only DV plasmepsins but non-digestive vacuole (non-DV) plasmepsins also play critical roles in the survival of the parasite. The overlapped functioning of DV plasmepsins is one of the challenges in the drug development. So we can’t rule out the possibility of non-digestive vacuole plasmepsins as drug targets [23].
Hydroxyethylamine-based molecules have been explored as strong antimalarial agents and also identified as inhibitors of malarial aspartic proteases [24–28]. In hydroxyethylamines, the secondary alcohol is the essential structural element, which plays a crucial role in inhibiting the proteolytic activity of aspartic proteases by mimicking the tetrahedral intermediate during peptide bond cleavage [27,29]. Encouraged by our previous results [30], a new series of functionalized phthalimides possessing chemical variability has been designed anticipating their strong antimalarial actions. In this paper, we report the rational design and synthesis of novel hydroxyethylamine derivatives and their in vitro antimalarial evaluation. The cytotoxic and hemolytic effects were also studied in order to correlate and satisfy the antimalarial activity and the inhibitory activity of the new compounds against plasmepsin 2 and 4.
Material and Methods
Chemistry
General procedure for regioselective ring opening of (2R, 3S)-3-(N-BOC-amino)-1-oxirane-4-phenylbutane by cyclic amines to afford compounds 5b-5e.
Synthesis of compound 5a and 5f was performed following the literature procedure [30]. To a solution of Boc-protected epoxide, 1 (500 mg, 1.9 mmol) in 50 mL of isopropanol, 1-methylpiperazine (0.315 mL, 2.85 mmol) was added and the contents were refluxed for 12 h (Fig 1). To obtain bis-hydroxyethylamine based molecules (5f, Fig 2), piperazine (82 mg, 0.95 mmol) was added to a solution of Boc-protected epoxide, 1 (500 mg, 1.9 mmol) in 50 mL of isopropanol and the contents were refluxed for 12 h. The resulting reaction mixture was concentrated under reduced pressure that afforded above mentioned compounds as white solid in good yield (Table 1) and used for further reactions without purification.
General procedure for deprotection and coupling reaction with N-phthaloyl-L-amino acids to afford 6a-6u.
A solution of Boc-protected hydroxyethylamines, 5a (101 mg, 0.28 mmol) in 5 mL of dichloromethane and 1 mL of trifluoroacetic acid (TFA) was stirred at room temperature for 3 h. The solvent was evaporated under reduced pressure and the oily substance obtained was used for further reactions without purification. In another round bottomed flask, N-phthaloyl-L-leucine (88 mg, 0.34 mmol) in dichloromethane (20 mL) along with triethylamine (TEA) (1 mL), N-(3-dimethylaminopropyl)-N'-ethylcarbodiimide hydrochloride (EDC·HCl) (97 mg, 0.51 mmol) and anhydrous 1-hydroxybenzotriazole (HOBt) (69 mg, 0.51 mmol) were stirred at 0°C for 1 h followed by addition of Boc-deprotected hydroxyethylamine diluted with dichloromethane (3 mL). This procedure was followed to prepare compounds 6a-q. The reaction mixture was stirred at 0°C for another 30 minutes and then at room temperature for 24 h. The reaction mixture was concentrated under reduced pressure. The semi-solid material obtained was diluted with 25 mL of ethylacetate and washed three times with distilled water. The organic layer obtained was dried over anhydrous Na2SO4 and solvent was evaporated under reduced pressure. The obtained compounds were purified by column chromatography using silica gel as stationary phase and chloroform-methanol (96:4) as eluent that afforded a white solid (6a-q) in fair to good yield (Table 2). The spectroscopic details are furnished in supporting information.
To obtain 6r, 6t and 6u, a solution of Boc-protected bis-hydroxyethylpiperazine, 5f (171 mg, 0.28 mmol) in 5 mL of dichloromethane and 1 mL of trifluoroacetic acid (TFA) was stirred at room temperature for 3 h. The solvent was evaporated under reduced pressure and the oily substance obtained was used for further reactions without purification. In another round bottomed flask, N-phthaloyl-L-amino acid (161 mg, 0.62 mmol) in dichloromethane (20 mL) along with triethylamine (TEA) (1 mL), N-(3-dimethylaminopropyl)-N'-ethylcarbodiimide hydrochloride (EDC·HCl) (177 mg, 0.93 mmol) and anhydrous1-hydroxybenzotriazole (HOBt) (125 mg, 0.93 mmol) were stirred at 0°C for 1 h followed by addition of Boc-deprotected hydroxyethylpiperazine diluted with dichloromethane (3 mL). The reaction mixture was stirred at 0°C for another 30 minutes and then at room temperature for 24 h. The reaction mixture was concentrated under reduced pressure. The semi-solid material obtained was diluted with 25 mL of ethylacetate and washed three times with distilled water. The organic layer obtained was dried over anhydrous Na2SO4 and solvent was evaporated under reduced pressure. The obtained compounds were purified by column chromatography using silica gel as stationary phase and chloroform-methanol (96:4) as eluent that afforded a white solid (6r-u) in fair to good yield (Table 3). The spectroscopic details are furnished in supporting information. Compound 6s was prepared following literature procedure [30].
Biology
Inhibition of plasmepsin 2 and 4.
Synthetic compounds were dissolved in DMSO to make stock solutions at concentrations around 100 μM. Dilutions from these stocks were made to obtain solutions to use for studying inhibition of the plasmepsins 2 and 4 from P. falciparum. Assays of enzymatic activity were typically done using 50 micromolar substrate and 50 nanomolar enzyme at pH 4.5. The protein samples were preincubated at the assay pH for 3 minutes to allow conversion from the proenzyme to the mature enzyme form, then immediately mixed with inhibitor solutions and incubated for a further 5 minutes to permit the enzyme-inhibitor complex to form. Following this, the samples were mixed with samples of the substrate to yield a 50 μM substrate concentration in the final assay solution and immediately placed in a recording spectrophotometer to follow the reaction. The cleavage of the chromogenic substrate yields a shift in absorbance that can be quantified by the decrease in absorbance at 300 nm. Compounds were mixed with the enzymes at concentrations ranging from 10 to 50 μM range and the effect on enzyme activity studied in duplicate assays.
Compounds that demonstrated inhibition of enzymatic activity greater than 50% at the higher concentrations were then studied at concentrations as low as 1 μM to determine the extent of binding. Those compounds that showed inhibition in a concentration dependent manner were further studied in a series of assays at varying substrate concentrations (10, 20, 30, 40, 50, 60 μM substrate) and at three concentrations of inhibitory compounds including zero inhibitor concentration. A second inhibitor concentration was chosen that would yield approximately 30–35% inhibition and a third inhibitor concentration was chosen to yield approximately 65–70% inhibition. The three resulting curves of initial velocity versus substrate concentration were fit simultaneously to the competitive inhibition scheme; this yields values of Vmax, Km, and Ki.
In vitro screening of synthesized compounds against P. falciparum 3D7.
For compound screening, SYBR green I-based fluorescence assay was setup as described in literature [31]. Sorbitol synchronized parasites were incubated under normal culture conditions at 2% hematocrit and 1% parasitemia in the absence or presence of increasing concentrations of the compounds. Chloroquine (CQ) was used as positive control, while 0.4% DMSO was used as the negative control. After 48 h of incubation, 100 μL of SYBR Green I solution (0.2 μL of 10,000 X SYBR Green I (Invitrogen)/mL) in lysis buffer (Tris (20 mM; pH 7.5), EDTA (5 mM), saponin (0.008%; w/v), and Triton X-100 (0.08%; v/v) was added to each well and mixed twice gently with multi-channel pipette and incubated in dark at 37°C for 1 h. Fluorescence was measured with a Victor fluorescence multi-well plate reader (Perkin Elmer) with excitation and emission wavelength bands centred at 485 and 530 nm, respectively. The fluorescence counts were plotted against the drug concentration and the 50% inhibitory concentration (IC50) was determined by analysis of dose-response curves.
In vitro measurement of cytotoxic activity against mammalian cell line (MCF7).
Phthalimides were serially diluted in DMSO in a final concentration of 100, 25 and 6.25 μg mL-1 and evaluated for their % cytotoxicity by using MTT-colorimetric assay on MCF-7 cell line [32]. Cell lines were maintained in RPMI-1640 (Hi-Media, Mumbai) medium supplemented with heat inactivated FCS (10% v/v) and 100 U/mL of streptomycin and were cultured in a humidified 5% CO2 atmosphere at 37°C. 100 μL (~2×105) of cells were seeded into wells (per well) on a 96-well plate (Nunc, Maxi Sorp) and incubated for 16 h for adherence. After 16 h, media was aspirated from the wells and the cells were washed with RPMI-1640 without FCS. 75 μL of each compound and control was added to separate wells containing cell lines and incubated at 37°C in humidified 5% CO2 atmosphere for 4 h. After 4 h, media containing the compounds was replaced with 200 μL of normal RPMI-1640 and cells were further incubated for 48 h under same conditions. After 48 h media was replaced by 200 μL of MTT (3-(4, 5-dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide, 0.5 mg/mL of RPMI-1640) and further incubated for 2 h. The formazan crystals were formed and suspended in 100 μL iso-propanol containing 0.06 M HCl and 0.5% SDS. Aliquots were drawn from each well and colour intensity was measured spectrophotometrically in an ELISA plate reader (Biotek, ELx800) at 540 nm. Untreated cells were taken as control with 100% viability and cells without addition of MTT were used as blank. The relative cell viability (%) compared to control cells was calculated by [abs] sample/ [abs] control × 100.
Hemolytic activity.
In vitro hemolytic activity of all synthesized compounds was performed as per reported protocol [33] with slight modifications. In brief, human erythrocytes were collected from Rotary Blood Bank (New Delhi, India) and suspended in phosphate buffer saline (PBS). The suspension (2% hemotocrit) was incubated with different concentrations (3.125 to 100 μg/mL) of the entire compound series at 37°C for 1 h. After 1 h, erythrocytes were pelleted at 3000 g for 10 min (Remi, R-8C/BL) at 4°C. The supernatant was collected and OD450 was determined spectrophotometrically (Biotek, ELx800 UV spectrophotometer) using PBS as a negative control and Triton X-100 for completely lysing the cells. The percent hemolysis was calculated using the following formula: Percent hemolysis = 100-(OD of test/OD of control) × 100. The percent hemolysis was plotted against concentration of compounds for determining the dose cytotoxic to erythrocytes.
Molecular docking
Preparation of Drug Targets and Binding site Analysis.
The crystal structure of plasmepsin 2 and plasmepsin 4 from Plasmodium falciparum was obtained from the RCSB Protein Data Bank (PDB ID: plasmepsin 2 (1LF3), plasmepsin 4 (1LS5) [34]. The binding sites of the both proteins were characterized by CASTp Server [35] and Q-Site Finder [36]. The interactions of the compounds (6p, 6r, 6s, 6t and 6u) with both selected protein residues in the active site were taken using Maestro (Schrödinger) [37]. The targets were prepared for docking by assigning of bond orders, adding hydrogen atoms, removal of all ions, water molecules and all heteroatom from the structures which were stored. The proteins were further optimized at neutral pH and minimized using OPLS_2005 (Optimized Potentials for Liquid Simulations) force field, by converging heavy atoms to RMSD of 0.3 Å.
Preparation of Ligands.
The three dimentional structures of compounds (6p, 6r, 6s, 6t and 6u) were drawn and minimized using Maestro (Schrödinger). Possible tautomers were generated using Ligprep [38] and their ADME properties were predicted using Qikprop [39].
Docking.
The docking study was completed for selected drug targets. All possible (322) tautomers were docked using GLIDE in extra precision mode (XP) with OPLS2005 force field [40,41]. The docked compounds were analysed for binding free energy using MM-GBSA (Molecular Mechanics Generalized Born Surface Area) [42]. All the docked compounds were ranked on the basis of energy estimation.
Results and Discussion
Molecular Design
The structure-based design of new bioactive molecules involves, as a crucial step, the development of a reaction strategy affording enough chemical diversity and structural variety of potentially bioactive products. In this way, not just a single molecule, but a significant part of the chemical space becomes available for the biological investigation making it possible to analyse the structure-activity relationship. Analysing the structure based design of some known potent inhibitors of aspartic proteases (plasmepsin 2) indicated that phenylalanine epoxide (1) moiety plays an important role. Two potential inhibitors, 3 [26] and 4 [26] are shown in Fig 3. Further, the structural analysis of these inhibitors and their interaction with plasmepsins revealed that most of the compounds possess some key structural features responsible for significant inhibition: i) the hydroxyethylamine moiety as a transition state mimicking scaffold [23]; ii) P1 and P2 positions as an aromatic system or a branched carbon side chain [43] and iii) P1’ position as an extended aromatic and aliphatic side chain. The above mentioned features were taken into consideration and a series of new hydroxyethylamine derivatives has been designed (Fig 3). The P1 position was kept constant and the P2 and P1’ moieties were varied. The P2 position was varied using different phthaloyl fragments with different natural amino-acid residues taking into account the high biological potentials of these scaffolds [44]. The medical efficacy of piperazines and piperidines has been explored by employing them as core isosteres [45–47]. Additionally, reports in the literature advocated the importance of piperazine and piperidine moieties in design of the potential candidates for antimalarial chemotherapy [48–50].
Chemistry
The synthetic route adopted to obtain various hydroxyethylamine derivatives is described in Figs 1 and 2. The latter one represents synthesis of hydroxyethylamine based phthalimides having C2 symmetry. Hydroxyethylamine moiety can be obtained by regioselective ring opening of phenylalanine epoxide (Table 1) [48]. The ring opening was performed by refluxing Boc-protected (2R,3S) phenylalanine epoxide (1) with moderately strong nucleophiles, piperazine, 1-ethylpiperazine, 1-methylpiperazine, 4-piperidinopiperidine, 1-(4-fluorobenzyl)piperazine or 1-(4-bromobenzyl)piperazine in isopropanol for 12 h resulting in Boc-protected hydroxyethylamines (5a-f). Compounds (5a-f, Table 2) were deprotected using 20% trifluoroacetic acid in dichloromethane stirred at room temperature for 3 h. Various N-phthaloyl-L-amino acids were prepared from the fusion of L-amino acids with phthalic anhydride following the literature procedure [51]. N-phthaloyl-L-amino acids were recrystallized using 20% ethanol and water (v/v) and their composition was confirmed by reported melting points [52]. N-phthaloyl-L-amino acids and deprotected amines were coupled using coupling reagent, EDC·HCl and HOBt in the presence of triethylamine at room temperature to give the final products (6a-u) in good yields (Tables 2 and 3). The newly prepared compounds were characterised by IR, NMR (1H and 13C) and HRMS which supported the proposed structures. The spectroscopic details of compounds 5a, 5f and 6s were confirmed with the literature [30].
Biology
Assays of inhibitory activity against plasmepsins 2 and 4.
Plasmepsins, a family of aspartic proteases of malaria parasite, are known to participate in a wide variety of cellular processes essential for parasite survival. Plasmepsin 2 from P. falciparum has proven to be an excellent enzyme system for testing the inhibitory potential of new synthetic compounds. It can be prepared by recombinant methods as the inactive zymogen, proplasmepsin 2, and it converts smoothly to the mature and catalytically active enzyme, plasmepsin 2, by simple incubation at pH 4.5. Similarly, plasmepsin 4 from P. falciparum can be prepared and activated. Catalytic activity of the plasmepsins can be studied by following the change in absorbance at 300 nm as cleavage occur between a Phe and a nitroPhe in the synthetic substrate, Lys-Pro-Ile-Glu-Phe-NO2Phe-Arg-Leu [53]. The compounds listed in Tables 1, 2 and 3 were assayed for inhibition by preincubating with differnt concentrations of the compound with equal aliquots of either plasmepsin 2 or 4 followed by mixing with the assay substrate. All 27 compounds described in this report were tested; those that showed inhibition at 50% or greater at 10 μM were studied by additional assays where the substrate concentration was varied at three different inhibitor concentrations. Fitting of the data (see Methods) yielded values of the inhibition constant, Ki, for the compounds and these values are reported in Table 4. It can be seen that for plasmepsin 2 compounds 6u and 6r have Ki values in the micromolar to submicromolar range. For compounds 6p and 6t tested against plasmepsin 4, Ki values were higher, in the micromolar range. These values are adequate to direct further drug optimization studies.
Ki values are reported in micromolar units.
Derivative 6r has a Leu-like side chain at the R position, which is expected to be the P2 position in the Schechter and Berger nomenclature [43] when binding to the active site of the plasmepsins. Compound 6u has an Ile-like side chain in the same position. Compounds 6n and 6s, with a benzyl-like side chain, apparently do not fit well into the pocket where the R group would fit, which would be the S2 pocket in the enzyme active site. The benzyl group next to the–OH would fit into the S1 pocket of the active site of enzyme. Furthermore, the compounds 8a-8u have the phthaloyl group on the N-terminal side of the structure and this should be accommodated in the S3 site of the plasmepsins. The phthaloyl group is much larger than the t-butyl group present on 5a-5f and this can account for why the compounds 5a-5f does not show inhibition when they are tested against plasmepsin 2 and plasmepsin 4.
Compound 6p in Table 2 has a combination of a sec-butyl group in the R position (equivalent to an Ile residue in the P2 position of the synthesised compound) and a 4-bromobenzyl group in the R1 position. This pairing produces a Ki value of 27.4 μM for plasmepsin 4. Apparently, the Ile residue in the P2 position is of optimal size, as smaller or larger substituents do not work as well and hydrophobicity is required to interact with the S2 pocket of the enzyme. Consideration of the substitution in the R1 position shows that the bromobenzyl group is optimal; the fluorobenzyl group does not provide as strong interaction as bromobenzyl.
Compound 6t may be compared to the other compounds in Table 3. It is seen that the isopropyl substitution (equivalent to a Valine side chain) at, possibly, the P2 position gives optimal binding with a measured Ki value of 3.3 μM. These results provide clues to the structure activity relationship for new inhibitors of the plasmepsins.
In vitro antimalarial activity.
We have studied the antimalarial effect of all the listed molecules in asynchronized culture conditions in order to identify a putative lead structure for pharmacological applications. Compounds, 5a-6u were screened for their ability to inhibit the growth of parasite activity, which was determined by the standardized SYBR green based assay [54]. The IC50 values for the inhibition of parasite growth by compounds (5a-6u) are listed in Table 5. The minimal inhibitory effect on the parasite growth activity was exhibited by 6q (IC50: 8.54 ± 0.04 μM). It was observed that four of the compounds, 5e (IC50: 1.16 ± 0.04 μM), 6j (IC50: 1.33 ± 0.04 μM), 6o (IC50: 1.25 ± 0.08 μM), and 6s (IC50: 1.30 ± 0.03 μM) inhibited the parasite growth in a concentration-dependent manner and showed the maximum inhibition as compared with other compounds. Further, these most active members were also evaluated for their antiplasmodial activity against CQ resistant strain (7GB) of P. falciparum and a significant inhibition of the parasite growth was observed.
IC50 values for most potent compounds (5e, 6j, 6o, 6s) against CQ resistant strain (7GB) over 48 h incubation are shown in the brackets.
Effect of most potent compound (5e) on P. falciparum developmental stages.
Further, in order to assess the effect of potent molecule 5e on the parasite developmental stages (asexual), we treated tightly synchronized ring-stage parasite cultures with 5e at 5 μM concentration and counted the parasites at different developmental stages (rings, trophozoites and schizonts) in Giemsa-stained smears at three time points (36, 44, and 52 h after treatment). The control (untreated) parasite progressed through the first cell cycle (0–48 h), increased in size and density, developed into trophozoites and schizonts and subsequently formed viable merozoites. These merozoites were able to invade fresh RBCs and form new ring-stage parasites. (Fig 4). While the drug-treated parasite set developed into the trophozoite and then continued to develop to initiate schizogony, however most of these parasites exhibited developmental arrest in the transition from early to late schizont and also showed swollen food vacuole with altered morphology at trophozoite stages (Fig 4).
Notably, drug-treated parasites were observed with abnormal morphologies mostly in trophozoite and schizont stages. It was observed that parasites developed into multiple daughter nuclei but failed to undergo cytokinesis to develop mature merozoites (Fig 4). Therefore development of new ring-stage parasites were significantly inhibited in the drug-treated sets and showed lower total parasitemia at 52 h as compared with the control set (Fig 5).
Cytotoxic and hemolytic assay.
Assaying all new compounds for cytotoxic effect, 13 compounds were found to display low cytotoxic effects on MCF7 cell line. Interestingly, four potent molecules with significant antimalarial activity (5e, 6j, 6o and 6s) showed low cytotoxic effect. Compound 5f showed maximal toxicity among all tested compounds followed by compound 5a, 5c, 6k, 6m and 6q with IC50 values above 8 μM. The compounds (5e and 6s) were found to be less toxic. Compound, 6s possessing C2 symmetry was identified as least cytotoxic as compared to other molecules.
In the course of the examination on the hemolysis of human erythrocytes, we found that 11 compounds showed lower hemotoxicity. Compound 5a, 6p and 6r possessed least hemotoxic character however, 5b possessing 4-ethylpiperazine showed highest hemotoxicity followed by 6k and 6m. Compounds 5e, 6d and 6j possessing significant antimalarial activity, were identified with lower hemolysis. The results indicate that inhibitory effect of potent molecules at their effective concentrations is independent of cytotoxicity and hemolysis.
Molecular Docking
Plasmepsin 2 and plasmepsin 4 [55] were selected as the drug target of P. falciparum whose structures were retrieved from PDB (3LF3 and 1LS5) and their binding sites were analyzed using CASTp, Q-Site Finder and also verified with literature survey. The top predicted sites were taken for grid formation in docking studies. The structures of compounds (6p, 6r, 6s, 6t, 6u) were prepared and all the possible tautomers were generated. These compounds were found to follow Lipinski's rule of 5 and show no reactive functional groups according to the ADME calculations (S1 Table).
Docking studies of all five compounds with plasmepsin 2 and plasmepsin 4 of Plasmodium falciparum were carried out using XP Glide (Schrödinger) followed by the estimation of binding free energy estimates (S2 Table and S3 Table). Based on Glide Scores and their binding free energy values, compounds (6r, 6u, 6s) were identified as top hits for plasmepsin 2 and 6r, 6u, 6t for plasmepsin 4 (Figs 6 and 7).
Compound 6u showed hydrogen bonding with plasmepsin 2 through residues Ser79 and Ser218 and hydrophobic interactions (In green, S49 Fig). Compound 6t also showed hydrogen bonding with plasmepsin 4 through residues Ser79, Gly78 and Tyr192 including π-π stacking with Phe120 and hydrophobic interactions (In green, S50 Fig). These results advocate 6u and 6t as potential starting molecules for the development of new inhibitors of plasmepsin 2 and 4.
Structure-Activity relationship
Functionalised hydroxyethylamine based molecules were synthesised and screened for in vitro antimalarial potency. During the design of molecules, scaffold 2 was kept constant and its flanking sides were diversified (Fig 3). It is evident that the substituent at P1’ position plays an important role regarding the potency of the compounds against P. falciparum. Bulkier the group, either aliphatic (5e) or aromatic (6o and 6j), the potency is increased. Compounds with smaller aliphatic groups viz., methyl and ethyl directly attached to fourth position of piperazine are less effective against parasitemia suppression. At the P2 position compounds that exhibit an extended aromatic system (R = benzyl group) enhanced the lethal effect on the parasite followed by sec-butyl and isobutyl groups. The above observation was further reinforced when 6l and 6p, possessing 4-(p-fluorobenzyl)piperazine and 4-(p-bromobenzyl)piperazine, respectively, displayed higher potency (6l, IC50: 3.36 ± 0.05 μM; 6p, IC50: 3.28 ± 0.04 μM) than their positional isomers 6i (IC50: 5.27 ± 0.06 μM) and 6m (IC50: 6.96 ± 0.10 μM), respectively as shown in Fig 8.
Conclusion
In conclusion, we have reported the antimalarial activity of hydroxyethylamine based phthalimides. Significantly, the obtained functionalized small molecules were found to provide a new avenue for the discovery of antimalarial agents against CQ resistant strain of P. falciparum. Almost all the molecules were active as anticipated owing to the presence of hydroxyethylamine and piperazine or piperidine moieties. The lower Ki values of all the screened compounds showed their efficacy as inhibitors of plasmepsin 2 (6u) and plasmepsin 4 (6t). The experiments demonstrated that the treatment of 5e inhibits the parasite to transit from early to late schizont stage, which was evidenced by absence of ring stage. Among all the compounds tested, 5e (an intermediate) was proved to be the most potent compound. However, compound 5e did not inhibit the activity of plasmepsins (2 and 4), which advocates the possibility of other drug targets like non-digestive vacuole plasmepsins which, in fact, might be the crucial drug targets according to the current understanding of plasmepsin biology. As there is lack of potent inhibitors of non-DV plasmepsins, we can’t rule out 5e and similar compounds as inhibitor of non-DV plasmepsin, which has been proved as good recipe for research in the area of plasmepsin medicinal chemistry.
Supporting Information
S49 Fig. Protein-ligand contacts of plasmepsin 2 (PDB ID: 1LF3) with top hit 6u.
https://doi.org/10.1371/journal.pone.0139347.s049
(TIF)
S50 Fig. Protein-ligand contacts of plasmepsin 4 (PDB ID: 1LS5) with top hit 6t.
https://doi.org/10.1371/journal.pone.0139347.s050
(TIF)
S1 Table. Physicochemical properties of compounds (5e, 6p, 6r, 6s, 6t, and 6u).
https://doi.org/10.1371/journal.pone.0139347.s051
(DOCX)
S2 Table. The XP GScore and binding free energy values of potent compounds docked to plasmepsin-2 (PDB ID: 2LF3).
https://doi.org/10.1371/journal.pone.0139347.s052
(DOCX)
S3 Table. The XP GScore and binding free energy values of potent compounds docked to plasmepsin-4 (PDB ID: 1LS5).
https://doi.org/10.1371/journal.pone.0139347.s053
(DOCX)
S2 Text. Spectroscopic data of synthesised compounds (5b-6r, 6t and 6u).
https://doi.org/10.1371/journal.pone.0139347.s055
(DOC)
Acknowledgments
The author Brijesh Rathi thanks Department of Science and Technology, Ministry of Science and Technology, India (SR/FT/CS-108/2010) for young scientist project award. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Author Contributions
Conceived and designed the experiments: BR. Performed the experiments: BR AKS SR YT NEG BMD VR NS. Analyzed the data: BR BMD PCG NL BKS MR. Contributed reagents/materials/analysis tools: BR AKS BKS BMD SR VR. Wrote the paper: BR BMD AKS.
References
- 1. Collins WE, Jeffery GM. Plasmodium malariae: parasite and disease. Clin Microbiol Rev. 2007; 20: 579–92. pmid:17934075
- 2. Good MF. Vaccine-induced immunity to malaria parasites and the need for novel strategies. Trends Parasitol. 2005; 21: 29–34. pmid:15639738
- 3. Greenwood BM, Bojang K, Whitty CJM, Targett GAT. Malaria. Lancet 2005; 365: 1487–98. pmid:15850634
- 4. Kappe SH, Vaughan AM, Boddey JA, Cowman AF. That was then but this is now: malaria research in the time of an eradication agenda. Science 2010; 328: 862–6. pmid:20466924
- 5. Snow RW, Guerra CA, Noor AM, Myint HY, Hay SI. The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature 2005; 434: 214–7. pmid:15759000
- 6. Trager W, Jensen JB. Human malaria parasites in continuous culture. Science 1976; 193: 673–5. pmid:781840
- 7.
World Health Organisation. Malaria. Fact sheet N°94 2012.
- 8.
World Health Organization. Guidelines for the treatment of malaria. 2006.
- 9.
Global Report on Antimalarial Drug Efficacy and Drug Resistance: 2000–2010. Geneva, Switzerland, WHO Press. 2010.
- 10. Alonso PL, Brown G, Arevalo-Herrera M, Binka F, Chitnis C, Collins F, et al. A research agenda to underpin malaria eradication. PLoS Med. 2011; 8: e1000406. pmid:21311579
- 11. Dondorp AM, Nosten F, Yi P, Das D, Phyo AP, Tarning J, et al. Artemisinin resistance in Plasmodium falciparum malaria. N Engl J Med. 2009; 361: 455–67. pmid:19641202
- 12.
Jacqueline, R. Halt called on single-drug antimalarial prescriptions. 2006.
- 13. Chugh M, Sundararaman V, Kumar S, Reddy VS, Siddiqui WA, Stuart KD, et al. Protein complex directs hemoglobin-to-hemozoin formation in Plasmodium falciparum. Proc Natl Acad Sci. 2013; 110: 5392–7. pmid:23471987
- 14. Francis SE, Sullivan DJ Jr, Goldberg DE. Hemoglobin metabolism in the malaria parasite Plasmodium falciparum. Annu Rev Microbiol. 1997; 51: 97–123. pmid:9343345
- 15. Moura PA, Dame JB, Fidock DA. Role of Plasmodium falciparum digestive vacuole plasmepsins in the specificity and antimalarial mode of action of cysteine and aspartic protease inhibitors. Antimicrob Agents Chemother. 2009; 53: 4968–78. pmid:19752273
- 16. Aureggi V, Ehmke V, Wieland J, Schweizer WB, Bernet B, Bur D, et al. Potent Inhibitors of Malarial Aspartic Proteases, the Plasmepsins, by Hydroformylation of Substituted 7-Azanorbornenes. Chem Eur J. 2013; 19: 155–64. pmid:23161835
- 17. Bhaumik P, Gustchina A, Wlodawer A. Structural studies of vacuolar plasmepsins. Biochim Biophys Acta 2012; 824: 207–23.
- 18. Miura T, Hidaka K, Azai Y, Kashimoto K, Kawasaki Y, Chen SE, et al. Optimization of plasmepsin inhibitor by focusing on similar structural feature with chloroquine to avoid drug-resistant mechanism of Plasmodium falciparum. Bioorg Med Chem Lett. 2014; 24: 1698–701. pmid:24631188
- 19. Valiente PA, Batista PR, Pupo A, Pons T, Valencia A, Pascutti PG. Predicting functional residues in Plasmodium falciparum plasmepsins by combining sequence and structural analysis with molecular dynamics simulations. Proteins 2008; 73: 440–457. pmid:18442137
- 20. Le Bonniec S, Deregnaucourt C, Redeker V, Banerjee R, Grellier P, Goldberg DE, et al. Plasmepsin II, an Acidic Hemoglobinase from the Plasmodium falciparum Food Vacuole, Is Active at Neutral pH on the Host Erythrocyte Membrane Skeleton. J Biol Chem. 1999; 274: 14218–23. pmid:10318841
- 21. Wyatt DM, Berry C. Activity and inhibition of plasmepsin IV, a new aspartic proteinase from the malaria parasite, Plasmodium falciparum. FEBS Lett. 2002; 513: 159–62. pmid:11904142
- 22. Janka L, Clemente J, Vaiana N, Sparatore A, Romeo S, Dunn BM. Targeting the Plasmepsin 4 Orthologs of Plasmodium sp. with “Double Drug” Inhibitors. Protein Pept Lett. 2008; 15: 868–73. pmid:18991760
- 23. Meyers MJ, Goldberg DE. Recent advances in plasmepsin medicinal chemistry and implications for future antimalarial drug discovery efforts. Curr Top Med Chem. 2012; 12: 445–55. pmid:22242846
- 24. Jaudzems K, Tars K, Maurops G, Ivdra N, Otikovs M, Leitans J, et al. Plasmepsin Inhibitory Activity and Structure-Guided Optimization of a Potent Hydroxyethylamine-Based Antimalarial Hit. ACS Med Chem Lett. 2014; 5: 373–7. pmid:24900843
- 25. Muthas D, Nöteberg D, Sabnis YA, Hamelink E, Vrang L, Samuelsson B, et al. Synthesis, biological evaluation, and modeling studies of inhibitors aimed at the malarial proteases plasmepsins I and II. Bioorg Med Chem. 2005; 13: 5371–90. pmid:16054370
- 26. Nöteberg D, Hamelink E, Hultèn J, Wahlgren M, Vrang L, Samuelsson B, et al. Design and Synthesis of Plasmepsin I and Plasmepsin II Inhibitors with Activity in Plasmodium falciparum-Infected Cultured Human Erythrocytes. J Med Chem. 2003; 46: 734–46. pmid:12593654
- 27.
Ogden RC, Flexner CW. Protease inhibitors in AIDS therapy; Marcel Dekker New York, USA; 2001.
- 28. Ciana CL, Siegrist R, Aissaoui H, Marx L, Racine S, Meyer S, et al. Novel in vivo active anti-malarials based on a hydroxy-ethyl-amine scaffold. Bioorg Med Chem Lett. 2013; 23: 658–62. pmid:23260352
- 29. Conceicao de Souza M, Goncalves-Silva T, Moreth M, RB Gomes C, Roland Kaiser C, de Oliveira Henriques Md, et al. Synthesis and In Vivo Antimalarial Evaluation of Novel Hydroxyethylamine Derivatives. Med Chem. 2012; 8: 266–72. pmid:22385178
- 30. Rathi B, Singh AK, Kishan R, Singh N, Latha N, Srinivasan S, et al. Functionalized hydroxyethylamine based peptide nanostructures as potential inhibitors of falcipain-3, an essential proteases of Plasmodium falciparum. Bioorg Med Chem. 2013; 21: 5503–9. pmid:23810423
- 31. Smilkstein M, Sriwilaijaroen N, Kelly JX, Wilairat P, Riscoe M. Simple and inexpensive fluorescence-based technique for high-throughput antimalarial drug screening. Antimicrob Agents Chem. 2004; 48: 1803–6.
- 32. Rajesh Sharma GL. Studies on antimycotic properties of Datura metel. J Ethnopharmacol. 2002; 80: 193–7. pmid:12007710
- 33. Latoud C, Peypoux F, Michel G, Genet R, Morgat JL. Interactions of antibiotics of the iturin group with human erythrocytes. Biochim Biophys Acta. 1986; 856: 526–35. pmid:3964695
- 34. Bernstein FC, Koetzle TF, Williams GJB, Meyer J, Brice MD, Rodgers JR, et al. The protein data bank: A computer-based archival file for macromolecular structures. J Mol Biol. 1977; 112: 535–42. pmid:875032
- 35. Dundas J, Ouyang Z, Tseng J, Binkowski A, Turpaz Y, Liang J. CASTp: computed atlas of surface topography of proteins with structural and topographical mapping of functionally annotated residues. Nucleic Acids Res. 2006; 34: W116–8. pmid:16844972
- 36. Laurie AT, Jackson RM Q-SiteFinder: an energy-based method for the prediction of protein-ligand binding sites. Bioinformatics 2005; 21: 1908–16. pmid:15701681
- 37.
Maestro, version 9.8 Schrödinger, LLC, New York, NY, 2014.
- 38.
LigPrep, version 3.0, Schrödinger, LLC, New York, NY, 2014.
- 39.
QikProp, version 4.0, Schrödinger, LLC, New York, NY, 2014.
- 40. Friesner RA, Banks JL, Murphy RB, Halgren TA, Klicic JJ, Mainz DT, et al. Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. J Med Chem. 2004; 47: 1739–49. pmid:15027865
- 41. Friesner RA, Murphy RB, Repasky MP, Frye LL, Greenwood JR, Halgren TA, et al. Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. J Med Chem. 2006; 49: 6177–96. pmid:17034125
- 42. Rastelli G, Rio AD, Degliesposti G, Sgobba M. Fast and accurate predictions of binding free energies using MM-PBSA and MM−GBSA. J Comput Chem. 2010; 31: 797–810. pmid:19569205
- 43. Schechter I, Berger A. On the size of the active site in proteases. I. Papain. Biochem Biophys Res Commun. 1967; 27: 157–62. pmid:6035483
- 44. Sharma U, Kumar P, Kumar N, Singh B Recent advances in the chemistry of phthalimide analogues and their therapeutic potential. Mini Rev Med Chem. 2010; 10: 678–704. pmid:20402635
- 45. Standen OD. Activity of piperazine, in vitro, against Ascaris lumbricoides. Br Med J. 1955; 2: 20–2. pmid:14378628
- 46. White RHR, Standen OD. Piperazine in the Treatment of Threadworms in Children. Br Med J. 1953; 2: 755–7. pmid:13082101
- 47. Alm B. Piperidine: effects on locomotor activity and brain monoamine turnover. Psychopharmacology 1976; 50: 301–4. pmid:826977
- 48. Cunico W, Gomes CR, Facchinetti V, Moreth M, Penido C, Henriques MG, et al. Synthesis, antimalarial evaluation and molecular modeling studies of hydroxyethylpiperazines, potential aspartyl protease inhibitors, part 2. Eur J Med Chem. 2009; 44: 3816–20. pmid:19403210
- 49. Mendoza A, Pérez-Silanes S, Quiliano M, Pabón A, Galiano S, González G, et al. Aryl piperazine and pyrrolidine as antimalarial agents. Synthesis and investigation of structure-activity relationships. Exp Parasitol. 2011; 128: 97–103. pmid:21354139
- 50. Ndakala AJ, Gessner RK, Gitari PW, October N, White KL, Hudson A, et al. Antimalarial Pyrido [1, 2-a] benzimidazoles. J Med Chem. 2011; 54: 4581–9. pmid:21644541
- 51.
Vogel AI, Tatchell AR, Furnis BS, Hannaford AJ, Smith PWG. 5th ed. Vogel's textbook of practical organic chemistry. Pearson Education India; 1996.
- 52. Al-Hazimi HM, El-Faham A, Ghazzali M, Al-Farhan K. Microwave irradiation: A facile, scalable and convenient method for synthesis of N-phthaloylamino acids. Arabian J Chem. 2012; 5: 285–9.
- 53. Westling J, Cipullo P, Hung SH, Saft H, Dame JB, Dunn BM. Active site specificity of plasmepsin II. Protein Sci. 1999; 8: 2001–9. pmid:10548045
- 54. Myers MA. Direct measurement of cell numbers in microtitre plate cultures using the fluorescent dye SYBR green I. J Immunol Methods 1998; 212: 99–103. pmid:9671157
- 55. Asojo OA, Gulnik SV, Afonina E, Yu B, Ellman JA, Silva AM. Novel Uncomplexed and Complexed Structures of Plasmepsin II, an Aspartic Protease from Plasmodium falciparum. J Mol Biol. 2003; 327: 173–81. pmid:12614616