Structural and evolutionary adaptation of rhoptry kinases and pseudokinases, a family of coccidian virulence factors

An atlas of bacterial serine-threonine kinases reveals functional diversity and key distinctions from eukaryotic kinases

Bacterial serine-threonine kinases (STKs) regulate diverse cellular processes associated with cell growth, virulence, and pathogenicity and are evolutionarily related to the druggable eukaryotic STKs. A deeper understanding of how bacterial STKs differ from their eukaryotic counterparts and how they have evolved to regulate diverse bacterial signaling functions is crucial for advancing the discovery and development of new antibiotic therapies. Here, we classified more than 300,000 bacterial STK sequences from the NCBI RefSeq nonredundant and UniProt protein databases into 35 canonical and seven pseudokinase families on the basis of the patterns of evolutionary constraints in the conserved catalytic domain and flanking regulatory domains. Through statistical comparisons, we identified features distinguishing bacterial STKs from eukaryotic STKs, including an arginine residue in a regulatory helix (C helix) that dynamically couples the ATP- and substrate-binding lobes of the kinase domain. Biochemical and peptide library screens demonstrated that evolutionarily constrained residues contributed to substrate specificity and kinase activation in the Mycobacterium tuberculosis kinase PknB. Together, these findings open previously unidentified avenues for investigating bacterial STK functions in cellular signaling and for developing selective bacterial STK inhibitors.

Atlas of the Bacterial Serine-Threonine Kinases expands the functional diversity of the kinome

Bacterial serine-threonine protein kinases (STKs) regulate diverse cellular processes associated with cell growth, virulence, and pathogenicity. They are evolutionarily related to the druggable eukaryotic STKs. However, an incomplete knowledge of how bacterial STKs differ from their eukaryotic counterparts and how they have diverged to regulate diverse bacterial signaling functions presents a bottleneck in targeting them for drug discovery efforts. Here, we classified over 300,000 bacterial STK sequences from the NCBI RefSeq non-redundant and UniProt protein databases into 35 canonical and seven non-canonical (pseudokinase) families based on the patterns of evolutionary constraints in the conserved catalytic domain and flanking regulatory domains. Through statistical comparisons, we identified distinguishing features of bacterial STKs, including a distinctive arginine residue in a regulatory helix (C-Helix) that dynamically couples ATP and substrate binding lobes of the kinase domain. Biochemical and peptide-library screens demonstrated that constrained residues contribute to substrate specificity and kinase activation in the Mycobacterium tuberculosis kinase PknB. Collectively, these findings open new avenues for investigating bacterial STK functions in cellular signaling and for the development of selective bacterial STK inhibitors.

In silico analysis and structural vaccinology prediction of Toxoplasma gondii ROP41 gene via immunoinformatics methods as a vaccine candidate

INTRODUCTION

Toxoplasma gondii (T. gondii) infects all warm-blooded animals, including humans. Currently, no effective treatments exist to prevent the generation of chronic tissue cysts in infected hosts. Therefore, developing a vaccine to protect to deal with toxoplasmosis is a promising strategy, as a single immunization could provide lifelong protective immunity. Rhoptry proteins (ROPs) play a vital role for the parasite's survival within host cells and perform critical functions during different phases of parasite invasion. Little is known about ROP41 gene. Nevertheless, Understanding the characteristics of ROP41 will enhance diagnostic and vaccine research.

MATERIALS AND METHODS

The current article provides a comprehensive analysis of the essential components of the ROP41 protein, including its transmembrane domain, physico-chemical properties, subcellular location, tertiary and secondary structures, and potential T- and B-cell epitopes. These features were determined by many bioinformatics approaches to identify possible epitopes for developing a highly effective vaccine.

RESULTS

ROP41 protein showed 36 possible post-translational modification regions. The ROP41 protein secondary structure contains 17.35 % extended strand, 33.47 % alpha-helix, and 49.18 % random coil. Also, ROP41 showed many possible B- and T-cell epitopes. According to the Ramachandran plot, 90.78 % of amino acid residues had been placed in favored, 3.28 % in outlier, and 5.94 % in allowed areas. Also, the allergenicity and antigenicity evaluation indicated that ROP41 is non-allergenic and immunogenic.

CONCLUSION

The current study offered critical basic and conceptual information on ROP41 to increase a successful vaccine in opposition to continual and acute toxoplasmosis for in addition in vivo assessments. Further research is necessary for the development of vaccines utilizing ROP41 alone or combined with various antigens.

A bioinformatic approach for the prediction and functional classification of Toxoplasma gondii long non-coding RNAs

Long non-coding RNAs (lncRNAs) have emerged as significant players in diverse cellular processes, including cell differentiation. Advancements in computational methodologies have facilitated the prediction of lncRNA functions, enabling insights even in non-model organisms like pathogenic parasites, in roles such as parasite development, antigenic variation, and epigenetics. In this work, we focus on the apicomplexan Toxoplasma gondii differentiation process, where the infective stage, tachyzoite, can develop into the cysted stage, bradyzoite, under stress conditions. Using a publicly available transcriptome dataset, we predicted putative lncRNA sequences associated with this differentiation process. Notably, a substantial proportion of these putative lncRNAs exhibited stage-specific expression, particularly at the bradyzoite stage. Furthermore, co-expression patterns between coding transcripts and putative TglncRNAs suggest their involvement in shared processes, such as bradyzoite development. Putative TglncRNA loci analysis revealed their potential influence on the expression of nearby coding genes, including subtelomeric genes unique to the T. gondii genome. Finally we propose a k-mer analysis approach to predict putative functional relationships between characterized lncRNAs from model organisms like Homo sapiens and the putative T. gondii lncRNAs. Our perspective led to predict putative T. gondii lncRNA that potentially could act mediating DNA damage repair pathways, opening a new study field to validate this kind of adaptive mechanisms of T. gondii in response to stress conditions.

A candidate virulence factor of Eimeria tenella (EtROP30) predicted by virulence enhancement of transgenic Toxoplasma gondii

Difficulties of in vitro culture and genetic manipulation of Eimeria tenella have hindered the screening of virulence factors in this parasite. In this study, the E. tenella rhoptry protein 30 (EtROP30) was expressed in Toxoplasma gondii (RH∆Ku80-EtROP30), and its effect on the proliferation and virulence of parasites was investigated. The results revealed that the expression of EtROP30 had no impact on the invasion and egress processes. However, the RH∆Ku80-EtROP30 strain formed larger plaques compared to the RH∆Ku80, indicating that the EtROP30 expression promotes T. gondii proliferation. Furthermore, the RH∆Ku80-EtROP30 strain exhibited greater pathogenicity, resulting in earlier mortality and shorter overall survival time compared to RH∆Ku80. These results imply that EtROP30 expression facilitates parasite intracellular proliferation and virulence in mice, suggesting that EtROP30 might be a candidate virulence factor of E. tenella.

The Hypervariable Tpr Multigene Family of Theileria Parasites, Defined by a Conserved, Membrane-Associated, C-Terminal Domain, Includes Several Copies with Defined Orthology Between Species

Multigene families often play an important role in host-parasite interactions. One of the largest multigene families in Theileria parva, the causative agent of East Coast fever, is the T. parva repeat (Tpr) gene family. The function of the putative Tpr proteins remains unknown. The initial publication of the T. parva reference genome identified 39 Tpr family open reading frames (ORFs) sharing a conserved C-terminal domain. Twenty-eight of these are clustered in a central region of chromosome 3, termed the "Tpr locus", while others are dispersed throughout all four nuclear chromosomes. The Tpr locus contains three of the four assembly gaps remaining in the genome, suggesting the presence of additional, as yet uncharacterized, Tpr gene copies. Here, we describe the use of long-read sequencing to attempt to close the gaps in the reference assembly of T. parva (located among multigene families clusters), characterize the full complement of Tpr family ORFs in the T. parva reference genome, and evaluate their evolutionary relationship with Tpr homologs in other Theileria species. We identify three new Tpr family genes in the T. parva reference genome and show that sequence similarity among paralogs in the Tpr locus is significantly higher than between genes outside the Tpr locus. We also identify sequences homologous to the conserved C-terminal domain in five additional Theileria species. Using these sequences, we show that the evolution of this gene family involves conservation of a few orthologs across species, combined with gene gains/losses, and species-specific expansions.

Overexpression of Eimeria tenella Rhoptry Kinase 2 Induces Early Production of Schizonts

Eimeria tenella is an obligate intracellular parasite responsible for avian coccidiosis. Like other apicomplexan parasites, such as Toxoplasma gondii, cell invasion and intracellular development rely on apical organelle content discharge, named micronemes and rhoptries. Some rhoptry (ROP) kinases (ROPK) are key virulence factors in T. gondii. To date, among the 28 ropk genes carried by E. tenella, only two to four were confirmed by proteomic analysis or immunostaining to be expressed at the sporozoite stage. We have previously shown that EtROP1 is implicated in the inhibition of host cell apoptosis by interacting with the cellular p53. This work functionally described the second ROP kinase expressed at the sporozoite stage in E. tenella. EtROP2 is an active kinase that phosphorylates cell substrates of approximately 50 kDa. Its overexpression leads to the shortening of the prepatent period and to the early development of first-generation schizonts. Conduction of RNA sequencing analysis and reverse transcriptase quantitative PCR (RT-qPCR) on the host cell allowed us to identify the mitogen-activated protein kinase (MAPK) pathway and the transcription factor cFos to be upregulated by EtROP2. We also showed by immunofluorescence assay that the active kinase EtROP2 is implicated in the p38 MAPK pathway activation. We established here that EtROP2 activates the p38 MAPK pathway through a direct or indirect phosphorylation, leading to the overexpression of the master transcription factor cFos known to be implicated in E. tenella development. IMPORTANCE Rhoptries are specialized secretory organelles found in zoite stages of apicomplexan parasites. In addition to well-conserved rhoptry neck proteins, their protein consists mostly of kinase proteins, highly divergent from eukaryotic kinases. Some of those kinases are described as major virulence factors in Toxoplasma gondii, secreted into the host cell to hijack signaling pathways. Most of those kinases remain to be characterized in Eimeria tenella. Deciphering their cellular function is a prerequisite to supporting their relevance as a druggable target in development of new means of Eimeria tenella control. Secreted divergent kinases that interact with host cell partners to modulate pathways are good candidates, as they coevolve with their host targets to ensure their function within the host and are less prone to mutations that would lead to drug resistance. The absence of any orthologous kinase in host cells makes these parasite kinases a promising drug target candidate.

Anti-Toxoplasma gondii Effects of Lipopeptide Derivatives of Lycosin-I

Toxoplasmosis, caused by Toxoplasma gondii (T. gondii), is a serious zoonotic parasitic disease. We previously found that Lycosin-I exhibited anti-T. gondii activity, but its serum stability was not good enough. In this study, we aimed to improve the stability and activity of Lycosin-I through fatty acid chain modification, so as to find a better anti-T. gondii drug candidate. The α/ε-amino residues of different lysine residues of Lycosin-I were covalently coupled with lauric acid to obtain eight lipopeptides, namely L-C12, L-C12-1, L-C12-2, L-C12-3, L-C12-4, L-C12-5, L-C12-6, and L-C12-7. Among these eight lipopeptides, L-C12 showed the best activity against T. gondii in vitro in a trypan blue assay. We then conjugated a shorter length fatty chain, aminocaproic acid, at the same modification site of L-C12, namely L-an. The anti-T. gondii effects of Lycosin-I, L-C12 and L-an were evaluated via an invasion assay, proliferation assay and plaque assay in vitro. A mouse model acutely infected with T. gondii tachyzoites was established to evaluate their efficacy in vivo. The serum stability of L-C12 and L-an was improved, and they showed comparable or even better activity than Lycosin-I did in inhibiting the invasion and proliferation of tachyzoites. L-an effectively prolonged the survival time of mice acutely infected with T. gondii. These results suggest that appropriate fatty acid chain modification can improve serum stability and enhance anti-T. gondii effect of Lycosin-I. The lipopeptide derivatives of Lycosin-I have potential as a novel anti-T. gondii drug candidate.

Bioinformatics-based prediction and screening of immunogenic epitopes of Toxoplasma gondii rhoptry proteins 7, 21 and 22 as candidate vaccine target

INTRODUCTION

Toxoplasmosis is a well-known zoonotic disease caused by Toxoplasma gondii. The main causes of the disease range from eating undercooked or contaminated meat and shellfish to cleaning litter trays into which cats that excreted toxoplasma via faeces. This pathogen can live for a very long time, possibly a lifetime, within the bodies of humans and other animals.

AIMS AND OBJECTIVES

This study aimed to predict and analyse candidate immunogenic epitopes for vaccine development by evaluating the physio-chemical properties, multiple sequence alignment, secondary and tertiary structures, phosphorylation sites, transmembrane domains, and signal peptides, of T. gondii rhoptry proteins ROP7, ROP21, and ROP22 using bioinformatics tools.

METHODS

To find immunogenic epitopes of rhoptry proteins, numerous bioinformatics web servers were used containing multiple sequence alignment, physiochemical properties, antigenicity and allergenicity, post-translational modification sites (PTMs), signal peptides, transmembrane domains, secondary and tertiary structures, and screening of predicted epitopes. We evaluated immunogenic linear B-cell epitopes as candidate proteins for vaccine development.

RESULTS

Nine epitopes were identified for each protein, and analysis of immunogenicity, revealed three candidate epitopes for ROP7, one for ROP21, and four for ROP22. Among all candidate epitopes, ROP22 contained the most immunogenic epitopes with immunogenicity score of 0.50575.

CONCLUSION

We acquired detailed information on predicted immunogenic epitopes using in-silico methods. The results provide a foundation for further experimental analysis of toxoplasmosis, and potential vaccine development.

ROP39 is an Irgb10-specific parasite effector that modulates acute Toxoplasma gondii virulence

Toxoplasma gondii (T. gondii) is a zoonotic apicomplexan parasite that is an important cause of clinical disability in humans. On a global scale, one third of the human population is infected with T. gondii. Mice and other small rodents are believed to be responsible for transmission of T. gondii to the domestic cat, its definitive host. Interferon-inducible Immunity-Related GTPases (IRG proteins) are important for control of murine T. gondii infections. Virulence differences between T. gondii strains are linked to polymorphic rhoptry proteins (ROPs) that cooperate to inactivate individual IRG family members. In particular, the pseudokinase ROP5 isoform B is critically important in laboratory strains of mice. We identified T. gondii ROP39 in complex with ROP5B and demonstrate its contribution to acute T. gondii virulence. ROP39 directly targets Irgb10 and inhibits homodimer formation of the GTPase leading to an overall reduction of IRG protein loading onto the parasitophorous vacuolar membrane (PVM). Maintenance of PVM integrity rescues the parasite from IRG protein-mediated clearance in vitro and in vivo. This study identifies a novel T. gondii effector that is important for specific inactivation of the IRG resistance system. Our data reveal that yet unknown T. gondii effectors can emerge from identification of direct interaction partners of ROP5B.

The protein phosphatase 2A holoenzyme is a key regulator of starch metabolism and bradyzoite differentiation in Toxoplasma gondii

Phenotypic switching between tachyzoite and bradyzoite is the fundamental mechanism underpinning the pathogenicity and adaptability of the protozoan parasite Toxoplasma gondii. Although accumulation of cytoplasmic starch granules is a hallmark of the quiescent bradyzoite stage, the regulatory factors and mechanisms contributing to amylopectin storage in bradyzoites are incompletely known. Here, we show that T. gondii protein phosphatase 2A (PP2A) holoenzyme is composed of a catalytic subunit PP2A-C, a scaffold subunit PP2A-A and a regulatory subunit PP2A-B. Disruption of any of these subunits increased starch accumulation and blocked the tachyzoite-to-bradyzoite differentiation. PP2A contributes to the regulation of amylopectin metabolism via dephosphorylation of calcium-dependent protein kinase 2 at S679. Phosphoproteomics identified several putative PP2A holoenzyme substrates that are involved in bradyzoite differentiation. Our findings provide novel insight into the role of PP2A as a key regulator of starch metabolism and bradyzoite differentiation in T. gondii.

Importance of aspartyl protease 5 in the establishment of the intracellular niche during acute and chronic infection of Toxoplasma gondii

Virulence and persistence of the obligate intracellular parasite Toxoplasma gondii involve the secretion of effector proteins belonging to the family of dense granule proteins (GRAs) that act notably as modulators of the host defense mechanisms and participate in cyst wall formation. The subset of GRAs residing in the parasitophorous vacuole (PV) or exported into the host cell, undergo proteolytic cleavage in the Golgi upon the action of the aspartyl protease 5 (ASP5). In tachyzoites, ASP5 substrates play central roles in the morphology of the PV and the export of effectors across the translocon complex MYR1/2/3. Here, we used N-terminal amine isotopic labeling of substrates to identify novel ASP5 cleavage products by comparing the N-terminome of wild-type and Δasp5 lines in tachyzoites and bradyzoites. Validated substrates reside within the PV or PVM in an ASP5-dependent manner. Remarkably, Δasp5 bradyzoites are impaired in the formation of the cyst wall in vitro and exhibit a considerably reduced cyst burden in chronically infected animals. More specifically two-photon serial tomography of infected mouse brains revealed a comparatively reduced number and size of the cysts throughout the establishment of persistence in the absence of ASP5.

Identification and Characterization of Eimeria tenella Rhoptry Protein 35 (EtROP35)

Simple Summary

Rhoptry proteins (ROPs) of phylum Apicomplexa parasites are important secretory virulence factors as well as candidate vaccines. However, studies on ROPs of Eimeria tenella are limited. In this study, the coding sequence of E. tenella rhoptry protein 35 (EtROP35) was cloned, and then its localization, expression in parasite, potential role within invasion, and protective efficacy were investigated. Sequence analysis and subcellular localization revealed that EtROP35 is a rhoptry protein of E. tenella. Sporozoite invasion-blocking assay and protective efficacy indicated that EtROP35 might be involved in the parasite invasion process and may be a potential vaccine candidate against E. tenella.

Abstract

Rhoptry proteins (ROPs) of Apicomplexa are crucial secreted virulence factors and sources of vaccine candidates. To date, Eimeria tenella ROPs are not well studied. This study identified and characterized a novel E. tenella ROP (EtROP35), which showed the highest levels among 28 putative ROPs in previous sporozoite and merozoite transcriptomes. Sequence analysis showed that EtROP35 contains an N-terminal secretory signal and a protein kinase domain including eight conserved ROP35-subfamily motifs. Subsequent experiments confirmed that it is a secretory protein. Subcellular localization revealed it localized at the apical end of the sporozoites and merozoites, which was consistent with the ROPs of other Apicomplexan parasites. To further understand the biological meaning of EtROP35, expression levels in different developmental stages and sporozoite invasion-blocking assay were investigated. EtROP35 showed significantly higher levels in sporozoites (6.23-fold) and merozoites (7.00-fold) than sporulated oocysts. Sporozoite invasion-blocking assay revealed that anti-EtROP35 polyclonal antibody significantly reduced the sporozoite invasion rate, suggesting it might participate in host cell invasion and be a viable choice as a vaccine candidate. The immunological protective assays showed that EtROP35 could induce a high level of serum IgY and higher mean body weight gain, and lower cecum lesion score and oocysts excretion than the challenged control group. These data indicated that EtROP35 had good immunogenicity and may be a promising vaccine candidate against E. tenella.

Characterization of the Eimeria tenella rhoptry protein with a nuclear localization sequence (EtROP30)

Rhoptry proteins (ROPs), secreted by specific rhoptry organelles of apicomplexan parasites, are determinants of parasite pathogenesis and sources of vaccine candidates. Twenty-eight ROPs of Eimeria tenella have been predicted by genomic approaches, and in the present study, E. tenella rhoptry protein 30 (EtROP30) was characterized. Subcellular localizations of EtROP30 in sporozoites and merozoites were in the apical complex and rhoptry-like bulb, suggesting that EtROP30 is a member of ROPs in E. tenella. Sequence analysis showed that EtROP30 contained an N-terminal secretory signal, a protein kinase domain with eight E. tenella-specific rhoptry kinase 1 subfamily (ROPK-Eten1) motifs, and a C-terminal nuclear localization sequence (NLS), making EtROP30 the only ROP that contains both a secretory signal and an NLS in E. tenella. Subsequent experiments showed that EtROP30 was a secreted protein in the sporozoite stage, relying on NLS for migration to the host nucleus. In addition, EtROP30 showed significantly higher expression levels in the parasite merozoite stage, indicating that EtROP30 plays a critical role during parasite reinvasion and development and may be a viable option as a vaccine candidate for anti-parasitic infection. The immunization protection efficacies of EtROP30 were evaluated. Significant improvements in mean body weight gain, reduction of cecum lesion score, and number of oocysts excreted were observed, indicating that EtROP30 has good immunogenicity against E. tenella. In the present study, a ROP of E. tenella with secretory and nuclear localization characteristics has been identified, and proved to be an effective vaccine candidate against this parasite.

Computational tools and resources for pseudokinase research

Pseudokinases regulate diverse cellular processes associated with normal cellular functions and disease. They are defined bioinformatically based on the absence of one or more catalytic residues that are required for canonical protein kinase functions. The ability to define pseudokinases based on primary sequence comparison has enabled the systematic mapping and cataloging of pseudokinase orthologs across the tree of life. While these sequences contain critical information regarding pseudokinase evolution and functional specialization, extracting this information and generating testable hypotheses based on integrative mining of sequence and structural data requires specialized computational tools and resources. In this chapter, we review recent advances in the development and application of open-source tools and resources for pseudokinase research. Specifically, we describe the application of an interactive data analytics framework, KinView, for visualizing the patterns of conservation and variation in the catalytic domain motifs of pseudokinases and evolutionarily related canonical kinases using a consistent set of curated alignments organized based on the widely used kinome evolutionary hierarchy. We also demonstrate the application of an integrated Protein Kinase Ontology (ProKinO) and an interactive viewer, ProtVista, for mapping and analyzing primary sequence motifs and annotations in the context of 3D structures and AlphaFold2 models. We provide examples and protocols for generating testable hypotheses on pseudokinase functions both for bench biologists and advanced users.

Dual RNA-seq transcriptome analysis of caecal tissue during primary Eimeria tenella infection in chickens

BACKGROUND

Coccidiosis is an infectious disease with large negative impact on the poultry industry worldwide. It is an enteric infection caused by unicellular Apicomplexan parasites of the genus Eimeria. The present study aimed to gain more knowledge about interactions between parasites and the host immune system during the early asexual replication phase of E. tenella in chicken caeca. For this purpose, chickens were experimentally infected with E. tenella oocysts, sacrificed on days 1-4 and 10 after infection and mRNA from caecal tissues was extracted and sequenced.

RESULTS

Dual RNA-seq analysis revealed time-dependent changes in both host and parasite gene expression during the course of the infection. Chicken immune activation was detected from day 3 and onwards with the highest number of differentially expressed immune genes recorded on day 10. Among early (days 3-4) responses up-regulation of genes for matrix metalloproteinases, several chemokines, interferon (IFN)-γ along with IFN-stimulated genes GBP, IRF1 and RSAD2 were noted. Increased expression of genes with immune suppressive/regulatory effects, e.g. IL10, SOCS1, SOCS3, was also observed among early responses. For E. tenella a general up-regulation of genes involved in protein expression and energy metabolism as well as a general down-regulation genes for DNA and RNA processing were observed during the infection. Specific E. tenella genes with altered expression during the experiment include those for proteins in rhoptry and microneme organelles.

CONCLUSIONS

The present study provides novel information on both the transcriptional activity of E. tenella during schizogony in ceacal tissue and of the local host responses to parasite invasion during this phase of infection. Results indicate a role for IFN-γ and IFN-stimulated genes in the innate defence against Eimeria replication.

Do All Coccidia Follow the Same Trafficking Rules?

The Coccidia are a subclass of the Apicomplexa and include several genera of protozoan parasites that cause important diseases in humans and animals, with Toxoplasma gondii becoming the ‘model organism’ for research into the coccidian molecular and cellular processes. The amenability to the cultivation of T. gondii tachyzoites and the wide availability of molecular tools for this parasite have revealed many mechanisms related to their cellular trafficking and roles of parasite secretory organelles, which are critical in parasite-host interaction. Nevertheless, the extrapolation of the T. gondii mechanisms described in tachyzoites to other coccidian parasites should be done carefully. In this review, we considered published data from Eimeria parasites, a coccidian genus comprising thousands of species whose infections have important consequences in livestock and poultry. These studies suggest that the Coccidia possess both shared and diversified mechanisms of protein trafficking and secretion potentially linked to their lifecycles. Whereas trafficking and secretion appear to be well conversed prior to and during host-cell invasion, important differences emerge once endogenous development commences. Therefore, further studies to validate the mechanisms described in T. gondii tachyzoites should be performed across a broader range of coccidians (including T. gondii sporozoites). In addition, further genus-specific research regarding important disease-causing Coccidia is needed to unveil the individual molecular mechanisms of pathogenesis related to their specific lifecycles and hosts.

Genome-Wide Expression Patterns of Rhoptry Kinases during the Eimeria tenella Life-Cycle

Kinome from apicomplexan parasites is composed of eukaryotic protein kinases and Apicomplexa specific kinases, such as rhoptry kinases (ROPK). Ropk is a gene family that is known to play important roles in host–pathogen interaction in Toxoplasma gondii but is still poorly described in Eimeria tenella, the parasite responsible for avian coccidiosis worldwide. In the E. tenella genome, 28 ropk genes are predicted and could be classified as active (n = 7), inactive (incomplete catalytic triad, n = 12), and non-canonical kinases (active kinase with a modified catalytic triad, n = 9). We characterized the ropk gene expression patterns by real-time quantitative RT-PCR, normalized by parasite housekeeping genes, during the E. tenella life-cycle. Analyzed stages were: non-sporulated oocysts, sporulated oocysts, extracellular and intracellular sporozoites, immature and mature schizonts I, first- and second-generation merozoites, and gametes. Transcription of all those predicted ropk was confirmed. The mean intensity of transcription was higher in extracellular stages and 7–9 ropk were specifically transcribed in merozoites in comparison with sporozoites. Transcriptional profiles of intracellular stages were closely related to each other, suggesting a probable common role of ROPKs in hijacking signaling pathways and immune responses in infected cells. These results provide a solid basis for future functional analysis of ROPK from E. tenella.

Structural insights into an atypical secretory pathway kinase crucial for Toxoplasma gondii invasion

Active host cell invasion by the obligate intracellular apicomplexan parasites relies on the formation of a moving junction, which connects parasite and host cell plasma membranes during entry. Invading Toxoplasma gondii tachyzoites secrete their rhoptry content and insert a complex of RON proteins on the cytoplasmic side of the host cell membrane providing an anchor to which the parasite tethers. Here we show that a rhoptry-resident kinase RON13 is a key virulence factor that plays a crucial role in host cell entry. Cryo-EM, kinase assays, phosphoproteomics and cellular analyses reveal that RON13 is a secretory pathway kinase of atypical structure that phosphorylates rhoptry proteins including the components of the RON complex. Ultimately, RON13 kinase activity controls host cell invasion by anchoring the moving junction at the parasite-host cell interface.

Biogenesis and discharge of the rhoptries: Key organelles for entry and hijack of host cells by the Apicomplexa

Rhoptries are specialized secretory organelles found in the Apicomplexa phylum, playing a central role in the establishment of parasitism. The rhoptry content includes membranous as well as proteinaceous materials that are discharged into the host cell in a regulated fashion during parasite entry. A set of rhoptry neck proteins form a RON complex that critically participates in the moving junction formation during invasion. Some of the rhoptry bulb proteins are associated with the membranous materials and contribute to the formation of the parasitophorous vacuole membrane while others are targeted into the host cell including the nucleus to subvert cellular functions. Here, we review the recent studies on Toxoplasma and Plasmodium parasites that shed light on the key steps leading to rhoptry biogenesis, trafficking, and discharge.

Dual RNA-Seq transcriptome analysis of chicken macrophage-like cells (HD11) infected in vitro with Eimeria tenella

The study aimed to monitor parasite and host gene expression during the early stages of Eimeria tenella infection of chicken cells using dual RNA-Seq analysis. For this, we used chicken macrophage-like cell line HD11 cultures infected in vitro with purified E. tenella sporozoites. Cultures were harvested between 2 and 72 h post-infection and mRNA was extracted and sequenced. Dual RNA-Seq analysis showed clear patterns of altered expression for both parasite and host genes during infection. For example, genes in the chicken immune system showed upregulation early (2–4 h), a strong downregulation of genes across the immune system at 24 h and a repetition of early patterns at 72 h, indicating that invasion by a second generation of parasites was occurring. The observed downregulation may be due to immune self-regulation or to immune evasive mechanisms exerted by E. tenella. Results also suggested pathogen recognition receptors involved in E. tenella innate recognition, MRC2, TLR15 and NLRC5 and showed distinct chemokine and cytokine induction patterns. Moreover, the expression of several functional categories of Eimeria genes, such as rhoptry kinase genes and microneme genes, were also examined, showing distinctive differences which were expressed in sporozoites and merozoites.

Deciphering Molecular Virulence Mechanism of Mycobacterium tuberculosis Dop isopeptidase Based on Its Sequence-Structure-Function Linkage

The pupylation pathway marks proteins for prokaryotic ubiquitin-like protein (Pup)-proteasomal degradation and survival strategy of mycobacteria inside of the host macrophages. Deamidase of Pup (Dop) plays a central role in the pupylation pathway. It is still a matter of investigation to know the function of Dop in virulence of mycobacterial lineage. Hence, the present study was intended to describe the sequence-structure-function-virulence link of Dop for understanding the molecular virulence mechanism of Mycobacterium tuberculosis H37Rv (Mtb). Phylogenetic analysis of this study indicated that Dop has extensively diverged across the proteasome-harboring bacteria. The functional part of Dop was converged across the pathogenic mycobacterial lineage. The genome-wide analysis pointed out that the pupylation gene locus was identical to each other, but its genome neighborhood differed from species to species. Molecular modeling and dynamic studies proved that the predicted structure of Mtb Dop was energetically stable and low conformational freedom. Moreover, evolutionary constraints in Mtb Dop were intensively analyzed for inferring its sequence-structure-function relationships for the full virulence of Mtb. It indicated that evolutionary optimization was extensively required to stabilize its local structural environment at the side chains of mutable residues. The sequence-structure-function-virulence link of Dop might have retained in Mtb by reordering hydrophobic and hydrogen bonding patterns in the local structural environment. Thus, the results of our study provide a quest to understand the molecular virulence and pathogenesis mechanisms of Mtb during the infection process.

A Novel Rhoptry Protein as Candidate Vaccine against Eimeria tenella Infection

Eimeria tenella (E. tenella) is a highly pathogenic and prevalent species of Eimeria that infects chickens, and it causes a considerable disease burden worldwide. The secreted proteins and surface antigens of E. tenella at the sporozoite stage play an essential role in the host–parasite interaction, which involves attachment and invasion, and these interactions are considered vaccine candidates based on the strategy of cutting off the invasion pathway to interrupt infection. We selected two highly expressed surface antigens (SAGs; Et-SAG13 and Et-SAG) and two highly expressed secreted antigens (rhoptry kinases Eten5-A, Et-ROPK-Eten5-A and dense granule 12, Et-GRA12) at the sporozoite stage. Et-ROPK-Eten5-A and Et-GRA12 were two unexplored proteins. Et-ROPK-Eten5-A was an E. tenella-specific rhoptry (ROP) protein and distributed in the apical pole of sporozoites and merozoites. Et-GRA12 was scattered in granular form at the sporozoite stage. To evaluate the potential of rEt-ROPK-Eten5-A, rEt-GRA12, rEt-SAG13 and rEt-SAG proteins as a coccidiosis vaccine, the protective efficacy was examined based on survival rate, lesion score, body weight gain, relative body weight gain and oocyst output. The survival rate was significantly improved in rEt-ROPK-Eten5-A (100%) and rEt-GRA12 (100%) immune chickens compared to the challenged control group (40%). The average body weight gains of rEt-ROPK-Eten5-A, rEt-GRA12, rEt-SAG13 and rEt-SAG immunized chickens were significantly higher than those of unimmunized chickens. The mean lesion score and oocyst output of the rEt-ROPK-Eten5-A immunized chickens were significantly reduced compared to unimmunized challenged chickens. These results suggest that the rEt-ROPK-Eten5-A protein effectively triggered protection against E. tenella in chickens and provides a useful foundation for future work developing anticoccidial vaccines.

Species diversity and genome evolution of the pathogenic protozoan parasite, Neospora caninum

Neospora caninum is a cyst-forming coccidian parasite of veterinary and economical significance, affecting dairy and beef cattle industries on a global scale. Comparative studies suggest that N. caninum consists of a globally dispersed, diverse population of lineages, distinguished by their geographical origin, broad host range, and phenotypic features. This viewpoint is however changing. While intraspecies diversity, and more specifically pathogenic variability, has been experimentally demonstrated in a myriad of studies, the underlying contributors and sources responsible for such diversity have remained nebulous. However, recent large-scale sequence and bioinformatics studies have aided in revealing intrinsic genetic differences distinguishing isolates of this species, that await further characterisation as causative links to virulence and pathogenicity. Furthermore, progress on N. caninum research as a non-model organism is hindered by a lack of robust, annotated genomic, transcriptomic, and proteomic data for the species, especially compared to other thoroughly studied Apicomplexa such as Toxoplasma gondii and Plasmodium species. This review explores the current body of knowledge on intra-species diversity within N. caninum. This includes the contribution of sequence variants in both coding and non-coding regions, the presence of genome polymorphic hotspots, and the identification of non-synonymous mutations. The implications of such diversity on important parasite phenotypes such as pathogenicity and population structure are also discussed. Lastly, the identification of potential virulence factors from both in-silico and next generation sequencing studies is examined, offering new insights into potential avenues for future research on neosporosis.

Emerging roles of the αC-β4 loop in protein kinase structure, function, evolution, and disease

The faithful propagation of cellular signals in most organisms relies on the coordinated functions of a large family of protein kinases that share a conserved catalytic domain. The catalytic domain is a dynamic scaffold that undergoes large conformational changes upon activation. Most of these conformational changes, such as movement of the regulatory αC-helix from an "out" to "in" conformation, hinge on a conserved, but understudied, loop termed the αC-β4 loop, which mediates conserved interactions to tether flexible structural elements to the kinase core. We previously showed that the αC-β4 loop is a unique feature of eukaryotic protein kinases. Here, we review the emerging roles of this loop in kinase structure, function, regulation, and diseases. Through a kinome-wide analysis, we define the boundaries of the loop for the first time and show that sequence and structural variation in the loop correlate with conformational and regulatory variation. Many recurrent disease mutations map to the αC-β4 loop and contribute to drug resistance and abnormal kinase activation by relieving key auto-inhibitory interactions associated with αC-helix and inter-lobe movement. The αC-β4 loop is a hotspot for post-translational modifications, protein-protein interaction, and Hsp90 mediated folding. Our kinome-wide analysis provides insights for hypothesis-driven characterization of understudied kinases and the development of allosteric protein kinase inhibitors.

Life cycle stages, specific organelles and invasion mechanisms of Eimeria species

Apicomplexans, including species of Eimeria, pose a real threat to the health and wellbeing of animals and humans. Eimeria parasites do not infect humans but cause an important economic impact on livestock, in particular on the poultry industry. Despite its high prevalence and financial costs, little is known about the cell biology of these 'cosmopolitan' parasites found all over the world. In this review, we discuss different aspects of the life cycle and stages of Eimeria species, focusing on cellular structures and organelles typical of the coccidian family as well as genus-specific features, complementing some 'unknowns' with what is described in the closely related coccidian Toxoplasma gondii.

Enzyme-Linked Aptamer Assay (ELAA) for Detection of Toxoplasma ROP18 Protein in Human Serum

Toxoplasma gondii engenders the common parasitic disease toxoplasmosis in almost all warm-blooded animals. Being a critical secretory protein, ROP18 is a major virulence factor of Toxoplasma. There are no reports about ROP18 detection in human serum samples with different clinical manifestations. New aptamers against ROP18 protein were developed through Systematic Evolution of Ligands by Exponential enrichment (SELEX). An Enzyme-Linked Aptamer Assay (ELAA) platform was developed using SELEX-derived aptamers, namely AP001 and AP002. The ELAA was used to evaluate total antigen from T. gondii RH strain (RH Ag) and recombinant protein of ROP18 (rROP18). The results showed that the ELAA presented higher affinity and specificity to RH Ag and rROP18, compared to negative controls. Detection limit of rROP18 protein in serum samples was measured by standard addition method, achieving a lower concentration of 1.56 μg/mL. Moreover, 62 seropositive samples with different clinical manifestations of toxoplasmosis and 20 seronegative samples were tested. A significant association between ELAA test positive for human serum samples and severe congenital toxoplasmosis was found (p = 0.006). Development and testing of aptamers-based assays opens a window for low-cost and rapid tests looking for biomarkers and improves our understanding about the role of ROP18 protein on the pathogenesis of human toxoplasmosis.

The secreted kinase ROP17 promotes Toxoplasma gondii dissemination by hijacking monocyte tissue migration

The protozoan parasite Toxoplasma gondii is thought to exploit monocyte trafficking to facilitate dissemination across endothelial barriers such as the blood-brain barrier. Here we analyzed the migration of parasitized monocytes in model endothelial and interstitial environments. We report that infection enhanced monocyte locomotion on the surface of endothelial cells, but profoundly inhibited monocyte transmigration across endothelial barriers. In contrast, infection robustly increased monocyte and macrophage migration through collagen-rich tissues in a Rho/ROCK-dependent manner consistent with integrin-independent interstitial migration. We further demonstrated that the secreted T. gondii protein kinase ROP17 was required for enhanced tissue migration. In vivo, ROP17-deficient parasites failed to upregulate monocyte tissue migration and exhibited an early dissemination delay, leading to prolonged mouse survival. Our findings indicate that the parasite-induced changes in monocyte motility likely facilitate the transport of T. gondii through tissues and promote systemic dissemination, rather than shuttle parasites across the blood-brain barrier via extravasation.

LF4/MOK and a CDK-related kinase regulate the number and length of cilia in Tetrahymena

The length of cilia is controlled by a poorly understood mechanism that involves members of the conserved RCK kinase group, and among them, the LF4/MOK kinases. The multiciliated protist model, Tetrahymena, carries two types of cilia (oral and locomotory) and the length of the locomotory cilia is dependent on their position with the cell. In Tetrahymena, loss of an LF4/MOK ortholog, LF4A, lengthened the locomotory cilia, but also reduced their number. Without LF4A, cilia assembled faster and showed signs of increased intraflagellar transport (IFT). Consistently, overproduced LF4A shortened cilia and downregulated IFT. GFP-tagged LF4A, expressed in the native locus and imaged by total internal reflection microscopy, was enriched at the basal bodies and distributed along the shafts of cilia. Within cilia, most LF4A-GFP particles were immobile and a few either diffused or moved by IFT. We suggest that the distribution of LF4/MOK along the cilium delivers a uniform dose of inhibition to IFT trains that travel from the base to the tip. In a longer cilium, the IFT machinery may experience a higher cumulative dose of inhibition by LF4/MOK. Thus, LF4/MOK activity could be a readout of cilium length that helps to balance the rate of IFT-driven assembly with the rate of disassembly at steady state. We used a forward genetic screen to identify a CDK-related kinase, CDKR1, whose loss-of-function suppressed the shortening of cilia caused by overexpression of LF4A, by reducing its kinase activity. Loss of CDKR1 alone lengthened both the locomotory and oral cilia. CDKR1 resembles other known ciliary CDK-related kinases: LF2 of Chlamydomonas, mammalian CCRK and DYF-18 of C. elegans, in lacking the cyclin-binding motif and acting upstream of RCKs. The new genetic tools we developed here for Tetrahymena have potential for further dissection of the principles of cilia length regulation in multiciliated cells.

Eimeria tenella ROP kinase EtROP1 induces G0/G1 cell cycle arrest and inhibits host cell apoptosis

Coccidia are obligate intracellular protozoan parasites responsible for human and veterinary diseases. Eimeria tenella, the aetiologic agent of caecal coccidiosis, is a major pathogen of chickens. In Toxoplasma gondii, some kinases from the rhoptry compartment (ROP) are key virulence factors. ROP kinases hijack and modulate many cellular functions and pathways, allowing T. gondii survival and development. E. tenella's kinome comprises 28 putative members of the ROP kinase family; most of them are predicted, as pseudokinases and their functions have never been characterised. One of the predicted kinase, EtROP1, was identified in the rhoptry proteome of E. tenella sporozoites. Here, we demonstrated that EtROP1 is active, and the N-terminal extension is necessary for its catalytic kinase activity. Ectopic expression of EtROP1 followed by co-immunoprecipitation identified cellular p53 as EtROP1 partner. Further characterisation confirmed the interaction and the phosphorylation of p53 by EtROP1. E. tenella infection or overexpression of EtROP1 resulted both in inhibition of host cell apoptosis and G0/G1 cell cycle arrest. This work functionally described the first ROP kinase from E. tenella and its noncanonical structure. Our study provides the first mechanistic insight into host cell apoptosis inhibition by E. tenella. EtROP1 appears as a new candidate for coccidiosis control.

Divergent kinase regulates membrane ultrastructure of the Toxoplasma parasitophorous vacuole

Apicomplexan parasites replicate within a protective organelle, called the parasitophorous vacuole (PV). The Toxoplasma gondii PV is filled with a network of tubulated membranes, which are thought to facilitate trafficking of effectors and nutrients. Despite being critical to parasite virulence, there is scant mechanistic understanding of the network's functions. Here, we identify the parasite-secreted kinase WNG1 (With-No-Gly-loop) as a critical regulator of tubular membrane biogenesis. WNG1 family members adopt an atypical protein kinase fold lacking the glycine rich ATP-binding loop that is required for catalysis in canonical kinases. Unexpectedly, we find that WNG1 is an active protein kinase that localizes to the PV lumen and phosphorylates PV-resident proteins, several of which are essential for the formation of a functional intravacuolar network. Moreover, we show that WNG1-dependent phosphorylation of these proteins is required for their membrane association, and thus their ability to tubulate membranes. Consequently, WNG1 knockout parasites have an aberrant PV membrane ultrastructure. Collectively, our results describe a unique family of Toxoplasma kinases and implicate phosphorylation of secreted proteins as a mechanism of regulating PV development during parasite infection.

Two Antagonistic Hippo Signaling Circuits Set the Division Plane at the Medial Position in the Ciliate Tetrahymena

In a single cell, ciliates maintain a complex pattern of cortical organelles that are arranged along the anteroposterior and circumferential axes. The underlying molecular mechanisms of intracellular pattern formation in ciliates are largely unknown. Ciliates divide by tandem duplication, a process that remodels the parental cell into two daughters aligned head-to-tail. In the elo1-1 mutant of Tetrahymena thermophila, the segmentation boundary/division plane forms too close to the posterior end of the parental cell, producing a large anterior and a small posterior daughter cell, respectively. We show that ELO1 encodes a Lats/NDR kinase that marks the posterior segment of the cell cortex, where the division plane does not form in the wild-type. Elo1 acts independently of CdaI, a Hippo/Mst kinase that marks the anterior half of the parental cell, and whose loss shifts the division plane anteriorly. We propose that, in Tetrahymena, two antagonistic Hippo circuits focus the segmentation boundary/division plane at the equatorial position, by excluding divisional morphogenesis from the cortical areas that are too close to cell ends.

Genome Wide Identification of Mutational Hotspots in the Apicomplexan Parasite Neospora caninum and the Implications for Virulence

Neospora caninum is an apicomplexan parasite responsible for neosporosis, a disease causing hind limb paralysis in dogs and abortion in cattle, resulting in substantial economic losses to beef and dairy industries. Marked differences in pathogenicity exist between N. caninum strains suggesting that intrinsic genetic differences exist between them. These differences likely exist in genes expressed during the tachyzoite lifecycle stage which is responsible for the pathogenesis of neosporosis. An improved understanding of these genetic differences is essential to understanding N. caninum virulence, though such knowledge is scarce. Using a variant detection workflow we compared the tachyzoite transcriptomes of two N. caninum strains with different virulence properties: NC-Liverpool (virulent) and NC-Nowra (avirulent). This workflow identified 3130 SNPs and 6123 indels between the strains, and nine markers capturing 30 variants were Sanger sequenced for both strains. Sequencing of these loci was extended to an additional eight strains and subsequent phylogenetic analysis supported a genetic population structure comprised of two major clades with no geographical segregation. Sequence polymorphisms within coding regions of tachyzoite-associated genes were concentrated on chromosomes XI and XII, with 19 distinct tachyzoite-associated SNP hotspot regions identified within coding regions of the N. caninum nuclear genome. The variants were predominantly located in loci associated with protein binding, protein-protein interactions, transcription, and translation. Furthermore, 468 nonsynonymous SNPs identified within protein-coding genes were associated with protein kinase activity, protein binding, protein phosphorylation, and proteolysis. This work may implicate these processes and the specific proteins involved as novel effectors of N. caninum tachyzoite virulence.

Characterization of Ligand Binding to Pseudokinases Using a Thermal Shift Assay

The protocol herein describes a robust and proven method for the measurement of pseudokinase-ligand interaction using a fluorescence-based thermal shift assay (TSA). Pseudokinases are kinase-like proteins that have recently emerged as crucial regulatory modules of signal transduction pathways and may well represent a novel class of drug targets. However, unlike kinases, the regulatory activity of pseudokinases is mainly conferred through protein-protein interactions. Understanding the mechanisms that underlie pseudokinase conformational changes through ligand binding and how such conformational changes can tune signaling pathways is a necessary step to unravel their biological functions.Thermal denaturation-based methods have proven to be a powerful method for determining the capacity of purified recombinant pseudokinases to bind ligands and can simultaneously inform on the potential druggability of the nucleotide-binding site. This assay takes advantage of a change in fluorescence arising when the dye, SYPRO Orange, binds to hydrophobic patches that become exposed when a protein undergoes thermal unfolding. Ligand binding to a protein is known to increase its thermal stability, which is reflected by a shift between the thermal denaturation curves of the unliganded protein and the liganded protein. Here, we illustrate the utility of the method with the pseudokinases, ErbB3/HER3, ILK, ROP5Bi, JAK1, JAK2, TYK2, MLKL, STRAD, TRIB1, VRK3, and ROR1. This method can also be used to determine optimal buffer conditions that may increase protein stability and can be tailored to other protein families.

Recognition of sites of functional specialisation in all known eukaryotic protein kinase families

The conserved function of protein phosphorylation, catalysed by members of protein kinase superfamily, is regulated in different ways in different kinase families. Further, differences in activating triggers, cellular localisation, domain architecture and substrate specificity between kinase families are also well known. While the transfer of γ-phosphate from ATP to the hydroxyl group of Ser/Thr/Tyr is mediated by a conserved Asp, the characteristic functional and regulatory sites are specialized at the level of families or sub-families. Such family-specific sites of functional specialization are unknown for most families of kinases. In this work, we systematically identify the family-specific residue features by comparing the extent of conservation of physicochemical properties, Shannon entropy and statistical probability of residue distributions between families of kinases. An integrated discriminatory score, which combines these three features, is developed to demarcate the functionally specialized sites in a kinase family from other sites. We achieved an area under ROC curve of 0.992 for the discrimination of kinase families. Our approach was extensively tested on well-studied families CDK and MAPK, wherein specific protein interaction sites and substrate recognition sites were successfully detected (p-value < 0.05). We also find that the known family-specific oncogenic driver mutation sites were scored high by our method. The method was applied to all known kinases encompassing 107 families from diverse eukaryotic organisms leading to a comprehensive list of family-specific functional sites. Apart from other uses, our method facilitates identification of specific protein interaction sites and drug target sites in a kinase family.

Dual RNA-seq reveals no plastic transcriptional response of the coccidian parasite Eimeria falciformis to host immune defenses

Background

Parasites can either respond to differences in immune defenses that exist between individual hosts plastically or, alternatively, follow a genetically canalized (“hard wired”) program of infection. Assuming that large-scale functional plasticity would be discernible in the parasite transcriptome we have performed a dual RNA-seq study of the lifecycle of Eimeria falciformis using infected mice with different immune status as models for coccidian infections.

Results

We compared parasite and host transcriptomes (dual transcriptome) between naïve and challenge infected mice, as well as between immune competent and immune deficient ones. Mice with different immune competence show transcriptional differences as well as differences in parasite reproduction (oocyst shedding). Broad gene categories represented by differently abundant host genes indicate enrichments for immune reaction and tissue repair functions. More specifically, TGF-beta, EGF, TNF and IL-1 and IL-6 are examples of functional annotations represented differently depending on host immune status. Much in contrast, parasite transcriptomes were neither different between Coccidia isolated from immune competent and immune deficient mice, nor between those harvested from naïve and challenge infected mice. Instead, parasite transcriptomes have distinct profiles early and late in infection, characterized largely by biosynthesis or motility associated functional gene groups, respectively. Extracellular sporozoite and oocyst stages showed distinct transcriptional profiles and sporozoite transcriptomes were found enriched for species specific genes and likely pathogenicity factors.

Conclusion

We propose that the niche and host-specific parasite E. falciformis uses a genetically canalized program of infection. This program is likely fixed in an evolutionary process rather than employing phenotypic plasticity to interact with its host. This in turn might limit the potential of the parasite to adapt to new host species or niches, forcing it to coevolve with its host.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4095-6) contains supplementary material, which is available to authorized users.

Metal coordination in kinases and pseudokinases

Protein phosphorylation, mediated by protein kinases, is a key event in the regulation of eukaryotic signal transduction. The majority of eukaryotic protein kinases perform phosphoryl transfer, assisted by two divalent metal ions. About 10% of all human protein kinases are, however, thought to be catalytically inactive. These kinases lack conserved residues of the kinase core and are classified as pseudokinases. Yet, it has been demonstrated that pseudokinases are critically involved in biological functions. Here, we show how pseudokinases have developed strategies by modifying amino acid residues in order to achieve stable, active-like conformations. This includes binding of the co-substrate ATP in a two metal-, one metal- or even no metal-binding mode. Examples of the respective pseudokinases are provided on a structural basis and compared with a canonical protein kinase, Protein Kinase A. Moreover, the functional roles of both independent metal-binding sites, Me1 and Me2, are discussed. Lack of phosphotransferase activity does not implicate a loss of function and can easily point to alternative roles of pseudokinases, i.e. acting as switches or scaffolds, and having evolved as components crucial for cellular cross-talk and signaling. Interestingly, pseudokinases are present in all kingdoms of life and their specific roles remain enigmatic. More studies are needed to unravel the crucial functions of those interesting proteins.

Genomics of apicomplexan parasites

The increasing prevalence of infections involving intracellular apicomplexan parasites such as Plasmodium, Toxoplasma, and Cryptosporidium (the causative agents of malaria, toxoplasmosis, and cryptosporidiosis, respectively) represent a significant global healthcare burden. Despite their significance, few treatments are available; a situation that is likely to deteriorate with the emergence of new resistant strains of parasites. To lay the foundation for programs of drug discovery and vaccine development, genome sequences for many of these organisms have been generated, together with large-scale expression and proteomic datasets. Comparative analyses of these datasets are beginning to identify the molecular innovations supporting both conserved processes mediating fundamental roles in parasite survival and persistence, as well as lineage-specific adaptations associated with divergent life-cycle strategies. The challenge is how best to exploit these data to derive insights into parasite virulence and identify those genes representing the most amenable targets. In this review, we outline genomic datasets currently available for apicomplexans and discuss biological insights that have emerged as a consequence of their analysis. Of particular interest are systems-based resources, focusing on areas of metabolism and host invasion that are opening up opportunities for discovering new therapeutic targets.

Toxoplasma Effectors Targeting Host Signaling and Transcription

Early electron microscopy studies revealed the elaborate cellular features that define the unique adaptations of apicomplexan parasites. Among these were bulbous rhoptry (ROP) organelles and small, dense granules (GRAs), both of which are secreted during invasion of host cells. These early morphological studies were followed by the exploration of the cellular contents of these secretory organelles, revealing them to be comprised of highly divergent protein families with few conserved domains or predicted functions. In parallel, studies on host-pathogen interactions identified many host signaling pathways that were mysteriously altered by infection. It was only with the advent of forward and reverse genetic strategies that the connections between individual parasite effectors and the specific host pathways that they targeted finally became clear. The current repertoire of parasite effectors includes ROP kinases and pseudokinases that are secreted during invasion and that block host immune pathways. Similarly, many secretory GRA proteins alter host gene expression by activating host transcription factors, through modification of chromatin, or by inducing small noncoding RNAs. These effectors highlight novel mechanisms by which T. gondii has learned to harness host signaling to favor intracellular survival and will guide future studies designed to uncover the additional complexity of this intricate host-pathogen interaction.

An in vitro model of intestinal infection reveals a developmentally regulated transcriptome of Toxoplasma sporozoites and a NF-κB-like signature in infected host cells

Toxoplasmosis is a zoonotic infection affecting approximately 30% of the world’s human population. After sexual reproduction in the definitive feline host, Toxoplasma oocysts, each containing 8 sporozoites, are shed into the environment where they can go on to infect humans and other warm-blooded intermediate hosts. Here, we use an in vitro model to assess host transcriptomic changes that occur in the earliest stages of such infections. We show that infection of rat intestinal epithelial cells with mature sporozoites primarily results in higher expression of genes associated with Tumor Necrosis Factor alpha (TNFα) signaling via NF-κB. Furthermore, we find that, consistent with their biology, these mature, invaded sporozoites display a transcriptome intermediate between the previously reported day 10 oocysts and that of their tachyzoite counterparts. Thus, this study uncovers novel host and pathogen factors that may be critical for the establishment of a successful intracellular niche following sporozoite-initiated infection.

Genome-Wide Identification and Evolutionary Analysis of Sarcocystis neurona Protein Kinases

The apicomplexan parasite Sarcocystis neurona causes equine protozoal myeloencephalitis (EPM), a degenerative neurological disease of horses. Due to its host range expansion, S. neurona is an emerging threat that requires close monitoring. In apicomplexans, protein kinases (PKs) have been implicated in a myriad of critical functions, such as host cell invasion, cell cycle progression and host immune response evasion. Here, we used various bioinformatics methods to define the kinome of S. neurona and phylogenetic relatedness of its PKs to other apicomplexans. We identified 97 putative PKs clustering within the various eukaryotic kinase groups. Although containing the universally-conserved PKA (AGC group), S. neurona kinome was devoid of PKB and PKC. Moreover, the kinome contains the six-conserved apicomplexan CDPKs (CAMK group). Several OPK atypical kinases, including ROPKs 19A, 27, 30, 33, 35 and 37 were identified. Notably, S. neurona is devoid of the virulence-associated ROPKs 5, 6, 18 and 38, as well as the Alpha and RIO kinases. Two out of the three S. neurona CK1 enzymes had high sequence similarities to Toxoplasma gondii TgCK1-α and TgCK1-β and the Plasmodium PfCK1. Further experimental studies on the S. neurona putative PKs identified in this study are required to validate the functional roles of the PKs and to understand their involvement in mechanisms that regulate various cellular processes and host-parasite interactions. Given the essentiality of apicomplexan PKs in the survival of apicomplexans, the current study offers a platform for future development of novel therapeutics for EPM, for instance via application of PK inhibitors to block parasite invasion and development in their host.

Classifying kinase conformations using a machine learning approach

Background

Signaling proteins such as protein kinases adopt a diverse array of conformations to respond to regulatory signals in signaling pathways. Perhaps the most fundamental conformational change of a kinase is the transition between active and inactive states, and defining the conformational features associated with kinase activation is critical for selectively targeting abnormally regulated kinases in diseases. While manual examination of crystal structures have led to the identification of key structural features associated with kinase activation, the large number of kinase crystal structures (~3,500) and extensive conformational diversity displayed by the protein kinase superfamily poses unique challenges in fully defining the conformational features associated with kinase activation. Although some computational approaches have been proposed, they are typically based on a small subset of crystal structures using measurements biased towards the active site geometry.

Results

We utilize an unbiased informatics based machine learning approach to classify all eukaryotic protein kinase conformations deposited in the PDB. We show that the orientation of the activation segment, measured by φ, ψ, χ1, and pseudo-dihedral angles more accurately classify kinase crystal conformations than existing methods. We show that the formation of the K-E salt bridge is statistically dependent upon the activation segment orientation and identify evolutionary differences between the activation segment conformation of tyrosine and serine/threonine kinases. We provide evidence that our method can identify conformational changes associated with the binding of allosteric regulatory proteins, and show that the greatest variation in inactive structures comes from kinase group and family specific side chain orientations.

Conclusion

We have provided the first comprehensive machine learning based classification of protein kinase active/inactive conformations, taking into account more structures and measurements than any previous classification effort. Further, our unbiased classification of inactive structures reveals residues associated with kinase functional specificity. To enable classification of new crystal structures, we have made our classifier publicly accessible through a stand-alone program housed at https://github.com/esbg/kinconform [DOI:10.5281/zenodo.249090].

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1506-2) contains supplementary material, which is available to authorized users.

Functional Characterization of Rhoptry Kinome in the Virulent Toxoplasma gondii RH Strain

Toxoplasma gondii is an obligatory intracellular apicomplexan protozoan which can infect any warm-blooded animal and causes severe diseases in immunocompromised individuals or infants infected in utero. The survival and success of this parasite require that it colonizes the host cell, avoids host immune defenses, replicates within an appropriate niche, and exits the infected host cell to spread to neighboring non-infected cells. All of these processes depend on the parasite ability to synthesis and export secreted proteins. Amongst the secreted proteins, rhoptry organelle proteins (ROPs) are essential for the parasite invasion and host cell manipulation. Even though the functions of most ROPs have been elucidated in the less virulent T. gondii (type II), the roles of ROPs in the highly virulent type I strain remain largely un-characterized. Herein, we investigated the contributions of 15 ROPs (ROP10, ROP11, ROP15, ROP20, ROP23, ROP31, ROP32, ROP33, ROP34, ROP35, ROP36, ROP40, ROP41, ROP46, and ROP47) to the infectivity of the high virulent type I T. gondii (RH strain). Using CRISPR-Cas9, these 15 ROPs genes were successfully disrupted and the effects of gene knockout on the parasite's ability to infect cells in vitro and BALB/c mice in vivo were investigated. These results showed that deletions of these ROPs did not interfere with the parasite ability to grow in cultured human foreskin fibroblast cells and did not significantly alter parasite pathogenicity for BALB/c mice. Although these ROPs did not seem to be essential for the acute infectious stage of type I T. gondii in the mouse model, they might have different functions in other intermediate hosts or play different roles in other life cycle forms of this parasite due to the different expression patterns; this warrants further investigations.

Secreted protein kinases regulate cyst burden during chronic toxoplasmosis

Toxoplasma gondii is an apicomplexan parasite that secretes a large number of protein kinases and pseudokinases from its rhoptry organelles. Although some rhoptry kinases (ROPKs) act as virulence factors, many remain uncharacterized. In this study, predicted ROPKs were assessed for bradyzoite expression then prioritized for a reverse genetic analysis in the type II strain Pru that is amenable to targeted disruption. Using CRISPR/Cas9, we engineered C-terminally epitope tagged ROP21 and ROP27 and demonstrated their localization to the parasitophorous vacuole and cyst matrix. ROP21 and ROP27 were not secreted from microneme, rhoptry, or dense granule organelles, but rather were located in small vesicles consistent with a constitutive pathway. Using CRISPR/Cas9, the genes for ROP21, ROP27, ROP28, and ROP30 were deleted individually and in combination, and the mutant parasites were assessed for growth and their ability to form tissue cysts in mice. All knockouts lines were normal for in vitro growth and bradyzoite differentiation, but a combined ∆rop21/∆rop17 knockout led to a 50% reduction in cyst burden in vivo. Our findings question the existing annotation of ROPKs based solely on bioinformatic techniques and yet highlight the importance of secreted kinases in determining the severity of chronic toxoplasmosis.

Functional Analysis of the Rhoptry Kinome during Chronic Toxoplasma gondii Infection

Toxoplasma gondii is one of the most common parasitic infections of humans worldwide. Once exposed, humans remain infected with T. gondii for life, and there are no therapeutics capable of eliminating a chronic infection. In the search for novel drug targets, T. gondii is known to contain several unique secretory organelles, one of which is called the rhoptries. Rhoptry organelles contain and secrete numerous proteins with kinase domains, but the roles of most of these kinases during infection remain unknown. In a recent mBio article, B. A. Fox et al. [mBio 7(3):e00193-16, 2016, http://dx.doi.org/10.1128/mBio.00193-16] performed a tour de force deletion analysis of 31 rhoptry kinases and examined their roles in the development of chronic infection. While rhoptry kinase deletion strains that displayed an acute infection defect also showed a reduction in chronic infection cyst burden, two rhoptry kinase deletion strains had decreased cyst burden without any change in acute virulence. These results indicate the necessity of the rhoptry kinases for the establishment and perhaps maintenance of chronic infection. They also highlight the potential of these kinases as drug targets to clear chronic infection or as candidates to generate a nonpersisting vaccine.

The Toxoplasma gondii Rhoptry Kinome Is Essential for Chronic Infection

Ingestion of the obligate intracellular protozoan parasite Toxoplasma gondii causes an acute infection that leads to chronic infection of the host. To facilitate the acute phase of the infection, T. gondii manipulates the host response by secreting rhoptry organelle proteins (ROPs) into host cells during its invasion. A few key ROP proteins with signatures of kinases or pseudokinases (ROPKs) act as virulence factors that enhance parasite survival against host gamma interferon-stimulated innate immunity. However, the roles of these and other ROPK proteins in establishing chronic infection have not been tested. Here, we deleted 26 ROPK gene loci encoding 31 unique ROPK proteins of type II T. gondii and show that numerous ROPK proteins influence the development of chronic infection. Cyst burdens were increased in the Δrop16 knockout strain or moderately reduced in 11 ROPK knockout strains. In contrast, deletion of ROP5, ROP17, ROP18, ROP35, or ROP38/29/19 (ROP38, ROP29, and ROP19) severely reduced cyst burdens. Δrop5 and Δrop18 knockout strains were less resistant to host immunity-related GTPases (IRGs) and exhibited >100-fold-reduced virulence. ROP18 kinase activity and association with the parasitophorous vacuole membrane were necessary for resistance to host IRGs. The Δrop17 strain exhibited a >12-fold defect in virulence; however, virulence was not affected in the Δrop35 or Δrop38/29/19 strain. Resistance to host IRGs was not affected in the Δrop17, Δrop35, or Δrop38/29/19 strain. Collectively, these findings provide the first definitive evidence that the type II T. gondii ROPK proteome functions as virulence factors and facilitates additional mechanisms of host manipulation that are essential for chronic infection and transmission of T. gondii.

Reactivation of chronic Toxoplasma gondii infection in individuals with weakened immune systems causes severe toxoplasmosis. Existing treatments for toxoplasmosis are complicated by adverse reactions to chemotherapy. Understanding key parasite molecules required for chronic infection provides new insights into potential mechanisms that can interrupt parasite survival or persistence in the host. This study reveals that key secreted rhoptry molecules are used by the parasite to establish chronic infection of the host. Certain rhoptry proteins were found to be critical virulence factors that resist innate immunity, while other rhoptry proteins were found to influence chronic infection without affecting virulence. This study reveals that rhoptry proteins utilize multiple mechanisms of host manipulation to establish chronic infection of the host. Targeted disruption of parasite rhoptry proteins involved in these biological processes opens new avenues to interfere with chronic infection with the goal to either eliminate chronic infection or to prevent recrudescent infections.

Comparative genomics reveals Cyclospora cayetanensis possesses coccidia-like metabolism and invasion components but unique surface antigens

Background

Cyclospora cayetanensis is an apicomplexan that causes diarrhea in humans. The investigation of foodborne outbreaks of cyclosporiasis has been hampered by a lack of genetic data and poor understanding of pathogen biology. In this study we sequenced the genome of C. cayetanensis and inferred its metabolism and invasion components based on comparative genomic analysis.

Results

The genome organization, metabolic capabilities and potential invasion mechanism of C. cayetanensis are very similar to those of Eimeria tenella. Propanoyl-CoA degradation, GPI anchor biosynthesis, and N-glycosylation are some apparent metabolic differences between C. cayetanensis and E. tenella. Unlike Eimeria spp., there are no active LTR-retrotransposons identified in C. cayetanensis. The similar repertoire of host cell invasion-related proteins possessed by all coccidia suggests that C. cayetanensis has an invasion process similar to the one in T. gondii and E. tenella. However, the significant reduction in the number of identifiable rhoptry protein kinases, phosphatases and serine protease inhibitors indicates that monoxenous coccidia, especially C. cayetanensis, have limited capabilities or use a different system to regulate host cell nuclear activities. C. cayetanensis does not possess any cluster of genes encoding the TA4-type SAG surface antigens seen in E. tenella, and may use a different family of surface antigens in initial host cell interactions.

Conclusions

Our findings indicate that C. cayetanensis possesses coccidia-like metabolism and invasion components but unique surface antigens. Amino acid metabolism and post-translation modifications of proteins are some major differences between C. cayetanensis and other apicomplexans. The whole genome sequence data of C. cayetanensis improve our understanding of the biology and evolution of this major foodborne pathogen and facilitate the development of intervention measures and advanced diagnostic tools.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2632-3) contains supplementary material, which is available to authorized users.

Analysis of Noncanonical Calcium-Dependent Protein Kinases in Toxoplasma gondii by Targeted Gene Deletion Using CRISPR/Cas9

Calcium-dependent protein kinases (CDPKs) are expanded in apicomplexan parasites, especially in Toxoplasma gondii where 14 separate genes encoding these enzymes are found. Although previous studies have shown that several CDPKs play a role in controlling invasion, egress, and cell division in T. gondii, the roles of most of these genes are unexplored. Here we developed a more efficient method for gene disruption using CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 (CRISPR-associated protein 9) that was modified to completely delete large, multiexonic genes from the genome and to allow serial replacement by recycling of the selectable marker using Cre-loxP. Using this system, we generated a total of 24 mutants in type 1 and 2 genetic backgrounds to ascertain the functions of noncanonical CDPKs. Remarkably, although we were able to confirm the essentiality of CDPK1 and CDPK7, the majority of CDPKs had no discernible phenotype for growth in vitro or infection in the mouse model. The exception to this was CDPK6, loss of which leads to reduced plaquing, fitness defect in a competition assay, and reduced tissue cyst formation in chronically infected mice. Our findings highlight the utility of CRISPR/Cas9 for rapid serial gene deletion and also suggest that additional models are needed to reveal the functions of many genes in T. gondii.