Multistage antiplasmodial activity of hydroxyethylamine compounds, in vitro and in vivo evaluations

Malaria, a global threat to the human population, remains a challenge partly due to the fast-growing drug-resistant strains of Plasmodium species. New therapeutics acting against the pathogenic asexual and sexual stages, including liver-stage malarial infection, have now attained more attention in achieving malaria eradication efforts. In this paper, two previously identified potent antiplasmodial hydroxyethylamine (HEA) compounds were investigated for their activity against the malaria parasite's multiple life stages. The compounds exhibited notable activity against the artemisinin-resistant strain of P. falciparum blood-stage culture with 50% inhibitory concentrations (IC50) in the low micromolar range. The compounds' cytotoxicity on HEK293, HepG2 and Huh-7 cells exhibited selective killing activity with IC50 values > 170 μM. The in vivo efficacy was studied in mice infected with P. berghei NK65, which showed a significant reduction in the blood parasite load. Notably, the compounds were active against liver-stage infection, mainly compound 1 with an IC50 value of 1.89 μM. Mice infected with P. berghei sporozoites treated with compound 1 at 50 mg kg−1 dose had markedly reduced liver stage infection. Moreover, both compounds prevented ookinete maturation and affected the developmental progression of gametocytes. Further, systematic in silico studies suggested both the compounds have a high affinity towards plasmepsin II with favorable pharmacological properties. Overall, the findings demonstrated that HEA and piperidine possessing compounds have immense potential in treating malarial infection by acting as multistage inhibitors.


Introduction
Malaria, the deadliest disease, continues to be a threat for the human population worldwide, inordinate in endemic tropical and Sub-Saharan Africa regions. Based on the World Malaria Report 2019, at least 228 million cases were noted in 2018, with an estimated number of 405 000 deaths globally. In 2018, 67% of all deaths reported were of children under the age of ve, making them the most vulnerable group affected by malaria. 1 This highlights the current burden of malaria and the advances achieved in disease control. Malaria is caused by eukaryotic protozoans belonging to Plasmodium species. Out of the ve species, Plasmodium falciparum, is accountable for pathogenic human malaria. 2 The complex life cycle of P. falciparum consists of multiple forms, i.e., asexual (liver and blood-stage) and sexual (gametocyte and transmission stages). All these life stages of the parasite development offer potential drug targets individually. 3 Although several approaches viz. insecticide spraying, insecticide-treated bed nets, and availability of effective chemotherapeutic drugs like artemisinin-based combination therapies (ACTs), 4 are available to control malaria. It is still considered a lethal infectious pathogen, particularly in the developing part of the world. 5 Therefore, ACTs remain the most effective treatments for P. falciparum malaria. 6,7 However, the emergence of increasing resistance to these front-line treatments has severely hindered current efforts to reduce disease burden and reduce mortality rates. [8][9][10] This scenario creates the possibility of not having even a single reliable effective drug capable of treating the P. falciparum malaria, particularly in endemic regions. 10 Therefore, generic ACTs are the only efficient treatments available in these regions with no personalized malarial drug resistance testing being done. 11 Moreover, the non-artemisinin drugs have also shown reduced efficiency against the P. falciparum infection. 9,12 More importantly, artemisinin resistance in P. falciparum is an unprecedented spread occurring in South-East Asia, eastern India and Central African countries. 13,14 It poses a higher risk to the control and eradication effort. Notably, the licensed RTS,S vaccine, Mosquirix, was introduced into clinical trials among the children aged 5-17 months with the efficacy of $35%. 15 Altogether, more efficacious therapeutics are in demand to overcome the issue of clinical resistance to the frontline drugs as well as develop new molecules that exhibit multistage activity.
Current therapeutics are selectively efficacious during the asexual blood-stage infection and fail to inhibit the sexual and liver stages. This eventually results in the development of clinical drug resistance in other life cycle stages of parasites. 16,17 The renewed guidelines for malaria eradication emphasize on identifying new compounds with multistage antimalarial activity i.e., able to inhibit liver, asexual and sexual blood stages. 18 This approach is highly essential as one of the critical paths for the elimination of malaria. Therefore, active drug molecules whose potency is not restricted to a spectrum of pathogenic asexual blood stages 19 are needed to overcome the issue of increasing resistance. The clinically approved antimalarial therapeutics, namely atovaquone (ATQ) and pyronaridine, Fig. 1 (A) Structure of compounds 1 and 2; (B) cytotoxicity data for Alamar blue assay (IC 50 ) were transformed, X ¼ log (10) , normalized between 0% and 100% and fitted with sigmoidal dose-response curves; (C) anti-plasmodial activity of compounds against ART-resistant strain, data (relative SYBR green-1 fluorescence) were transformed, X ¼ log (10) , normalized between 0% and 100% and fitted with sigmoidal dose-response curves. DHAdihydroartemisinin, CQchloroquine. displayed signicant multistage activity to discover bioactive entities. 16,20 On a particular note, the discovery of drugs active against the liver stage and the sexual blood stage (gametocyte) was hampered partly due to the limited availability of established screening assays. 21,22 Precisely, the liver stage is a necessary progressive step that leads to the symptomatic asexual blood-stage, representing an attractive molecular drug-targets. Previously, we reported that hydroxyethylamine (HEA) analogs inhibited cultured blood-stage parasites at sub-micromolar range, 23 explicitly supporting the two compounds, 1 and 2 ( Fig. 1). Encouraged by these initial results, we intended to evaluate antiplasmodial potency against various life-stages of P. falciparum in culture condition and mice model of P. berghei infection to establish their multistage activity. Therefore, we assessed the efficacy of both the compounds against the blood stage of the artemisinin (ART) resistant eld isolate (in vitro), liver stage, ookinete, and gametocyte stages. We have also systematically performed in silico studies to gain insight into the Quantitative-Structure-Activity Relationship (QSAR) and identify the possible targets.

Assessment of blood stage antiplasmodial activity and cytotoxicity
Antimalarial drug discovery is now partly shied towards the therapeutics with multistage activity, probably due to a part of the malaria eradication agenda. 16,24 Encouraged by the blood stage antiplasmodial activity of compounds 23 (1 and 2; Fig. 1A), we decided to explore their potency against various life forms of P. falciparum. These two compounds were synthesized following the reported procedures. 23 As the rst part of the investigation, we evaluated the cytotoxicity of the compounds by exposing HEK293 (Kidney) and HepG2/Huh 7.1 (liver) cells to varying drug concentrations corresponding to exponentially higher than the blood-stage IC 50 values. This was to establish if the observed antiplasmodial effects of the compounds are due to non-specic inhibitory activity. As shown in Fig. 1B, both compounds were non-toxic at pharmacological concentration ranges. Testing these compounds on kidney and liver cells showed no observable toxic effects when used up to 100 times higher concentration required to kill 50% of blood-stage parasites, indicating the specic killing ability of the compounds.
Interestingly, compounds are required in high micromolar concentrations (>170 mM) to cause toxicity to mammalian cells, whereas the IC 50 for P. falciparum strains found to be in the submicromolar range ($1 mM). The compounds (Fig. 1) possess HEA scaffold 23,25 that is crucial to achieving the potent antiplasmodial activity. 24 Both compounds 1 and 2 displayed IC 50 values of 1.16 mM and 1.25 mM, respectively against P. falciparum (3D7) strain and $3 mM against P. falciparum (7GB) strain aer 48 h of treatment as reported earlier. 23 While testing against the ART-resistant P. falciparum [MRA-1241 (IPC4912, ARTresistant), BEI-MR4], 26 both compounds had substantial inhibitory activity (Fig. 1C). The IC 50 testing was performed with ART resistant strain to determine cross-resistance with compounds tested as shis in IC 50 s. The compounds markedly inhibited the growth of ART resistant strain with IC 50 , 785.23 AE 20.9 nM for compound 1 and 1097.82 AE 17.6 nM for compound 2 aer 72 h of treatment. Our results clearly showed that the compounds have anti-plasmodial activities against P. falciparum (CQ and ART sensitive and resistant strains) in culture. Drug resistance is a complex phenotype comprising all over changes in the mutant parasites such as overexpressing efflux pumps and various cellular components affecting a range of physiological characteristics. The consistent efficacy of compounds against resistant parasites shows no cross-resistance with the inherent resistance of CQ and ART resistant strains and against non-specic resistant phenotype. This also negates the possibility of observed efficacy to be reference strain 3D7 specic anomaly.
To further evaluate the in vivo efficacy of compounds, a mouse model of malaria infected with Plasmodium berghei NK-65 was developed. The results of these studies are presented in Fig. 2A and B. There was a dose-dependent suppression in blood-stage parasites with no signs of toxicity. It was observed that subcutaneous administration of four consecutive doses of compound 1 at 100 mg kg À1 signicantly reduced the bloodstage parasites on day 5 (P < 0.05), 8 (P < 0.01), and 15 (P < 0.01) as compared to the untreated control. We found that mice treated with 1 exhibited a signicant decline in the parasite load compared to 2 ( Fig. 2A). The mean survival time of the treated group was slightly higher than in untreated ones (Fig. 2B). Notably, compound 1 could not eliminate the blood-stage infection at both the dosages. We reasoned that this effect is probably due to the difference in the drug distribution and varied uptake towards infected RBC. Our ndings also indicate that partial clearance of blood parasites is due to shorter halflife (less than 4 h) and the rapid elimination of 1 from the circulation. Our results demonstrated that 1 exhibited inhibitory action on the growth of P. berghei infection in mice.
Assessment of inhibition potential of compounds during P. berghei liver-stage infection Based on the results from blood-stage inhibition assays and cytotoxicity assays, the inhibitory activity of compounds was tested against the liver stage P. berghei infection under in vitro and in vivo condition. Both compounds displayed an inhibitory effect against liver-stage infection. Since the parasite density is very low in the liver stage as compared to other developmental stages, this represents one of the potential targets for the antimalarial drugs that prohibit blood-stage infection and subsequent transmission to mosquitoes can also be quenched. [27][28][29] Therefore, HepG2 cells infected P. berghei sporozoites were treated at 1, 10 and 100 mM concentrations of both the compounds in culture. The parasite load was systematically quantied in control and the treated cells. The moderate efficacy (44% inhibition) was noted at the 1 mM concentration of compound 1. Compound 1 showed a signicant inhibition ($99%) upon treating cells with 10 mM. However, the efficacy was reduced when the cells were treated at 100 mM, possibly due Fig. 3 (A) Each scatters plot column represents different sets of samples from different time points and treatments. A 100 mL 3 rd , 4 th and 5 th -day gametocyte rich culture was seeded in microtiter plates in triplicates for each treatment (0.5 mM) and control and three wells each was fixed and stained by the Giemsa staining. Both compounds had an apparent reduction of intact gametocytes compared to control; (B) the most representative morphologies were documented. Both compounds had an evident detrimental effect on the normal morphology to the extent that none of the treated samples had any mature gametocytes even though with 50% viability in SYBR green screenings.
to compounds precipitation at higher concentration. The IC 50 value of compound 1 was 1.89 mM, whereas compound 2 did not show signicant activity at a lower concentration, with an IC 50 value of 30 mM.
We also tested the liver stage activity of both compounds in mice infected with P. berghei sporozoites through mosquito bites, which result in the maturation of merozoites in $55 hours. The infected mice were treated with a xed-dose, 50 mg kg À1 of compounds 1 and 2 for three days. Interestingly,  compound 1 showed 67% inhibition, as shown in Fig. 2C. For compound 2, the starting of the blood-stage parasite infection was monitored aer three days post-challenge by Giemsa staining. One day delay was observed at the onset of the prepatent period. Single-day delay in blood-stage infection represents $10-fold decrease in liver parasite burden. 24 Also, survival and parasitemia for the infected mice were monitored aer the treatment. Control mice showed a mean survival time (MST) of $10 days, whereas the treated mice displayed an improved MST of $16.5 days in the case of compound 2 (Fig. 2D). One day delay decreases the load of the parasite in the liver that could probably contribute to the survival benets. The parasitemia percentage remained low until the day 7 post-sporozoite challenge. The percent of parasitemia increased aer day 7 until the  death of all the animals. The results indicated that both the molecules exhibited comparable activity against P. berghei liver stage infection in mice models.

Assessment of plasma clearance of compound 1 in mice
It is evident from our ndings that compound 1 showed better efficacy than 2 on the parasite growth in P. berghei infection. Therefore, we investigated the half-life of compound 1 in plasma aer a single, subcutaneous dose (20 mg kg À1 ) at different time intervals using a validated HPLC method. We found that compound 1 showed the maximum plasma concentration (C max ) of 0.2 mg mL À1 , the area under the curve (AUC) was 0.5 mg mL À1 h À1 . These results indicated that compound 1 remained in the blood circulation with less than 4 h, which signicantly correlates with a therapeutic effect in suppressing the blood-stage and liver-stage infection in mice.
The developed compound retained in the blood circulation for a shorter time resulting in less exposure to the blood parasites causing the suppression in parasitemia.

Anti-gametocyte activity
Next, we evaluated the effect of both the compounds on the development of sexual erythrocytic stages of Pf3D7 cultures. Stage I and II gametocytes were incubated for 48 h with drug concentrations (0.5 mM) corresponding to the median IC 50 values of compounds at asexual stage assays. The counts of dead or deformed gametocytes, were recorded at stipulated time points. Treatment with the compounds resulted in a decreased proportion of stage III to V gametocytes (Fig. 3A), and morphologically distorted gametocytes compared to DMSO control (Fig. 3B). Enumeration of the counts showed survival of only 30-45% gametocytes aer 24 h of treatment; however, at later time points, there were no observable mature gametocytes (Fig. 3B). Aside from this, we observed that the untreated gametocytes displayed a distinct morphological appearance like ghost cells (remnant RBCs) due to the total consumption of hemoglobin by parasites, whereas the treated ones failed to show.

Activity against ookinete
Next, efficacy for both compounds was tested against the transmission stage of the parasite because widely available antimalarial agents selectively target blood-stage infection.
The drug molecules with the capability to intervene in the transmission stages of the parasites are considered as highvalue partners for the combination treatments due to their potential role in achieving malaria eradication. 30 Plasmodium transmission relies upon gametocytes capacity to convert into ookinetes, the prominent type of the parasite that attacks and develops a malarial infection in the mosquito midgut. By blocking this transmission stage, further malarial infection rates can be minimized to a greater extent. 31 We observed that both the compounds, 1 and 2 showed signicant inhibition against ookinete with the IC 50 values of 0.704 AE 0.054 mM and 0.7025 AE 0.004 mM, respectively. An increase in the concentration of the compounds reduced the parasites growth, leading to deformed parasites (Fig. 4).
Computational studies 3D QSAR and ADME calculations. Compounds were investigated for QSAR and ADME proling. The ADME results are illustrated in Table 1. Compound 1 showed optimal log P value i.e., 3.29. However, compound 2 found to be slightly above lipophilic/hydrophilic balance, i.e., 6.11. Compound 1 follows the Lipinski's rule of ve more stringently, whereas compound 2 showed one violation. 32 The oral absorption prognosis stated that both compounds have maximum tendency to absorb through the gastrointestinal tract. It is also noticeable that 1 tends to cross the blood-brain barrier. However, compound 2 lacks these characteristics. The metabolism study of the compounds on various sites of cytochrome P450 revealed that compound 1 has lower metabolism score at CYP450-3A4 and CYP450-2D6 (0.43 and 0.22 probability score, respectively), however, it does not have a metabolic tendency at the site of cytochrome CYP450 1A2, CYP450 2C19, and CYP450 2C9. Compound 2 is more likely to metabolize on cytochrome CYP450 3A4 (i.e. 0.91 probability score). However, this was reported a nearly negligible metabolic tendency on the site of cytochrome CYPP450 2D6, CYP450 1A2, and CYP450 2C9. At the same time low cytotoxicity proles with Huh-7 cells were observed, a liver carcinoma line with a xenobiotic degrading enzyme, negates the possibility of a toxic metabolite formed by these degradations. The overall ADME result showed that the compounds have substantial efficacy for further investigations.
The distribution of pharmacophore, anti-pharmacophore (an atom or group of atoms that reduce biological activity, providing, for example, steric hindrances), and ballast part in 1 and 2 is shown in the Fig. 5A and B. The result showed that the '-OH' group present in 1 have an anti-pharmacophore property and need to be replaced by an amide group. However, compound 2 doesn't show any pharmacophoric and antipharmacophoric part. The toxicity studies showed that the oxygen atom present in side chains 1 and 2 has a role in toxicity and needs to be replaced by an atom or group that reduces the toxicity of the compound, such as an amino group shown in Fig. 5C and D. Such a replacement, according to the 3D QSAR prediction, will lead to a decrease in the probability of toxicity by 0.28 and 0.32 for compounds 1 and 2, correspondingly. At large, 3D QSAR study indicated that oxygen atom is connected to the basic skeleton of the compounds act as an anti-pharmacophoric skeleton. To increase the anti-malarial activity of both the compounds, the hydroxyl group (bridged chain) should be replaced by suitable electrophile.

Proteinsinhibitors interaction analysis
Both compounds were docked into the active site of the studied proteins, and the resulted outcomes are shown in Table 2. The highest binding affinity of both compounds was observed for plasmepsin-II [PLM II], malarial aspartyl protease. Compound 2 showed the highest combination of free energy of binding and inhibition constant of À8.49 kcal mol À1 and 600.03 nM, respectively. On the other hand, compound 1 showed the interaction energy of À7.56 with an inhibition constant to be predicted around 10.55 mM. Apart from PLM II, both the compounds showed relatively higher binding energy for CDK as compared to the other proteins ( Table 2).
The analyses of the active site of the crystal structure of PLM (PDB ID -1LEE) showed that the residues Val78 and Ser79 are the ap residues while the Asp34 and Asp214 form the catalytic dyad, which is surrounded by Ser218 and Gly36 residues. 33 Compound 2 was observed to have interactions with Asp34, Asn39, Val78, Ser79, Ile123, Leu131 and Gly216 (Fig. 6A). Moreover, the interaction pocket of 1 residue contained Ile32, Asp34, Gly36, Tyr77, Val78, Ile123, Asp214, Thr217, Leu292 and Ile300 (Fig. 6B). These analyses indicated that 2 could inhibit the functions of the multiple protein targets with a relatively higher affinity than compound 1.
MD simulations (100 ns each) were performed on generated PLM II docked complexes with both the compounds for further conformational assessments and validations. The interaction patterns between the PLM II and compounds in the generated trajectories were analyzed in the form of the distance deviations from initial docking and hydrogen bonding dynamic patterns during MD simulations. The computed distances between the PLM II and the studied inhibitors were found to be uctuating between 0.15À0.25 nm, indicating that all the inhibitors were present in the interacting pocket throughout during MD simulations (Fig. 7). The different patterns of the hydrogen bonds (H-bonds) were observed for the studied system. Compound 1 showed the presence of up to ve hydrogen bonds with PLM (Fig. 7), while with compound 2, the protein only forms up to four hydrogen bonds during MD simulations. While observing the generated pattern of H-bonds, compound 1 even showed relatively higher repeat in the H-bonds followed by compound 2.
The compactness of the studied compounds was further assessed in terms of R g values, which showed that the docked complex of 1 showed a relatively lower R g value of around 1.95 nm (Fig. 8). Furthermore, both the systems showed comparatively similar stability as indicated from the RMSD values, which are uctuating between 0.25-4.2 nm (Fig. 8). Moreover, the MM/PBSA based algorithm was used to calculate the interaction energies between the PLM II and the studied inhibitors (Fig. 9). The energy analysis showed that the primary contributor for the total energy is electrostatic energy instead of Van der Waals forces in all the studied systems. The highest average total energy of the interaction of À1780.61 kJ mol À1 seen between PLM II and compound 1. While compound 2 showed similar patterns in which total energy was observed around À1474.49 kJ mol À1 . These outcomes provided a more indepth insight into the interaction of both compounds with PLM II and validated that they bind to the protein efficiently. Still, the interactions for compound 1 were more favorable as assessed from the MD simulation analyses.

Conclusion
Regardless of signicant advances in controlling and ultimately eradicating malaria, it remains a global health problem, particularly in endemic countries. Hydroxyethylamine based compounds were selected for their potency against drugresistant strains, including multistage activity in culture and in animal models. Both compounds showed potent blood-stage inhibitory activity against ART-resistant strain, anti-gametocyte action probably through the reduced capability to digest hemoglobin (due to inhibition of PLMs) and thereby failed to mature. The compounds were effective against the blood stage of the P. berghei infected mice and mainly compound 1 reduced the parasitemia load more efficiently. For better efficacy, compound 1 could be administered in conjunction with standard partner drugs with greater half-life and that acts synergistically. Our ndings also suggested that compound 1 is equally potent against blood-stage and liver-stage infection in mice and inhibited the progression of the gametocyte stage. The in silico studies suggested PLM II as a potential target for both compounds based on their high binding affinities. However, the preliminary screening resulted in partial inhibitory activity against PLM as reported earlier, 23 suggesting a possibility of other drug targets through which compound inhibits parasite growth. Notably, neither of these compounds were found toxic to human cells at therapeutic dosages. Taken together, our data provide new information on compound 1 as a multistage inhibitor against all life forms of the Plasmodium. This could be a potential choice for further chemical modications and the generation of derivatives with the improved biochemical and pharmacokinetic prole.

Experimental
Compounds were synthesized and puried following the previously reported procedures. 23  In vitro anti-plasmodial activity against ART-resistant P. falciparum Determination of blood stages and percent parasitemia was performed by Giemsa-stained thin smears from parasite cultures. Antimalarial inhibition activity (IC 50 ) was measured using the SYBR Green I assay, as reported by our group 24,25 early ring-stage parasites (>80% of rings) from ART resistant (IPC5202, MRA-1240, MR4) cultures were incubated in the presence of 1 : 2 serial dilution of each compound tested for 72 h, in concentrations ranging from 100 000-195.31 nM. Fluorescence intensity was determined using the SpectraMax® M5 Multi-Mode Microplate Reader, with excitation and emission wavelengths of 485 and 530 nm, respectively, and analyzed by nonlinear regression with Normalization & Baselinecorrections using soware GraphPad Prism to determine IC 50 values.

The in vivo antimalarial study
The mouse model of malaria infection was developed by intraperitoneal administration of 1 Â 10 7 Plasmodium berghei-NK65 infected RBC to healthy C57BL6 mice. The in vivo efficacy of compounds was performed in accordance to a slightly modied version of the 'Peters four day suppressive test'. 34 For experimental infections, each animal was injected with 1 Â 10 7 infected RBC through intraperitoneal injection. Aer 48 h of post-infection, mice were treated with 50 or 100 mg kg À1 body weight of compounds suspended in vehicle control (10% DMSO) for four consecutive days through subcutaneous route. The control group was treated with only vehicle control. The efficacy of the treatment was determined by measuring the parasitemia and survival of mice on days 5, 8, and 15 by making thin smears of blood, from the tail vein of mice, stained with 10% Giemsa and examined at 1000Â magnication. The percent parasitemia was measured by counting infected and non-infected RBCs from randomly selected optical elds and expressed as the number of infected RBCs per 100 RBCs. The survival rate of mice (mean survival time, MST) was determined aer treatment until 30 days post-infection. The reduction in the number of parasitemia was taken as the index for the therapeutic activity of the test compounds. The percentage of parasitemia was calculated manually using Cell Counting Aid soware 35 using the formula %parasitemia ¼ total number of parasitized RBCs total number of RBCs Â 100 Pharmacokinetic study A pharmacokinetic study was performed for compound 1 using Swiss albino mice weighing 25 to 30 g. The experimental animals were separated into two groups consisting of three animals each. The mice of the rst group were given subcutaneous administration of 1 (20 mg kg À1 ) suspended in DMSO : PEG-400 : 5% glucose (2 : 10 : 88) (vehicle control) while the second group was administrated with vehicle control (without test compound). At different time intervals (0.25, 0.5, 1, 2, 4, 8, and 24 h), a volume of 0.2 mL of blood was withdrawn from retro-orbital plexus and collected in tubes containing heparin and subjected to centrifugation at 5000 rpm for 15 min for separating the plasma. The concentrations in plasma were determined using Gilson (USA) HPLC instrument. Various pharmacokinetic parameters were determined by plotting drug concentration versus time.

In vitro liver-stage infection in HepG2 cells
Salivary glands of Anopheles stephensi mosquitoes infected with GFP-expressing P. berghei ANKA (Insectaries at National Institute of Immunology, New Delhi, India) were dissected out. Sporozoites were perched in inoculation medium (DMEM with 10% FBS, 3% P/S), and allowed to infect HepG2 cells. The cells were treated with compounds before the infection at various concentrations. The plate was centrifuged at 3000 rpm for 5 minutes to ensure the sporozoites settled down in the plate and incubated at 37 C for 3 h. The medium was replaced with pre-warmed medium (DMEM, 10% FBS, 3% P/S, and 0.1% fungizone). Aer two days, the cells were harvested in trizol reagent and freezed at À80 C. RNA isolation followed by quantitative RT-PCR was performed to estimate the parasite load in the infected cells.
Liver-stage infection in mice, treatment and survival assay GFP-expressing P. berghei ANKA was maintained in Anopheles stephensi mosquitoes (Insectaries at the National Institute of Immunology (New Delhi, India)). Sporozoite count per salivary gland and infectivity was measured before the experimental infection of mice. Briey, anesthetized mice were exposed to a cage of starved mosquitoes having sporozoites (100/5 mouse) for 15 minutes with intermittent disturbing of mosquito every 2-3 minutes to encourage the biting process. Biting was done at 22 + 1 C and further maintained in dark conditions in a controlled condition. The rst dose (50 mg kg À1 ) was administered 24 h before the infection, second dose was given 2 h post sporozoite challenge (biting) through the intraperitoneal route. For calculating parasite load, the liver of experimental mice was isolated 50 h post-infection. The mice were anesthetized, sterilized the abdominal area with 70% ethanol and dissected the liver and homogenized it in a precold denaturating solution. RNAs were isolated from the homogenized sample followed by real-time PCR. Survival rate was monitored on a daily basis with a subsequent monitoring on blood stage infection on day 3 onwards. Blood smears were obtained from the tail vein of mice and stained with 10% Giemsa and the parasites were counted using a 100Â oil objective.

Cytotoxicity assays in human liver and kidney cells in culture
Cytotoxicity of test compounds was evaluated by assessing the cell viability of Huh 7.1 and HEK293 cells using the Alamar Blue assay. 36 Cells (1 Â 10 3 per well) were exposed to different concentrations of test compounds for 24 h in 96-well plate format along with 10/100 mL Alamar Blue solution (10X ¼ 4.8 mM in phosphate buffer) per well. Cell viability was measured using the SpectraMax® M5 Multi-Mode Microplate Reader, at 530 nm excitation wavelength and 590 nm emission wavelength for metabolically reduced resazurin dye. The percent cell viability were expressed as CC 50 at which 50% of the cells died/metabolically slowed down due to drug exposure compared to the non-treated and vehicular control. The average minima and maxima were added for normalizations with negative control and blanked respectively to retain their relative global trend. The toxicities observed in part were contributed by DMSO vehicle's cytostatic effect as seen in vehicular controls, although the volume was 1 : 100 (Vehicle to Media).
In vitro gametocyte culture of P. falciparum Gametocyte cultures were initiated as described previously. 24 Briey, on day 0, cultures were synchronized at the ring stage by lysis of the ask pellet with 5 volumes of 5% sorbitol for 10 min at 37 C and maintained at 0.2% parasitemia and 4% hematocrit. The culture was grown in RPMI 1640 medium (Gibco) supplemented with hypoxanthine (Sigma-Aldrich), sodium bicarbonate and 1.5% (wt/vol) AlbuMAX II. Cultures with >8% gametocytes (on day 12 normally) were treated with 50 mM N-acetyl-D-glucosamine (GLcNA, Sigma-Aldrich) and 50 ng mL À1 of bistratene A for 3-5 days to remove asexual forms. The percentage of asexual forms and gametocytes were calculated on Giemsa-stained smears for performing anti-gametocyte assays.

In vitro gametocyte parasite growth inhibition assay
Aliquots 100 mL of cultures containing gametocytes at the indicated stages (preferably stage I) were transferred to 96-well plate in duplicates and incubated at 37 C in the presence of 5% CO 2 . The culture medium was changed aer 48 h and the appropriate concentration of test compounds was maintained. Blood smears of parasites were prepared and stained with Giemsa and the parasitemia was monitored using oil immersion microscopic analysis at 1000Â magnication.

Activity against ookinete
To study whether the test compound blocks the transmission stage, the sexual stage of the malaria parasite (PbGFPcon) was used to check its ability to inhibit ookinete development in vitro. For this, 1.2 mg phenylhydrazine (Sigma-Aldrich, St. Louis, USA) in 0.9% NaCl was intraperitoneally (i.p.) injected into BALB/c mice 3 days before P. berghei infection. The treated mice were then inoculated with 2 Â 10 6 P. berghei-infected erythrocytes intraperitoneally to initiate the blood-stage infection. Parasitemia was allowed to reach $5-10% at three days postinfection, and the gametocytemia was observed to be 1-2%. The mice were anesthetized and the blood was withdrawn in heparin containing tubes and added to the ookinete culture medium. The parasites were grown in the medium with a blood/ medium volume ratio of 1 : 10 at 19 C for 24 h. Aer 24 h, the culture was Giemsa-stained and analyzed for ookinetes morphology. The number of ookinetes of both normal and abnormal morphology was calculated.

Theoretical calculations
3D QSAR studies were performed via the "Cinderella Shoe" 37 methods and within the CoMIn algorithm. 38 The pharmacophore part (in red), antipharmacophore skeleton (blue), and inactive parts (grey color) of the synthesized compounds (1 and 2) are displayed in Fig. 5A and B. The grey part of the ring does not participate in pharmacological volubility but liable for structure conguration and matrix function to provide the necessary mutual arrangement of pharmacophoric fragments. However, the toxicity prognosis was shown in Fig. 5C and D (toxicophoric red and non-toxicophoric blue). The ADME prole like log P, toxicity, metabolism at CYP450:3A4, CYP450:2D6, and predicted anti-inammatory activity was carried out using the website www.chemosophia.com. 39 However, TPSA, Lipinski violation, H donor, H acceptor, water solubility score, GI absorption score, and the possibility of metabolism at CYP1A2, CYP2C19, and CYP2C9 measured by http://www.swissadme.ch and Cresset soware. 40,41 Target selection and molecular modelling To identify the possible receptors for the designed inhibitors, a target based search was performed in the DrugBank, 42 which resulted in the identication of 10 potential protein molecules in P. falciparum. These are listed as "approved" drug targets in DrugBank and their structural information was obtained from the PDB and Uniprot databases. Moreover, the structures of calcium-dependent protein kinase 1 (CDK), cytochrome B (CYB), dihydropteroate synthetase (DPS) and malaria protein EXP-1 (MPEXP) were predicted by satisfying the spatial restraints using homology modelling. The available crystal structures of the proteins such as bifunctional dihydrofolate reductase thymidylate synthase (DRTS), dihydroorotate dehydrogenase quinone mitochondrial (DOD), ferredoxin NADP reductase apicoplast (FR), glutathione reductase (GR), glutathione-S-transferase (GST), plasmepsin II (PLM II), were further optimized using the "Prepare protein" module of Schrodinger suite. The structures of studied inhibitors were generated using the "2-D sketcher" of Schrodinger and then converted to three dimensional (3-D) coordinates, which were further optimized utilizing the "LigPrep" module.

Molecular docking studies
Consequently, the generated optimized structures of the respective proteins and ligands were subjected to molecular docking based interaction analyses using AutoDock 4 package. 43 The available protocol resulted in the generation of bound protein-ligand conformations, based on free energy calculations, which were performed based on of the empirical force eld as well as the Lamarckian genetic algorithms. 44 Primarily, the grid boxes of varied dimensions along the XYZ directions were created using the AutoGrid module with a spacing of 0.375Å. The efficiency of the predictions was further enhanced by setting the docking parameters coupled with the Lamarckian genetic algorithm to their maximum values, which involved setting the number of individuals in the population to 250 as well as the maximum number of energy assessments to a "longer" interval. The 100 bound conformations were generated for each protein and ligand system (Total 60 systems) which were further grouped according to the RMSD tolerance of 2.0Å. The free energy of binding was used as a parameter for assessing the efficiency of interactions, was calculated on the basis of the given equation Here, DG vdW symbolizes the van der Waals, DG H-bond represents hydrogen bonding while DG elec and DG sol stand for the free energy of electrostatics and solvation, respectively. The rescoring of the generated docked complexes was performed using DrugScoreX server, 45 which increases the efficiency of the predictions by incorporating knowledgebase scoring functions.

Molecular dynamics simulations
The docked complexes of PLM were subjected to MD simulations using GROMACS 46 (version -2018-2) for the validation of the interaction studies. Primarily, the topologies of the protein structure in different complexes were generated using the GROMOS96-53a6 force eld 46 and for the ligands topology generation, the PRODRG server 47 was used, as GROMACS does not possess the topology generation ability for natural ligands. The PRODRG server lacks the functionality of generating the partial charges, therefore, the correction was performed using the DFT method implemented in GAUSSIAN which utilized the B3LYP6-31G(d,p) basis set and the CHELPG program. 48 Aer successful topology generation, all the docked complexes were solvated by immersion in the SPC/E water model 49 and in order to generate the stable systems, neutralization was performed by adding the counterions. The Na + and Cl À ions were added to neutralize the systems. Consequently, the solvated and neutralized systems were subjected to energy minimization step by using combined steepest descent and conjugate gradient algorithms, with a convergence criterion of 0.005 kcal mol À1 . The position restraints were applied to the structure of the minimized system ligands before the equilibration phase. The equilibration step was performed in a combination of NVT (constant volume) and NPT (constant pressure) ensemble conditions, each at 100 ps time scale. The temperature of 300 K was maintained for the system using Berendsen weak coupling method and pressure of 1 bar was maintained utilizing Parrinello-Rahman barostat in the equilibration stage. The nal conformational production stage was carried out using the LINCS algorithm for 100 ns timescale and trajectories were generated, which were analyzed in order to understand the behavior of each complex in the explicit water environment. The changes in the H-bonds, as well as distances between the PLM II and designed inhibitors in the complex systems, were analyzed. Furthermore, the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) protocols implemented in g_mmpbsa package 50 were used for the calculation of free energy of binding, which provides deeper insights into the interaction mechanisms of the protein and the ligand molecules.

Author contributions
BR, PK, APS and Poonam jointly conceptualized the idea, analyzed data and executed the study; NS, MB and SS performed synthesis and characterization; VR tested in vivo blood stage efficacy; YG and RM performed cytotoxicity and antiplasmodial efficacy against resistant strain and gametocyte activity; PP and MS performed the in silico analysis. NS, JS, MK undertook testing experiments on liver stage (in vitro and in vivo), and ookinete stage. All the authors contributed to writing, independent reviewing, and approved the manuscript.