SCIENTIFIC REPORTS

<figure><img src="image_1.png" /></figure>

OPEN

# Biochemical analysis of leishmanial and human GDP-Mannose Pyrophosphorylases and selection of inhibitors as new leads

Received: 10 October 2016

Accepted: 16 March 2017

Published online: 07 April 2017

Wei Mao<fnref n="1" />, Pierre Daligaux<fnref n="1" />, Noureddine Lazar², Tâp Ha-Duong³, Christian Cavé<fnref n="1" />, Herman van Tilbeurgh², Philippe M. Loiseau<fnref n="1" /> & Sébastien Pomel<fnref n="1" />

<figure><img src="image_2.png" /></figure>

Leishmaniases are an ensemble of diseases caused by the protozoan parasite of the genus *Leishmania*. Current antileishmanial treatments are limited and present main issues of toxicity and drug resistance emergence. Therefore, the generation of new inhibitors specifically directed against a leishmanial target is an attractive strategy to expand the chemotherapeutic arsenal. GDP-Mannose Pyrophosphorylase (GDP-MP) is a prominent therapeutic target involved in host-parasite recognition which has been described to be essential for parasite survival. In this work, we produced and purified GDP-MPs from *L. mexicana* (*LmGDP-MP*), *L. donovani* (*LdGDP-MP*), and human (*hGDP-MP*), and compared their enzymatic properties. From a rationale design of 100 potential inhibitors, four compounds were identified having a promising and specific inhibitory effect on parasite GDP-MP and antileishmanial activities, one of them exhibits a competitive inhibition on *LdGDP-MP* and belongs to the 2-substituted quinoline series.

Leishmaniases are Neglected Tropical Diseases (NTD) caused by protozoan parasites belonging to the genus *Leishmania*, and transmitted by the female sandfly. Currently, this NTD threatens 310 million people in 98 countries with 1.3 million new cases per year<fnref n="1" />. There are three main clinical manifestations of leishmaniasis: cutaneous (CL), which causes skin sores and chronic ulcers, muco-cutaneous (MCL), affecting naso-oro-pharyngeal mucosa and causing facial disfigurement, and visceral (VL) which is lethal without treatment. *L. mexicana* and *L. donovani* are two leishmanial species distributed in the New and Old World and responsible of CL and VL, respectively.

The existing treatments suggested by World Health Organization are mainly based on the use of antimonials, liposomal amphotericin B, miltefosine and paromomycin². All of them having main issues of toxicity and drug resistance emergence³, the development of new inhibitors specifically directed against a parasite target is an attractive approach for identifying new drugs. The GDP-mannose Pyrophosphorylase (GDP-MP) is an attractive therapeutic target due to its essential role in amastigote survival in macrophages both *in vitro* and *in vivo* and its involvement in the biosynthesis of diverse glycoconjugates necessary for host cell recognition⁴, ⁵. Furthermore, by molecular modeling, we identified several structural differences between the parasite and the human orthologs, including the presence of a specific motif in the catalytic pocket, making this enzyme a target of choice for the development of new specific antileishmanial agents⁶, ⁷. A high throughput screening performed on *L. major* GDP-MP allowed the identification of a quinoline derivative, presenting an IC₅₀ value on the enzyme at the sub-micromolar range and on intramacrophage parasites at about 20 µM⁸. Various quinoline derivatives have been identified exhibiting antileishmanial activity⁹. Our laboratory has previously revealed the 2-substituted quinoline series as antileishmanial lead¹⁰, ¹¹. Thus, in the present study, we decided to integrate this series in the design of GDP-MP competitive inhibitors. Preliminary molecular modeling studies allowed us to hypothesize that the

<footnote n="1">¹Chimiothérapie Antiparasitaire, UMR CNRS 8076 BioCIS, Université Paris-Sud, Université Paris-Saclay, 5 rue Jean-Baptiste Clément, 92290, Châtenay-Malabry, France. ²Institut de Biologie Intégrative de la Cellule, UMR 9198, CNRS, Université Paris-Sud, Université Paris-Saclay, Bâtiment 430, Orsay, 91405, France. ³Molécules Fluorées et Chimie Médicinale, UMR CNRS 8076 BioCIS, Université Paris-Sud, Université Paris-Saclay, 5 rue Jean-Baptiste Clément, 92290, Châtenay-Malabry, France. Correspondence and requests for materials should be addressed to S.P. (email: sebastien.pomel@u-psud.fr)</footnote>
quinoline motif could replace the guanine group of GDP-mannose within the GDP-MP catalytic site. Therefore, GDP-MP competitive inhibitors could be designed by including such quinoline group in the inhibitor scaffold. Various pharmacomodulations were also carried out without quinoline, supplying an in-house library of 100 compounds that have been studied in the present work.

Starting from mannose-1-phosphate (Man-1-P) and GTP, GDP-MP catalyzes the formation of GDP-mannose. This activated form of mannose is a key substrate for different glycosylation processes such as N-glycosylation or O-mannosylation which are essential for post-translational modifications in eukaryotes¹². In *Leishmania* sp., GDP-mannose is a precursor for the synthesis of a range of mannose containing glycoconjugates involved in host cell recognition¹³. Many GDP-MPs appear to be dimeric, like in bacteria and in some protozoa such as *Trypanosoma brucei*, and in some bacterial species, the protein presents both GDP-MP and phosphomannose isomerase activities, leading to bifunctional enzymes¹⁴⁻¹⁷. In mammalian, GDP-MP is composed of two subunits α and β, the latter containing the GDP-MP catalytic activity *per se*, while the former would have a regulatory role allowing allosteric inhibition of the β subunit by GDP-mannose¹⁸⁻²⁰. In human, the gene encoding for GDP-MP β subunit is transcribed as two isoforms β1 and β2. In opposition to β1, β2 isoform is strongly expressed within a wide tissue distribution. Unlike in other organisms, in *L. mexicana*, GDP-MP is a hexamer which can dissociate in trimers and monomers in low ionic strength conditions leading to a mixture of the three forms²¹. Although the GDP-MP activity was detected in the size range of a hexamer (240 kDa), the trimer and monomer forms were not described to be inactive.

In this study, we compared enzymatic parameters of GDP-MPs from two geographically distant leishmanial species, *L. donovani* (*Ld*GDP-MP), found in Africa and Asia, and *L. mexicana* (*Lm*GDP-MP), encountered in South America, with those from the human ortholog (*h*GDP-MP), with the objective to select specific inhibitors of leishmanial GDP-MPs. These compounds rationally designed from docking analysis were evaluated on the purified enzymes as well as on both parasite species *in vitro*.

## Results

**Purification of GDP-MPs.** Although *L. mexicana* and *L. donovani* are geographically distant parasite species, *Lm*GDP-MP and *Ld*GDP-MP share 94% identity in sequence. The human GDP-MP isoform analyzed in this work is the β2 subunit (*h*GDP-MP), which presents a wide tissue distribution, and has a higher identity score for the leishmanial counterparts (~50%), compared to the β1 subunit (~46%). While amino acids involved in substrate binding are conserved between *Lm*GDP-MP and *Ld*GDP-MP, several differences have been identified in the active site of *h*GDP-MP⁶,⁷. In order to compare their enzymatic properties, His6-tagged *Ld*GDP-MP, *Lm*GDP-MP and *h*GDP-MP were produced and purified using nickel affinity chromatography (Fig. 1a,d and h). *Ld*GDP-MP was further purified by size exclusion chromatography (Fig. 1b). The profile of gel filtration shows one main peak at 50 mL and three smaller peaks at 67, 76 and 86 mL. As revealed in SDS-PAGE, the fractions from peaks 1 to 4 showed the presence of a markedly major band at 42 kDa, corresponding to the predicted molecular weight of *Ld*GDP-MP and showing 63% sequence coverage with *Ld*GDP-MP in mass spectrometry (MS; Fig. 1c, Supplementary Fig. S1). The estimated molecular weight of the peak 1 at around 2,000 kDa (Fig. 1b and Supplementary Table S1) indicates that *Ld*GDP-MP is an aggregated soluble form in these fractions. Major and comparable GDP-MP activities were measured in peaks 2, 3 and 4 whose elution volumes correspond to hexameric (240 kDa), trimeric (120 kDa) and monomeric (40 kDa) forms of the enzyme, respectively. The fractions of peaks 2, 3 and 4 were further pooled and concentrated for Dynamic Light Scattering analysis (Table S1). The size of the objects detected was estimated at 240 kDa with a low polydispersity, showing that *Ld*GDP-MP is mainly in the hexameric form in this concentrated pool of fractions.

Following purification on a nickel affinity column, His6-MBP tagged *Lm*GDP-MP was further purified by anion exchange chromatography. The purified fusion protein of 84 kDa was digested by TEV to remove the tag and loaded on a gel filtration column. The elution profile shows two peaks containing either two bands of 41 and 75 kDa or one main band of 43 kDa, as revealed by SDS-PAGE (Fig. 1e and f). The 43 kDa band corresponding to the predicted molecular weight of His6-MBP tag, and the MS analysis of the 41 kDa band revealing 41% sequence coverage with *Lm*GDP-MP (Supplementary Fig. S1), gel filtration allowed to separate efficiently *Lm*GDP-MP from His6-MBP tag. Thereafter, the contaminant at 75 kDa was mostly removed by anion exchange chromatography, as revealed in Fig. 1g. As opposed to *Ld*GDP-MP, *Lm*GDP-MP was eluted in a single gel filtration peak (peak 1 in Fig. 1e), with an estimated molecular weight of 240 kDa corresponding with a hexameric form of the protein.

Three steps of purification were used to purify *h*GDP-MP: nickel affinity, anion exchange and gel filtration chromatographies. The elution profile of the last step shows that *h*GDP-MP is eluted in one main peak corresponding to the dimeric form of *h*GDP-MP (Fig. 1i). The analysis of the fractions of this peak in SDS-PAGE revealed a single band at 39 kDa which showed 63% sequence coverage with *h*GDP-MP in MS (Fig. 1j and Supplementary Fig. S1).

Fractions corresponding to purified *Ld*GDP-MP (peaks 2, 3 and 4 in Fig. 1c), *Lm*GDP-MP (Fig. 1g) and *h*GDP-MP (Fig. 1j) were pooled and concentrated for further enzymatic studies.

**Enzymatic properties.** As shown in Fig. 2a and Supplementary Fig. S3, the best reaction conditions were obtained at 37 °C with a 20 min reaction time for the three purified GDP-MPs. From Michaelis-Menten plots, the optimal concentrations of Man-1-P and GTP have been determined at 100 µM for the three GDP-MPs, although enzyme saturation was already reached at 50 µM of both substrates (Supplementary Fig. S5), and a significant inhibition was observed with a concentration higher than 150 µM of GTP, but not Man-1-P (data not shown). Enzyme concentration was reduced by a factor 8 (0.25 ng/µl) for *Lm*GDP-MP compared to the two other GDP-MPs in order to stay in the linear regime of the detection assay (Supplementary Fig. S4). The pH optima for *Ld*GDP-MP and *Lm*GDP-MP are similar (pH 7.5 and 7.6), and slightly more acidic for *h*GDP-MP (pH 6.5; Fig. 2b). In these optimal conditions, the activities of *Ld*GDP-MP and *h*GDP-MP are only 15–20% relative to that
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Purification of *LdGDP-MP*, *LmGDP-MP* and *hGDP-MP* (*a*,*d*,*h*) Coomassie blue stained SDS-PAGE of *LdGDP-MP* (*a*), *LmGDP-MP* (*d*), and *hGDP-MP* (*h*) purified from Ni-NTA column. Lane TOT: total, Lane SN: supernatant, Lane P: Pellet, Lane FT: flow through, Lane W: wash, Lanes E1, E2 and E4: elution fraction at 100 µM, 200 µM and 400 µM imidazole, respectively. The arrowheads show the proteins of interest. The protein ladders on the left are indicated in kDa. Full length gels are presented in Supplementary Fig. S7. (*b*,*e*,*i*) Size exclusion chromatography profiles of *LdGDP-MP* (*b*), *LmGDP-MP* (*e*), and *hGDP-MP* (*i*). In the inset, the filled triangles and open circles represent standard proteins and the peaks of proteins of interest, respectively. (*c*,*f*,*j*) Coomassie blue stained SDS-PAGE of fraction volumes of peaks 1 (49–51 mL), 2 (66–68 mL), 3 (75–78 mL) and 4 (84–86 mL) of *LdGDP-MP* (*c*), peaks 1 (12–13 mL) and 2 (16 mL) of *LmGDP-MP* (*f*) and peak 1 (77–82 mL) of *hGDP-MP* (*j*). Lo is the loading control. The arrowheads show the proteins of interest. The protein ladder on the left is indicated in kDa. Full length gels are presented in Supplementary Fig. S7. (*g*) Coomassie blue stained SDS-PAGE of a concentrated pool of fractions containing *LmGDP-MP* eluted from anion exchange chromatography (AEC). The arrowhead shows the protein of interest. The protein ladder on the left is indicated in kDa. Full length gel is presented in Supplementary Fig. S7.</figcaption></figure>
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

d.

<table><thead><tr><th></th><th>[GDP-MP]<br>(ng/µl)</th><th>[GTP]<br>(µM)</th><th>[Man-1-P]<br>(µM)</th><th>Temperature<br>(°C)</th><th>pH</th><th>[Mg<sup>2+</sup>]<br>(mM)</th></tr></thead><tbody><tr><td>LdGDP-MP</td><td>2</td><td>100</td><td>100</td><td>37</td><td>7.6</td><td>5</td></tr><tr><td>LmGDP-MP</td><td>0,25</td><td>100</td><td>100</td><td>37</td><td>7.5</td><td>2.5</td></tr><tr><td>hGDP-MP</td><td>2</td><td>100</td><td>100</td><td>37</td><td>6.5</td><td>1</td></tr></tbody></table>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

**Figure 2.** Determination of *LdGDP-MP*, *LmGDP-MP* and *hGDP-MP* optimal conditions of enzymatic reaction and specificity. (a,b,c) Determination of *LdGDP-MP*, *LmGDP-MP* and *hGDP-MP* optimal temperatures (a), pH (b) and Mg²⁺ concentrations (c). For the determination of the optimal pH, the buffers used are 50 mM MES for pH 5.6 and 6.5, and 50 mM Tris-HCl for pH between 6.8 and 10. In each graph, enzyme activities were calculated relative to the maximal specific activity obtained with *LmGDP-MP* and are expressed in percent. The results correspond to the mean of three independent experiments ± SD. (d) Optimal parameters of *LdGDP-MP*, *LmGDP-MP* and *hGDP-MP* enzymatic reactions. All reactions were performed with the optimal reaction time of 20 min. (e,f,g) Determination of *LdGDP-MP*, *LmGDP-MP* and *hGDP-MP* specificities for cofactors (e), nucleotide triphosphates (f) or sugar monophosphates (g). All cofactors (Mg²⁺, Mn²⁺, Ca²⁺, Ni²⁺, Cu²⁺, Zn²⁺) were used at 5 mM. The NTPs (GTP, ATP, CTP, TTP, UTP and ITP) and sugar monophosphates (Man-1-P, Glu-1-P and Gal-1-P) were used at 100 µM. Enzyme activities were calculated relative to each enzyme specific activity obtained with Mg²⁺ (e), GTP (f) or Man-1-P (g), and are expressed in percent. The results correspond to the mean of three independent experiments ± SD.

of *LmGDP-MP* (Fig. 2a-c). In absence of Mg²⁺, no enzyme activity was observed, showing this cofactor is necessary for the activity of the three purified GDP-MPs (Fig. 2e). Optimal concentrations of Mg²⁺ were determined at 1 mM, 2.5 mM and 5 mM for *hGDP-MP*, *LmGDP-MP* and *LdGDP-MP*, respectively (Fig. 2c). No GDP-MP activity was observed with any other divalent ions, such as Mn²⁺, Ca²⁺, Ni²⁺, Cu²⁺ and Zn²⁺, except for Mn²⁺ which can be used as a cofactor by *LdGDP-MP* to generate 33% of activity compared to Mg²⁺ (Fig. 2e). Likewise, no GDP-MP activity was observed with substrate analogues, except for ITP which can be used by *LmGDP-MP* or *hGDP-MP* to generate 22-24% of activity compared to GTP (Fig. 2f and g). These results show that the three purified GDP-MPs present high cofactor and substrate specificities. All further enzymatic analyses were performed using the optimal conditions described in Fig. 2d.

Kinetic constants of each purified GDP-MP for both substrates, GTP and Man-1-P, were determined by using non-linear regression analysis (Fig. 3g) and graphical representations were depicted using Lineweaver-Burk double reciprocal plots (Fig. 3a-f). The maximal velocities ($V_m$) of *LdGDP-MP* and *hGDP-MP* have comparable
<figure><img src="image_10.png" /><figcaption>b.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>c.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>c.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>e.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>f.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>g.</figcaption></figure>

<table><thead><tr><th></th><th>Substrate</th><th>V<sub>m</sub> (µM/min) ± SD</th><th>K<sub>m</sub> (µM) ± SD</th><th>k<sub>cat</sub> (min<sup>-1</sup>) ± SD</th><th>k<sub>cat</sub>/K<sub>m</sub> (min<sup>-1</sup>µM<sup>-1</sup>) ± SD</th></tr></thead><tbody><tr><td rowspan="2">LdGDP-MP</td><td>Man-1-P</td><td>0.30 ± 0.06</td><td>9.98 ± 6.15</td><td>6.23 ± 1.24</td><td>0.75 ± 0.33</td></tr><tr><td>GTP</td><td>0.39 ± 0.03</td><td>16.84 ± 3.29</td><td>8.05 ± 0.71</td><td>0.48 ± 0.05</td></tr><tr><td rowspan="2">LmGDP-MP</td><td>Man-1-P</td><td>1.05 ± 0.13</td><td>11.58 ± 1.76</td><td>87.31 ± 10.59</td><td>7.56 ± 0.22</td></tr><tr><td>GTP</td><td>0.95 ± 0.17</td><td>7.67 ± 2.28</td><td>78.92 ± 13.86</td><td>10.52 ± 1.14</td></tr><tr><td rowspan="2">hGDP-MP</td><td>Man-1-P</td><td>0.47 ± 0.01</td><td>6.98 ± 0.64</td><td>9.30 ± 0.19</td><td>1.34 ± 0.11</td></tr><tr><td>GTP</td><td>0.44 ± 0.01</td><td>7.08 ± 0.57</td><td>8.74 ± 0.14</td><td>1.24 ± 0.10</td></tr></tbody></table>

**Figure 3.** Determination of *LdGDP-MP*, *LmGDP-MP* and *hGDP-MP* kinetic constants. (a,c,e) Lineweaver Burk double reciprocal plots 1/V = f(1/[Man-1-P]) of *LdGDP-MP* (a), *LmGDP-MP* (c) and *hGDP-MP* (e). GTP concentration was held constant at 150 µM for *LdGDP-MP* and 80 µM for *LmGDP-MP* and *hGDP-MP*. The results expressed correspond to the mean of three independent experiments ± SD. (b,d,f) Lineweaver Burk double reciprocal plots 1/V = f(1/[GTP]) of *LdGDP-MP* (b), *LmGDP-MP* (d) and *hGDP-MP* (f). Man-1-P concentration was held constant at 150 µM for *LdGDP-MP* and 80 µM for *LmGDP-MP* and *hGDP-MP*. The results expressed correspond to the mean of three independent experiments ± SD. (g) Kinetic constants (*V*<sub>m</sub>, *K*<sub>m</sub>, *k*<sub>cat</sub>, *k*<sub>cat</sub>/*K*<sub>m</sub>) of *LdGDP-MP*, *LmGDP-MP* and *hGDP-MP* for both substrates Man-1-P and GTP.

values for GTP and Man-1-P ranging from 0.30 to 0.47 µM/min, while the *LmGDP-MP* *V*<sub>m</sub> value is two to three times higher for both substrates. Accordingly, *LmGDP-MP* exhibits a turnover rate (*k*<sub>cat</sub>) 10 times higher than the two other purified enzymes, with a value at around 80 min<sup>-1</sup> for both substrates. The *K*<sub>m</sub> values of *hGDP-MP* is at around 7 µM for both substrates and is comparable or slightly lower than the *K*<sub>m</sub> values obtained with leishmanial GDP-MPs, reflecting a moderate higher affinity of *hGDP-MP* for both substrates compared to leishmanial GDP-MPs. The *K*<sub>m</sub> values of *LmGDP-MP* are 11.58 ± 1.76 µM for Man-1-P and 7.67 ± 2.28 µM for GTP
suggesting it has a higher affinity for GTP than for Man-1-P. This result is in opposition to LdGDP-MP which presents a smaller K<sub>m</sub> for Man-1-P compared to GTP. The calculated catalytic efficiencies of LmGDP-MP (k<sub>cat</sub>/K<sub>m</sub>) at 7.56 ± 0.22 min<sup>-1</sup> µM<sup>-1</sup> and 10.52 ± 1.14 min<sup>-1</sup> µM<sup>-1</sup> are also 10 to 20 times higher than for the two other GDP-MPs. In addition, the catalytic efficiency of hGDP-MP is similar for Man-1-P and GTP at around 1.3 min<sup>-1</sup> µM<sup>-1</sup>, and is twice higher than the k<sub>cat</sub>/K<sub>m</sub> calculated for LdGDP-MP for both substrates.

In order to investigate the mechanism of reaction of each purified GDP-MPs and to assess the kinetics of substrates binding, enzyme activities were analyzed with a range of Man-1-P concentrations and a fixed GTP concentration (Fig. 4a,c and e), and inversely with a range of GTP concentrations and a fixed Man-1-P concentration (Fig. 4b,d and f). In these two series of enzyme assays the double inverse plots intersect on the left of the 1/y axis for the three GDP-MPs (Fig. 4), which is consistent with the formation of a ternary complex, and thus a sequential mechanism²². Furthermore, for both substrates, the intersection of double inverse plots is above the 1/x axis for LdGDP-MP and hGDP-MP, and on the 1/x axis for LmGDP-MP (Fig. 4a,b,e and f). These analyses support a sequential random mechanism for the three enzymes. Moreover, for LdGDP-MP and hGDP-MP, the fixation of one substrate would be favored by the binding of the other. For LmGDP-MP, both substrates would bind sequentially to the enzyme, but independently from one another.

**Selection of leishmanial GDP-MP inhibitors.** Based on 3D-models of leishmanial GDP-MPs, an in-house library of 100 compounds was designed and synthesized⁶, ⁷, ²³. The activity of these compounds was firstly evaluated at a concentration of 100 µM on the three purified GDP-MPs to select the molecules that specifically inhibit the enzyme of the parasite (Fig. 5a). Among the 100 products evaluated, 11 showed an enzyme inhibition above 50% on a leishmanial GDP-MP, namely compounds **46, 47, 71, 76, 81, 82, 83, 92, 98, 99** and **100**. In order to determine the affinity of these compounds for the GDP-MPs, their Kᵢ was investigated on the three purified enzymes. Out of the 11 molecules, only 5 (**46, 83, 92, 99** and **100**) exhibit a significant Kᵢ on a leishmanial GDP-MP (Fig. 5b), indicating that the 6 remaining products have a low affinity for the enzyme of the parasite. With a promising Kᵢ at 7.00 ± 3.39 µM, compound **99** inhibits specifically and competitively LdGDP-MP, but not the two other purified GDP-MPs, for which no Kᵢ could be determined (Fig. 5b). A competitive inhibition was also observed with compound **100** on both LdGDP-MP and hGDP-MP with Kᵢ values at 61.79 ± 16.32 µM and 19.74 ± 3.87 µM, respectively (Fig. 5b). These results show that despite a modest activity on LdGDP-MP, compound **100** could inhibit hGDP-MP with a higher affinity. Compounds **46, 83** and **92** exhibited Kᵢ values between 15 to 25 µM on LmGDP-MP. No Kᵢ was obtained with these products on LdGDP-MP or hGDP-MP indicating a specificity of action on LmGDP-MP. Nonetheless, these compounds exert a non-competitive (compounds **46** and **92**) or uncompetitive (compound **83**) inhibition on LmGDP-MP (Fig. 5b). The types of inhibition were also graphically determined and confirmed by Dixon plots for competitive or non-competitive inhibitors or Cornish-Bowden plots for uncompetitive inhibitors (Supplementary Fig. S6).

**Docking analysis of competitive inhibitors.** In order to further study and compare the interactions of competitive inhibitors identified in this work (compounds **99** and **100**) in the catalytic site of LdGDP-MP and hGDP-MP, a docking analysis of these products was performed on structural models of the enzymes⁷, with the natural ligand GDP-mannose used as a reference (Fig. 6a). While compound **100** binds to both LdGDP-MP and hGDP-MP with similar potency as for GDP-mannose, compound **99** binds preferentially and more potently to LdGDP-MP. These results are in good agreement with the Kᵢ measurements of these two competitive inhibitors on LdGDP-MP and hGDP-MP (Fig. 5b). Furthermore, the preferential position of the quinoline group of compound **99** is similar in both enzyme species, while its diisopropylphosphonate groups are located more deeply in the LdGDP-MP binding site (Fig. 6b) than in the hGDP-MP (Fig. 6c), in line with the more potent activity of compound **99** on LdGDP-MP than on hGDP-MP. More specifically, the quinoline moiety is surrounded by the three residues L14, V15 and E88 in LdGDP-MP, and, similarly, by the four residues L6, S54, E80 and P89 in hGDP-MP. In both enzyme species, the magnesium ion is chelated by two aspartate residues (D188 and D226 in LdGDP-MP, and D110 and D218 in hGDP-MP) and makes favorable electrostatic interactions with the triazole nitrogens. Nevertheless, the triazole cycle is more tightly held by the three residues K31, T94 and N116 in LdGDP-MP than by the only residue K23 in hGDP-MP. But the most striking difference is that the diisopropylphosphonate groups interact with the residues Y152, E169, K170 and N180 in LdGDP-MP, whereas they only interact with Y144 and R13 (which is located slightly outside of the active site) in hGDP-MP. When comparing the pose of compound **99** in the protein catalytic site with the one of GDP-Mannose⁷, it can be noted that one of its two diisopropylphosphonate groups is located in the cavity that is naturally occupied by the mannose group in LdGDP-MP, but this cavity is not filled by compound **99** in hGDP-MP. Altogether, these docking results confirm the selective and competitive inhibition observed with compound **99** on LdGDP-MP (Fig. 5b).

**In vitro antileishmanial activity and cytotoxicity of compounds 46, 83, 92, 99 and 100.** The antileishmanial activity of compounds **46, 83, 92, 99** and **100** was investigated on both axenic and intramacrophage amastigotes of *L. donovani* and *L. mexicana* on two cell host models, RAW264.7 macrophages and bone marrow derived macrophages (BMDM), the latter being closer to clinical conditions (Table 1). Concerning the RAW264.7 model, compound **99** shows a promising IC₅₀ on both axenic and intramacrophage amastigotes of *L. donovani* at 1.06 ± 0.10 µM and 0.63 ± 0.14 µM, respectively. However, this compound has a CC₅₀ at 1.53 ± 0.17 µM resulting in a modest but noticeable SI value of 2.4 on *L. donovani*. Interestingly, compound **99** did not show any activity on axenic amastigotes of *L. mexicana*, showing thus a specificity of action on *L. donovani* which is in agreement with the enzymatic results (Fig. 5b). Compound **100** showed similar antileishmanial activities of both parasite species with IC₅₀ between approximately 30 to 50 µM and 20 to 27.5 µM on axenic and intramacrophage amastigotes, respectively. With a cytotoxicity at 62.06 ± 7.39 µM, the SI of this compound is comparable to compound **99** on both *L. donovani* and *L. mexicana*. The IC₅₀ of compounds **46, 83** and **92** is above
<caption>a.</caption>

LdGDP-MP

<figure><img src="image_16.png" /><figcaption>c.</figcaption></figure>

LmGDP-MP

<figure><img src="image_17.png" /><figcaption>e.</figcaption></figure>

hGDP-MP

<figure><img src="image_18.png" /><figcaption>b.</figcaption></figure>

LdGDP-MP

<figure><img src="image_19.png" /><figcaption>d.</figcaption></figure>

LmGDP-MP

<figure><img src="image_20.png" /><figcaption>f.</figcaption></figure>

hGDP-MP

<figure><img src="image_21.png" /></figure>

**Figure 4.** Determination of LdGDP-MP, LmGDP-MP and hGDP-MP mechanisms of reaction. (a and b) Lineweaver Burk plots 1/V = f(1/[Man-1-P]) (a) and 1/V = f(1/[GTP]) (b) of LdGDP-MP. Constant concentrations of GTP (a) and Man-1-P (b) were held at 10 µM, 20 µM, 40 µM and 150 µM. The results expressed correspond to the mean of three independent experiments ± SD. (c and d) Lineweaver Burk plots 1/V = f(1/[Man-1-P]) (c) and 1/V = f(1/[GTP]) (d) of LmGDP-MP. Constant concentrations of GTP (c) and Man-1-P (d) were held at 5 µM, 10 µM, 20 µM and 50 µM. The results expressed correspond to the mean of three independent experiments ± SD. (e and f) Lineweaver Burk plots 1/V = f(1/[Man-1-P]) (e) and 1/V = f(1/[GTP]) (f) of hGDP-MP. Constants concentrations of GTP (e) and Man-1-P (f) were held at 10 µM, 20 µM, 40 µM and 80 µM. The results expressed correspond to the mean of three independent experiments ± SD.

100 µM on axenic amastigotes of *L. donovani*, showing a specific activity on *L. mexicana* which is consistent with the GDP-MP inhibiton assays (Fig. 5b). Compound **46** presents an IC₅₀ at 7.69 ± 0.56 µM and 11.72 ± 1.13 µM on axenic and intramacrophage amastigotes of *L. mexicana*, respectively. Interestingly, compound **92** shows an IC₅₀ on intramacrophage amastigotes of *L. mexicana* below 10 µM. With an absence of cytotoxicity at 100 µM, compounds **46** and **92** exhibit attractive SI above 8.5 and 12.1 on *L. mexicana*, respectively. Compound **83** exhibits a moderate antileishmanial activity with IC₅₀ between 30 and 45 µM on both amastigote forms of *L. mexicana*. On
<figure><img src="image_22.png" /><figcaption>**Figure 5.** Evaluation of compound activities on *Ld*GDP-MP, *Lm*GDP-MP and *h*GDP-MP. (a) Percentage of *Ld*GDP-MP (O), *Lm*GDP-MP (Δ) and *h*GDP-MP (■) inhibition as a function of compounds numerical ID (NID). Each compound was used at 100 µM. Compounds showing inhibition above 50% are represented in green and red for leishmanial and human GDP-MPs, respectively. For a question of readability, all values at 0% inhibition were removed. The results expressed correspond to the mean of three independent experiments. (b) Lineweaver plots double reciprocal plots 1/V = f(1/[Man-1-P]) of compounds **99**, **100**, **46**, **92** and **83** on *Ld*GDP-MP, *Lm*GDP-MP and *h*GDP-MP. Each compound was evaluated with a range of concentrations: 0–67.5 µM for compound **99** on *Ld*GDP-MP, 0–50 µM for compound **46** on *Lm*GDP-MP and compound **100** on *h*GDP-MP, and 0–100 µM for compounds **83** and **92** on *Lm*GDP-MP and compound **100** on *Ld*GDP-MP. Man-1-P was used at 10–100 µM for *Ld*GDP-MP and *h*GDP-MP and 2.5–50 µM for *Lm*GDP-MP. GTP was held constant at 100 µM for *Ld*GDP-MP and *h*GDP-MP, and 50 µM for *Lm*GDP-MP. The formula of each compound is represented on the left of the figure. The compounds *K*ᵢ are indicated on each plot, as well as the type of inhibition. ND: no *K*ᵢ could be determined since the double reciprocal plots did not fit with any (competitive, non-competitive, uncompetitive or mixed) inhibition model. The results expressed correspond to the mean of three independent experiments ± SD.</figcaption></figure>
<table><thead><tr><th>a.</th><th>GDP-Ma</th><th>99</th><th>100</th></tr></thead><tbody><tr><th>LdGDP-MP</th><td>23 %</td><td>40 %</td><td>24 %</td></tr><tr><th>hGDP-MP</th><td>22 %</td><td>6 %</td><td>26 %</td></tr></tbody></table>

<figure><img src="image_23.png" /><figcaption>**Figure 6.** Docking analysis of compounds **99** and **100** on *LdGDP-MP* and *hGDP-MP*. (a) Percentage of “good” poses of GDP-Mannose, compound **99** and **100** over the hundred conformations generated by the ten docking calculations onto *LdGDP-MP* and *hGDP-MP*. For GDP-Mannose, a “good” pose is the one observed in the PDB crystallographic structure 2X5Z. For the two compounds **99** and **100**, a “good” pose has the quinoline group located on the GDP-Mannose guanine. (b,c) Position and orientation of the lowest energy “good” pose of compound **99** in the catalytic site of *LdGDP-MP* (b) and *hGDP-MP* (c). The protein surface is colored as a function of the charge density (red, white and blue colors indicating negative, neutral and positive area, respectively). The amino-acids that make contacts with compound **99** are indicated by their one-letter code and number in the sequence.</figcaption></figure>

the BMDM model, compound **99** appears to be the most promising one as it is active both on *L. donovani* and *L. mexicana* intramacrophage amastigotes at the 1–10 micromolar range with the best selectivity index. On both host models, compound **99** showed the highest selectivity index, rendering this compound as the most promising one.

## Discussion

In this work, two leishmanial GDP-MPs from two geographically distant species (*LdGDP-MP* and *LmGDP-MP*), as well as the most homologous human counterpart (*hGDP-MP*) were expressed and purified in order to compare their enzymatic properties and to select new specific antileishmanial agents. Following expression of recombinant GDP-MPs, several successive chromatography methods were used for proteins purification. Size exclusion chromatography showed that *LdGDP-MP* adopts different oligomerization states: one at very high molecular weight presumably corresponding to soluble aggregates, and three others matching with hexamer, trimer and monomer forms of the enzyme. Several approaches were used to dissociate *LdGDP-MP* soluble aggregates into lower molecular weight oligomeric forms. Co-expression with chaperone proteins, changing buffer conditions (pH, salt concentration, glycerol) or different types of chromatography were unfruitful. An attempt to fuse *LdGDP-MP* with MBP tag was also performed, but without success. Nonetheless, according to DLS analysis, the concentrated pool of fractions of peaks 2, 3 and 4 is predominantly in the hexamer form and could be used for enzyme assays, showing that *LdGDP-MP* has the propensity to hexamerize when the concentration of enzyme is increased. Similar phenomena have been observed for *LmGDP-MP*, whose hexamerization was dependent on enzyme concentration²¹. Moreover, *LmGDP-MP* gel filtration profile showed a single peak corresponding to a hexamer, which is in agreement to our study (Fig. 1e). These results indicate that the *LmGDP-MP* quaternary structure would be more stable than *LdGDP-MP*. Although *LmGDP-MP* and *LdGDP-MP* share 94% identity and amino acids directly involved in substrate binding are conserved in both enzymes, the difference observed in catalytic efficiency between both GDP-MPs (Fig. 3g) could be explained by residues that would influence protein quaternary structure and thus enzyme activity. Indeed, in other NTP transferases, such as UDP-glucose pyrophosphorylase, single amino acid replacements have been reported to alter the quaternary structure as well as the activity of the enzyme²⁴. Furthermore, the hexameric form of *LmGDP-MP* has been shown to dissociate in trimers and monomers at low ionic strength²¹. Thus, once diluted in low ionic strength conditions for enzymatic assays, both *LmGDP-MP* and *LdGDP-MP* hexamers could dissociate in a mixture of oligomeric forms that would be able to catalyze the enzyme reaction. In gel filtration, *hGDP-MP* is eluted as a dimer showing that this subunit is able to self-associate. In mammalian, native GDP-MP is an oligomer of about 450 kDa composed of subunits α and β, confirming a self-association of both subunits to form a high molecular weight complex¹⁸, ²⁵. This native complex has been described to synthesize both GDP-mannose and GDP-glucose. However, only the β subunit has a high activity in synthesizing GDP-mannose. Likewise, in our study, no significant activity was observed with *hGDP-MP* in the presence of substrate analogues such as Glu-1-P sugar monophosphate, or nucleotide triphosphates. Only ITP exhibited 24% of activity on *hGDP-MP* compared to GTP. Similar results were obtained with
<table><caption>**Table 1.** Antileishmanial and cytotoxic activities of compounds 46, 83, 92, 99 and 100. The results expressed correspond to the mean of three independent experiments (±SD). (a) Miltefosine: reference compound. (b) BMDM: Bone Marrow Derived Macrophages.</caption><thead><tr><th rowspan="3">Compounds</th><th colspan="3">IC<sub>50</sub> (µM) ± SD</th><th colspan="3">L. mexicana</th><th colspan="3">Cytotoxicity CC<sub>50</sub> (µM) ± SD</th><th colspan="3">Selectivity index (CC<sub>50</sub>/IC<sub>50</sub>)</th></tr><tr><th colspan="3">L. donovani</th><th colspan="3">L. mexicana</th><th colspan="3">L. donovani</th><th colspan="3">L. mexicana</th></tr><tr><th>axenic amastigotes</th><th>infected RAW.264.7 macrophages</th><th>infected BMDM<sup>b</sup></th><th>axenic amastigotes</th><th>infected RAW.264.7 macrophages</th><th>infected BMDM<sup>b</sup></th><th>RAW264.7 macrophages</th><th>BMDM<sup>b</sup></th><th>RAW264.7 macrophages</th><th>BMDM<sup>b</sup></th><th>RAW264.7 macrophages</th><th>BMDM<sup>b</sup></th></tr></thead><tbody><tr><td>46</td><td>&gt;100</td><td>&gt;100</td><td>12.54 ± 0.93</td><td>7.69 ± 0.56</td><td>11.72 ± 1.13</td><td>15.71 ± 3.48</td><td>&gt;100</td><td>13.10 ± 1.67</td><td>—</td><td>1.0</td><td>&gt;8.5</td><td>0.8</td></tr><tr><td>83</td><td>&gt;100</td><td>51.88 ± 7.49</td><td>27.14 ± 2.04</td><td>30.64 ± 3.29</td><td>45.41 ± 6.21</td><td>25.92 ± 5.37</td><td>&gt;100</td><td>23.76 ± 1.98</td><td>&gt;1.9</td><td>0.9</td><td>&gt;2.2</td><td>0.9</td></tr><tr><td>92</td><td>&gt;100</td><td>41.60 ± 13.99</td><td>82.45 ± 7.61</td><td>43.15 ± 3.42</td><td>8.25 ± 1.07</td><td>79.29 ± 9.76</td><td>&gt;100</td><td>77.54 ± 5.34</td><td>&gt;2.4</td><td>0.9</td><td>&gt;12.1</td><td>1.0</td></tr><tr><td>99</td><td>1.06 ± 0.10</td><td>0.63 ± 0.14</td><td>1.06 ± 0.41</td><td>&gt;100</td><td>1.49 ± 0.26</td><td>8.59 ± 2.44</td><td>1.53 ± 0.17</td><td>&gt;100</td><td>2.4</td><td>&gt;94.3</td><td>1.0</td><td>&gt;11.6</td></tr><tr><td>100</td><td>30.68 ± 6.62</td><td>19.52 ± 4.53</td><td>12.18 ± 4.74</td><td>49.25 ± 0.26</td><td>27.52 ± 2.23</td><td>12.05 ± 1.27</td><td>62.06 ± 7.39</td><td>51.36 ± 3.45</td><td>3.2</td><td>4.2</td><td>2.3</td><td>4.3</td></tr><tr><td>Miltefosine<sup>a</sup></td><td>2.08 ± 0.24</td><td>1.83 ± 0.22</td><td>0.83 ± 0.12</td><td>23.74 ± 2.81</td><td>52.62 ± 4.98</td><td>51.39 ± 5.34</td><td>&gt;25</td><td>&gt;25</td><td>13.7</td><td>&gt;30.1</td><td>&gt;0.5</td><td>&gt;0.5</td></tr></tbody></table>

leishmanial GDP-MPs, showing that the three purified enzymes present high substrate specificity. These data are consistent with the high substrate specificity described in the *T. brucei* GDP-MP, except that this enzyme is also able to use ATP as a substrate with a high $K_m$ at 290 µM¹⁵. In addition, a high cofactor specificity was observed with the three purified enzymes, except that LdGDP-MP is also able to generate a moderate activity in the presence of $Mn^{2+}$. Surprisingly, hGDP-MP has been previously reported to be active in the presence of $Mn^{2+}$, but at a concentration 10 times lower than in our conditions¹⁸.

Our enzymatic assays showed that the three GDP-MPs have very similar optimal reaction conditions with the concentration of $Mg^{2+}$ ranging from 1 mM to 5 mM and a slightly more acidic pH for hGDP-MP compared to the leishmanial enzymes. Except for LmGDP-MP which presents a $K_m$ at 7.67 ± 2.28 µM for GTP, a moderately higher affinity was observed with hGDP-MP for both substrates ($K_m$ ≈ 7 µM) compared to LdGDP-MP and LmGDP-MP ($K_m$ ≈ 10–17 µM). These micromolar $K_m$ values are consistent with those obtained with bacterial or trypanosomal GDP-MPs¹⁴,¹⁵. The investigation of substrate binding order showed that the three purified GDP-MPs exhibit a sequential random mechanism. Although several NTP transferases, such a bacterial GDP-MP or diverse ADP-glucose or TDP-glucose pyrophosphorylases, have been described to exhibit a sequential ordered mechanism, with the NTP binding prior to sugar-1-P¹⁴,²⁶⁻²⁸, a sequential random mechanism was observed in a mammalian UDP-glucose pyrophosphorylase²⁹ which is consistent with our findings on human and leishmanial GDP-MPs.

From an in-house library of 100 compounds designed and synthesized from molecular models of leishmanial GDP-MPs⁶,⁷,²³, 11 molecules presented LdGDP-MP and/or LmGDP-MP inhibition above 50%, and 5 exhibit a $K_i$ on a leishmanial GDP-MP at the micromolar range. Among these 5 compounds, three are derivatives of quinoline substituted in position 2 (compounds **99** and **92**) or 4 (compound **100**). Compound **100**, a 4-substituted quinoline with a chloride in position 7, has been previously described to be the most potent L. major GDP-MP inhibitor identified by high throughput screening with an IC₅₀ at 0.58 µM⁸. This compound presents an IC₅₀ on L. donovani and L. mexicana intramacrophage amastigotes between 20 and 30 µM, which is consistent with the IC₅₀ previously reported at 21.9 µM in L. major⁸. Nonetheless, compound **100** exhibits only moderate or no inhibition on the leishmanial enzymes LdGDP-MP and LmGDP-MP, respectively, while hGDP-MP is inhibited with a noticeable $K_i$ at around 20 µM. Thus, this inhibitor could probably not be used against L. donovani or L. mexicana for the development of new antileishmanials. Indeed, 2-substituted quinolines have been previously shown to present promising antileishmanial activities both *in vitro* and *in vivo*¹¹,³⁰. Thus, the presence of a 2-substituted quinoline in leishmanial GDP-MP inhibitors could potentiate their antileishmanial activities through GDP-MP inhibition. The most active and specific competitive inhibitor on LdGDP-MP is the compound **99** ($K_i$ of 7 µM) which does not inhibit the other GDP-MPs studied. This result is in agreement with the docking analyses showing a significant higher potency of binding of the compound **99** in the catalytic pocket of LdGDP-MP compared to hGDP-MP. Indeed, this 2-substituted quinoline presents a diisopropylphosphonate group located in place of the mannose group of GDP-mannose which could account for the competitive inhibition activity of the molecule. Compound **99** has been designed from a structural analysis of GDP-MP active site molecular models with the aim to get inhibitors with molecular volumes suitable for filling the catalytic pocket. Consequently, this compound initially designed as a substrate analogue has a chemical structure with substituents able to achieve these criteria. Regarding the $K_i$ values of all evaluated inhibitors, compound **99** is the inhibitor presenting the highest affinity for a leishmanial GDP-MP. Although compound **99** exhibits a very interesting IC₅₀ on intramacrophage parasites at the sub-micromolar range and is specifically active on axenic amastigotes of L. donovani, it presents a substantial cytotoxicity leading to a SI value of 2.4 on the RAW264.7 model. However, in the BMDM model, it is active both on L. donovani and L. mexicana intramacrophage amastigotes without any cytotoxicity. These discrepancies can be ascribed to different mechanisms of uptake and accumulation within the host cell. Nevertheless, both BMDM and RAW264.7 cell models are recognized and routinely used in classical drug screening and to decipher cell mechanisms in *Leishmania*³¹⁻³⁵. It is not surprising to observe discrepancies from one model to another as a function of an evaluated chemical series. In our case, the most promising inhibitor, compound **99**, exhibited this tendency with CC₅₀ values at 1.53 µM and >100 µM on RAW264.7 and BMDM cell models, respectively. However, it is the most active compound on both L. donovani and L. mexicana intramacrophage
amastigotes. These promising results have now to be confirmed by a further *in vivo* evaluation. Moreover, in a subsequent fundamental study focused on GDP-MP cell biology, we will extend the investigation by comparing the results in diverse cell models, including human ones, in order to understand the mechanism of action of compound 99. Additionally, this compound could be further used as a pharmacological lead compound to develop new antileishmanial drugs that would be specifically active on the parasite, without toxicity on mammalian host. Compound 92, which is also a 2-substituted quinoline with a chloride in position 6 and a diphenyl group in place of the 2 diisopropylphosphonates of compound 99, exhibits a specific inhibition on *Lm*GDP-MP, and not on the other GDP-MPs studied. This specific inhibition of *Lm*GDP-MP is also observed with compounds 46 and 83, which are quinoxaline and guanine derivatives, respectively. Interestingly, no activity was observed on *L. donovani* axenic amastigotes with compounds 46, 83 and 92, which is consistent with the specific activity observed on *Lm*GDP-MP (Fig. 5b). Interestingly, these three inhibitors do not present any cytotoxicity on RAW264.7 macrophages, leading to an attractive SI above 8.5 and 12.5 for compounds 46 and 92, respectively. Compound 83 exhibits an uncompetitive inhibition on *Lm*GDP-MP and thus binds only to the complex formed between the enzyme and the substrate, while compounds 46 and 92 present a non-competitive inhibition on the same leishmanial enzyme and can bind to the enzyme indifferently of the presence of the substrate. These results indicate that these three compounds are potential allosteric inhibitors of *Lm*GDP-MP. Indeed, targeting allosteric sites has been previously shown to provide new opportunities for the development of specific enzyme inhibitors³⁶.

Altogether, these results allowed identifying compounds 46, 83, 92 and 99 as pharmacological tools which could be further pharmacomodulated to improve their specificity and affinity for leishmanial GDP-MPs as well as their antileishmanial activities, and to generate new specific antileishmanial drugs and expand the current chemotherapeutic arsenal for the treatment of leishmaniasis. This study also shows that a conservation of amino acids at the level of the catalytic site does not prevent to get specific inhibitors since probable extra-catalytic site interactions could occur between the inhibitor and the target, as probably happening with compounds 46, 83 and 92. More interestingly, compound 99 clearly appears to be the most promising one for the following reasons: (i) it is the single competitive inhibitor of *Ld*GDP-MP, (ii) it is active both on *L. donovani* and *L. mexicana* intramacrophage amastigotes at the micromolar range in the BMDM model, (iii) its selectivity index values are higher than 10 on the BMDM host cell model. As perspectives, compound 99 is worth of further *in vivo* evaluation on *L. donovani* Balb/c mice model. In addition, next pharmacomodulations will be set up to understand the structure-activity relationship of derivatives of this lead compound towards GDP-MP and intracellular parasites through deciphering their uptake mechanism.

## Methods

**Chemical synthesis.** Detailed description of synthetic methods and compounds characterization are provided in the Supplementary Methods. For biological evaluations, all compounds were dissolved in DMSO at 10 mM stock solution, aliquoted, and stored at -80 °C before use.

**Cloning, production and purification of recombinants GDP-MPs.** Detailed materials and methods of cloning, production and purification of *Ld*GDP-MP, *Lm*GDP-MP and *h*GDP-MP are provided in Supplementary Methods. The band stained by Coomassie blue on SDS-PAGE corresponding to the each purified *Ld*GDP-MP, *Lm*GDP-MP or *h*GDP-MP was excised, digested and analyzed by peptide mass fingerprint using MALDI-TOF mass spectrometry as previously described³⁷.

**Enzymatic assays.** GDP-MP catalyzes the following reaction:

$$ \text{Man} - \text{1-P} + \text{GTP} \rightarrow \text{GDP} - \text{Mannose} + \text{Pyrophosphate} $$

The reaction media were buffered at pH optima and contained optimal concentrations of Man-1-P, GTP, GDP-MP and Mg²⁺ determined for each purified GDP-MP (Fig. 2d), supplemented with DTT (1 mM). Inorganic pyrophosphatase, which hydrolyses the inorganic pyrophosphates generated by the GDP-MPs in inorganic phosphates, was added at 0.1 U/ml (Sigma-Aldrich). Inorganic pyrophosphatase was added in excess to the reaction medium in order to minimize the probability that the putative inhibitors of GDP-MPs interfere with this enzyme. The reaction was initiated by the addition of GDP-MP, and was stopped after adding 100 µl of revelation buffer containing malachite green (0.03% w/v), ammonium molybdate (0.2% w/v) and Triton X-100 (0.05% v/v) in HCl (0.7 M) for 5 min at the optimal temperature³⁸. The quantity of inorganic phosphate generated, representing GDP-MP activity, was determined by measuring OD₆₅₀nm (Lab systems Multiskan MS). The molar extinction coefficient ε of the end product detected was determined by measuring OD₆₅₀nm with a range of inorganic phosphate concentrations diluted in the revelation buffer and was calculated at 27207 M⁻¹ cm⁻¹ (Supplementary Fig. S2). The kinetic constants Kₘₚ, Vₘ and kₐt were determined by non linear regression using GraphPad Prism 6.0.

**Evaluation of compounds on purified GDP-MPs.** Enzyme activity was measured as described above in the presence of 100 µM of compounds and the percentage of inhibition was calculated in reference to untreated enzyme. As a control, compounds were evaluated in parallel in the same conditions on inorganic pyrophosphatase alone in order to discard false-positives. No significant inhibition (<5%) of inorganic pyrophosphatase was observed with the positive compounds (46, 83, 92, 99 and 100) selected in this study. Furthermore, no effect was observed on the activity of the three purified GDP-MPs with 1% DMSO, corresponding to the maximal solvent concentration used at 100 µM of compound. Compounds presenting a minimum of 50% inhibition on a leishmanial GDP-MP were subsequently selected to determine their type of inhibition and their Kᵢ values. Inhibition data were fitted to competitive, non-competitive or uncompetitive models of enzyme inhibition by non linear least squares regression analysis using GraphPad Prism 6.0.
**Cell Cultures.** All details of bacterial, parasite, macrophage cultures are described in Supplementary Methods.

**In vitro antileishmanial evaluation of compounds on axenic and intramacrophage amastigotes.** The evaluations of activity on axenic and intramacrophage amastigotes were adapted from the protocols previously described³⁹. Briefly, for the evaluation on axenic amastigotes, two fold serial dilutions of the compounds from a maximal concentration of 100 µM were performed in 100 µl of complete medium (see above) in 96-well microplates. Axenic amastigotes were then added to each well at a density of 10⁶/ml in a 200 µl final volume. After 72 h of treatment at 37 °C for *L. donovani* or at 32 °C for *L. mexicana* with 5% CO₂, 20 µl of resazurin (450 µM) was added to each well and further incubated in the dark for 24 h at 37 °C for *L. donovani* or at 32 °C for *L. mexicana* with 5% CO₂. In living cells, resazurin is reduced in resorufin and this conversion is monitored by measuring OD₅₇₀nm (resorufin) and OD₆₀₀nm (resazurin; Lab systems Multiskan MS). The activity of the compounds was expressed as IC₅₀. Miltefosine was used as the reference drug.

Concerning the evaluation on intramacrophage amastigotes, RAW 264.7 macrophages (from ATCC) and bone marrow derived macrophages (BMDM) were plated in 96-well microplates at a density of 2 × 10⁴ cells per well and incubated for 24 h at 37 °C with 5% CO₂ for RAW264.7 cells or at 37 °C with 7% CO₂ for BMDM. Axenic amastigotes were differenciated as described above, centrifuged at 2,000 g for 10 min, resuspended in DMEM complete medium, and added to each well to reach a 16:1 parasite to macrophage ratio. After 24 h of infection at 37 °C with 5% CO₂ for RAW264.7 cells or at 37 °C with 7% CO₂ for BMDM, extracellular parasites were removed, and DMEM complete medium (100 µl) containing two fold serial dilutions of the compounds from a maximal concentration of 100 µM was added to each well. After 48 h of treatment, the medium was removed and replaced by DirectPCR Lysis Reagent (100 µl; Euromedex) before 3 freeze-thaw cycles at room temperature, addition of 50 µg/ml proteinase K, and a final incubation at 55 °C overnight to allow cell lysis. 10 µl of each cell extract was then added to 40 µl of DirectPCR Lysis reagent containing Sybr Green I (0.05%; Invitrogen). DNA fluorescence was monitored using Mastercycler® realplex (Eppendorf). The activity of the compounds was expressed as IC₅₀. Miltefosine was used as the reference drug.

**Evaluation of compounds cytotoxicity.** Cytotoxicity was evaluated on RAW 264.7 macrophages and BMDM. Cells were plated in 96-well microplates at a density of 2 × 10⁴ cells per well. After an incubation of 24 h at 37 °C with 5% CO₂ for RAW264.7 cells or at 37 °C with 7% CO₂ for BMDM, the medium was removed in each well, and 100 µl of DMEM complete medium containing two fold serial dilutions of the compounds was added to each well. After 48 h of incubation at 37 °C with 5% CO₂ for RAW264.7 cells or at 37 °C with 7% CO₂ for BMDM, 10 µl of resazurin (450 µM) was added to each well, and further incubated in the dark for 4 h at 37 °C with 5% CO₂. Cell viability was then monitored as described above. The cytotoxicity of the compounds was expressed as CC₅₀ (Cytotoxic Concentration 50%: concentration inhibiting the macrophages growth by 50%).

**Molecular docking.** To provide a structural basis for the activity of competitive inhibitors on GDP-MPs, three-dimensional structures of enzyme-ligand complexes were predicted using docking calculations. All details of docking calculations are described in Supplementary Methods.

**References**

1. World Health Organization. Investing to overcome the global impact of neglected tropical diseases. Third WHO report on neglected tropical diseases. WHO/HTM/NTD/2015.1 (2015).

2. Sundar, S. & Singh, A. Recent developments and future prospects in the treatment of visceral leishmaniasis. *Ther. Adv. Infect. Dis.* **3**, 98–109, doi:10.1177/2049936116646063 (2016).

3. Gutierrez, V., Seabra, A. B., Reguera, R. M., Khandare, J. & Calderon, M. New approaches from nanomedicine for treating leishmaniasis. *Chem. Soc. Rev.* **45**, 152–168, doi:10.1039/c5cs00674k (2016).

4. Garami, A. & Ilg, T. Disruption of mannose activation in *Leishmania mexicana*: GDP-mannose pyrophosphorylase is required for virulence, but not for viability. *EMBO J* **20**, 3657–3666, doi:10.1093/emboj/20.14.3657 (2001).

5. Stewart, J. et al. Characterization of *Leishmania mexicana* knockout lacking guanosine diphosphate-mannose pyrophosphorylase. *Int. J. Parasitol.* **35**, 861–873, doi:10.1016/j.ijpara.2005.03.008 (2005).

6. Pomel, S., Rodrigo, J., Hendra, F., Cavé, C. & Loiseau, P. M. *In silico* analysis of a therapeutic target in *Leishmania infantum*: the guanosine-diphospho-D-mannose pyrophosphorylase. *Parasite* **19**, 63–70, doi:10.1051/parasite/2012191063 (2012).

7. Daligaux, P. et al. Comparative study of structural models of *Leishmania donovani* and human GDP-mannose pyrophosphorylases. *Eur. J. Med. Chem.* **107**, 109–118, doi:10.1016/j.ejmech.2015.10.037 (2016).

8. Lackovic, K. et al. Inhibitors of *Leishmania* GDP-mannose pyrophosphorylase identified by high-throughput screening of small-molecule chemical library. *Antimicrob. Agents Chemother.* **54**, 1712–1719, doi:10.1128/AAC.01634-09 (2010).

9. Reynolds, K. A., Loughlin, W. A. & Young, D. J. Quinolines as chemotherapeutic agents for leishmaniasis. *Mini Rev. Med. Chem* **13**, 730–743, doi:10.2174/1389557511313050010 (2013).

10. Fournet, A. et al. 2-substituted quinoline alkaloids as potential antileishmanial drugs. *Antimicrob. Agents Chemother.* **37**, 859–863, doi:10.1128/AAC.37.4.859 (1993).

11. Nakayama, H. et al. Efficacy of orally administered 2-substituted quinolines in experimental murine cutaneous and visceral leishmaniasis. *Antimicrob. Agents Chemother.* **49**, 4950–4956, doi:10.1128/AAC.49.12.4950-4956.2005 (2005).

12. Freeze, H. H. & Elbein, A. D. Glycosylation precursors in *Essentials of Glycobiology 2nd edition* (eds Varki, A., Cummings, R.D., Esko, J. D., Freeze, H. H., Stanley, P., Bertozzi, C. R., Hart, G. W. & Etzler, M. E.) 47–61 (Cold Spring Harbor Laboratory Press, 2009).

13. Descoteaux, A. & Turco, S. J. Glycoconjugates in *Leishmania* infectivity. *Biochim. Biophys. Acta* **1455**, 341–352, doi:10.1016/S0925-4439(99)00065-4 (1999).

14. Pelissier, M. C., Lesley, S. A., Kuhn, P. & Bourne, Y. Structural insights into the catalytic mechanism of bacterial guanosine-diphospho-D-mannose pyrophosphorylase and its regulation by divalent ions. *J. Biol. Chem.* **285**, 27468–27476, doi:10.1074/jbc. M109.095182 (2010).

15. Denton, H., Fyffe, S. & Smith, T. K. GDP-mannose pyrophosphorylase is essential in the bloodstream form of *Trypanosoma brucei*. *Biochem. J.* **425**, 603–614, doi:10.1042/BJ20090896 (2010).

16. Shinabarger, D. et al. Purification and characterization of phosphomannose isomerase-guanosine diphospho-D-mannose pyrophosphorylase. A bifunctional enzyme in the alginate biosynthetic pathway of *Pseudomonas aeruginosa*. *J. Biol. Chem.* **266**, 2080–2088 (1991).
17. Wu, B., Zhang, Y., Zheng, R., Guo, C. & Wang, P. G. Bifunctional phosphomannose isomerase/GDP-D-mannose pyrophosphorylase is the point of control for GDP-mannose biosynthesis in *Helicobacter pylori*. *FEBS Lett* **519**, 87–92, doi:10.1016/S0014-5793(02)02717-5 (2002).

18. Ning, B. & Elbein, A. D. Cloning, expression and characterization of the pig liver GDP-mannose pyrophosphorylase. Evidence that GDP-mannose and GDP-Glc pyrophosphorylases are different proteins. *Eur. J. Biochem* **267**, 6866–6874, doi:10.1046/j.1432-1033.2000.01781.x (2000).

19. Carss, K. J. et al. Mutations in GDP-mannose pyrophosphorylase B cause congenital and limb-girdle muscular dystrophies associated with hypoglycosylation of α-dystroglycan. *Am. J. Hum. Genet.* **93**, 29–41, doi:10.1016/j.ajhg.2013.05.009 (2013).

20. Koehler, K. et al. Mutations in GMPPA cause a glycosylation disorder characterized by intellectual disability and autonomic dysfunction. *Am. J. Hum. Genet.* **93**, 727–734, doi:10.1016/j.ajhg.2013.08.002 (2013).

21. Davis, A. J. et al. Properties of GDP-mannose pyrophosphorylase, a critical and drug target in *Leishmania mexicana*. *J. Biol. Chem.* **279**, 12462–12468, doi:10.1074/jbc.M312365200 (2004).

22. Rudolf, F. B. & Fromm, H. J. Plotting methods for analyzing enzyme rate data in *Contemporary enzyme kinetics and mechanism* (ed Purich, D. L.) 53–74 (Academic Press, 1983).

23. Daligaux, P., Pomel, S., Leblanc, K., Loiseau, P. M. & Cavé, C. Simple and efficient synthesis of 5'-aryl-5'-deoxyguanosine analogs by azide-alkyne click reaction and their antileishmanial activities. *Mol. Divers.* **20**, 507–519, doi:10.1007/s11030-015-9652-9 (2016).

24. Führing, J. *et al*. Octamerization is essential for enzymatic function of human UDP-glucose pyrophosphorylase. *Glycobiology* **23**, 426–437, doi:10.1093/glycob/cws217 (2013).

25. Szumilo, T., Drake, R. R., York, J. L. & Elbein, A. D. GDP-mannose pyrophosphorylase. *Purification to homogeneity, properties, and utilization to prepare phosphoaffinity analogs*. *J. Biol. Chem* **268**, 17943–17950 (1993).

26. Boehlein, S. K., Shaw, J. R., Hwang, S. K., Stewart, J. D. & Hannah, L. C. Deciphering the kinetic mechanisms controlling selected plant ADP-glucose pyrophosphorylases. *Arch. Biochem. Biophys.* **535**, 215–226, doi:10.1016/j.abb.2013.04.003 (2013).

27. Barton, W. A. et al. Structure, mechanism and engineering of a nucleotidylyltransferase as a first step toward glycorandomization. *Nat. Struct. Biol* **8**, 545–551, doi:10.1038/88618 (2001).

28. Zuccotti, S. *et al*. Kinetic and crystallographic analyses support a sequential-ordered bi bi mechanism for *Escherichia coli* glucose-1-phosphate thymidylyltransferase. *J. Mol. Biol.* **313**, 831–843, doi:10.1006/jmbi.2001.5073 (2001).

29. Persat, F., Azzar, G., Martel, M. B. & Got, R. Properties of uridine diphosphate glucose pyrophosphorylase from Golgi apparatus of liver. *Biochim. Biophys. Acta* **479**, 329–332, doi:10.1016/0167-4838(83)90243-1 (1983).

30. Loiseau, P. M. *et al*. *In vitro* activities of new 2-substituted quinolines against *Leishmania donovani*. *Antimicrob. Agents Chemother.* **55**, 1777–1780, doi:10.1128/AAC.01299-10 (2011).

31. Aulner, N. *et al*. High content analysis of primary macrophages hosting proliferating *Leishmania* amastigotes: application and antileishmanial drug discovery. *PLoS Neglect. Trop. Dis* **7**, e2154, doi:10.1371/journal.pntd.0002154 (2013).

32. Khraiwesh, M. *et al*. Antileishmanial activity of compounds derived from the medicines for malaria venture open access box against intracellular *Leishmania major* amastigotes. *Am. J. Trop. Med. Hyg.* **94**, 340–347, doi:10.4269/ajtmh.15-0448 (2016).

33. Bastos, M. S. *et al*. Achievement of constitutive fluorescent pLEXSY-egfp *Leishmania braziliensis* and its application as an alternative method for drug screening *in vitro*. *Mem. Inst. Oswaldo Cruz* **112**, 155–159, doi:10.1590/0074-02760160237 (2017).

34. Ricardo-Carter, C. *et al*. *Leishmania major* inhibits IL-12 in macrophages by signaling through CR3 (CD11b/CD18) and down regulation of ETS-mediated transcription. *Parasite Immunol.* **35**, 409–420, doi:10.1111/pim.12049 (2013).

35. Banerjee, S. *et al*. Attenuated *Leishmania* induce pro-inflammatory mediators and influence leishmanicidal activity by p38 MAPK dependent phagosome maturation in *Leishmania donovani* co-infected macrophages. *Sci. Rep* **6**, 22335, doi:10.1038/srep22335 (2016).

36. Wu, P., Clausen, M. H. & Nielsen, T. E. Allosteric small-molecule kinase inhibitors. *Pharmacol. Ther.* **156**, 59–68, doi:10.1016/j.pharmthera.2015.10.002 (2015).

37. Webster, J. & Oxley, D. Protein identification by MALDI-TOF mass spectrometry. *Methods Mol. Biol* **800**, 227–240, doi:10.1007/978-1-61779-349-3_15 (2012).

38. Kodama, T., Fukui, K. & Kometani, K. The initial phosphate burst in ATP hydrolysis by myosin and subfragment-1 as studied by a modified malachite green method for determination of inorganic phosphate. *J. Biochem* **99**, 1465–1472, doi:10.1093/oxfordjournals.jbchem.a135616 (1986).

39. Cheikh-Ali, Z. *et al*. "Squalenoylcurcumin" nanoassemblies as water-dispersible drug candidates with antileishmanial activity. *ChemMedChem* **10**, 411–418, doi:10.1002/cmdc.201402449 (2015).

## Acknowledgements

We thank Eric Prina for providing bone marrow derived macrophages, Dominique Durand for her assistance in DLS experiment, Didier Busso for the construction of the plasmid pnEAvMH, and Stéphanie Yen-Nicolaý and David Cornu for their expertise in MALDI-TOF analyses. This work was funded by DIM Malinf from Région Ile de France.

## Author Contributions

W.M. performed the experiments, analyzed the data and wrote the initial manuscript; P.D. synthesized the compounds; T.H.D. performed docking analyses and participated in writing and editing of the article; C.C. supervised chemical synthesis experiments and participated in writing and editing of the article; N.L. and H.v.T. conceived purification experiments and edited the article; P.M.L. conceived the project and edited the article; SP conceived the project and experiments, analyzed the data and wrote the article.

## Additional Information

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-00848-8

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_24.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons license and your intended use is not per-
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017