# Leishmania Donovani 90 kD Heat Shock Protein - Impact of Phosphosites on Parasite Fitness, Infectivity and Casein Kinase Affinity

Antje Hombach-Barrigah, Katharina Bartsch, Despina Smirlis, Heidi

Rosenqvist, Andrea Macdonald, Florent Dingli, Damarys Loew, Gerald Späth, Najma Rachidi, Martin Wiese, et al.

► To cite this version:

Antje Hombach-Barrigah, Katharina Bartsch, Despina Smirlis, Heidi Rosenqvist, Andrea Macdonald, et al.. Leishmania Donovani 90 kD Heat Shock Protein - Impact of Phosphosites on Parasite Fitness, Infectivity and Casein Kinase Affinity. Scientific Reports, 2019, 9 (1), pp.5074. 10.1038/s41598-019-41640-0. pasteur-02874587

HAL Id: pasteur-02874587

https://pasteur.hal.science/pasteur-02874587

Submitted on 19 Jun 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

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

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

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

OPEN

# *Leishmania donovani* 90 kD Heat Shock Protein – Impact of Phosphosites on Parasite Fitness, Infectivity and Casein Kinase Affinity

Antje Hombach-Barrigah<fnref n="1" />, Katharina Bartsch<fnref n="1" />, Despina Smirlis²,⁴, Heidi Rosenqvist³,⁶, Andrea MacDonald<fnref n="1" />, Florent Dingli⁵, Damarys Loew⁵, Gerald F. Späth², Najma Rachidi², Martin Wiese³ & Joachim Clos<fnref n="1" />

Received: 19 November 2018

Accepted: 11 March 2019

Published online: 25 March 2019

*Leishmania* parasites are thought to control protein activity at the post-translational level, e.g. by protein phosphorylation. In the pathogenic amastigote, the mammalian stage of *Leishmania* parasites, heat shock proteins show increased phosphorylation, indicating a role in stage-specific signal transduction. Here we investigate the impact of phosphosites in the *L. donovani* heat shock protein 90. Using a chemical knock-down/genetic complementation approach, we mutated 11 confirmed or presumed phosphorylation sites and assessed the impact on overall fitness, morphology and *in vitro* infectivity. Most phosphosite mutations affected the growth and morphology of promastigotes *in vitro*, but with one exception, none of the phosphorylation site mutants had a selective impact on the *in vitro* infection of macrophages. Surprisingly, aspartate replacements mimicking the negative charge of phosphorylated serines or threonines had mostly negative impacts on viability and infectivity. HSP90 is a substrate for casein kinase 1.2-catalysed phosphorylation *in vitro*. While several putative phosphosite mutations abrogated casein kinase 1.2 activity on HSP90, only Ser<sub>289</sub> could be identified as casein kinase target by mass spectrometry. In summary, our data show HSP90 as a downstream client of phosphorylation-mediated signalling in an organism that depends on post-transcriptional gene regulation.

*Leishmania donovani* is the causative agent of visceral leishmaniasis, also known as *Kala-Azar*, a generalised infection that causes persistent fever, weight loss, anaemia, and hepato-spleno-megaly in afflicted humans and has a fatal outcome in untreated and treatment failure cases. The annual global incidence is estimated as between 300,000 and 400,000 cases with up to 30,000 fatalities<fnref n="1" />.

The *L. donovani* life cycle is bipartite, with a flagellated, elongated promastigote stage that proliferates in the digestive tract of the transmitting sand flies, and an ovoid, aflagellated intracellular stage, the amastigote, which proliferates inside mammalian phagocytic cells. The conversion between the life cycle stages is crucial for parasite survival in sand flies and mammalian hosts and triggered by changes in the environment, i.e. temperature and pH. In keeping with this, *L. donovani* promastigotes can be triggered into amastigote conversion by elevated culture temperatures combined with acidic growth medium².

An elevated cultivation temperature also causes an increased synthesis of several parasite heat shock proteins³⁻⁸ suggesting their involvement in intracellular survival and life cycle control. Heat shock proteins also play other decisive roles in the intracellular amastigote stage, being a prominent part of the protein payload

<footnote n="1">¹Bernhard Nocht Institute for Tropical Medicine, Hamburg, Germany. ²Institut Pasteur and Institut National de Santé et Recherche Médicale INSERM U1201, Unité de Parasitologie Moléculaire et Signalisation, Paris, France. ³Strathclyde Institute of Pharmacy and Biomedical Sciences (SIPBS) University of Strathclyde, Glasgow, Scotland, UK. ⁴Hellenic Pasteur Institute, Athens, Greece. ⁵Laboratoire de Spectrométrie de Masse Protéomique, Centre de Recherche, Institut Curie, PSL Research University, Paris, France. ⁶Present address: Novo Nordisk A/S, Gentofte, Denmark. Correspondence and requests for materials should be addressed to J.C. (email: clos@bnitm.de)</footnote>
of immune modulatory secreted vesicles known as exosomes. Nearly all major chaperones, including HSP90, HSP70, CPN60.2 and CPN60.3, STIP1 and HSP100 are released into the host cell cytoplasm, with the latter playing a crucial role in the sorting of exosome content. Lack of chaperones in the exosomes abrogates their immune modulatory potential⁹,¹⁰.

Chaperone proteins of the HSP90 family play well known roles in the maturation and activation of a multitude of regulatory proteins, such as transcription factors, protein kinases, hormone receptors and cytoskeletal proteins¹¹,¹² and can even silence mutations in regulatory proteins¹³.

Consequently, HSP90 is an essential protein in all eukaryotic organisms precluding the viability of HSP90 null mutants in any known system. Moreover, HSP90 is a highly abundant protein species, accounting for up to 3% of a cell's total protein³,¹⁴. Therefore, reverse genetics approaches, i.e. creating null mutants, to unravel the role and function of HSP90 are likely to fail. The discovery of the HSP90-specific inhibitor geldanamycin (GA)¹⁵ allowed the pharmacological inactivation of this protein family and the study of the phenotypical consequences in tumor cell lines¹⁶⁻¹⁹, but also in parasitic microorganisms such as *Leishmania donovani*, *Trypanosoma cruzi*, and *Plasmodium falciparum*²⁰⁻²². Geldanamycin has a specific affinity for the HSP90 ATP-binding pocket, an affinity it shares with another specific inhibitor, radicicol (RAD)²³,²⁴. Both inhibitor classes cause a cell cycle arrest in logarithmically growing culture cells and indirectly inhibit the maturation and function of HSP90-dependent client proteins, initially offering great promise as anti-cancer drugs²⁵. In *Leishmania donovani*, pharmacological inhibition of HSP90 can induce morphological and biochemical promastigote-to-amastigote differentiation²⁰,²⁶,²⁷, mimicking the known environmental triggers, i.e. heat shock and low pH, and indicating a pivotal role for HSP90 in environmental sensing and life cycle control.

HSP90 is usually found as a dimer in large complexes known as foldosomes. In their function they depend on ATP hydrolysis and on a cohort of so-called co-chaperones that assist in the recruitment of client proteins and in the regulation of the HSP90 reaction cycle²⁸,²⁹. The HSP90 foldosome can vary in composition depending on the client proteins, but usually includes subunits such as HSP70/40, P23, AHA1 and STIP1 (HOP)²⁹. HSP90, HSP70 and STIP1 are also known as substrates of protein kinases³⁰⁻³², affecting function and subcellular localisation of these chaperones.

*Leishmania* spp also possess an almost complete set of co-chaperones, with the notable absence of CDC-37, but including the essential co-chaperones STIP1, SGT²⁷,²⁹,³³⁻³⁵ and others with less impact on the viability, such as HIP, P23, and AHA1³⁶⁻³⁸. The HSP90 and HSP70 chaperones, but also co-chaperones such as STIP1 and cyclophilin 40, have been shown to be phosphorylated during amastigote stage conversion³⁴,³⁹. Mutagenesis of two phosphorylation sites (P-sites) in *L. donovani* STIP1 caused either general lethality or no phenotype at all³⁴. Mutation of the single P-site in cyclophilin 40 had no phenotypic consequences⁴⁰.

HSP90 is involved in a variety of diverse cellular processes and its complex functions are highly regulated by post-translational modifications, such as phosphorylation, acetylation, SUMOylation, methylation, ubiquitination and s-nitrosylation³¹,⁴¹. Phosphorylation of HSP90 does not only regulate its activity directly, but also its interaction with other chaperones, co-chaperones, nucleotides and client proteins³²,⁴²⁻⁴⁵. The kinases responsible for phosphorylation of HSP90 include double-stranded DNA protein kinase, B-raf, Akt, c-Src, protein kinase A, SweWe1, GSK-3beta, casein kinase 1 and casein kinase 2³¹,³².

Essential protein kinases are discussed as attractive drug targets against cancer but also against infectious diseases, such as leishmaniasis. In *Leishmania*, a variety of kinases have already been identified as potential drug targets, such as mitogen-activated protein kinases (MAPKs), cdc-related kinase 3 (CRK3) and GSK3⁴⁶⁻⁵⁰. The members of the *Leishmania* casein kinase 1 (CK1) family have also emerged as potential drug targets⁵¹. The CK1 family consists of multifunctional Ser/Thr protein kinases, characterised by a highly conserved kinase domain and a specific C-terminal domain responsible for kinase regulation and localisation⁵². CK1 isoform 2 (CK1.2) of *L. donovani* was identified as exokinase released via immune-modulatory exosomes into the host cell cytosol, where it may phosphorylate and modulate host cell proteins⁹,¹⁰,⁵³⁻⁵⁵. This kinase was shown to be essential for intracellular parasite survival.

Thus far, *Leishmania* kinases and phosphatases involved in HSP90 modification are poorly understood. Nevertheless, a recent study demonstrated that the *L. donovani* MAPK1, whose *L. mexicana* orthologue is involved in parasite viability and drug resistance⁵⁶, interacts with HSP90 and HSP70, affects the expression of HSP90 and HSP70 and phosphorylates both chaperones *in vitro*⁵⁷. Given the co-localisation of HSP90 and other major chaperones with CK1.2 in the exosome export pathway¹⁰, this kinase may also act as upstream kinase for HSP90 regulation.

Thus far, three phosphorylation sites have been identified in the *L. donovani* HSP90, Thr₂²³, Ser₅²⁶ and Thr₂¹⁶⁴,⁵⁸, all showing increased modification during promastigote-to-amastigote differentiation. However, an analysis of the impact of phosphorylation on HSP90, a key regulator of the *Leishmania* life cycle²⁰, was hindered so far by the high copy numbers of HSP90-coding genes in *Leishmania*⁵⁹ and the essential nature of HSP90. However, the introduction of an inhibitor-resistant sequence variant of *L. donovani* HSP90, allowing a conditional mutagenesis of this protein²⁷, has paved the way for an analysis of HSP90 phosphorylation sites and their impact on fitness, morphology and infectivity. Also, the ability of CK1.2 to phosphorylate HSP90 *in vitro* was tested on a panel of phosphorylation site mutants, identifying at least one CK1.2-specific site.

## Results

**Identification of HSP90 P-sites and mutagenesis.** The P-sites at Thr₂²³ and Ser₅²⁶ were identified in *L. donovani* HSP90 by phosphoproteome analysis³⁴.

To identify further phosphorylation sites within the *Leishmania* HSP90, phosphoproteomics analyses were performed on cultured promastigotes, *in vitro* differentiated axenic amastigotes and mouse-lesion-derived amastigotes of *L. mexicana*. This led to the identification of seven previously unknown phosphorylation sites at Thr₂¹¹, Thr₂¹⁶, Ser₂⁸⁹, Ser₅²⁶, Ser₅⁹⁴ and Ser₅⁹⁵. Fragmentation spectra can be found in the supplementary material
A

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

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

B

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Schematic representation of the targeted phosphorylation sites and domain organisation of HSP90. (A) HSP90 sequence alignment for *L. donovani* (LdBPK_330360.1) and *L. mexicana* (LmxM.32.0314). Phosphorylation sites (see Table 1) are highlighted. (B) Putative domain architecture of HSP90 and distribution of phosphorylation sites over N-terminal (ND), charged linker (CLD), middle (MD) and C-terminal domain (CD). Standard one-letter amino acid codes apply.</figcaption></figure>

<table><caption>**Table 1.** List of the Ser and Thr residues in HSP90 that were subjected to mutagenesis. Three-letter amino acid codes apply.</caption><thead><tr><td>Position</td><td>Reference</td><td>Constitutive de-phosphorylation</td><td>Phosphomimetic exchange</td><td>Surrounding Sequence</td></tr></thead><tbody><tr><td>Thr<sub>21</sub></td><td>(Mollapour et al., 2011)</td><td>Thr<sub>21</sub> - Ile<sub>21</sub></td><td>Thr<sub>21</sub> - Asp<sub>21</sub></td><td>Ser-Leu-Ile-Ile-Asp-<b>Thr</b>-Phe-Tyr-Ser-Asn-Lys</td></tr><tr><td>Thr<sub>100</sub></td><td>(Hawle et al., 2006)</td><td>Thr<sub>100</sub> - Ile<sub>100</sub></td><td>Thr<sub>100</sub> - Asp<sub>100</sub></td><td>Ile-Ala-Arg-Ser-Gly-<b>Thr</b>-Lys-Ala-Phe-Met-Glu</td></tr><tr><td>Thr<sub>211</sub></td><td>this paper</td><td>Thr<sub>211</sub> - Ala<sub>211</sub></td><td>Thr<sub>211</sub> - Asp<sub>211</sub></td><td>Met-Val-Glu-Lys-Thr-<b>Thr</b>-Glu-Lys-Glu-Val-Thr</td></tr><tr><td>Thr<sub>216</sub></td><td>this paper</td><td>Thr<sub>216</sub> - Ala<sub>216</sub></td><td>Thr<sub>216</sub> - Asp<sub>216</sub></td><td>Thr-Glu-Lys-Glu-Val-<b>Thr</b>-Asp-Glu-Asp-Glu-Glu</td></tr><tr><td>Thr<sub>223</sub></td><td>(Morales et al., 2010)</td><td>Thr<sub>223</sub> - Ala<sub>223</sub></td><td>Thr<sub>223</sub> - Asp<sub>223</sub></td><td>Glu-Asp-Glu-Glu-Asp-<b>Thr</b>-Lys-Kys-Ala-Ala-Glu</td></tr><tr><td>Ser<sub>289</sub></td><td>this paper</td><td>Ser<sub>289</sub> - Ala<sub>289</sub></td><td>Not done</td><td>Phe-Tyr-Lys-Ala-Ile-<b>Ser</b>-Asn-Asp-Trp-Glu-Asp</td></tr><tr><td>Ser<sub>526</sub></td><td>(Morales et al., 2010)</td><td>Ser<sub>526</sub> - Ala<sub>526</sub></td><td>Not done</td><td>Val-His-Phe-Glu-Glu-<b>Ser</b>-Glu-Glu-Glu-Lys-Gln</td></tr><tr><td>Ser<sub>594</sub></td><td>this paper</td><td>Ser<sub>594</sub> - Ala<sub>594</sub></td><td>Ser<sub>594</sub> - Asp<sub>594</sub></td><td>Gln-Ala-Lys-Arg-Asp-<b>Ser</b>-Ser-Met-Ala-Gln-Tyr-Met</td></tr><tr><td>Ser<sub>595</sub></td><td>this paper</td><td>Ser<sub>595</sub> - Ala<sub>595</sub></td><td>Ser<sub>595</sub> - Asp<sub>595</sub></td><td>Gln-Ala-Leu-Arg-Asp-Ser-<b>Ser</b>-Met-Ala-Gln-Tyr-Met</td></tr><tr><td>Thr<sub>693</sub></td><td>(Muller et al., 2013)</td><td>Thr<sub>693</sub> - Ala<sub>693</sub></td><td>Thr<sub>693</sub> - Asp<sub>693</sub></td><td>Glu-Val-Thr-Ala-Gly-<b>Thr</b>-Ser-Ser-Met-Glu-Gln-Val-Asp</td></tr><tr><td>Ser<sub>694</sub></td><td>(Muller et al., 2013)</td><td>Ser<sub>694</sub> - Ala<sub>694</sub></td><td>Ser<sub>694</sub> - Asp<sub>694</sub></td><td>Glu-Val-Thr-Ala-Gly-Thr-<b>Ser</b>-Ser-Met-Glu-Gln-Val-Asp</td></tr></tbody></table>

(Fig. S1). All sites are conserved in *L. donovani* (Fig. 1A). The Thr<sub>223</sub> P-site of *L. donovani* is not conserved in *L. mexicana*.

We also identified putative P-sites by literature search and sequence alignments. Two P-sites at Thr<sub>22</sub> and Thr<sub>101</sub> are located in the ATP-binding pocket of the *S. cerevisiae* HSP90, and Thr to Ile exchanges affect the ATPase activity of ScHSP90<sup>43,60</sup>. In *L. donovani* HSP90, these residues are conserved at positions Thr<sub>21</sub> and Thr<sub>100</sub> and were mutated identically, namely Thr<sub>21</sub> to Ile<sub>21</sub> and Thr<sub>100</sub> to Ile<sub>100</sub>.

The putative P-sites at Thr<sub>693</sub> and Ser<sub>694</sub> were inferred from the known P-sites Thr<sub>725</sub> and Ser<sub>726</sub> in the human cytosolic αHSP90 32. A list of the serine and threonine residues targeted for replacement with non-modifiable alanine and/or phosphomimetic aspartate is shown in Table 1 and schematically in Fig. 1.

We have previously established a system for a conditional phenotype analysis of HSP90 mutants *in vitro*. Briefly, a Leu33Ile mutant of HSP90, dubbed HSP90rr, was created which is resistant to the HSP90-specific inhibitor RAD. Ectopic expression of HSP90rr from episomal transgenes has no effect under normal culture conditions, but under RAD challenge, HSP90rr facilitates normal growth, morphology and infectivity. Any mutation introduced into HSP90rr only shows a phenotype under RAD challenge when the endogenously coded HSP90 is inactivated<sup>27</sup>.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Impact of phosphorylation site mutations on the *in vitro* proliferation. *L. donovani* promastigotes were seeded at 1 × 10⁶ cells ml⁻¹ in Medium199 containing radicicol at IC₉₀ (batch-dependent; 0,5-1,25 µg/ml) and allowed to proliferate for 96 h. Median (±range) cell density was recorded and plotted as [%] relative to the cell density of promastigotes expressing HSP90rr. Standard one-letter amino acid codes apply. (n.s.): not significant, (***): p < 0.001 by U-test; *n* = 8. (A) ATP-binding pocket mutations. (B) Charged linker domain mutations. (C) Middle domain mutations. (D) C-terminal domain mutations.</figcaption></figure>

Using a non-commercial procedure for targeted mutagenesis²⁷ and specific oligonucleotide primer pairs for each target sequence (Table S1), we introduced the desired codon changes into the HSP90rr coding sequence and placed the mutated and verified genes in the shuttle expression plasmid pTLv6²⁷. The various mutated HSP90rr transgenes were then stably transfected under G418 selection in *L. donovani* strain 1SR. Faithful transcription from the transgenes was verified by transgene-specific RT qPCR (not shown).

The growth kinetics of parasites expressing the transgenes for HSP90wt, HSP90rr and mutants thereof were monitored in the absence of RAD and found to be indistinguishable (not shown), indicating that the expression of HSP90wt from the chromosomal gene copies masks the mutant phenotypes, leading us to conclude that none of the mutations has a dominant negative effect.

**Effect of P-site mutations on *in vitro* growth.** To assess the effect of the P-site mutations in HSP90rr on *in vitro* growth, we seeded promastigotes in medium without selective antibiotic and added RAD at 90% growth-inhibitory concentration (IC₉₀). Cell density was then monitored over 4 days. The density of HSP90rr cells after 96 h was then defined as 100% growth and parasites expressing the non-RAD-resistant HSP90wt transgene were taken as standard for the loss of HSP90 function.

Mutations in the N-terminal ATPase domain of HSP90rr at Thr₂₁ and Thr₁₀₀ both had deleterious effects on the proliferation under RAD inhibition (Fig. 2A). Interestingly, both Thr to Ile (non-phosphorylatable) and Thr to Asp (phosphomimetic exchange) had the same effect, suggesting that the phosphomimetic exchange was not functional.

All three known or suspected P-sites in the charged linker domain (CLD), Thr₂₁₁, Thr₂₁₆ and Thr₂₂₃, were mutated into alanine or aspartate, respectively. All mutations caused a significant ~30% reduction of the *in vitro* growth rates compared to the parental HSP90rr transgene, but no differences between alanine and aspartate exchanges were observed. Combined mutations of both Thr₂₁₁ and Thr₂₁₆ had significant additive effects, reducing growth by up to 65% (Fig. 2B).

Both mutations in the middle domain (MD) of HSP90rr caused a significant growth retardation of 55% for Ser₂₈₉Ala and 25% for Ser₅₂₆Ala (Fig. 2C). Since the MD is thought to specify client protein affinity, phosphorylation-mediated changes may affect client protein recognition.

Two pairs of adjacent P-sites were targeted in the C-terminal domain which is important for the interaction with co-chaperones 27 and for HSP90 dimerisation. For this analysis, the P-sites Ser₅₉₄, Ser₅₉₅, Thr₆₉₃ and Ser₆₉₄ were mutated into Ala or Asp, respectively. While the Ser₅₉₄Ala exchange caused a significant (~50%) reduction of *in vitro* growth, the Ser₅₉₄Asp exchange abolishes parasite growth entirely. The Ser₅₉₅Ala exchange also caused a significant but weaker (~25%) reduction of the *in vitro* growth, but the Ser₅₉₅Asp exchange facilitates parasite proliferation at ~100%. The Thr₆₉₃Ala but not the Thr₆₉₃Asp mutation reduced growth by ~20%. A strong effect was observed for Ser₆₉₄Asp where a permanent negative charge completely abrogates promastigote growth under RAD challenge. By contrast, the Ser₆₉₄Ala exchange had no significant effect (Fig. 2D).

**Effect of phosphosite mutations on promastigote morphology.** HSP90 plays a pivotal role in *L. donovani* life cycle control. Inhibition of HSP90 by geldanamycin or radicicol induces promastigote-to-amastigote conversion in the absence of elevated temperature and/or acidic pH²⁰,²⁶,²⁷,⁶¹. Morphologically, this leads to a change of the cell shape from spindle-shaped to ovoid and to the shortening of the flagellum²⁷. We therefore analysed length and width of the cell body and length of the flagellum of *L. donovani* transfected with the P-site-mutated HSP90rr variants 72 h after RAD treatment, using either scanning electron microscopy (SEM, Fig. 3A)²⁷ or anti-tubulin labelling and fluorescence microscopy (Fig. 3B)³⁸. The lengths of the mutants were compared to the lengths of HSP90rr and HSP90wt over expressing parasites (Fig. 3A,B).

While the parasites bearing the HSP90rr transgene maintain a spindle shape with a long flagellum (~9 µm cell body length, ~1 µm cell body width, ~4 µm flagella length) after a 72 h RAD treatment, the cells over expressing
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Morphological analysis of recombinant *L. donovani* populations. *L. donovani* promastigotes were seeded at 1 × 10⁶ cells ml⁻¹ in Medium199 containing radicicol at IC₉₀ (batch-dependent; 0,5-1,25 µg/ml) and allowed to proliferate for 72 h. The cells were fixed and then visualised by SEM or fluorescence microscopy using anti-α-tubulin mAb. By using image analysis, the cell body lengths of 50 (in µm, upper panels), the cell body widths of 25 (in µm, middle panels) and the flagella lengths (in µm, lower panels) of 25 randomly selected cells were analysed and mean values were compiled. Bars represent the mean measurements with SEM, n = 25. (n.d.): not detectable, (n.s.): not significant, *p* ≤ 0.05, **p** ≤ 0.01, ***p** ≤ 0.001, respectively, by U-test. (A) analysis performed by using SEM images. (B) analysis performed by using fluorescence microscopy images. (C) SEM images of HSP90rr and HSP90rr-T100D expressing cells -RAD and + RAD.</figcaption></figure>

the non-RAD-resistant HSP90wt are shorter, with a slightly but significantly shorter flagellum (~6 µm cell body length, ~2 µm cell body width, ~3.5 µm flagellar length) (Fig. 3A).

Changing the two N-terminal P-sites Thr₂₁ and Thr₁₀₀ to Ile or to Asp abrogates the ability of HSP90rr to promote the promastigote shape under RAD challenge. For all these mutants, we observed a significant reduction of the cell body length to ~5–6 µm, an ovoid cell shape with an increase of the cell body width to ~3 µm and a significant reduction of the flagellum to ~2 µm. The Thr₁₀₀Asp exchange in particular resulted in damaged cells under RAD treatment (Fig. 3C) and non-viable promastigotes. Therefore, it was impossible to obtain measurements of this mutant.

Due to the closeness of the positions Thr₂₁₁ and Thr₂₁₆ we created double Ala-mutants (Thr₂₁₁/Thr₂₁₆Ala) and analysed only the morphology of the double mutant due to its indistinguishable growth phenotype from the single mutants (Fig. 2A). The simultaneous Ala-exchange of these two Thr-residues (Thr₂₁₁/Thr₂₁₆Asp) in the CLD reduces the cell body length significantly to ~6 µm and the flagellum length to ~2.5 µm under RAD. The cell body width is also significantly increased.

Peculiar morphologies correlate with the mutations at Thr₂₂₃. Thr₂₂₃Ala causes a shortening of the cell body (~6 µm), but normal cell body width and flagellar length (~4 µm), reminiscent of metacyclic promastigotes⁶². The cell body length phenotype is reversed in the phosphomimetic Thr₂₂₃Asp mutant (8–9 µm), which however shows a shortened flagellum. This result hints at an involvement of Thr₂₂₃ phosphorylation in the maintenance of the procyclic promastigote.

The middle domain P-sites Ser₂₈₉ and Ser₅₂₆ are also critical for maintaining spindle shape. The Ala-exchanges at these positions render the cells shorter (Ser₂₈₉Ala: ~6 µm, Ser₅₂₆Ala: ~7 µm) with a reduced flagellum (Ser₂₈₉Ala: ~2 µm, Ser₅₂₆Ala: ~2 µm).

The impact of the mutations at the C-terminus of HSP90rr were analysed by fluorescence microscopy using anti-tubulin antibodies. This showed identical effects when comparing HSP90wt and HSP90rr over expressing cells (Fig. 3B). The Ser₅₉₄Ala mutation slightly reduces cell body length (~8 µm) and the length of the flagellum (~5.5 µm). The cell body width is slightly increased (~3 µm). The mutation Ser₅₉₄Asp in contrast has a strong and significant impact on the promastigotes morphology, with a shortening of the cell body length to ~7 µm and the flagellum to ~4 µm. The cell body width of these mutants is significantly enlarged to ~4 µm. The mutations Ser₅₉₅Ala and Ser₅₉₅Asp have no significant effects on the cell body lengths, but increased the cell body widths slightly and decreased the length of the flagella to ~6 µm. Also for Thr₆₉₃Ala and Thr₆₉₃Asp we observed an intermediate impact of the mutations on the cell morphology. Ser₆₉₄Ala has no impact on the cell body length
~10 µm) with a slight increase of the cell body width. The length of the flagellum is slightly but not signifi-
cantly reduced. These morphological promastigote-like features are reversed by the phosphomimetic exchange
Ser<sub>694</sub>Asp. The mutation Ser<sub>694</sub>Asp reduces the cell body length (~7 µm) and the length of the flagellum (5 µm).
The cell body width is significantly enlarged to ~5 µm.

The observed impact on the morphologies of the HSP90 variant-expressing parasites largely reflects their
growth rates under RAD. One exemption to this correlation is the Thr<sub>223</sub> site where lack of phosphorylation
appears to drive the mutant towards metacyclogenesis.

**Effect of phosphosite mutations on *in vitro* infectivity.** Several P-sites in HSP90 were identified pre-
viously<sup>34</sup> or in this paper as being preferentially phosphorylated in amastigotes by comparative phosphopro-
teomics of promastigote and axenic amastigote lysates (see Results, 1st section). While axenic amastigotes are
convenient models for stage conversion, we were interested in the impact of phosphorylation on the intracel-
lular survival of true *L. donovani* amastigotes. Testing the sensitivity of primary mouse bone marrow-derived
macrophages against RAD, we found that the concentrations required to render *L. donovani* dependent on the
episomal HSP90rr gene variants already affected the ability of the macrophages to support a *Leishmania* infection
(not shown). We therefore decided to pre-treat the parasites bearing the HSP90rr transgenes with RAD for 48 h
prior to macrophage infection to avoid exposure of the host cells to RAD, as described previously<sup>27</sup>. The effects
of the various HSP90rr mutants on the *in vitro* infectivity were analysed 48 h post-infection and compared with
the recombinant cells bearing the resistant HSP90rr transgene and the HSP90wt over expressing cells (Fig. 4). We
applied a duplex qPCR approach to quantify the ratio of parasite actin to macrophage actin DNA as described
previously<sup>63</sup>.

Under these conditions, over expression of HSP90rr prevents a loss of infectivity due to the RAD treatment,
resulting in a significantly higher median infection ratio compared with HSP90wt over expressing parasites, with
differences slightly varying between individual experiments (Fig. 4A–F).

The Thr<sub>21</sub>Ile mutation resulted in ~50% reduced parasite load, albeit twice as high as the HSP90wt control, while the phosphomimetic exchange Thr<sub>21</sub>Asp reduced the parasite load even below HSP90wt levels. The mutations Thr<sub>100</sub>Ile and Thr<sub>100</sub>Asp both lower the parasite load by ~60% (Fig. 4A). Obviously, the putative P-sites Thr<sub>21</sub> and Thr<sub>100</sub> that are located in the ATP-binding domain are not only critical for promastigote growth but also for infection of macrophages. Surprisingly, the Thr<sub>100</sub> mutants support increased intracellular persistence compared with the HSP90rr control, in spite of their inability to support growth of the promastigote (Fig. 2).

The Ala-exchanges at the P-sites in the charged linker region, Thr<sub>211</sub>Ala and Thr<sub>216</sub>Ala, resulted in non-significantly reduced parasite loads (Fig. 4B). The phosphomimetic exchange Thr<sub>211</sub>Asp further lowered the parasite to host cell ratio, arguing against a simple ionic effect of phosphorylated Thr<sub>211</sub>. Thr<sub>216</sub>Asp by contrast restored infectivity significantly over the Thr<sub>216</sub>Ala mutant (Fig. 4B). Phosphorylation of Thr<sub>216</sub> therefore appears to boost intracellular persistence of *L. donovani*.

The Thr<sub>223</sub>Ala exchange also affected infectivity, reducing parasite loads in macrophages to 57% (Fig. 4C,
Table 2). However, the phosphomimetic Thr<sub>223</sub>Asp exchange cannot restore infectivity significantly, also arguing
against simple ionic effects of phosphorylated Thr<sub>223</sub>.

The phenotypes of Ser<sub>289</sub>Ala and Ser<sub>526</sub>Ala mutations, which are located in the middle domain, are not signif-
icant and suggest a rather limited effect on virulence (Fig. 4D).

The Ser<sub>594</sub>Ala, Ser<sub>595</sub>Ala and Ser<sub>595</sub>Asp mutations in the C-terminal region of HSP90rr reduce infectivity slightly but significantly by ~20–30% (Fig. 4E). By contrast, the cells expressing the phosphomimetic Ser<sub>594</sub>Asp mutant were unable to establish an infection (Fig. 4E). This is presumably caused by the inability of these transgenic parasites to grow under RAD *in vitro* (Fig. 2D).

We also tested the effects of the mutated putative P-sites Thr<sub>693</sub> and Ser<sub>694</sub> located directly in front of the STIP-1
recognition motif 27 at the extreme C-terminus of HSP90. The exchanges Thr<sub>693</sub>Ala, Thr<sub>693</sub>Asp and Ser<sub>694</sub>Ala do
not affect infection ratios. However, the phosphomimetic Ser<sub>694</sub>Asp exchange abrogates intracellular survival
completely (Fig. 4F). Like Ser<sub>594</sub>Asp, the Ser<sub>694</sub>Asp mutant of HSP90rr is unable to promote promastigote prolif-
eration under RAD (Fig. 2D), suggesting an inability to tolerate the pre-treatment before infection.

**Identification of HSP90 as *in vitro* target of casein kinase 1.2.** To gain insight into the mechanisms involved in HSP90 regulation, we investigated whether kinases that have been described to target mammalian HSP90 may also target *Leishmania* HSP90. In particular, we focused on three protein kinases, GSK-3, CK1.2 and DYRK1. Of these, GSK-3 and CK1.2 have been shown to phosphorylate the C-terminal domain of mammalian HSP90. This phosphorylation is particularly important to prevent the binding of HSP90 to CHIP, a co-chaperone containing an ubiquitin ligase activity and thus to prevent the client protein from being degraded³². We also selected DYRK1 as it has been shown to be a priming kinase for GSK-3; it phosphorylates the substrates on a specific site which can be recognised and phosphorylated by GSK-3⁶⁴. Therefore, we expressed and purified recombinant *L. major* CK1.2 51, *L. major* GSK-3⁶⁵ and *L. infantum* DYRK1 to perform *in vitro* kinase assays using recombinant *L. donovani* HSP90wt.

As shown in Fig. 5B lane 4, we detected a signal at 80 kDa revealing the transfer of ³²P onto HSP90 by CK1.2 but not by GSK-3 or DYRK1 (lanes 5 and 6). In contrast, all three kinases were able to phosphorylate myelin basic protein (MBP), a canonical kinase substrate (lanes 1–3). Since we can exclude MBP autophosphorylation⁵¹, these data suggest that HSP90 is a substrate of CK1.2, which is consistent with previous findings on mammalian HSP90, but not a substrate of GSK-3 or DYRK1. This is different from mammalian HSP90 where it has been shown that mammalian GSK-3 phosphorylates HSP90 at its C-terminus³². This may be due to the use of the *L. major* (GSK3 and CK1.2) and *L. infantum* (DYRK1) orthologues in our *in vitro* assay; however, the sequence identity between the GSK3 of *L. major* and *L. donovani* is 98%, and the sequence identity between the DYRK1 of *L. infantum* and *L. donovani* is 99%, making a complete change of client specificity very unlikely.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Effects of HSP90 phosphorylation site mutations on the *in vitro* infectivity (A–F). Recombinant parasites were first pre-treated with RAD for 48 h at IC~90~ and then used to infect bone marrow-derived macrophages (MOI 10:1). After 4 h, free parasites were removed and the infected macrophage cultures were further incubated at 37°C and 5% CO~2~ for 44 h. Genomic DNAs were isolated from free and attached macrophages and subjected to a duplex qPCR against mouse and *Leishmania* actin DNA (Bifeld *et al.*, 2016). The median parasite:host DNA ratio was then compiled. Abbreviations: (n.s.): not significant, \*p ≤ 0.05, \*\*p ≤ 0.01 in a one-tailed, paired Student's t test (A), or in a two-tailed U-test (B–F); *n* = 6.</figcaption></figure>

**Identification of the P-site(s) targeted by CK1.2.** Among the P-sites that we identified, the residues Thr~211~, Thr~216~, Thr~223~, Thr~693~ and Ser~694~, match the canonical consensus site for CK1 family members: S/T(X)~2~_3S/T or D/E(X)~2~_3S/T⁶⁶. To identify which of the previously studied P-sites may be targeted by CK1.2, we expressed HSP90 with serine or threonine to alanine mutations to test which mutation abrogates HSP90 phosphorylation by CK1.2. We tested mutants for all known or suspected P-sites since members of the CK1 family can phosphorylate canonical as well as non-canonical sites such as SLS-Xn-(E/D)n⁶⁷ or K/R-X-K/R-(X)₂-S/T⁶⁸,⁶⁹. We performed *in vitro* kinase assays using recombinant CK1.2 with HSP90wt or HSP90 phosphosite mutants as substrates. As shown in Fig. 5D, HSP90-Thr~693~Ala, -Ser~694~Ala and -Thr~693~Ala/Ser~694~Ala mutants are still phosphorylated by CK1.2 at levels comparable with the wild type control, although we observe that the levels of Ser~693~Ala and
<table><caption>**Table 2.** Compilation of results. Relative (in % of HSP90rr values) *in vitro* growth, morphology, infectivity and substrate potential for CK1.2 of the tested mutants.</caption><thead><tr><th>domains</th><th>AA</th><th>to</th><th>growth [%]</th><th>length [%]</th><th>width [%]</th><th>flagellum [%]</th><th>infectivity [%]</th><th>CK1.2 Substrate</th></tr></thead><tbody><tr><td></td><td>HSP90rr</td><td></td><td>100.0</td><td>100</td><td>100</td><td>100</td><td>100</td><td>n.d.</td></tr><tr><td></td><td>HSP90</td><td></td><td>7.4</td><td>71</td><td>189</td><td>78</td><td>25</td><td>+</td></tr><tr><td rowspan="4">N domain</td><td>Thr<sub>21</sub></td><td>I</td><td>7.3</td><td>70</td><td>229</td><td>75</td><td>29</td><td>±</td></tr><tr><td></td><td>D</td><td>5.1</td><td>59</td><td>247</td><td>61</td><td>11</td><td>n.d.</td></tr><tr><td>Thr<sub>100</sub></td><td>I</td><td>6.3</td><td>69</td><td>221</td><td>53</td><td>38</td><td>±</td></tr><tr><td></td><td>D</td><td>3.9</td><td>n.d.</td><td>n.d</td><td>nd</td><td>10</td><td>n.d.</td></tr><tr><td rowspan="6">CL domain</td><td>Thr<sub>211</sub></td><td>A</td><td>72.9</td><td>n.d.</td><td>n.d</td><td>nd</td><td>79</td><td>n.d.</td></tr><tr><td></td><td>D</td><td>69.2</td><td>n.d.</td><td>n.d</td><td>nd</td><td>54</td><td>n.d.</td></tr><tr><td>Thr<sub>216</sub></td><td>A</td><td>69.6</td><td>n.d.</td><td>n.d</td><td>nd</td><td>61</td><td>n.d.</td></tr><tr><td></td><td>D</td><td>69.4</td><td>n.d.</td><td>n.d</td><td>nd</td><td>114</td><td>n.d.</td></tr><tr><td>Thr<sub>211,216</sub></td><td>A</td><td>52.9</td><td>76</td><td>120</td><td>64</td><td>—</td><td>—</td></tr><tr><td>Thr<sub>211,216</sub></td><td>D</td><td>37.4</td><td>n.d.</td><td>n.d.</td><td>n.d.</td><td>n.d.</td><td>n.d.</td></tr><tr><td rowspan="2">M domain</td><td>Thr<sub>223</sub></td><td>A</td><td>76.0</td><td>74</td><td>111</td><td>85</td><td>57</td><td>—</td></tr><tr><td></td><td>D</td><td>62.0</td><td>94</td><td>134</td><td>51</td><td>61</td><td>n.d.</td></tr><tr><td rowspan="6">C domain</td><td>Ser<sub>289</sub></td><td>A</td><td>47.5</td><td>75</td><td>123</td><td>54</td><td>66</td><td>+</td></tr><tr><td>Ser<sub>526</sub></td><td>A</td><td>71.6</td><td>80</td><td>100</td><td>54</td><td>76</td><td>—</td></tr><tr><td>Ser<sub>594</sub></td><td>A</td><td>58.1</td><td>91</td><td>111</td><td>80</td><td>65</td><td>n.d.</td></tr><tr><td></td><td>D</td><td>6.3</td><td>76</td><td>156</td><td>57</td><td>3</td><td>n.d.</td></tr><tr><td>Ser<sub>595</sub></td><td>A</td><td>74.0</td><td>98</td><td>117</td><td>87</td><td>70</td><td>n.d.</td></tr><tr><td></td><td>D</td><td>98.6</td><td>95</td><td>125</td><td>79</td><td>80</td><td>n.d.</td></tr><tr><td rowspan="2">Ser<sub>594/5</sub></td><td>A</td><td></td><td>n.d.</td><td>n.d.</td><td>n.d.</td><td>n.d.</td><td>—</td></tr><tr><td></td><td>D</td><td>n.d.</td><td>n.d.</td><td>n.d.</td><td>n.d.</td><td>n.d.</td></tr><tr><td rowspan="3">Thr<sub>693</sub></td><td>A</td><td>77.9</td><td>94</td><td>125</td><td>80</td><td>91</td><td>±</td></tr><tr><td></td><td>D</td><td>93.1</td><td>94</td><td>132</td><td>92</td><td>81</td><td>n.d.</td></tr><tr><td>Ser<sub>694</sub></td><td>A</td><td>88.2</td><td>97</td><td>112</td><td>85</td><td>93</td><td>±</td></tr><tr><td></td><td></td><td>D</td><td>6.5</td><td>65</td><td>217</td><td>60</td><td>4</td><td>n.d.</td></tr></tbody></table>

Ser<sub>693</sub>Ala/Thr<sub>694</sub>Ala are slightly lower than those of Thr<sub>694</sub>Ala. By contrast, phosphorylation by CK1.2 was abrogated in Thr<sub>21</sub>Ala, Thr<sub>100</sub>Ala, Thr<sub>211/216</sub>Ala, Thr<sub>223</sub>Ala, Thr<sub>289</sub>Ala, Ser<sub>526</sub>Ala and Ser<sub>594</sub>Ala and Ser<sub>595</sub>Ala mutants of HSP90 (Fig. 5D,F).

We also performed a “cold”, i.e. non-radioactive *in vitro* kinase assay with recombinant HSP90wt and CK1.2 in the presence or absence of 10 µM D4476, a specific CK1 inhibitor<sup>51</sup>. The samples of three independent kinase assays were analysed by mass spectrometry to identify the sites phosphorylated by CK1.2. We found only one site, Ser<sub>289</sub> (Fig. 5G), which is phosphorylated in the presence of CK1.2 but not in the presence of CK1.2 + D4476, suggesting that CK1.2 targets Ser<sub>289</sub>. No phosphorylation was detected at Thr<sub>211</sub>, Thr<sub>216</sub>, Ser<sub>526</sub>, Ser<sub>594</sub>, Ser<sub>595</sub> and Thr<sub>223</sub> in our dataset (data not shown), suggesting that they are not targeted by CK1.2. We cannot exclude that the sites Thr<sub>21</sub>, Thr<sub>100</sub>, Thr<sub>693</sub> and Ser<sub>694</sub> may be true CK1.2 sites as the peptides containing those sites were not detected in our mass spectrometric analyses. However, those sites are only presumed phosphorylation sites based on analogies with model HSP90s. Moreover, only Thr<sub>693</sub> and Ser<sub>694</sub> resemble canonical CK1 sites.

In summary (Table 2), mutations at the established and putative phosphorylation sites in the *L. donovani* HSP90 result in a variety of phenotypic manifestations. The most drastic and general effects are observed after introducing mutations at sites within the highly conserved ATP-binding domain. These HSP90rr variants apparently cannot support growth or infectivity under RAD challenge. Mutation at sites in the charged linker domain show moderate effects on growth and infectivity. Middle domain site mutations result in moderate to strong growth inhibition, morphologic changes and a moderate loss of infectivity, even though one of them is established as the only CK1.2 phosphorylation site. In the C-terminal domain we observe the strongest effects for two of the phosphomimetic Ser-Asp exchanges with no or moderate impact of Ser/Thr to Ala mutations. Clearly, permanent negative charges at these sites are not conducive to parasite fitness. At other sites, Ser/Thr to Asp restore growth, morphology and/or infectivity to a various extent compared with the Ser/Thr to Ala mutations.

## Discussion

Trypanosomatida and perhaps the entire phylum Euglenozoa lack inducible transcription as a means for gene expression control<sup>70,71</sup>. The existing genome databases contain no confirmed genes for transcriptional regulators and cis-regulatory DNA sequence elements such as gene promoters and enhancer elements known from crown group Eukaryota. Moreover, the established mode of multicistronic transcription precludes individual transcriptional control of gene expression<sup>71,72</sup>. Cis-splicing is very rare but trans-splicing has been discussed as point of regulation<sup>73</sup>. This leaves mRNA processing and stability, translation, protein folding and post-translational protein modifications as conceivable targets of gene expression control.
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Protein kinase of HSP90. (A) Three protein kinases, rLmCK1.2, rDYRK 1, and rGSK 1, were incubated with γ⁻³²P-ATP with MBP or with HSP90wt as substrate as indicated below and analysed by SDS-PAGE and Coomassie Brilliant Blue staining. The positions of protein size markers are indicated to the left, the positions of the kinases and of HSP90 are indicated on the right. (B) As in (A), but autoradiograph. (C,D) *In vitro* protein kinase assay of rLmCK1.2 with HSP90 and 5 phosphosite mutants thereof, using γ⁻³²P-ATP. (C) shows a Coomassie Brilliant Blue stain, (D) shows the corresponding autoradiograph. (E,F) As in (C,D), but with 4 other phosphosite mutants of HSP90. (G) Mass spectrometry-based quantification of the phosphoserine/serine ratio at Ser₂₈₉ without (−D4476) or with (+D4476) a CK1.2-specific inhibitor. The y-axis is in log₁₀ scale. All greyscale images were cropped from contiguous parts of the original image. Digital enhancements, using Adobe Photoshop CS3, were performed over the entire greyscale images and restricted to tonality optimisation and size adjustments.</figcaption></figure>
The major chaperone protein HSP90 is subject to various post-translational modifications such as phosphorylation, acetylation, s-nitrosylation, oxidation, SUMOylation or ubiquitination³¹,⁷⁴. Post-translational modifications of HSP90 are known to be involved in the fine-tuning of the chaperone cycle which includes conformational changes, ATPase activity, interaction with co-chaperones and the recruitment and activation of various client proteins³¹,³²,⁴²,⁷⁵. Differences in this fine-tuning processes also adapt HSP90 regulation and function to the unique intracellular environments of different organisms or to various physiological states within different tissues of metazoans²⁹,⁷⁶.

Previous work showed that HSP90, HSP70 and the co-chaperone STIP1 are subject to increased phosphorylation during the promastigote-to-amastigote differentiation of *L. donovani*³⁴. These data suggested that the function of HSP90 is adaptable to the different environments of the distinct life cycle stages. It was previously described that a heat shock of yeast cells decreases the phosphorylation level of HSP90, while an opposite effect was observed looking at the turnover of phosphate groups in heat-shocked HeLa cells⁷⁷,⁷⁸. Phosphorylation of this essential chaperone is one possible mechanism to adapt its function to various environmental conditions and in a species-specific manner.

This paper focuses on the impact of HSP90 phosphorylation by introducing amino acid exchange mutations into the radicicol-resistant HSP90rr variant²⁷ and monitoring their impact on *in vitro* promastigote proliferation, morphology, *in vitro* infectivity and casein kinase 1.2 substrate properties.

Phosphorylation sites are distributed over the whole polypeptide sequence of HSP90 and all putative domains. In our study, we included sites shown as phosphorylation sites for *L. mexicana* HSP90 (Thr₂₁, Thr₂₁₆, Ser₂₈₉, Ser₅₹₄, Ser₅₹₅), *L. donovani* HSP90 (Thr₂₂₃, Ser₅₂₆) and putative sites identified by analogy (Thr₂₁, Thr₁₀₀, Thr₆₉₃, Ser₆₉₄).

We observed that amino acid exchanges at the two putative phosphorylation sites Thr₂₁ and Thr₁₀₀ which are implicated in the activity of the ATPase domain from previous work in *Saccharomyces cerevisae*⁴⁴,⁶⁰,⁷⁹, severely inhibit growth of the promastigote stage. Even phosphomimetic exchanges to Asp in both positions cannot restore the proliferation capacity. We can offer three explanations for this observation, (i) the exchanges may render HSP90 non-functional as chaperone, (ii) the structural changes caused by the amino acid exchanges in the nucleotide binding pocket may impair ATP turnover, or (iii) the changes may render HSP90rr sensitive to RAD again. The latter view is supported by the finding that expression of these four mutant HSP90rr variants results in phenotypes indistinguishable from the cells that are over expressing the RAD sensitive HSP90wt (Table 2). While phosphorylation by casein kinase 2 at the residue Thr₂₂, equivalent to Thr₂₁ in *Leishmania* HSP90, was established for the yeast HSP90⁴³, it was also shown that Thr₂₂Ala and Thr₂₂Glu amino acid exchanges impact on ATPase activity and interaction of HSP90 with co-chaperones such as AHA1 and increase the sensitivity against RAD and GA⁴⁴. However, we could recently demonstrate that lack of AHA1 in *Leishmania* has no impact on RAD sensitivity³⁸, arguing against an impaired HSP90-AHA1 interaction as the cause for the observed, severe phenotypes. Testing the mutants *in vitro* for ATPase activity, RAD sensitivity and AHA1 binding may shed light on this phenomenon in the future.

Ala and Asp exchanges at Thr₂₁ and Thr₂₁₆ in the highly diverged charged linker domain have opposing effects on the parasites' infectivity. A stable negative charge at position 211 reduces parasite loads by almost 50%, while the effect of a non-phosphorylatable Ala at 216 is reverted by the phosphomimetic Asp exchange. Although neither site could be identified as CK1.2 phosphorylation site, a double Thr-Ala exchange at 211 and 216 abrogates CK1.2 phosphorylation of HSP90 and causes a reduction of the flagellum. Another site, Thr₂₂₃, also in the CL domain and specific for Trypanosomatida HSP90s, moderately affects growth and infectivity, as well as CK1.2-dependent phosphorylation, but is no CK1.2 site itself. Of note are the effects on the morphology where abrogation of phosphorylation (Thr to Ala) reduces cell body length, but not flagellar length, thereby resembling the shift from procyclic to metacyclic promastigotes⁸⁰. This is a first hint at a possible role of HSP90 in the transition from proliferative procyclic promastigotes to the infective metacyclic form. The Thr₂₂₃Asp exchange restores cell body length, but increases diameter and reduces flagellar length, resembling early amastigote development. This fits with the identification of Thr₂₂₃ as an amastigote-specific phosphorylation site³⁴.

In the middle domain, which in model HSPs has a function in client protein recognition⁶⁰,⁸¹, we analysed two sites, Ser₂₈₉ and Ser₅₂₆. Ser₂₈₉ has an impact on proliferation and the maintenance of promastigote morphology. It was also identified as the only CK1.2 site known so far in HSP90, and its replacement with Ala reduces *in vitro* infectivity by 24%. Ser₅₂₆, an amastigote-specific site³⁴ moderately affects growth, virulence, flagellar length and CK1.2 phosphorylation, either at Ser₂₈₉ or at Ser₅₂₆, since this site was not detected in mass spectrometry. It is conceivable that the effects of Thr₂₁₁/2₁₆, Thr₂₂₃ and Ser₅₂₆ are at least in part due to their impact on Ser₂₈₉ phosphorylation by CK1.2. The question remains open how mutations at Thr₂₁, Thr₁₀₀, Thr₂₁₁/2₁₆, Thr₂₂₃, but also Ser₅₂₆, Ser₅₹₄, and Ser₅₹₅ can affect phosphorylation at Ser₂₈₉.

The C-terminal domain of HSP90 chaperones is known to facilitate dimerisation and interaction with co-chaperones STIP1²⁷ and prolyl-peptidyl isomerases such as cyclophilin 40⁴⁰,⁸² and others. Four phosphorylation sites were detected or inferred by analogy. Replacement of Ser₅₹₄ caused a ~40% reduction of *in vitro* growth and infectivity while the Ser₆₉₄ to Ala change had little impact. However, the phosphomimetic Ser-to-Asp exchange in these positions completely abrogated HSP90-dependent growth and infectivity. We conclude that permanent negative charges in these positions have a deleterious effect on the HSP90 functionality. The HSP90 reaction cycle is known to include the ordered assembly and disassembly of HSP90-co-chaperone complexes at the C-terminal domain⁸³, possibly requiring frequent reversal of post-translational modifications, such as phosphorylation. A permanent negative charge may therefore impede the ordered transitions of the HSP90 foldosome complex. The effects observed for mutation in two other sites, Ser₅₹₅ and Thr₆₉₃, are moderate by comparison.

One lesson from our data is the strong deleterious effects of presumed phosphomimetic Asp mutations. Unlike MAP kinases where Asp replacements in the activation sites render downstream kinases permanently active, the situation in HSP90 with its multiplicity of active phosphorylation sites and a complex reaction cycle with multiple partnering chaperones and co-chaperones may be too complicated for simplistic approaches. Another problem
resides in our experimental set-up where we express RAD-resistant, mutated HSP90rr before a strong back-
ground of RAD-sensitive, endogenous HSP90. A treatment with RAD is necessary to express the mutant pheno-
types, yet one cannot exclude interactions and even heterologous complexes of HSP90rr and HSP90 monomers.
For example, the HSP90-STIP1 interaction is not affected by RAD treatment²⁷, opening up the possibility that
phenotypes of mutations outside the ATP binding pocket are “watered down” by interaction with endogenously
coded HSP90wt. To rectify this problem, we hope for a future implementation of advanced reverse genetics such
as CRISPR/cas9 based gene editing⁸⁴,⁸⁵ to eliminate the entire cluster of ~17 tandemly arranged HSP90 gene cop-
ies and to express the mutant HSP90s before a null mutant background.

So far, not much is known about the upstream and downstream players in the presumed signal transduction pathways. In the eukaryotic model organisms, from yeasts to humans, signal transduction pathways govern the activity of many gene expression regulators through reversible phosphorylation⁸⁶. Many of these regulators are trans-acting transcription factors, a protein class that is almost entirely absent from the *Leishmania proteome*²¹. Therefore, post-transcriptional regulation, in particular on the protein synthesis and protein modification levels, are important for Trypanosomatida⁴⁹,⁸⁷. How these processes are linked to the signal transduction pathways remains one of the most important questions in the field. Since HSP90 and its co-chaperones are integral parts of the signal transduction pathways governing the *Leishmania* life cycle²⁰,²⁷,⁶¹ and are targets of stage-specific phosphorylation events³⁴, the post-translational modification of HSP90 is likely to have multiple effects on *Leishmania* viability and fitness, some of which were uncovered in our analysis.

**Materials and Methods**

**L. donovani culture, genetic complementation and in vitro infection model.** *L. donovani* 1SR⁸⁸ were cultivated as described²⁷. Batch-specific IC₉₀ were determined for radicicol (Sigma-Aldrich, München, Germany) prior to each experiment. Concentration range was 0.5–1.25 µg/ml. Electrotransfection and antibiotic selection of *L. donovani* promastigotes was carried out as described²⁷.

**In vitro infection model.** In vitro infection of murine bone marrow-derived macrophages with RAD-treated *L. donovani* promastigotes was performed as described earlier²⁷ with two modifications: (i) bone marrow macrophage precursor cells were differentiated using 10–30% supernatant of LADMAC cells⁸⁹, and (ii) quantification of infections were performed by quantitative PCR⁶³. Adherent BMMs were infected at a multiplicity of infection of 10 parasites per macrophage. After 4 hours of incubation at 37 °C in DMEM/F-12, supplemented with GlutaMAX, free parasites were removed by 3 washing steps with PBS, and the infected cells were incubated at 37 °C and 9% CO₂ for another 44 hours. At 48 h post-infection, free cells in the culture supernatant and attached cells were pooled by sedimentation and lysed. Genomic DNA (gDNA) was isolated from the infected cells using the ISOLATE II Genomic DNA Kit (Bioline, Luckenwalde, Germany). Parasites were then quantified by semi-quantitative real-time PCR (qPCR) targeting host cell and parasite actin-coding genes with double labeled probes and using total parasite and host gDNA as the template⁶³. The relative parasite load was defined as the ratio of parasite actin DNA against mouse actin DNA.

**Targeted mutagenesis of HSP90-coding sequences.** The generation of the RAD-resistant HSP90 variant (HSP90rr) was described earlier²⁷. Plasmid pUC:HSP90rr served as template for site-directed mutagenesis. Primers (Supplementary data, Table S1) were phosphorylated using ATP and polynucleotide kinase (PNK, New England Biolabs)⁹⁰ and then used to prime a PCR amplification of pUC:HSP90rr using the iProof-PCR kit for GC-rich DNA (#172–5320) from Bio-Rad Laboratories (München, Germany). Following amplification (30 cycles), the linear PCR product was subjected to ligation (3 h, RT) to form circular plasmids. These were then used to transform competent *E. coli* DH5-α cells (#18265–017, Invitrogen, Karlsruhe, Germany). After amplification and purification⁹⁰ via caesium-chloride density gradient ultracentrifugation (50% w/v CsCl, 6 h, 90,000 rpm, 20 °C, rotor NVT90, Beckman, Krefeld, Germany), the pUC:HSP90rr-derived mutants were verified by DNA sequencing (LGC Genomics, Berlin, Germany).

The derivative of the plasmid pJC45⁹¹, pJC45:HSP90, was also used for site-directed mutagenesis. The same set of primers (Table S1) was used. PCR, ligation, transformation of *E. coli* and CsCl purification of plasmids were performed the same as for the site-directed mutagenesis of pUC:HSP90rr.

**Construction and preparation of recombinant DNA.** The expression plasmid pTLv6 has been described²⁷. Mutated HSP90 coding sequences derived from pUC:HSP90rr were excised with enzymes *Kpn*I and *Bam*HI (compatible with *Bgl*II sticky ends) and ligated into linearised pTLv6 to create the pTLv6:HSP90rr expression plasmids coding for the various HSP90rr variants with P-site codon exchanges.

**Expression profiling.** Semi-quantitative real-time RT-PCR was performed essentially as described⁹². HSP90 transgene-specific primers were pTLv6.HSP90.F2 (CGACGAGGAGGAGGAGGCAG) and pTLv6.HSP90.B1 (GCCAGTACATCACAAGACTC ATAGATCC). HSP90 mRNA abundance was calculated relative to the actin signal.

**Scanning electron microscopy (SEM).** *Leishmania* cells were cultivated for 72 h without RAD or with the IC₉₀ of RAD. The treated and non-treated cells were then washed twice in PBS, fixed in 2% glutaraldehyde in sodium cacodylate buffer and postfixed with 1% osmium. Samples were dehydrated at increasing ethanol concentrations (30–100%). After critical point drying, samples were treated with gold and analysed on a Philips SEM 500 electron microscope. Images were taken using a conventional 35 mm camera, and the developed black-and-white films were digitalised using a HAMA 35 mm film scanner.
**Fluorescence microscopy.** Promastigote cells were incubated for 72 h without RAD or with the IC₉₀ of RAD. After the incubation period, cells (1 × 10⁷ cells) were sedimented, washed twice with PBS and suspended in 1 ml of PBS. Aliquots of the suspension (2 × 10⁵ cells) were applied on microscopic slides and air-dried. After fixing the cells for 2 min in ice-cold methanol, the slides were air-dried for 20 min. Non-adherent cells were removed by gentle washing (0.1% Triton X-100 in PBS) followed by incubation in blocking solution (2% BSA, 0.1% Triton-X 100 in PBS). Slides were then incubated for 1 h with monoclonal mouse anti-tubulin (Sigma-Aldrich, München, Germany, 1:4000), washed thrice and then incubated for 1 h with anti-mouse Alexa 594 (Dianova, Hamburg, Germany, 1:250) and DAPI (Sigma-Aldrich, München, Germany, 1:25). After washing the slides thrice, Mowiol and coverslips were applied and the slides were left to dry for 24 h at 4°C. Fluorescence microscopy was carried out on an EVOS FL Auto epifluorescence microscope (Life Technologies, Darmstadt, Germany).

**Phosphoproteomics.** *Leishmania mexicana* MNYC/BZ/62/M379 promastigotes were cultured in SDM-79 medium⁹³ supplemented with 10% heat-inactivated FCS (Sigma, Steinheim, Germany), 7.5 µg ml⁻¹ hemin (Sigma, Steinheim, Germany) and 100 U ml-1 penicillin/100 µg ml-1 streptomycin (Pen/Strep; Gibco, UK). Cultures were incubated at 27°C until they reached a density of 4–5 × 10⁷ parasites ml-1 (late-log phase). Cultures were either harvested or differentiated into amastigotes by inoculation in Schneider's Drosophila medium (PAN Biotech, Aidenbach, Germany) supplemented with 20% FCS (PAN Biotech), 2 mM L-glutamine, 100 U ml-1 Pen/Strep, and 20 mM 2-morpholinoethane sulfonic acid monohydrate [MES, pH 5.5] (Serva, Heidelberg, Germany). Cultures were incubated at 34 °C/5% CO₂ for 72 h. Parasites were harvested by centrifugation at 2000 × g, 4 °C, and washed in 1) ice-cold wash buffer (21 mM HEPES pH 7.5, 137 mM NaCl, 5 mM KCl) and 2) ice-cold wash buffer with protease and phosphatase inhibitors (1 mM Na-orthovanadate, 0.1 µM okadaic acid, 10 mM NaF, 10 mM o-phenanthroline, EDTA-free protease inhibitors). For lysis, parasites were resuspended in a solution of 7 M urea, 2 M thiourea, 40 mM Tris, 1% n-octyl-β-D-glycopyranoside, 1 mM MgCl₂, 1 mM o-phenanthroline, 300 U benzonase, 1 mM Na-pervanadate (Na-orthovanadate activated in 18% H₂O₂), protease inhibitors (Roche EDTA-free protease inhibitor tablets), and phosphatase inhibitor cocktails (P2850 and P5726 from Sigma), and sonicated for 3 × 15 s on ice. Lysates were incubated at –80 °C for 30 min prior to reduction (50 mM DTT) and alkylation with 50 mM iodoacetamide. Proteins were precipitated in an 8-fold excess of ice-cold acetone-ethanol, 1:1, v/v, by overnight incubation at –20 °C.

Protein precipitates were reconstituted in 6 M urea/2 M thiourea and diluted in 50 mM NH₄HCO₃ for digestion with trypsin at a 1:75 enzyme-substrate ratio overnight. Digestion was quenched by addition of formic acid (FA) to a final concentration of 1%. Aliquots of 330 µg digested proteins were used for phosphopeptide enrichment by TiO₂ (Titansphere, 5 µm; GL Sciences, Tokyo, Japan). Enrichment was performed with home-made columns packed in p200 pipette tips, essentially as described before⁹⁴. Samples for LC-MS/MS were prepared by desalting on *Stop And Go Extraction* (STAGE)-tips⁹⁵. A plug of C8 (phosphopeptides) or C18 (non-phosphorylated peptides) Empore disc (3 M Bioanalytical Technologies, St. Paul, MN, USA) was placed in the thin end of an ordinary p200 pipette tip and washed in 100% acetonitrile (MeCN). The reversed-phase (RP) material was equilibrated with 5% FA and the sample loaded in a 1:1 sample-to-5% FA volume ratio. Peptides were eluted into a clean microcentrifuge tube by 70% MeCN/5% FA and vacuum-dried. Upon preparation for LC-MS/MS analysis, peptides were reconstituted in 0.4 µl 100% FA and diluted with a suitable volume of 0.1% FA. Samples were analysed on LTQ Orbitrap XL mass spectrometers (Thermo Fisher Scientific, Bremen, Germany). For peptide separation, an Easy-nLC™-system (Thermo Fisher Scientific (previously Proxeon), Odense, Denmark) was coupled directly to the mass spectrometer. The Easy-nLC™ was fitted with a home-made analytical column (50 µm i.d.; 100 µm o.d., 18 cm) (ReproSil-Pur C18-AQ, 3 µm; Dr. Maisch, GmbH, Ammerbuch-Entringen, Germany). Phosphopeptides were analysed using collision induced dissociation (CID) MultiStage Activation (MSA) in an MS/MS Top10 set-up for fragmentation of the 10 most intense peptide peaks in each MS scan with 45 seconds exclusion time.

**Phosphoproteomics data analysis and phosphorylation validation.** The LCMS data were processed with DTASuperCharge⁹⁶ and searched against the TriTrypDB *L. mexicana* database (version 4.0) (TriTrypDB, http://tritrypdb.org/tritrypdb) using an in-house Mascot server (version 2.2.06, Matrix Science, London, UK). Mascot searches were conducted with 5 ppm peptide tolerance and 0.6 Da fragment ion tolerance. For database searches, cysteine carbamidomethylation was included as a fixed modification, while oxidation of methionine, Gln → pyro-Glu (N-term Q), Glu → pyro-Glu (N-term E), pyrocarbamidomethyl (N-term C) and phospho (S, T, and Y) were added as variable modifications. Data sets were compiled in Scaffold (version 3.3.1, Proteome Software Inc., Portland, OR, USA) with a protein threshold of minimum 95% and a peptide threshold of minimum 95% for at least 1 peptide. The phosphorylated proteins fulfilling these criteria were further evaluated in Scaffold PTM (version ScaffoldPTM_1.1.3, Proteome Software Inc., Portland, OR, USA).

Phosphorylation sites were validated by a semi-automatic approach, applying Ascore⁹⁷ in combination with Mascot Delta Score⁹⁸ and manual inspection of selected spectra. Phosphorylation sites with Ascores >19 were accepted unconditionally.

**Phosphoproteomic identification of CK1.2 phosphosites in HSP90.** For identification of CK1.2 phosphosites, we performed three independent cold kinase assays that were then transferred into fresh tubes to perform in-solution digestion, according to standard protocols. Briefly, after reduction (to a final concentration of 5 mM DTT) and alkylation (to a final concentration of 10 mM iodoacetamide) samples were diluted 10-fold with ammonium bicarbonate and incubated overnight with 0.2 µg trypsin/LysC (Promega) at 37°C. Samples were then loaded onto a homemade C₁₈ stage tip for desalting. Desalted samples were reconstituted in 2% MeCN/0.3% TFA and analysed by nano-LC-MS/MS by using an RSLCnano system (Ultimate 3000, Thermo Scientific) coupled to an Orbitrap Fusion Tribrid mass spectrometer (Thermo Fisher Scientific). Peptides were
analysed using HCD fragmentation with normalised collision energy of 30, in a top speed mode with 45 seconds exclusion time. Spectra were interrogated by Sequest™ through Proteome Discoverer™ 2.1 (Thermo Scientific) with the in-house database containing the sequences of the *Leishmania* HSP90, the *Leishmania* casein kinase 1 (E9AHM9, E9AHM8 and A4IAZ8) and the most abundant contaminants (244 protein sequences). Enzyme specificity was set to trypsin and a maximum of two missed cleavage site were allowed. Oxidised methionine, N-terminal acetylation, phosphorylation of Ser, Thr and Tyr and carbamidomethyl cysteine were set as variable modifications. Maximum allowed mass deviation was set to 10 ppm for monoisotopic precursor ions and 0.6 Da for MS/MS peaks. The resulting files were further processed using myProMS⁹⁹. FDR calculation used Percolator and was set to 1% at the peptide level for the whole study. We validated phosphorylated peptides by combining the phosphoRS information and by manually inspecting the peak assignment.

To quantify the phosphorylated S289 peptide, we extracted from the MS survey of the nanoLC-MS/MS raw files the extracted ion chromatogram (XIC) signal by using the retention time and m/z values of the well characterised tryptic peptide ions using the Xcalibur softwares (manually). XIC areas were integrated in Xcalibur under the QualBrowser interface using the ICIS algorithm. Areas were normalised by using the non-phosphorylated ion's signal¹⁰⁰. Mean values and standard deviation were calculated from three independent experiments.

**Recombinant expression of the mutated HSP90 variants.** The constructs for the expression of recombinant wild type or mutated HSP90 were introduced into *E. coli* BL21(DE3) [pAPlacIQ] and expression was induced using 0.4 mM IPTG at 37 °C for 2 h¹⁰¹. The His-tagged proteins were purified using a Ni-NTA column (Qiagen), as per the manufacturer's instructions. Purity of the proteins was verified by SDS-PAGE and Coomassie Brilliant Blue staining. The purified proteins were used for the kinase assays.

**Expression and purification of recombinant LmCK1.2-V5-His<sub>6</sub>, GSK3 and DYRK.** The *L. major* Friedlin GSK-3 short open reading frame (LmjF18.0270, LmaGSK-3s) was amplified from genomic DNA, cloned into the pTriEx-1.1 vector (Novagen) as previously described⁴⁸,¹⁰², while a synthetic ORF of the *L. infantum* DYRK1 (LinJ.15.0180), optimised for codon usage in *E. coli*, was synthesised in pET14b (EMBL, BamHI/XhoI sites) to encode a His<sub>6</sub>-LinDYRK1 protein (gift from Dr Vinicius Pinto Rocha and Dr Milena Soares). LmCK1.2-V5-His6 (LmjF.35.1010)⁵¹, GSK-3 and DYRK1 were introduced into *E. coli* BL21 Rosetta. The expression of recombinant kinases was induced for 4 hr, with 1 mM IPTG at 37 °C or RT for GSK-3 and DYRK1 respectively, while CK1.2 was expressed with 0.02% (w/v) arabinose at RT as previously described⁵¹. The His-tag proteins were purified using Co-NTA agarose (Pierce) as previously described⁵¹.

**In vitro kinase assay.** In vitro kinase assays were performed as described⁵¹. Briefly, recombinant CK1.2 (rCK1.2), GSK-3 and DYRK1 were assayed in buffer C (60 mM β-glycerophosphate, 30 mM p-nitrophenyl phosphate, 25 mM MOPS, 5 mM EGTA, 15 mM MgCl₂, 1 mM dithiothreitol, 0.1 mM sodium vanadate, pH 7.0) with 15 µM γ-³²P-ATP in 30 µl final volume. For the identification of the CK1.2 phosphosites in HSP90, a similar experiment was performed using CK1.2, recombinant HSP90 and cold ATP at 15 µM. We used the following substrates: 2 µg of either recombinant HSP90, recombinant HSP90 mutants or 4 µg of myelin basic protein (MBP, Sigma-Aldrich), except if stated otherwise. The MBP used was batch-tested previously for absence of autophosphorylation⁵¹. After a 30 min incubation at 30 °C in a 20 µl reaction volume, the reaction was stopped by adding 1 volume 2 × loading buffer¹⁰³. The incorporation of γ-³²P-ATP was monitored by SDS-PAGE and autoradiography. Samples for mass spectrometry analysis were prepared identically, but with non-labelled ATP.

**In silico procedures.** DNA and protein sequence analysis was performed using the MacVector® software (versions 12 to 16). Numerical data were analysed using the Prism® software (version 5). Greyscale and colour images were optimised for contrast using Photoshop® CS3 (Adobe). Composite figures were assembled using the Intaglio® software (Purgatory). For the length measurements of cells, the digital images were imported into the ImageJ 1.42q software (Wayne Rasband, National Institute for Health, USA) analysed using the measurement line tool. Measurements in centimetres were then normalised using the integrated size bars. For the analysis of the cell body lengths 50 randomly selected cells were measured. For the analysis of the cell body widths and flagella length 25 cells were analysed. Significance was determined using the two-sided U-test¹⁰⁴. All greyscale images were cropped from contiguous parts of the original image; no recombination of lanes was performed. Digital enhancements, using Adobe Photoshop CS3, were performed over the entire greyscale images and restricted to tonality (curve) optimisation and size adjustments.

**Data Availability**

All data and materials generated for this study may be obtained from the corresponding author upon written request with the stipulation that any work derived from the data and materials will cite the source.

**References**

1. Alvar, J. et al. Leishmaniasis worldwide and global estimates of its incidence. *PloS one* 7, e35671, https://doi.org/10.1371/journal.pone.0035671 (2012).

2. Barak, E. et al. Differentiation of Leishmania donovani in host-free system: analysis of signal perception and response. *Mol Biochem Parasitol* 141, 99–108 (2005).

3. Brandau, S., Dresel, A. & Clos, J. High constitutive levels of heat-shock proteins in human-pathogenic parasites of the genus *Leishmania*. *Biochem J* 310, 225–232 (1995).

4. Hübel, A., Brandau, S., Dresel, A. & Clos, J. A member of the ClpB family of stress proteins is expressed during heat shock in *Leishmania* spp. *Mol Biochem Parasitol* 70, 107–118 (1995).

5. Krobitsch, S. et al. *Leishmania donovani* heat shock protein 100: characterization and function in amastigote stage differentiation. *J. Biol. Chem.* 273, 6488–6494 (1998).
6. Schlüter, A. et al. Expression and Subcellular Localization of Cpn60 Protein Family Members in *Leishmania donovani*. *Biochim. Biophys. Acta* **1491**, 65–74 (2000).

7. Zamora-Veyl, F. B., Kroemer, M., Zander, D. & Clos, J. Stage-specific expression of the mitochondrial co-chaperonin of Leishmania donovani, CPN10. *Kinetoplastid Biol Dis* **4**, 3 (2005).

8. Hombach, A., Ommen, G., MacDonald, A. & Clos, J. A small heat shock protein is essential for thermotolerance and intracellular survival of Leishmania donovani. *J Cell Sci*, https://doi.org/10.1242/jcs.157297 (2014).

9. Silverman, J. M. et al. Leishmania exosomes modulate innate and adaptive immune responses through effects on monocytes and dendritic cells. *J Immunol* **185**, 5011–5022, https://doi.org/10.4049/jimmunol.1000541 (2010).

10. Silverman, J. M. et al. An exosome-based secretion pathway is responsible for protein export from Leishmania and communication with macrophages. *J Cell Sci* **123**, 842–852, https://doi.org/10.1242/jcs.056465 (2010).

11. Rutherford, S. L. & Zuker, C. S. Protein folding and the regulation of signaling pathways. *Cell* **79**, 1129–1132 (1994).

12. Scheibel, T. & Buchner, J. The Hsp90 complex—a super-chaperone machine as a novel drug target. *Biochem Pharmacol* **56**, 675–682 (1998).

13. Rutherford, S. L. & Lindquist, S. Hsp90 as a capacitor for morphological evolution. *Nature* **396**, 336–342 (1998).

14. Lindquist, S. & Craig, E. A. The heat-shock proteins. *Annual review of genetics* **22**, 631–677 (1988).

15. Whitesell, L., Mimnaugh, E. G., De Costa, B., Myers, C. E. & Neckers, L. M. Inhibition of heat shock protein HSP90-pp60v-src heteroprotein complex formation by benzoquinone ansamycins: essential role for stress proteins in oncogenic transformation. *Proc Natl Acad Sci USA* **91**, 8324–8328 (1994).

16. Whitesell, L. & Cook, P. Stable and specific binding of heat shock protein 90 by geldanamycin disrupts glucocorticoid receptor function in intact cells. *Mol Endocrinol* **10**, 705–712 (1996).

17. Ali, A., Bharadwaj, S., O'Carroll, R. & Ovsenek, N. HSP90 interacts with and regulates the activity of heat shock factor 1 in Xenopus oocytes. *Mol Cell Biol* **18**, 4949–4960 (1998).

18. Kim, H. R., Lee, C. H., Choi, Y. H., Kang, H. S. & Kim, H. D. Geldanamycin induces cell cycle arrest in K562 erythroleukemic cells. *IUBMB Life* **48**, 425–428 (1999).

19. Srethapakdi, M., Liu, F., Tavorath, R. & Rosen, N. Inhibition of Hsp90 function by ansamycins causes retinoblastoma gene product-dependent G1 arrest. *Cancer Res* **60**, 3940–3946 (2000).

20. Wiesgigl, M. & Clos, J. Heat Shock Protein 90 Homeostasis Controls Stage Differentiation in *Leishmania donovani*. *Mol Biol Cell* **12**, 3307–3316 (2001).

21. Graefe, S. E., Wiesgigl, M., Gaworski, I., Macdonald, A. & Clos, J. Inhibition of HSP90 in Trypanosoma cruzi Induces a Stress Response but No Stage Differentiation. *Eukaryot Cell* **1**, 936–943 (2002).

22. Banumathy, G., Singh, V., Pavithra, S. R. & Tatu, U. Heat shock protein 90 function is essential for Plasmodium falciparum growth in human erythrocytes. *The Journal of biological chemistry* **278**, 18336–18345, https://doi.org/10.1074/jbc.M211309200 (2003).

23. Schulte, T. W. et al. Antibiotic radicicol binds to the N-terminal domain of Hsp90 and shares important biologic activities with geldanamycin. *Cell Stress Chaperones* **3**, 100–108 (1998).

24. Schulte, T. W. et al. Interaction of radicicol with members of the heat shock protein 90 family of molecular chaperones. *Mol Endocrinol* **13**, 1435–1448 (1999).

25. Whitesell, L. & Lindquist, S. L. HSP90 and the chaperoning of cancer. *Nature reviews* **5**, 761–772 (2005).

26. Bente, M. et al. Developmentally induced changes of the proteome in the protozoan parasite Leishmania donovani. *Proteomics* **3**, 1811–1829 (2003).

27. Hombach, A., Ommen, G., Chrobak, M. & Clos, J. The Hsp90-StI1 Interaction is Critical for *Leishmania donovani* Proliferation in Both Life Cycle Stages. *Cell Microbiol* **15**, 585–600, https://doi.org/10.1111/cmi.12057 (2013).

28. Prodromou, C. et al. Regulation of Hsp90 ATPase activity by tetratricopeptide repeat (TPR)-domain co-chaperones. *The EMBO journal* **18**, 754–762, https://doi.org/10.1093/emboj/18.3.754 (1999).

29. Johnson, J. L. & Brown, C. Plasticity of the Hsp90 chaperone machine in divergent eukaryotic organisms. *Cell Stress Chaperones* **14**, 83–94 (2009).

30. Longshaw, V. M., Chapple, J. P., Balda, M. S., Cheetham, M. E. & Blatch, G. L. Nuclear translocation of the Hsp70/Hsp90 organizing protein mSTI1 is regulated by cell cycle kinases. *Journal of cell science* **117**, 701–710 (2004).

31. Mollapour, M. & Neckers, L. Post-translational modifications of Hsp90 and their contributions to chaperone regulation. *Biochimica et biophysica acta* **1823**, 648–655, https://doi.org/10.1016/j.bbamcr.2011.07.018 (2012).

32. Muller, P. et al. C-terminal phosphorylation of Hsp70 and Hsp90 regulates alternate binding to co-chaperones CHIP and HOP to determine cellular protein folding/degradation balances. *Oncogene* **32**, 3101–3110, https://doi.org/10.1038/onc.2012.314 (2013).

33. Webb, J. R., Campos-Neto, A., Skeiky, Y. A. & Reed, S. G. Molecular characterization of the heat-inducible LmSTI1 protein of Leishmania major [In Process Citation]. *Mol Biochem Parasitol* **89**, 179–193 (1997).

34. Morales, M. et al. Phosphoproteome dynamics reveals heat shock protein complexes specific to the *Leishmania* infectious stage. *Proc Natl Acad Sci USA* **107**, 8381–8386 (2010).

35. Ommen, G., Chrobak, M. & Clos, J. The co-chaperone SGT of *Leishmania donovani* is essential for the parasite's viability. *Cell Stress and Chaperones* **39**, 541–546, https://doi.org/10.1007/s12192-009-0160-7 (2010).

36. Ommen, G., Lorenz, S. & Clos, J. One-step generation of double-allele gene replacement mutants in Leishmania donovani. *Int J Parasitol* **39**, 541–546 (2009).

37. Hombach, A., Ommen, G., Sattler, V. & Clos, J. Leishmania donovani P23 protects parasites against HSP90 inhibitor-mediated growth arrest. *Cell Stress Chaperones* **20**, 673–685, https://doi.org/10.1007/s12192-015-0595-y (2015).

38. Bartsch, K., Hombach-Barrigah, A. & Clos, J. Hsp90 inhibitors radicicol and geldanamycin have opposing effects on Leishmania Aha1-dependent proliferation. *Cell Stress Chaperones* **22**, 729–742, https://doi.org/10.1007/s12192-017-0800-2 (2017).

39. Morales, M. A. et al. Phosphoproteomic analysis of *Leishmania donovani* pro- and amastigote stages. *Proteomics* **8**, 350–363, https://doi.org/10.1002/pmic.200700697 (2008).

40. Yau, W.-L. et al. Cyclophilin 40-deficient Leishmania donovani fail to undergo stress-induced development of the infectious metacyclic stage. *Cell. Microbiol.* **93**, 80–97, https://doi.org/10.1111/mmi.12639 (2014).

41. Zuehlke, A. D., Beebe, K., Neckers, L. & Prince, T. Regulation and function of the human HSP90AA1 gene. *Gene* **570**, 8–16, https://doi.org/10.1016/j.gene.2015.06.018 (2015).

42. Zhao, Y. G. et al. Hsp90 phosphorylation is linked to its chaperoning function. *Assembly of the reovirus cell attachment protein. J Biol Chem* **276**, 32822–32827, https://doi.org/10.1074/jbc.M105562200 (2001).

43. Mollapour, M., Tsutsumi, S., Kim, Y. S., Trepel, J. & Neckers, L. Casein kinase 2 phosphorylation of Hsp90 threonine 22 modulates chaperone function and drug sensitivity. *Oncotarget* **2**, 407–417 (2011).

44. Mollapour, M. et al. Threonine 22 phosphorylation attenuates Hsp90 interaction with cochaperones and affects its chaperone activity. *Molecular cell* **41**, 672–681, https://doi.org/10.1016/j.molcel.2011.02.011 (2011).

45. Xu, W. et al. Dynamic tyrosine phosphorylation modulates cycling of the HSP90-P50(CDC37)-AHA1 chaperone machine. *Molecular cell* **47**, 434–443, https://doi.org/10.1016/j.molcel.2012.05.015 (2012).

46. Grant, K. M. et al. Inhibitors of *Leishmania mexicana* CRK3 cyclin-dependent kinase: chemical library screen and antileishmanial activity. *Antimicrob Agents Chemother* **48**, 3033–3042, https://doi.org/10.1128/AAC.48.8.3033-3042.2004 (2004).

47. Wiese, M. *Leishmania* MAP kinases—familiar proteins in an unusual context. *Int J Parasitol* **37**, 1053–1062, https://doi.org/10.1016/j.ijpara.2007.04.008 (2007).
48. Xingi, E. et al. 6-Br-5methylindirubin-3'oxime (5-Me-6-BIO) targeting the leishmanial glycogen synthase kinase-3 (GSK-3) short form affects cell-cycle progression and induces apoptosis-like death: exploitation of GSK-3 for treating leishmaniasis. *Int J Parasitol* **39**, 1289–1303, https://doi.org/10.1016/j.ijpara.2009.04.005 (2009).

49. Spath, G. F., Drini, S. & Rachidi, N. A touch of Zen: post-translational regulation of the Leishmania stress response. *Cell Microbiol* **17**, 632–638, https://doi.org/10.1111/cmi.12440 (2015).

50. Duncan, S. M. et al. Conditional gene deletion with DiCre demonstrates an essential role for CRK3 in Leishmania mexicana cell cycle regulation. *Mol Microbiol* **100**, 931–944, https://doi.org/10.1111/mmi.13375 (2016).

51. Rachidi, N. et al. Pharmacological assessment defines Leishmania donovani casein kinase 1 as a drug target and reveals important functions in parasite viability and intracellular infection. *Antimicrob Agents Chemother* **58**, 1501–1515, https://doi.org/10.1128/AAC.02022-13 (2014).

52. Budini, M. et al. Autophosphorylation of carboxy-terminal residues inhibits the activity of protein kinase CK1alpha. *J Cell Biochem* **106**, 399–408, https://doi.org/10.1002/jcb.22019 (2009).

53. Vieira, L. L., Sacerdoti-Sierra, N. & Jaffe, C. L. Effect of pH and temperature on protein kinase release by Leishmania donovani. *Int J Parasitol* **32**, 1085–1093 (2002).

54. Silverman, J. M. et al. Proteomic analysis of the secretome of Leishmania donovani. *Genome biology* **9**, R35 (2008).

55. Liu, J. et al. Mammalian casein kinase 1alpha and its leishmanial ortholog regulate stability of IFNAR1 and type I interferon signaling. *Mol Cell Biol* **29**, 6401–6412, https://doi.org/10.1128/MCB.00478-09 (2009).

56. Wiese, M. A mitogen-activated protein (MAP) kinase homologue of Leishmania mexicana is essential for parasite survival in the infected host. *Embo J* **17**, 2619–2628 (1998).

57. Kaur, P., Garg, M., Hombach-Barrigah, A., Clos, J. & Goyal, N. MAPK1 of Leishmania donovani interacts and phosphorylates HSP70 and HSP90 subunits of foldosome complex. *Sci Rep* **7**, 10202, https://doi.org/10.1038/s41598-017-09725-w (2017).

58. Tsigankov, P. et al. Regulation dynamics of Leishmania differentiation: deconvoluting signals and identifying phosphorylation trends. *Mol Cell Proteomics* **13**, 1787–1799, https://doi.org/10.1074/mcp.M114.037705 (2014).

59. Zilka, A., Garlapati, S., Dahan, E., Yaolsky, V. & Shapira, M. Developmental regulation of HSP83 in Leishmania: transcript levels are controlled by the efficiency of 3' RNA processing and preferential translation is directed by a determinant in the 3'UTR. *J Biol Chem* **111**, 11 (2001).

60. Hawle, P. et al. The middle domain of Hsp90 acts as a discriminator between different types of client proteins. *Mol Cell Biol* **26**, 8385–8395, https://doi.org/10.1128/MCB.02188-05 (2006).

61. Bifeld, E. et al. Ribosome Profiling Reveals HSP90 Inhibitor Effects on Stage-specific Protein Synthesis in Leishmania donovani. *mSystems* accepted.

62. Sacks, D. L. & Perkins, P. V. Identification of an infective stage of Leishmania promastigotes. *Science* **223**, 1417–1419 (1984).

63. Bifeld, E., Tejera Nevado, P., Bartsch, J., Eick, J. & Clos, J. A versatile qPCR assay to quantify trypanosomatidic infections of host cells and tissues. *Med Microbiol Immunol* **205**, 449–458, https://doi.org/10.1007/s00430-016-0460-3 (2016).

64. Woods, Y. L. et al. The kinase DYRK phosphorylates protein-synthesis initiation factor eIF2Bepsilon at Ser539 and the microtubule-associated protein tau at Thr212: potential role for DYRK as a glycogen synthase kinase 3-priming kinase. *Biochem J* **355**, 609–615 (2001).

65. Ojo, K. K. et al. Structure determination of glycogen synthase kinase-3 from Leishmania major and comparative inhibitor structure-activity relationships with Trypanosoma brucei GSK-3. *Mol Biochem Parasitol* **176**, 98–108, https://doi.org/10.1016/j. molbiopara.2010.12.009 (2011).

66. Knippschild, U. et al. The CK1 Family: Contribution to Cellular Stress Response and Its Role in Carcinogenesis. *Front Oncol* **4**, 96, https://doi.org/10.3389/fonc.2014.00096 (2014).

67. Marin, O. et al. A noncanonical sequence phosphorylated by casein kinase 1 in beta-catenin may play a role in casein kinase 1 targeting of important signaling proteins. *Proc Natl Acad Sci USA* **100**, 10193–10200, https://doi.org/10.1073/pnas.1733909100 (2003).

68. Kawakami, F., Suzuki, K. & Ohtsuki, K. A novel consensus phosphorylation motif in sulfatide- and cholesterol-3-sulfate-binding protein substrates for CK1 *in vitro*. *Biol Pharm Bull* **31**, 193–200 (2008).

69. Amanchy, R. et al. Identification of Novel Phosphorylation Motifs Through an Integrative Computational and Experimental Analysis of the Human Phosphoproteome. *J Proteomics Bioinform* **4**, 22–35, https://doi.org/10.4172/jpb.1000163 (2011).

70. Keller, M., Chan, R. L., Tessier, L. H., Weil, J. H. & Imbault, P. Post-transcriptional regulation by light of the biosynthesis of Euglena ribulose-1,5-bisphosphate carboxylase/oxygenase small subunit. *Plant Mol Biol* **17**, 73–82 (1991).

71. Clayton, C. E. Life without transcriptional control? From fly to man and back again. *Embo J* **21**, 1881–1888 (2002).

72. Martinez-Calvillo, S., Vizuet-de-Rueda, J. C., Florencio-Martinez, L. E., Manning-Cela, R. G. & Figueroa-Angulo, E. E. Gene expression in trypanosomatid parasites. *Journal of biomedicine & biotechnology* **2010**, 525241, https://doi.org/10.1155/2010/525241 (2010).

73. Liang, X. H., Haritan, A., Uliel, S. & Michaeli, S. trans and cis splicing in trypanosomatids: mechanism, factors, and regulation. *Eukaryot Cell* **2**, 830–840 (2003).

74. Mollapour, M. et al. Asymmetric Hsp90 N domain SUMOylation recruits Aha1 and ATP-competitive inhibitors. *Mol Cell* **53**, 317–329, https://doi.org/10.1016/j.molcel.2013.12.007 (2014).

75. Xu, Y., Singer, M. A. & Lindquist, S. Maturation of the tyrosine kinase c-src as a kinase and as a substrate depends on the molecular chaperone Hsp90. *Proc Natl Acad Sci USA* **96**, 109–114 (1999).

76. Zuehlke, A. D. et al. An Hsp90 co-chaperone protein in yeast is functionally replaced by site-specific posttranslational modification in humans. *Nat Commun* **8**, 15328, https://doi.org/10.1038/ncomms15328 (2017).

77. Legagneux, V., Morange, M. & Bensaude, O. Heat shock increases turnover of 90 kDa heat shock protein phosphate groups in HeLa cells. *FEBS Lett* **291**, 359–362 (1991).

78. Mollapour, M., Tsutsumi, S. & Neckers, L. Hsp90 phosphorylation, Wee1 and the cell cycle. *Cell Cycle* **9**, 2310–2316, https://doi.org/10.4161/cc.9.12.12054 (2010).

79. Nathan, D. F. & Lindquist, S. Mutational analysis of Hsp90 function: interactions with a steroid receptor and a protein kinase. *Mol Cell Biol* **15**, 3917–3925 (1995).

80. Howard, M. K., Sayers, G. & Miles, M. A. Leishmania donovani metacyclic promastigotes: transformation *in vitro*, lectin agglutination, complement resistance, and infectivity. *Exp Parasitol* **64**, 147–156 (1987).

81. Pearl, L. H. & Prodromou, C. Structure and mechanism of the Hsp90 molecular chaperone machinery. *Annu Rev Biochem* **75**, 271–294, https://doi.org/10.1146/annurev.biochem.75.103004.142738 (2006).

82. Yau, W. L. et al. Phenotypic Characterization of a Leishmania donovani Cyclophilin 40 Null Mutant. *J Eukaryot Microbiol* **63**, 823–833, https://doi.org/10.1111/jeu.12329 (2016).

83. Li, J., Richter, K. & Buchner, J. Mixed Hsp90-cochaperone complexes are important for the progression of the reaction cycle. *Nat Struct Mol Biol* **18**, 61–66, https://doi.org/10.1038/nsmb.1965 (2011).

84. Sollelis, L. et al. First efficient CRISPR-Cas9-mediated genome editing in Leishmania parasites. *Cell Microbiol* **17**, 1405–1412, https://doi.org/10.1111/cmi.12456 (2015).

85. Zhang, W. W. & Matlashewski, G. CRISPR-Cas9-Mediated Genome Editing inLeishmania donovani. *MBio* **6**, e00861, https://doi.org/10.1128/mBio.00861-15 (2015).

86. Cohen, P. The regulation of protein function by multisite phosphorylation-a 25 year update. *Trends Biochem Sci* **25**, 596–601 (2000).
87. Lahav, T. et al. Multiple levels of gene regulation mediate differentiation of the intracellular pathogen Leishmania. *FASEB journal: official publication of the Federation of American Societies for Experimental Biology* **25**, 515–525, https://doi.org/10.1096/fj.10-157529 (2011).

88. Rosenzweig, D. et al. Retooling Leishmania metabolism: from sand fly gut to human macrophage. *FASEB J* **22**, 590–602, https://doi.org/10.1096/fj.07-9254com (2008).

89. Bifeld, E. In *Leishmania: Methods and Protocols Methods in Molecular Biology* (ed Joachim Clos) Ch. 12, in press (Springer Nature, 2019).

90. Sambrook, J. & Russell, D. W. *Molecular Cloning*. 3rd edn, 1.55 ff (Cold Spring Harbor Laboratory Press, 2001).

91. Schlüter, A. et al. Expression and Subcellular Localization of Cpn60 Protein Family Members in *Leishmania donovani*. *Biochim. Biophys. Acta* **1491**, 65–74 (2000).

92. Choudhury, K., Zander, D., Kube, M., Reinhardt, R. & Clos, J. Identification of a Leishmania infantum gene mediating resistance to miltefosine and SbIII. *Int J Parasitol* **38**, 1411–1423 (2008).

93. Brun, R. & Schonenberger. Cultivation and *in vitro* cloning or procyclic culture forms of Trypanosoma brucei in a semi-defined medium. *Short communication. Acta Trop* **36**, 289–292 (1979).

94. Jensen, S. S. & Larsen, M. R. Evaluation of the impact of some experimental procedures on different phosphopeptide enrichment techniques. *Rapid Commun Mass Spectrom* **21**, 3635–3645, https://doi.org/10.1002/rcm.3254 (2007).

95. Rappsilber, J., Ryder, U., Lamond, A. I. & Mann, M. Large-scale proteomic analysis of the human spliceosome. *Genome Res* **12**, 1231–1245, https://doi.org/10.1101/gr.473902 (2002).

96. Mortensen, P. et al. MSQuant, an open source platform for mass spectrometry-based quantitative proteomics. *Journal of proteome research* **9**, 393–403, https://doi.org/10.1021/pr900721e (2010).

97. Beausoleil, S. A., Villen, J., Gerber, S. A., Rush, J. & Gygi, S. P. A probability-based approach for high-throughput protein phosphorylation analysis and site localization. *Nat Biotechnol* **24**, 1285–1292, https://doi.org/10.1038/nbt1240 (2006).

98. Savitsky, M. et al. Recording information on protein complexes in an information management system. *J Struct Biol* **175**, 224–229, https://doi.org/10.1016/j.jsb.2011.05.009 (2011).

99. Poullet, P., Carpentier, S. & Barillot, E. myProMS, a web server for management and validation of mass spectrometry-based proteomic data. *Proteomics* **7**, 2553–2556, https://doi.org/10.1002/pmic.200600784 (2007).

100. Valot, B., Langella, O., Nano, E. & Zivy, M. MassChroQ: a versatile tool for mass spectrometry quantification. *Proteomics* **11**, 3572–3577, https://doi.org/10.1002/pmic.201100120 (2011).

101. Clos, J. & Brandau, S. pJC20 and pJC40— two high-copy-number vectors for T7 RNA polymerase-dependent expression of recombinant genes in *Escherichia coli*. *Prot. Expression Purif.* **5**, 133–137 (1994).

102. Efstathiou, A. et al. An inhibitor-driven study for enhancing the selectivity of indirubin derivatives towards leishmanial Glycogen Synthase Kinase-3 over leishmanial cdc2-related protein kinase 3. *Parasites & vectors* **7**, 234, https://doi.org/10.1186/1756-3305-7-234 (2014).

103. Laemmli, U. K. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. *Nature* **227**, 680–685 (1970).

104. Mann, H. B. & Whitney, D. R. On a test of whether one of two random variables is stochastically larger than the other. *Ann Math Statistics* **18**, 50–60 (1947).

## Acknowledgements

We acknowledge Sandra Arriens for preparing the kinase substrate proteins and Anna-Lena Thiel for a first analysis of two phosphosite mutants. We are indebted to Luke Whitesell for crucial advice at the start of the project. We thank former and current members of our laboratories for helpful discussions and suggestions. We thank Dr Vinicius Pinto Costa Rocha and Dr Milena Botelho Pereira Soares [Instituto Gonçalo Moniz, Fundação Oswaldo Cruz (FIOCRUZ), Hospital São Rafael, Salvador, Bahia, Brazil] for providing the DYRK1 plasmid. A.H.-B. and K.B. were supported by a grant from the Deutsche Forschungsgemeinschaft (CL120/8-1). Despina Smirlis and Najma Rachidi were supported by the ANR-13-ISV3-0009 grant. This work was supported by “Région Ile-de-France” and FRM grants (D.L.).

## Author Contributions

A.H.-B. Creation of point mutants, expression and phenotype analysis in *Leishmania*, experimental planning, manuscript preparation. K.B. Creation of point mutants, expression and phenotype analysis in *Leishmania*, manuscript preparation. D.S. Production of Dyrk1 and GSK3, manuscript preparation. A.M. Scanning electron microscopy. H.R. Design and execution of phosphoproteomics. F.D. Mass spectrometry of the label free kinase quantification. D.L. Mass spectrometry of the label free Kinase quantification and study design. G.S. Design of study, supervision, manuscript preparation. N.R. Protein kinase assays, design and supervision, manuscript preparation. M.W. Design and supervision of phosphoproteomics, manuscript preparation. J.C. Design of study, experimental planning, overall supervision, manuscript preparation.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-41640-0.

**Competing Interests:** The authors declare 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_10.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 Creative 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 permitted 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) 2019