Quantitative Analyses Reveal Calcium-dependent Phosphorylation Sites and Identifies a Novel Component of the Invasion Motor Complex

Summary

Apicomplexan parasites depend on the invasion of host cells for survival and proliferation. Calcium-dependent signaling pathways appear to be essential for micronemal release and gliding motility, yet the target of activated kinases remains largely unknown. We have characterized calcium-dependent phosphorylation events during Toxoplasma host cell invasion. Stimulation of live tachyzoites with Ca²⁺-mobilizing drugs leads to phosphorylation of numerous parasite proteins, as shown by differential 2-DE display of ³²[P]-labeled protein extracts. Multi-dimensional Protein Identification Technology (MudPIT) identified ∼546 phosphorylation sites on over 300 Toxoplasma proteins, including 10 sites on the actomyosin invasion motor. Using a Stable Isotope of Amino Acids in Culture (SILAC)-based quantitative LC-MS/MS analyses we monitored changes in the abundance and phosphorylation of the invasion motor complex and defined Ca²⁺-dependent phosphorylation patterns on three of its components - GAP45, MLC1 and MyoA. Furthermore, calcium-dependent phosphorylation of six residues across GAP45, MLC1 and MyoA is correlated with invasion motor activity. By analyzing proteins that appear to associate more strongly with the invasion motor upon calcium stimulation we have also identified a novel 15-kDa Calmodulin-like protein that likely represents the MyoA Essential Light Chain of the Toxoplasma invasion motor. This suggests that invasion motor activity could be regulated not only by phosphorylation but also by the direct binding of calcium ions to this new component.

Introduction

The phylum Apicomplexa is a large group of obligate intracellular parasites of wide medical and agricultural significance. Alone, Toxoplasma gondii infects between 30- 80% of people worldwide and is one the most common infectious agents of humans. Toxoplasma transmission occurs by exposure to the feces of an infected cat, eating undercooked meat or ingestion of contaminated water harboring oocysts [1]. An ocular infection route is also common and is a leading cause of blindness in some countries [2]. Primary exposure or reactivation of tissue cysts in pregnant women can lead to congenital birth defects and spontaneous abortion, while toxoplasmosis is a common secondary infection of AIDS patients and other immuno-compromised individuals and can lead to death if untreated.

Invasion of host cells by apicomplexan parasites is an obligatory step for their survival and proliferation. Despite each species having a range of host and cell types that they target the intracellular processes governing invasion appear to be largely conserved. Toxoplasma’s ease of growth in the laboratory and high genetic tractability has provided researchers with an excellent model for investigating the molecular basis of invasion in related species such as Plasmodium spp - the causative agent of malaria. For example, a highly conserved actomyosin-based invasion motor that drives parasite motility and invasion was first identified and has been largely characterized in Toxoplasma. At its core the invasion motor consists of a novel class XIV Myosin, MyoA [3] and Myosin Light Chain 1, MLC1 [4] which is anchored into the outer side of the inner membrane complex (IMC) by the Glideosome-Associated Proteins GAP40, GAP45, GAP50 and GAPM’s [5], [6], [7], [8]. Recently, the architecture of the invasion motor has been mapped and the central role and absolute requirement of GAP45 has been demonstrated [7]. Furthermore, it was shown that GAP45 spans the supra-alveolar space providing cohesion between the IMC and plasma membrane [7]. The current model of invasion suggests that upon host cell contact GAP40/45/50, GAPM’s, MyoA and MLC1 complex and filamentous actin forms [9]. Transmembrane host cell adhesins linked to actin filaments through the glycolytic enzyme aldolase are then pulled rearwards by the action of MyoA, thus driving the parasite forward into the host cell [10].

The release of apical organelles is also required for successful invasion. Micronemes are first to be secreted and contain high affinity host cell ligands that are necessary for strong host cell attachment and invasion [11]. Rhoptries, large club-shaped organelles, rich in lipids are released after micronemes upon contact of the apical end of the parasite with the host cell. Rhoptry release is correlated with host membrane penetration and the creation of the parasitophorous vacuole [12]. The molecular machinery that drives exocytosis of these apical organelles is unknown.

Activation of the invasion motor and regulated release of apical organelles occurs after a change in extracellular environment and contact with the host cell [13]. Pharmacologically, calcium-dependent signal transduction pathways have been identified in Toxoplasma and P. falciparum and appear to be essential for micronemal release, motility and invasion [14], [15]. Calcium mobilizing agents potently induce exocytosis of micronemes and gliding motility, whereas membrane permeable calcium chelators such as BAPTA-AM inhibit these processes [11], [15]. In addition, ethanol induces a transient [Ca²⁺]_i increase and micronemal release via a putative signaling pathway involving the activation of phospholipase C (PLC) [16]. Further supporting a role of calcium signaling in apicomplexan parasite invasion is the visualization of intracellular calcium ([Ca²⁺]_i) oscillations that accompany invasion processes in Toxoplasma [14], [17] and Plasmodium [15]. Together, this data indicates that [Ca²⁺]_i signals mediate multiple signaling pathways, including micronemal release and gliding motility and these are essential for regulating the processes involved in host cell invasion.

In mammalian cells, spatiotemporally distinct [Ca²⁺]_i signals are translated into enzymatic activity using a variety of serine/threonine-specific protein kinases and phosphatases, which control the function of downstream effectors by promoting local changes in protein conformation, interactions and activity. In Apicomplexa, a conserved family of plant-like Calcium-Dependent Protein Kinases (CDPKs) has been strongly implicated in mediating several critical Ca²⁺-dependent signal transduction pathways during the complex parasite life cycles [18]. Recent chemical genetic and conditional gene knockout approaches have demonstrated that Toxoplasma CDPK1 is an essential regulator of Ca²⁺-dependent micronemal exocytosis [19], [20], whereas P. falciparum CDPK1 (PfCDPK1) appears to play a role in regulating red blood cell invasion [21]. Furthermore, PfCDPK1 phosphorylates P. falciparum Myosin A Tail Interacting Protein (PfMTIP) and Glideosome-Associated Protein 45 (PfGAP45) in vitro [21], [22], which suggests a role in regulating invasion motor function during host cell invasion. Together, this data suggests that CDPK’s play crucial roles in regulating invasion processes, yet the in vivo substrates of these kinases remain largely unexplored.

To identify substrates of calcium-dependent kinases we have employed proteomics approaches to obtain a global snapshot of the phosphorylation pattern of Toxoplasma proteins following stimulation of Ca²⁺ signaling pathways. This identified targets of Ca²⁺-dependent phosphorylation pathways potentially involved in regulating invasion processes. Phosphorylation sites on Toxoplasma GAP45, MLC1 and MyoA were identified and given that their is calcium-dependent phosphorylation deposition on only some sites and the presence of a phospho-tyrosine it appears that the phosphorylation status of the invasion motor is likely regulated by multiple kinases. Further, a novel protein was identified as a new component of the Toxoplasma invasion motor that associates more tightly upon calcium signaling suggesting that the invasion motor may also be directly regulated by Ca²⁺ binding.

Results

In vivo stimulation of intracellular Ca²⁺ signaling pathways triggers phosphorylation of multiple proteins

To identify Ca²⁺-dependent phospho-substrates during Toxoplasma invasion we conducted a 2-DE-based screen by labeling parasites with ³²[P] orthophosphate and stimulating with either ethanol (thought to activate PLC) or the calcium ionophore ionomycin. Both calcium pathway agonists have previously been shown to induce micronemal secretion and motility after short periods of stimulation [16]. As a control parasite samples were either mock treated (DMSO alone) or pretreated with the membrane permeable calcium chelator BAPTA-AM followed by stimulation with ethanol. Extracts were separated on 2-DE gels and over 50 ³²[P]-labeled protein spots were identified that were absent or barely detectable in BAPTA-AM or mock-treated negative controls (Figure 1A, C) and exclusively present in samples stimulated with either ethanol or ionomycin (Figure 1B, D – marked with red arrows). False colour 2-DE image overlays of corresponding autoradiographs highlight ³²[P]-labeled protein spots that were specifically phosphorylated following stimulation with Ca²⁺ mobilizing drugs (Figure 1E). Spots that appear yellow likely represent Ca²⁺-insensitive phosphoproteins or GPI-anchored parasite proteins which are unaffected by stimulation, whereas spots that appear red represent potential substrates of calcium-dependent phosphorylation (Figure 1E). Interestingly, the same proteins appeared to be phosphorylated when stimulated by either calcium ionophore or ethanol (Figure 1B, D), strongly suggesting that both agonists stimulate the same signal transduction pathway(s) in Toxoplasma tachyzoites.

A total of 50 2-DE spots that were consistently phosphorylated following Ca²⁺ pathway stimulation were subjected to LC-MS/MS-based protein identification (Figures 1E, Supplementary Figure S1 and Table S1). In some instances proteins from separate 2-DE spots shared the same identity, as indicated by bounding boxes (Figure 1E, Table S1). Assuming that phosphorylated 2-DE spots are the product of a single gene the present data identified a total of 63 Toxoplasma proteins - some of which represent molecules potentially involved in mediating intracellular signaling cascades, regulating exocytosis of invasion organelles, or controlling parasite motility (Table S1). These included amongst others, a cAMP kinase regulatory subunit (PKA-R, TGME49_042070), a casein kinase subunit II beta subunit (CKIIß, TGME49_072400), a soluble NSF-attachment protein (SNAP, TGME49_018760), an armadillo repeat-containing protein (ARM1, TGME49_061440), a serine/threonine phosphatase 2C (PP2C, TGME49_054770), a Rab GTPase (RAB5B, TGME49_007460), two highly acidic proteins similar to I1PP2A/ANP32 (ANP32; TGME49_071810) or I2PP2A/SET (SET; TGME49_044110), as well as a number of putative uncharacterized proteins (ie. Hypo) (Figure 1E). Importantly, previously identified phosphoproteins were also identified in our analysis including Toxofilin (TGME49_014080), an actin binding protein shown to be phosphorylated by parasitic type 2C phosphatase and CKII activities [23] (Figure 1E, Table S1). Furthermore, we identified ³²[P]-labeled spots representing GAP45 and MLC1 (Figure 1E, Table S1), indicating that these components of the invasion motor are also likely targets of calcium-dependent phosphorylation in vivo.

A limitation in our protein identification methodology is that highly abundant proteins that co-migrate with our calcium-dependent phosphoproteins will be preferentially identified. We therefore sought to confirm the Ca²⁺-dependent phosphorylation of some of the interesting putative phosphoproteins that we identified. To do this we tagged candidates with a C-terminal triple HA tag. Upon radiolabelling parasites and immunoprecipitation with HA antibodies we assessed the phosphorylation status of proteins during calcium stimulation, using ethanol alone or in the presence of BAPTA-AM (Figure 1F). This confirmed that HA-tagged Toxoplasma PKA-R, PP2C, ARM1 and MLC1 are in-vivo phospho-substrates (Figure 1F). As a control we compared the signal intensity of immunoreactive bands detected on Western blots probed with anti-HA antibody. The ³²[P]-labeling intensity of the corresponding HA-PKA-R PP2C-HA, ARM1-HA or MLC1-HA protein bands was significantly lower in samples treated with BAPTA-AM prior to stimulation with ethanol (Figure 1F). This is consistent with the increased number and intensity of ³²[P]-labeled 2-DE spots for native PKA-R, PP2C, ARM1 and MLC1 observed by 2-DE display (Figure 1E) and validates that the phosphorylation of these proteins is dependent on intracellular Ca²⁺ flux.

Focusing in on MLC1 and GAP45 we could see in our 2-DE autoradiographs that both motor components are basally phosphorylated and upon calcium stimulation additional phosphorylated isoforms are detectable (Figure 1E). Multiple ³²P-labeled protein bands seen in immunoprecipitates from parasite lines expressing HA-tagged MLC1 are consistent with this protein co-purifying with other phosphorylated components of the native invasion motor complex (Figure 1F). Overall, our ³²[P]/2DE analysis has identified a suite of validated and potential phosphoproteins that respond to calcium signaling and further, identified the Toxoplasma invasion motor as a substrate of both calcium-dependent and independent phosphorylation events.

Identification of phosphorylation sites for Toxoplasma calcium-dependent phosphoproteins

Protein phosphorylation often occurs at low stoichiometry and/or involves proteins with low expression levels. Indeed, some of the proteins identified from ³²[P]/2DE spots could be more abundant co-migrating proteins that are not responsible for the signal in the autoradiographs. To overcome this issue we employed a Multi-dimensional liquid chromatography Protein Identification Technology (MudPIT)-based strategy for large-scale analysis of phosphorylation sites following in vivo Ca²⁺ pathway stimulation of Toxoplasma tachyzoites. Proteins were extracted from ethanol-stimulated parasites and digested with trypsin, fractionated via hydrophilic interaction chromatography (HILIC) and partitioned via titanium dioxide (TiO₂) affinity chromatography. The resulting unbound or phosphopeptide-enriched (bound) fractions were subjected to MudPIT analyses on a LTQ linear ion trap mass spectrometer. Detailed analyses of TiO₂-bound LC-MS/MS spectra using both Sequest and Mascot-based search algorithms resulted in a total of 305 non-redundant Toxoplasma phosphoproteins at 0.9% false discovery rate (FDR) supported by 496 manually approved unique phosphopeptide identifications (Figure 2, Supplementary Tables S2–S4 and Text S1). This data has been integrated into ToxoDB (www.toxodb.org) and will be available in the September 2011 release.

Combining 2-DE analyses of intact ³²[P]-labeled proteins with global MudPIT phosphoproteomics enabled us to cross-correlate our analysis for the targeted study of parasite proteins modified by phosphorylation following in vivo stimulation of Ca²⁺ signaling pathways. This identified 10 major Ca²⁺-dependent tachyzoite phosphoproteins represented by 21 valid phosphoprotein spots and 86 phosphopeptide spectra (Figure 2A). The table in Figure 2B summarises the 2-DE and MudPIT LC-MS/MS data for this subset of the parasite phosphoproteome, which includes 17 putative Ca²⁺-dependent phosphorylation sites. In agreement with differential 2-DE display data (Figure 1E), the existing phosphopeptide MS evidence supports the presence of multiple phosphorylation sites on some phosphoproteins in more than one ³²[P]-labeled 2-DE spot (Figure 2B and Supplementary Table S1). This includes a cluster of five residues localized on GAP45 (S153, Y158, S163, S167, S184) (Figure S2, A–E) and at least two phosphorylation sites each for PKA-R (S27, S94) (Figure S2, G–H), CKIIß (S22, S309) (Figure S2, I–J), and a hypothetical protein of unknown function TGME49_048700 (S249, S250) (Figure S2, N-O). By contrast, we only detected MS evidence for a single phosphorylation site on MLC1 (S98) (Supplementary Figure S2F), ARM1 (S33) (Supplementary Figure S2M), ANP32 (S187), SET (S103) and SNAP (S6) (Supplementary Figure S2, K–L), most of which were authenticated in multiple ³²[P]-labeled 2-DE spots (Figures 1E, 2B). We failed to observe any phosphopeptide evidence for PP2C by LC-MS/MS even though our IP analyses of ³²[P]-labeled parasites expressing transgenic HA-tagged suggested that this protein is a bone fide phosphoprotein (Figures 1F, 2B).

The Toxoplasma invasion motor is heavily phosphorylated

Complementary 2-DE and MudPIT analyses of whole parasite extracts provided direct evidence that at least 2 of the 5 known myosin motor complex components, namely GAP45 and MLC1, are phosphorylated in vivo in a Ca²⁺-dependent manner. To further understand this we employed a Stable Isotope Labeling by Amino Acids in Culture (SILAC)-based proteomics approach and quantitatively assess the in vivo phosphorylation of tachyzoite invasion motor complexes following Ca²⁺ pathway stimulation (Figure 3A). First, we demonstrated that over the labeling period we achieved >95% incorporation of ‘heavy’ stably isotopically labeled arginine and lysine into tachyzoite protein, as detailed in Materials & Methods and Supplementary Figure S3. Then we grew SILAC-labeled bulk cultures and prepared protein lysates of un-stimulated (light) or ethanol-stimulated (heavy) tachyzoites. Samples were quantified using Sypro Ruby protein stain (Figure 3B, lanes 1&2), total protein was mixed at a ratio of 1∶1 and used to immunoaffinity-purify Toxoplasma myosin motor complexes using anti-GAP45 specific antibodies (Figure 3B, lane 3). As expected anti-GAP45 column eluates were specifically enriched in the five major components of the Toxoplasma glideosome complex – MyoA, GAP50, GAP45, GAP40 and MLC1 (Figure 3B) [5], [7]. Intact protein complexes were then digested with trypsin and partitioned using TiO₂-based phosphopeptide enrichment [24]. TiO₂ eluates and flow through were then run through the high mass accuracy Orbitrap instrument. (Phospho)peptides were then identified and relative abundances calculated using the MaxQuant computational platform [25] (Figure 3A). From this analysis we were able to achieve 70.8%, 36.9%, 52.2%, 3.7% and 73.7% sequence coverage across MyoA, GAP50, GAP45, GAP40 and MLC1, respectively (Supplementary Figure S4).

A total of 61 Toxoplasma phosphopeptide sequences were identified in digests of immunoaffinity-purified invasion motor complexes (Supplementary Table S5). They were exclusively singly phosphorylated and recovered in the TiO₂-bound fraction consistent with a specific enrichment of mono-phosphorylated peptides by TiO₂ beads [26]. Computational analyses of Orbitrap LC-MS/MS data identified up to 13 potential phosphorylation sites on the invasion motor complex proteins GAP45, MLC1 and MyoA, but none on GAP50 (Table S5). Evaluation of evidence spectra resulted in manual approval of ten unambiguous phosphorylation sites with highly confident Mascot (>25), ΔMascot (>10) and MaxQuant Phospho S/T/Y site probability scores (>90%) (Table 1 and Supplementary Figure S5, A–J). The phosphorylation of neighboring GAP45 residues S184/S185 was the only uncertainty and was narrowed to either site with 50% probability (Table 1 and Supplementary Figure S5E).

The OrbiTrap SILAC data was further analyzed to quantitate the response of calcium signaling on invasion motor phosphorylation. We calculated the ratio of Heavy:Light phosphopeptides normalized for particular components using the MaxQuant algorithm (Table 1, Supplementary Table S5). This showed that phosphorylation of 4 of the identified residues on GAP45 (S153, Y158, S163, S169) did not change upon calcium stimulation, whereas S184/185 and T189 had H:L ratios of 3.9 and 3.7 respectively (Table 1), strongly suggesting that phosphorylation on these sites is deposited by a calcium-dependent protein kinase (the quantitation of Ca²⁺-dependent GAP45 phosphorylation was manually validated as shown in Supplementary Figure S6, and a detailed discussion of results is provided in Supplementary Text S2). For MLC1 we could determine that while one site (S55) did not respond to calcium signaling, T98 and S132 had mild increases in phosphorylation deposition, again suggesting a role of a calcium-dependent kinase in modulating the post-translationally modified state of this protein (the quantitation of Ca²⁺-dependent MLC1 phosphorylation was manually validated as shown in Supplementary Figure S7 and is discussed in Supplementary Text S3). The only MyoA phosphorylation site that we were able to identify appeared to change negligibly upon calcium stimulation and further work will be needed to identify if S21 phosphorylation truly responds to calcium signaling (Table 1). Given that we identified both serine/threonine and tyrosine phosphorylation and that some residues respond to calcium signaling, while others do not, this suggests that at least three kinases are responsible for the phosphorylation profile of the Toxoplasma invasion motor.

Quantitative proteomics suggests that calcium signaling induces changes in invasion motor assembly

Given that the invasion motor must form a complex to promote parasite motility [7] we wanted to see if we could use our SILAC approach to see what effect calcium signaling has on motor assembly. To do this we analyzed the abundance of non-phosphopeptides from motor complex components and in their different labeling states from anti-GAP45 column eluates. Using our LC-MS/MS analyses on the Orbitrap we identified abundant MyoA, GAP50, GAP45, GAP40 and MLC1 peptides, with a total of 186, 58, 54 or 43 assigned spectra from six independent experiments (including different modifications and/or labeling states) (Table 2 and Supplementary Table S6). Analysis of Heavy:Light ratios of non-phosphopeptides indicated a 2–3 fold increase in the relative abundance of heavy-labeled invasion motor complex components from calcium-stimulated parasites (Table 2 and Supplementary Table S7). By contrast, the average H/L ratio of total Toxoplasma protein in these samples was estimated to be 1.1 ± 0.9, as expected given that a 1∶1 mixture of H/L-labeled parasite lysates was used in our pull down experiments (Table 2). This strongly suggests that more invasion motor is available for immunoprecipitation after heavy labeled parasites underwent calcium pathway stimulation. Two possible interpretations of this result are that the invasion motor complex forms more readily upon calcium signaling or the invasion motor could potentially form higher order structures after calcium signaling. This is the first evidence to suggest that changes in the structure or assembly of the invasion motor occurs upon calcium signaling, and until a more in-depth analysis is performed these results should be treated with caution.

Identification of a Calmodulin-like protein as new component of the Toxoplasma MyoA motor complex

By interrogating our SILAC data we were able to show that there are potentially changes in the assembly of the invasion motor upon calcium stimulation. We therefore reasoned that we might also be able to identify other unknown components of this complex by this virtue. The most highly represented additional protein was a hypothetical Calmodulin-like Toxoplasma protein with a predicted mass of 15 kDa (TGME49_069440) (Tables 2 and Supplementary Tables S6 and S7), that we have termed Essential Light Chain 1 (ELC1)(see below).

To investigate the localization and potential association of this 15 kDa Calmodulin-like protein with the motor complex, we fused this protein with a triple HA epitope tag at the endogenous locus and performed immunofluorescence (IFA) and co-immunoprecipitation/Western blot analyses (CoIP) using anti-HA or anti-GAP45 antibodies (Figure 4). Immunofluorescence microscopy of stably expressing parasites demonstrates this protein resides both at the apical end (likely the conoid) in intracellular parasites and at the parasite periphery as demonstrated by co-localization with GAP45 (Figure 4A). Interestingly, upon egress from host cells extracellular parasites appeared to lose conoid localization and adopt a more typical pattern of an IMC protein (Figure 4A). The significance of this is not yet understood. Clostridium septicum alpha-toxin-induced swelling of extracellular parasites revealed this protein is more firmly anchored to the IMC than the plasma membrane (PM), providing further evidence that this protein could be tightly associated with the motor complex (Figure 4B). CoIP’s performed with anti-HA antibody on tagged and wildtype lines (as a control) followed by Western blot analyses formally established the association of this protein with the invasion motor complex (Figure 4C). In addition, pull-downs using antibodies against MLC1-HA (Figure 4D) or GAP45 (Figure 4E, F) precipitated the endogenous or HA-tagged ELC1. The endogenous ELC1 protein was readily detectable as a 15-kDa Sypro-Ruby stained protein band present in anti-HA IP’s in a MLC1-HA expressing line and proven by LC-MS/MS (Figure 4D and Supplementary Figure S8). The co-purification of ELC1-HA by anti-GAP45 antibodies (Figure 4E and Supplementary Tables S6 and 7) was also shown to be highly specific, as demonstrated by IP/Western blot analysis of immunoprecipitates prepared from ELC1-HA expressing parasite lines using anti-GAP45-specific or non-specific rabbit IgG (Figure 4F). Reproducibly, we noticed that this protein could be seen as two bands (Figure 4C, F) but further work will be needed to determine if this is a degredation product or has physiological significance. Toxoplasma ELC1 displays limited homology with apicomplexan calmodulins, but has significant structural similarity to the essential light chains of more conventional myosins. We found that the sequences of MLC1 and ELC1 to be highly compatible with the structure of the essential light chain of Scallop myosin II complex (PDB 2BL0, chains C and B, respectively) and therefore we constructed a model of MyoA, MLC1 and ELC1 based on this structure (Figure 4G, protein alignments used are presented in Supplementary ). We therefore suggest that TGME49_069440 is the Essential Light Chain of MyoA (here named ELC1) and given its predicted calcium-binding ability this new component suggests that Ca²⁺ binding could also directly regulate invasion motor activity.

Discussion

Invasion of apicomplexan parasites is a complex, multistep process critical for the survival of this group of pathogens. In Toxoplasma calcium signaling is required for the release of apical organelles and the activation of the invasion motor [18], yet the molecular pathways mediating these processes are unclear. Pharmacological evidence suggests that upon host cell recognition by an unknown ligand, PLC is activated to produce soluble Inositol 1,4,5 triphosphate (IP₃), which then promotes the release of Ca²⁺ from ER stores. An increase in cytoplasmic Ca²⁺ concentration activates a range of calmodulin and calcium-dependent protein kinases, allowing for phosphorylation of specific targets. This then is thought to change the cellular activity of substrates, thus activating the invasion process [18].

Recent work has identified several kinases that potentially play a role in invasion [18], [19], [21], yet little is known about the substrates that they target. To address this problem we have used global proteomics approaches to identify these substrates and understand the patterns of calcium-dependent phosphorylation upon motility and invasion. We have identified over 50 potential calcium-dependent phosphorylation substrates and detected phosphorylation sites on at least 10 on these proteins. Further, by analyzing the 2-DE pattern of phosphoproteins after stimulation with either ethanol (thought to activate PLC) or the calcium-mobilizing drug Ionomycin we have shown that both these agonists stimulate the phosphorylation of a largely overlapping set of proteins, strongly suggesting that both these agents act to stimulate the same signal transduction cascade. Our analysis therefore provides further evidence that PLC acts in the same pathway as calcium signaling to activate host cell invasion [16], [18].

Using a combination of molecular and biochemical approaches we identified a wide range of Toxoplasma phosphoproteins with different predicted cellular functions and validated at least 10 calcium-dependent protein kinase substrates that are phosphorylated. Our work therefore presents a snapshot of the complexity associated with calcium-mediated signal transduction in apicomplexan parasites. It will now be interesting to understand the function of these individual molecules in mediating host cell invasion. For example, does the phosphorylation of the vesicle fusion components αSNAP and ARM1 (homologous to Vac8 from yeast) activate apical organelle release? The calcium-dependent phosphorylation of the kinase regulatory subunits CKIIß and PKA-R and the protein phosphatase PP2C are also interesting observations as they suggest modulation of the activity of these molecules by calcium signal transduction pathways. Dissecting the role of individual calcium-dependent phosphoproteins is clearly the next important step in revealing how this group of pathogens regulates host cell invasion.

We have also shown that upon calcium-pathway stimulation the invasion motor is a major in vivo target of calcium-dependent phosphorylation. This strongly suggests that calcium-dependent phosphorylation events stimulate invasion motor activity. GAP45 is the major target of phosphorylation with seven identified sites. GAP45 is a multifunctional protein essential for both maintaining cohesion between the IMC and the PM and also for recruiting MyoA-MLC1 to the parasite periphery for productive movement [7]. GAP45 is tethered to the cytoplasmic face of the plasma membrane via it’s N-terminal lipid anchor and to extend across the ∼300 Å supra-alveolar space to the outer face of the inner membrane complex, where it interacts with invasion motor components via it’s C-terminal globular domain [7]. All seven GAP45 phosphorylation sites that we identified are clustered together in a semi-conserved region (aa 152-192) between the predicted coiled-coil region and the highly conserved globular domain (Figure 5A, B). Although it has been demonstrated that the coiled-coil and globular domain are both important for parasite motility [7] nothing is yet known about the role of this highly phosphorylated region. Bioinformatic analysis of the pairwise energy content of the GAP45 amino acid sequence using the IUPred algorithm [27] predicts this region to be intrinsically unstructured (Figure 5B, blue line). Furthermore, when the S/T/Y phosphorylation sites are replaced by glutamic acid (mimicking phosphorylated residues), IUPred predicts an additional ∼25% increase in the disorder tendency of this domain (Figure 5B, red line). Interestingly, clusters of phosphorylation are often located in regions predicted to be intrinsically unstructured [28] and are commonly protein-protein interaction domains [29]. This suggests that GAP45’s hyperphosphorylated region may represent a protein-protein interaction domain regulated by phosphorylation. Indeed, recent studies have implicated phosphorylation of S163 and S167, two calcium-independent phosphorylation sites within this region, in modulating an interaction of GAP45 with GAP50 [30], further supporting the notion that this region is involved in motor regulation through protein interaction.

It is also apparent that only the calcium-dependent phosphorylation sites on Toxoplasma GAP45 are conserved amongst apicomplexan species (Figure 5C). Of the three Ca²⁺-dependent residues authenticated in the present study, S184 is replaced by a glutamic acid in sequences of other Plasmodium species, Toxoplasma GAP45 residue S185 maps to in vitro phosphorylation site S149 on PfGAP45 [31], and TgGAP45 residue T189 did not map to any identified PfCDPK1 in vitro phosphorylation site (Figure 5C). Homologous Toxoplasma residues corresponding to other PfGAP45 phosphorylation sites identified as in vitro substrates of PfCDPK1 (S156, S173) [31] were not detected in our present study, putting some doubt on the relevance of the sites identified in vitro (Figure 5C). This data suggests that essential aspects of Ca²⁺-dependent phosphorylation of GAP45 might be highly conserved amongst apicomplexan species (Figure 5C).

In Toxoplasma MLC1 anchors MyoA to GAP45 through a long N-terminal extension (Figure 6A)[7], separate from the degenerate EF hands (Figure 6B). Alignments of MLC1/MTIP’s from other Apicomplexa suggest that this domain can be further broken down into a relatively well-conserved N-terminal portion and a C-terminal region, which, like the hyperphosphorylated segment of GAP45, is also predicted to be intrinsically unstructured (Figure 6A, Supplementary Table S8). The calcium-independent phosphorylation site S55 that we have identified is found within this disordered segment of the N-terminal extension (Figure 6A). Within this region two other phosphorylation sites have been identified in vitro in P. falciparum [21], [22], and like S55 also appear to be unconserved across the phylum. Another post-translational modification has been identified neighboring phospho-S55 and although the nature of this modification is currently unknown it has been implicated in regulating the activity of the invasion motor [32]. This region of MLC1 therefore appears to be a hotspot of posttranslational modifications, suggesting that this could also be a regulatory domain controlling the activity and/or assembly of the invasion motor. Given the role of the N-terminal segment of MLC1 in anchoring MLC1-MyoA to GAP45 [7] it is possible that S55 and other posttranslational modifications found within this region could regulate this interaction.

The two calcium-dependent phosphorylation sites that we identified on MLC1 map to the C-terminal region, containing the degenerate EF hands responsible for binding the MyoA tail (Figure 6B). The two Ca²⁺-dependent MLC1 phosphorylation sites that were located in this work (T98, S132) do not correspond to residues critical for MyoA binding [33], nor are they conserved amongst Apicomplexa (Figure 6B). To understand a potential role of Ca²⁺-dependent phosphorylation we mapped these sites onto a modeled structure of MLC1-MyoA (Figure 6C). The Calmodulin-like C-terminal region of MLC1 shows significant structural homology to the regulatory domain of scallop (Physarum polycephalum) myosin, which is made up of two highly ordered degenerate EF-hand domains (N- and C-lobe), connected by a central helical ‘hinge’ region (Figure 6B, C)[34]. Interestingly, the calcium-dependent phosphorylation sites that we identified map to the N-lobe and on the surface of one particular side of MLC1’s predicted three-dimensional structure (Figure 6C). This poses the interesting possibility that this particular region could be an interaction face regulating the attachment of MLC1 with another invasion motor component.

The only phosphorylation site to be identified on MyoA was at S21. We saw that this site changes mildly in abundance upon calcium signaling. This region is found in the head of the MyoA a region that binds to the actin filament. Although this residue appears conserved across Apicomplexan species (Supplementary Figure S10) we at this stage are unable to speculate on its function.

Upon calcium pathway stimulation we also noticed that another protein more strongly associates with the invasion motor and we show that this Calmodulin-like protein is a bona fide component of the Toxoplasma motor complex. This Calmodulin-like Toxoplasma protein is phylogenetically distinct from both canonical calmodulins and myosin-light chains of apicomplexans [35]. Based on the presence of the intact EF hands of this component and the predictions from protein threading we suggest that this new component could represent the Essential Light Chain of MyoA.

Myosin activity in muscle tissue is regulated by one of two mechanisms, either by the calcium-dependent phosphorylation of the regulatory light chain (RLC) or by the direct binding of calcium ions to an EF hand within ELC. Both RLC and ELC bind to the lever arm of myosin and either mode of activation is thought to rigidify the myosin motor and allow it to function more efficiently [36]. We found that the sequences of MLC1 and ELC1 to be highly compatible with the structure of the essential light chain of scallop myosin II complex (PDB 2BL0, chains C and B, respectively) and therefore we constructed a model of MyoA, MLC1 and ELC1 based on this structure (Figure 4G and Supplementary Figure S9). The model of Toxoplasma ELC1 includes a calcium-binding site juxtaposed near the region of closest approach between the ELC1 (residues T16-D17) and MLC1 (residues Y177-G178-E179) – homologues to G178 that is highly conserved amongst RLC’s and other calmodulins. Therefore the model predicts that an interaction between Toxoplasma MLC1 and ELC1 is also possible (Figure 4G). This opens the way to investigating the theory that invasion motor activity may be regulated by direct binding of calcium ions as well as calcium-dependent phosphorylation.

What kinases are responsible for invasion motor phosphorylation? We hypothesized that enzymes engaged in these events may be associated with the invasion motor complex. We therefore interrogated our proteomic data set to identify possible interacting kinases, but unfortunately we were unable to detect any significant hits which we attribute to be likely due to the transitory nature of kinases with their substrates. Appraisal of the literature and the motif signature of calcium-dependent phosphorylation sites suggest that members of the CPDK family are the most likely candidates (data not shown). Given that PfCDPK1 has been implicated in motor phosphorylation in P. falciparum the most likely candidate in Toxoplasma is the phylogenetically related Toxoplasma kinase CDPK3 [18]. TgCDPK3, like PfCDPK1, has a predicted motif for N-terminal acylation, which could potentially anchor it to the PM, putting this enzyme in the vicinity to modulate invasion motor phosphorylation in this parasite.

There are several possible kinases that could modulate the calcium-independent phosphorylation status of the motor; Protein Kinase B (PKB) has been shown to phosphorylate PfGAP45 in vitro, whereas Toxoplasma Protein Kinase G (PKG) is found at the periphery and a specific inhibitor to this enzyme prevents parasite motility and invasion [37]. Another candidate is Casein Kinase II (CKII). The beta regulatory subunit of CKII in Toxoplasma also localizes to the periphery of parasite [38] and two calcium-independent phosphorylation sites on GAP45 conform to a CKII-like substrate motif (data not shown).

A tyrosine-based kinase also likely modulates the phosphorylation status of the invasion motor. Tyrosine 158 on GAP45 is phosphorylated in a calcium-independent manner. Given that tyrosine kinases form a unique class of enzyme this suggests a third pathway regulates that invasion motor phosphorylation. Tyrosine kinase and tyrosine phosphorylation have not yet been reported for any apicomplexan parasite, and therefore, it is not possible to speculate on what the identity of this kinase might be. Overall, our data suggests that there are at least three different kinases, and therefore, at least three potential pathways that modulate the phosphorylation status of the invasion motor. Applying specific conditional mutants to our quantitative proteomics approach has the potential to reveal the in vivo kinase for these important phosphorylation events.

We have also shown that upon induction of calcium signaling more invasion motor components are found associated with GAP45. This suggests that although the invasion motor can be found pre-formed in resting parasites, calcium signaling induces changes in the association between GAP45 and the other components. Mechanistically such a change could be one of the following: 1. Inducing more invasion motor complexes to form from individual components, 2. The invasion motor complex is rendered more stable after calcium signaling or 3. A higher order structure, such as multimerization is induced. To distinguish from these possibilities and to confirm our finding is physiologically relevant biochemical analyses of isolated invasion motor complexes will need to be performed.

Materials and Methods

Cloning of DNA constructs

Predicted open reading frames of phosphoproteins (as described at www.toxodb.org) were amplified and inserted into Toxoplasma transfection vectors. Toxoplasma MLC1, ARM1 and PP2C were amplified with primers 1&2, 3&4, 5&6 (Table 3) respectively and inserted into BglII/AflII sites of pCT3H (Tonkin, unpublished), which enabled the tagging of these genes with at the C-terminus with a 3xHA tag. PKA-R was amplified with primers 7&8 and a HA tag introduced by adding this sequence ending this tag onto the 5’ primer (Table 3). PKA-R PCR product was then inserted behind the fkbp-derived destabilization domain (DD) at AvrII/PstI of pCTDD (Tonkin, unpublished).

The invasion motor associated Calmodulin-like protein (TGME49_069440) was tagged at the 3’ end of the endogenous gene using ΔKu80 parasite line [39]. A 3’ flank of this gene was amplified with primers 9 and 10 and inserted into pLIC-HA3/HX [40] using the ligation independent cloning strategy. Plasmids were linearized within the gene flank for efficient homologous integration after transfection [39].

Toxoplasma culture and transfection

Toxoplasma parasites were grown using standard procedures. Briefly, tachyzoites were grown in confluent human foreskin fibroblasts (HFF) or Vero cells maintained in Dulbecco’s Modified Eagle’s Medium (DMEM; GIBCO, Invitrogen) supplemented with 10% cosmic calf serum (CCS) and an additional 2 mM glutamine (Gibco). If HFF were used DMEM/10%CCS was replaced with DMEM supplemented with 1% fetal calf serum (DMEM/1%FCS) with additional 2 mM glutamine (Gibco) upon infection with parasites.

Toxoplasma was transfected using standard procedures [41], [42]. Electroporated tachyzoites were inoculated onto confluent HFF cells and selected on 6 ug/ml of chloramphenicol [41] or Mycophenolic Acid (25 µg/ml)/Xanthine (50 µg/ml) for stable transfectants [43].

Immunofluorescence analyses and microscopy

IFA was performed using standard procedures. Briefly, Parasites, either intracellular or extracellular, were fixed in 4% paraformaldehyde, permeabilized in 0.1% Triton-X/PBS and blocked in 3% BSA/PBS. Anti-HA mAb (clone 3F10;Roche), anti-GAP45 [5] and anti-SAG1 (a kind gift from L.D Sibley) were decorated onto blocked slides and primary antibodies were detected with antibodies conjugated to AlexaFluor-594 and 488 (Molecular Probes). Extracellular parasites were stuck to slides using 0.1% polyethylenimine. Parasites were treated with dialyzed C. septicum culture supernatant containing alpha-Toxin for ∼3 hours. Parasites were then imaged on a Zeiss Inverted Axioscope equipped with AxioCam MRM and Axiovision with Deconvolution software.

Radioactive or stable-isotope labeling and Ca²⁺ pathway stimulation of Toxoplasma

For radioactive or heavy isotope labeling and Ca²⁺ pathway stimulation studies, specialized cell culture conditions were applied. For radioactive labeling, Vero cells were infected with T. gondii tachyzoites and metabolically labeled in 25 ml of custom sodium phosphate-free DMEM (GIBCO, Invitrogen) supplemented with 10% fetal calf serum, 2 mM glutamine and 20 µCi/ml of ³²[P]-labeled monosodium phosphate, overnight.

For stable isotope labeling studies, HFF were infected with Toxoplasma and grown in custom DMEM medium devoid of lysine, arginine and leucine. This media was then supplemented with 1% fetal calf serum (dialyzed against sterile PBS), 2 mM glutamine, 0.802 mM L-leucine, and either 0.398 mM L-Arginine + 0.798 mM L-Lysine for the „light“ condition (R0/K0 media), or 0.398 mM (¹⁵N₄) L-Arginine + 0.798 mM (¹³C₆, ¹⁵N₂) L-Lysine (98% isotopic purity; Sigma-Isotec) for the “heavy” condition (R4/K8 media). Heavy label incorporation was applied to parasites for ∼6 days to get maximal stable isotope incorporation before bulking up and subsequent parasite harvest. Analysis demonstrated >95% incorporation of R4 or K8 labels after this time (Supplementary Figure S3).

The preparation of free, viable tachyzoites was done as follows; Upon the >80% disruption of HFF cells, parasites were scraped, needle passed and pushed through a 5 µM pore-size membrane filters (Millipore) and then repeatedly washed in ice cold Invasion Medium (DMEM/HEPES/1%FCS) and counted using a hemocytometer. Aliquots of radioactive or heavy-isotope-labeled tachyzoites were resuspended in invasion medium and incubated at 37°C for 10 min prior to Ca²⁺ pathway stimulation. In vivo Ca²⁺ pathway stimulation was achieved by adding 1% (171 mM) ethanol or 1 µM ionomycin. As a control, samples of parasites were either mock treated (equal volume DMSO) or were pretreated for 10 minutes with the membrane-permeable calcium chelator BAPTA-AM, followed by stimulation with 1% ethanol. Samples were stimulated for 60 seconds, immediately mixed with 0.5 ml ice cold 2 × lysis buffer (40 mM HEPES pH 7.4, 300 mM NaCl, 2% Triton X-100, 1% NP-40) containing protease and phosphatase inhibitors (HALT protease and phosphatase inhibitor solution, PIERCE) and snap frozen in liquid nitrogen.

Preparation of parasite protein extracts

Parasites were disrupted by freeze thawing and ultrasonication in ice cold lysis buffer (20 mM HEPES pH 7.4, 150 mM NaCl, 1% Triton X-100, 0.5% NP-40) including protease and phosphatase inhibitors (Pierce), and proteins extracted for 30 min on ice with vortexing. Parasite extracts were clarified by ultracentrifugation at 75,000 g for 30 min and the detergent-soluble protein fraction was precipitated using 2-D Clean-Up kit (GE-healthcare).

2-D gel electrophoresis

For 2-DE gel electrophoresis, precipitated Toxoplasma protein preparations (∼50 µg protein for analytical gels/∼500 µg protein for preparative gels) were redissolved in 300 µl rehydration/sample buffer (7 M Urea, 2 M Thiourea, 2% ASB-14, 1% DTT, 1% ampholytes), loaded onto 13 cm pI 4–7 IPG strips by passive rehydration and focused at a current limit of 50 µA/IPG strip using a fast voltage gradient (8000 V max, 24,000 Vh) at 15°C. The second dimension was carried out on 10% polyacrylamide gels (18 cm×16 cm×1.5 mm) using a Hoefer SE 600 system (GE Healthcare) at 75 V constant voltage and 10°C overnight. Analytical 2-DE gels were electrophoretically transferred to Immobilon-PSQ PVDF membranes (Millipore). Protein spots on PVDF membranes were visualized using Deep Purple protein stain (GE Healthcare), and protein spots in preparative 2-DE gels were stained with Sypro Ruby Protein stain (Molecular Probes), according to manufacturer’s protocols. Imaging of ³²[P]-labeled protein spots was achieved by direct autoradiography (7 day exposure) of dry PVDF membrane blots using FUJIFILM BAS-TR2040 tritium imaging plates. Fluorescent or autoradiographic 2-DE images were digitized on a FLA-3000 laser-scanning detection system (Fuji) and manually matched and annotated using Image Master Platinum v7.0 2D image analysis software (GE Healthcare). Preparative 2-DE gels were counter-stained using Colloidal Coomassie Brilliant Blue (Sigma) and regions matching ³²[P]-labeled protein spots manually excised and subjected to LC-MS/MS analysis.

Gel excision, in-gel digestion and Nano-LC-MS/MS

For 1-DE or 2-DE samples, protein spots or bands were manually excised from preparative SDS-PAGE gels and subjected to automated in-gel reduction, alkylation, and tryptic digestion using a MassPREP Station (Micromass, UK). All gel samples were reduced with 10 mM DTT (SIGMA) for 30 min, alkylated for 30 min with 50 mM iodoacetic acid (SIGMA) and digested with 375 ng trypsin (Promega) for 16 hrs at 37°C. The extracted peptide solutions were then acidified (0.1% formic acid) and concentrated to approximately 10 ul by centrifugal lyophilisation using a SpeedVac AES 1010 (Savant). Briefly, extracted peptides were injected and fractionated by nanoflow reversed-phase liquid chromatography on a nano LC system (1200 series, Agilent) using a nanoAcquity C18 150 mm×0.15 mm I.D. column (Waters) developed with a linear 60-min gradient with a flow rate of 0.5 µl/min at 45°C from 100% solvent A (0.1% Formic acid in Milli-Q water) to 100% solvent B (0.1% Formic acid, 60% acetonitrile, 40% Milli-Q water). The nano HPLC was coupled on-line to an LTQ-Orbitrap mass spectrometer equipped with a nanoelectrospray ion source (Thermo Fisher Scientific) for automated MS/MS. Up to five most intense ions per cycle were fragmented and analysed in the linear trap, with target ions already selected for MS/MS being dynamically excluded for 3 min.

Large-scale phosphopeptide purification and MudPIT analysis

For proteome-wide phosphopeptide analyses using MudPIT, precipitated Toxoplasma protein preparations (2.5 mg of protein of ionomycin-stimulated tachyzoites) were dissolved in a total volume of 1.25 ml of 6 M Urea, 100 mM Tris-HCl (pH 8.0) and sequentially reduced with 5 mM DTT and alkylated with 10 mM iodoacetamide. Samples were diluted 1∶5 (v/v) in digestion buffer (40 mM ammonium bicarbonate, pH 8, 10% acetonitrile, 1 mM Ca₂Cl) and digested with LysC (Roche) and trypsin (Sigma). Approximately 1 mg tryptic digest was desalted using Sep-Pak C18 6 cc/1 g cartridges (Waters) and fractionated by hydrophilic interaction chromatography on a TSKgel Amide 80 column (TOSOH Biosciences) using an optimized phosphopeptide gradient, as per published protocol [44]. Ten phosphopeptide fractions were collected and lyophilized. Fractions 1-4 were pooled and the resulting 7 samples were resuspended in 100 µl TiO₂ loading/wash buffer (1 M glycolic acid, 80% ACN, 5% TFA) and purified on self-packed TiO₂ micro-columns using a highly efficient method for the selective enrichment of phosphorylated peptides, as detailed elsewhere [45]. The flow-through and all washings were combined and constituted the phosphopeptide-depleted (TiO₂-unbound) fraction. Individual TiO₂-bound phosphopeptide fractions were analyzed by multi-dimensional LC-MS/MS on an LTQ linear ion trap mass spectrometer (Thermo Scientific), according to published protocols [46].

Mass spectra database searching

MudPIT LC-MS/MS spectra were analyzed with SEQUEST 2.7 [47] using a non-redundant protein decoy database (Ludwig NR_Q309_con reverse; 9870917 entries for all species). The SEQUEST outputs were analyzed by DTASelect 2.0.37 [48]. DTASelect 2.0 uses a quadratic discriminant analysis to dynamically set XCorr and DeltaCN thresholds for the entire data set to achieve a user-specified false positive rate (1% in this analysis). The false positive rates were estimated by the program from the number and quality of spectral matches to the decoy database [49]. After filtering the results from SEQUEST using DTASelect, MS/MS spectra were analysed using the DeBunker algorithm for automatic validation of phosphopeptide identifications from tandem mass spectra [50]. This software package uses a support vector machine binary classifier to assess the correctness of phosphopeptide/spectrum matches.

For protein identification of gel samples and the identification of additional MudPIT phosphopeptide sequences (ie. not included in the LudwigNR database) LC-MS/MS data were searched against a redundant protein decoy database comprising sequences from the latest version of Swiss-Prot (Human, Bovine, Plasmodium, Toxoplasma species), Trembl (Toxoplasma entries), PlasmoDB/ToxoDB (Toxoplasma entries), as well as their reverse sequences (Toxoplasma_decoy; 117496 entries). Mass spectra peak lists were extracted using extract-msn as part of Bioworks 3.3.1 (Thermo Fisher Scientific) linked into Mascot Daemon (Matrix Science, UK). The parameters used to generate the peak lists for the LTQ Orbitrap were as follows: minimum mass 400; maximum mass 5000; grouping tolerance 0.01 Da; intermediate scans 1; minimum group count 1; 10 peaks minimum and total ion current of 100. Peak lists for each nano-LC-MS/MS run were used to search MASCOT v2.2.04 search algorithm (Matrix Science, UK) provided by the Australian Proteomics Computational Facility (www.apcf.edu.au). The search parameters consisted of carboxymethylation of cysteine as a fixed modification (+58 Da, for gel samples only), with variable modifications set for NH₂-terminal acetylation (+42 Da) and oxidation of methionine (+16 Da), phosphorylation of serine, threonine or tyrosine (+80 Da). A precursor mass tolerance of ±3 Da (LTQ spectra) or 20 ppm (Orbitrap spectra), #13C defined as 1, fragment ion mass tolerance of ±0.8 Da, and an allowance for up to three missed cleavages for tryptic searches was used.

Generation of anti-GAP45 immunoaffinity resin

To generate immunoaffinity resin anti-GAP45 IgG was purified using a 1 ml HiTrap Protein-A HP column (GE Healthcare) using an ÄKTA prime FPLC chromatography system (Pharmacia). Purified IgG was desalted using a PD-10 buffer exchange column and chemically cross-linked to CN-Br-activated Sepharose 4B resin (GE Healthcare) in 0.1 M NaHCO₃, 0.5 M NaCl pH 8.5 for 2 h at room temperature. Un-reacted sites were blocked with 1 M triethanolamine pH 8.0 for 2 h at room temperature and the resin stringently washed with 50 mM glycine-HCl, 0.5 M NaCl pH 3.5 and 50 mM Tris-HCl, 0.5 M NaCl pH 8.0 and resuspended in PBS.

Immunoaffinity purification of Toxoplasma invasion motor complexes

Large-scale immunoaffinity purification of intact MyoA motor complexes from Toxoplasma parasites was carried out using anti-GAP45 rabbit serum, as previously described [5]. For our SILAC-based quantitative LC-MS/MS analyses, Toxoplasma tachyzoites were grown in confluent HFF’s in “light” R0/K0 media or “heavy” R4/K8 media for 48 h. We found that 4×T150 cm² flasks each of light- or heavy-labeled Toxoplasma cultures yielded sufficient material to obtain >60% sequence coverage for most of the complex components. Free, viable tachyzoites were prepared and the heavy isotope-labeled parasites were stimulated with 1% ethanol to measure the effect of Ca²⁺-pathways stimulation, as detailed above. Freshly prepared parasite lysates of un-stimulated (light) or 1% ethanol-stimulated (heavy) tachyzoites were prepared as detailed above, mixed at a protein ratio of 1∶1 in 2 ml lysis buffer with protease and phosphatase inhibitors and incubated with 0.25 ml of the anti-GAP45 resin for 1 h at 4°C. Unbound protein was removed and the resin washed 5 times with 1 ml lysis buffer and 2 times with 0.5 ml water by centrifugation at 2000 g for 5 min in 1 ml microcentrifuge spin columns (Pierce). Invasion motor complexes were eluted with 200 µl 0.1% TFA in water and the eluted material dried in a SpeedVac.

In-solution digestion of Toxoplasma invasion motor complexes and phosphopeptide enrichment

Dried samples of anti-GAP45 column eluates were dissolved in a total volume of 250 µl of 6 M Urea, 50 mM ammonium bicarbonate (pH 8.0) and sequentially reduced with 10 mM DTT and alkylated with 55 mM iodoacetamide. Samples were diluted 1∶5 (v/v) in digestion buffer (40 mM ammonium bicarbonate, pH 8, 10% acetonitrile, 1 mM Ca₂Cl) and digested with proteomics-grade trypsin (Sigma) at 37°C overnight. The resulting peptide solutions were vacuum-dried, peptides resolved in 0.2%TFA in water, desalted on MacroSpin C18 columns (The Nest Group Inc.), eluted with 0.2%TFA/60% ACN in water and again dried in a SpeedVac. Desalted peptide samples were resuspended in 100 µl TiO₂ loading buffer (1 M glycolic acid in 80% ACN, 5% TFA) and affinity purified using self-packed TiO₂ microcolumns containing ∼10 µg of 5 µm Titansphere TiO₂ beads (GL Sciences), using a published method [24]. The flow-through and all washings were combined and constituted the phosphopeptide-depleted (TiO₂-unbound) fraction. The TiO₂ eluate comprising the phosphopeptide enriched fraction and the phosphopeptide-depleted fraction were vacuum-dried. Both fractions were redissolved in 0.2%TFA in water, desalted on MacroSpin C18 columns (The Nest Group Inc.) and again dried in a SpeedVac. For LC-MS/MS analyses, dried (phospho)peptide fractions were redissolved in 1 µl 100% formic acid and diluted to 20 µl in MilliQ water. Digested (phospho)peptides were then subjected to nano-LC-MS/MS on an LTQ-Orbitrap instrument (Thermo Scientific), as described above.

SILAC-based phosphopeptide identification and quantification

LTQ-Orbitrap LC-MS/MS data were searched against the Toxoplasma_decoy database using the MASCOT search engine, as detailed above. Search parameters were identical, except that additional variable modifications were set for heavy-isotope labeling of arginine (+4 Da) or lysine (+8 Da) and the precursor mass tolerance was 20 ppm for the SILAC/Orbitrap datasets. Mascot search results were loaded into MaxQuant (v1.1.08) [51] for peptide and protein quantification as well as Scaffold (v3.0) (www.proteomesoftware.com) for manual validation of peptide spectral matches. Manual validation of all peptide spectral matches was done irrespective of peptide scores or expectation values (E-values) to ensure that all major fragment ions were annotated in accordance with known rules of peptide fragmentation. The Xcalibur program (Thermo Scientific) was used for generating XIC's of selected peptides in order to validate and verify the MaxQuant results.

Co-immunoprecipitation of invasion motor complexes and Western blot analyses

For anti-HA CoIP’s, freshly released tachyzoites expressing wild-type (negative control) or HA-tagged parasite proteins were harvested as above, washed in PBS, and lysed in ice cold lysis buffer (20 mM HEPES pH 7.4, 150 mM NaCl, 1% Triton X-100, 0.5% NP-40) containing protease and phosphatase inhibitors (Pierce). Protein was extracted for 30 min on ice with vortexing and centrifuged at 50,000 g for 30 min at 4°C. Supernatants were mixed with monoclonal anti-HA antibody coupled to paramagnetic MicroBeads (Miltenyi Biotec) and incubated for 1 h at 4°C. Supernatants were then subjected to CoIP, according to a manufacturer’s protocol for the specific isolation of HA-tagged proteins utilizing a MultiMACS separator (Miltenyi Biotech). For anti-GAP45 CoIP’s, 1 ml protein extracted from parasites expressing ELC1-HA were mixed with 100 µl of ProteinG-coupled paramagnetic microbeads (Miltenyi Biotech), mixed with 5 µl of anti-GAP45-specific rabbit IgG or non-specific rabbit IgG (negative control), and processed as above. For Western blot analyses eluted protein samples were separated by SDS-PAGE on precast 10% or 4-12% NuPAGE Bis-Tris gels, transferred onto Nylon membranes using an iBlot blotting system (Invitrogen), and then probed with specific primary or HRP-conjugated secondary antibody as indicated.

Protein fold recognition and comparative modeling

Protein fold recognition was conducted using the WURST protein-threading web server [52]. Homology models were constructed using the sequence alignments predicted from WURST with the MODELLER (9v7) comparative modeling software [53].

List of accession numbers of genes mentioned in the text

PKA-R, TGME49_042070; CKIIß, TGME49_072400; SNAP, TGME49_018760; ARM1, TGME49_061440; PP2C, TGME49_054770; Rab5B, TGME49_007460; ANP32, TGME49_071810; SET, TGME49_044110; Toxofilin, TGME49_014080; Hypothetical protein, TGME49_048700; Conserved hypothetical protein, TGME49_032440; CaM-like (ELC1), TGME49_069440; MLC1, TGME49_057680; GAP45, TGME49_023940; MyoA, TGME49_035470; GAP50, TGME49_019320; GAP40, TGME49_049850.

Supporting Information

Attachment 1

Attachment 2

Attachment 3

Attachment 4

Attachment 5

Attachment 6

Attachment 7

Attachment 8

Attachment 9

Attachment 10

Attachment 11

Attachment 12

Attachment 13

Attachment 14

Attachment 15

Attachment 16

Attachment 17

Attachment 18

Attachment 19

Attachment 20

Attachment 21