Genome-Wide Bimolecular Fluorescence Complementation-Based Proteomic Analysis of Toxoplasma gondii ROP18's Human Interactome Shows Its Key Role in Regulation of Cell Immunity and Apoptosis

Proteomics Applications in Toxoplasma gondii: Unveiling the Host-Parasite Interactions and Therapeutic Target Discovery

Toxoplasma gondii, a protozoan parasite with the ability to infect various warm-blooded vertebrates, including humans, is the causative agent of toxoplasmosis. This infection poses significant risks, leading to severe complications in immunocompromised individuals and potentially affecting the fetus through congenital transmission. A comprehensive understanding of the intricate molecular interactions between T. gondii and its host is pivotal for the development of effective therapeutic strategies. This review emphasizes the crucial role of proteomics in T. gondii research, with a specific focus on host-parasite interactions, post-translational modifications (PTMs), PTM crosstalk, and ongoing efforts in drug discovery. Additionally, we provide an overview of recent advancements in proteomics techniques, encompassing interactome sample preparation methods such as BioID (BirA*-mediated proximity-dependent biotin identification), APEX (ascorbate peroxidase-mediated proximity labeling), and Y2H (yeast two hybrid), as well as various proteomics approaches, including single-cell analysis, DIA (data-independent acquisition), targeted, top-down, and plasma proteomics. Furthermore, we discuss bioinformatics and the integration of proteomics with other omics technologies, highlighting its potential in unraveling the intricate mechanisms of T. gondii pathogenesis and identifying novel therapeutic targets.

Improved understanding of biorisk for research involving microbial modification using annotated sequences of concern

Regulation of research on microbes that cause disease in humans has historically been focused on taxonomic lists of 'bad bugs'. However, given our increased knowledge of these pathogens through inexpensive genome sequencing, 5 decades of research in microbial pathogenesis, and the burgeoning capacity of synthetic biologists, the limitations of this approach are apparent. With heightened scientific and public attention focused on biosafety and biosecurity, and an ongoing review by US authorities of dual-use research oversight, this article proposes the incorporation of sequences of concern (SoCs) into the biorisk management regime governing genetic engineering of pathogens. SoCs enable pathogenesis in all microbes infecting hosts that are 'of concern' to human civilization. Here we review the functions of SoCs (FunSoCs) and discuss how they might bring clarity to potentially problematic research outcomes involving infectious agents. We believe that annotation of SoCs with FunSoCs has the potential to improve the likelihood that dual use research of concern is recognized by both scientists and regulators before it occurs.

Preparation and Preliminary Application of Epitope Peptide-Based Antibody against Toxoplasma gondii GRA3

Toxoplasma gondii dense granule protein GRA3 has been shown to promote Toxoplasma gondii transmission and proliferation by interacting with the host cell endoplasmic reticulum (ER) through calcium-regulated cyclophilin ligands (CAMLG). Although many studies have focused on the interaction between the host cell endoplasmic reticulum and GRA3, no polyclonal antibodies (PcAbs) against GRA3 have been reported to date. According to the antigenicity prediction and exposure site analysis, three antigen peptide sequences were selected to prepare polyclonal antibodies targeting GRA3. Peptide scans revealed that the major antigenic epitope sequences were 125ELYDRTDRPGLK136, 202FFRRRPKDGGAG213, and 68NEAGESYSSATSG80, respectively. The GRA3 PcAb specifically recognized the GRA3 of T. gondii type Ⅱ ME49. The development of PcAbs against GRA3 is expected to elucidate the molecular mechanisms by which GRA3 regulates host cell function and contribute to the development of diagnostic and therapeutic strategies for toxoplasmosis.

Overview of Apoptosis, Autophagy, and Inflammatory Processes in Toxoplasma gondii Infected Cells

Toxoplasma gondii (T. gondii) is an obligate intracellular parasite. During the parasitic invasion, T. gondii creates a parasitophorous vacuole, which enables the modulation of cell functions, allowing its replication and host infection. It has effective strategies to escape the immune response and reach privileged immune sites and remain inactive in a controlled environment in tissue cysts. This current review presents the factors that affect host cells and the parasite, as well as changes in the immune system during host cell infection. The secretory organelles of T. gondii (dense granules, micronemes, and rhoptries) are responsible for these processes. They are involved with proteins secreted by micronemes and rhoptries (MIC, AMA, and RONs) that mediate the recognition and entry into host cells. Effector proteins (ROP and GRA) that modify the STAT signal or GTPases in immune cells determine their toxicity. Interference byhost autonomous cells during parasitic infection, gene expression, and production of microbicidal molecules such as reactive oxygen species (ROS) and nitric oxide (NO), result in the regulation of cell death. The high level of complexity in host cell mechanisms prevents cell death in its various pathways. Many of these abilities play an important role in escaping host immune responses, particularly by manipulating the expression of genes involved in apoptosis, necrosis, autophagy, and inflammation. Here we present recent works that define the mechanisms by which T. gondii interacts with these processes in infected host cells.

Host cell proteins modulated upon Toxoplasma infection identified using proteomic approaches: a molecular rationale

Toxoplasma gondii is a pathogenic protozoan parasite belonging to the apicomplexan phylum that infects the nucleated cells of warm-blooded hosts leading to an infectious disease known as toxoplasmosis. Apicomplexan parasites such as T. gondii can display different mechanisms to control or manipulate host cells signaling at different levels altering the host subcellular genome and proteome. Indeed, Toxoplasma is able to modulate host cell responses (especially immune responses) during infection to its advantage through both structural and functional changes in the proteome of different infected cells. Consequently, parasites can transform the invaded cells into a suitable environment for its own replication and the induction of infection. Proteomics as an applicable tool can identify such critical proteins involved in pathogen (Toxoplasma)-host cell interactions and consequently clarify the cellular mechanisms that facilitate the entry of pathogens into host cells, and their replication and transmission, as well as the central mechanisms of host defense against pathogens. Accordingly, the current paper reviews several proteins (identified using proteomic approaches) differentially expressed in the proteome of Toxoplasma-infected host cells (macrophages and human foreskin fibroblasts) and tissues (brain and liver) and highlights their plausible functions in the cellular biology of the infected cells. The identification of such modulated proteins and their related cell impact (cell responses/signaling) can provide further information regarding parasite pathogenesis and biology that might lead to a better understanding of therapeutic strategies and novel drug targets.

An integrated analysis of hypoxic-ischemic encephalopathy-related cell sequencing outcomes via genes network construction

Hypoxic-ischemic encephalopathy (HIE) is one of the main causes of morbidity and severe neurological deficits in neonates. This study aimed to find core genes and their potential roles in HIE with the help of single-cell sequencing (SCS) technology and genes network construction. We collected and screened an HIE genes data set from the Pubmed database to analyze differential expression, and the differential values of genes were ≥3 or ≤-3 in gene expression. We constructed a protein-protein interaction (PPI) network by the string, which was also verified by Cytoscape 3.8.2. Functional enrichment analysis was performed to determine the characteristics and pathways of the core genes. We examined two meaningful papers and integrated all genes by SCS, which were classified into 12,093 genes without duplicates, 217 shared genes, and 11,876 distinct genes. Among 217 genes, the signal transducer and activator of transcription (STAT) family was the most targeted gene in the PPI network. Moreover, Gene Ontology and Kyoto encyclopedia of genes and genome analysis showed that the process in response to virus and the JAK-STAT signaling pathway play significant roles in HIE. We also found that 54 screened genes were highly expressed, while three genes (B2M, VIM, and MRPS30) were different in the heat map and differential genes expression exhibition. VIM, as an essential portion of the brain's cytoskeleton, is closely linked to STAT and neurologic development. From the findings of SCS and bioinformatics predictive analytics model, our outcomes provided a better understanding of the roles of STAT, the JAK-STAT signaling pathway, and VIM, which can pave an alternative avenue for further studies on HIE progression.

Toxoplasma gondii Extends the Life Span of Infected Human Neutrophils by Inducing Cytosolic PCNA and Blocking Activation of Apoptotic Caspases

Toxoplasma gondii is an intracellular protozoan parasite that has the remarkable ability to infect and replicate in neutrophils, immune cells with an arsenal of antimicrobial effector mechanisms. We report that T. gondii infection extends the life span of primary human peripheral blood neutrophils by delaying spontaneous apoptosis, serum starvation-induced apoptosis, and tumor necrosis alpha (TNF-α)-mediated apoptosis. T. gondii blockade of apoptosis was associated with an inhibition of processing and activation of the apoptotic caspases caspase-8 and -3, decreased phosphatidylserine exposure on the plasma membrane, and reduced cell death. We performed a global transcriptome analysis of T. gondii-infected peripheral blood neutrophils using RNA sequencing (RNA-Seq) and identified gene expression changes associated with DNA replication and DNA repair pathways, which in mature neutrophils are indicative of changes in regulators of cell survival. Consistent with the RNA-Seq data, T. gondii infection upregulated transcript and protein expression of PCNA, which is found in the cytosol of human neutrophils, where it functions as a key inhibitor of apoptotic pro-caspases. Infection of neutrophils resulted in increased interaction of PCNA with pro-caspase-3. Inhibition of this interaction with an AlkB homologue 2 PCNA-interacting motif (APIM) peptide reversed the infection-induced delay in cell death. Taken together, these findings indicate a novel strategy by which T. gondii manipulates cell life span in primary human neutrophils, which may allow the parasite to maintain an intracellular replicative niche and avoid immune clearance.IMPORTANCEToxoplasma gondii is an obligate intracellular parasite that can cause life-threatening disease in immunocompromised individuals and in the developing fetus. Interestingly, T. gondii has evolved strategies to successfully manipulate the host immune system to establish a productive infection and evade host defense mechanisms. Although it is well documented that neutrophils are mobilized during acute T. gondii infection and infiltrate the site of infection, these cells can also be actively infected by T. gondii and serve as a replicative niche for the parasite. However, there is a limited understanding of the molecular processes occurring within T. gondii-infected neutrophils. This study reveals that T. gondii extends the life span of human neutrophils by inducing the expression of PCNA, which prevents activation of apoptotic caspases, thus delaying apoptosis. This strategy may allow the parasite to preserve its replicative intracellular niche.

Bimolecular Fluorescence Complementation (BiFC) and Multiplexed Imaging of Protein-Protein Interactions in Human Living Cells

Deciphering protein-protein interactions (PPIs) in vivo is crucial to understand protein function. Bimolecular fluorescence complementation (BiFC) makes applicable the analysis of PPIs in many different native contexts, including human live cells. It relies on the property of monomeric fluorescent proteins to be reconstituted from two separate subfragments upon spatial proximity. Candidate partners fused to such complementary subfragments can form a fluorescent protein complex upon interaction, allowing visualization of weak and transient PPIs. It can also be applied for investigation of distinct PPIs at the same time using a multicolor setup. In this chapter, we provide a detailed protocol for analyzing PPIs by doing BiFC in cultured cells. Proof-of-principle experiments rely on the complementation property between the N-terminal fragment of mVenus (designated VN173) and the C-terminal fragment of mCerulean (designated CC155) and the partnership between HOXA7 and PBX1 proteins. This protocol is compatible with any other fluorescent complementation pair fragments and any type of candidate interacting proteins.

ROP18-Mediated Transcriptional Reprogramming of HEK293T Cell Reveals New Roles of ROP18 in the Interplay Between Toxoplasma gondii and the Host Cell

Toxoplasma gondii secretes a number of virulence-related effector proteins, such as the rhoptry protein 18 (ROP18). To further broaden our understanding of the molecular functions of ROP18, we examined the transcriptional response of human embryonic kidney cells (HEK293T) to ROP18 of type I T. gondii RH strain. Using RNA-sequencing, we compared the transcriptome of ROP18-expressing HEK293T cells to control HEK293T cells. Our analysis revealed that ROP18 altered the expression of 750 genes (467 upregulated genes and 283 downregulated genes) in HEK293T cells. Gene ontology (GO) and pathway enrichment analyses showed that differentially expressed genes (DEGs) were significantly enriched in extracellular matrix– and immune–related GO terms and pathways. KEGG pathway enrichment analysis revealed that DEGs were involved in several disease-related pathways, such as nervous system diseases and eye disease. ROP18 significantly increased the alternative splicing pattern “retained intron” and altered the expression of 144 transcription factors (TFs). These results provide new insight into how ROP18 may influence biological processes in the host cells via altering the expression of genes, TFs, and pathways. More in vitro and in vivo studies are required to substantiate these findings.

TgROP18 targets IL20RB for host-defense-related-STAT3 activation during Toxoplasma gondii infection

Background

Toxoplasma gondii is an opportunistic protozoan infecting almost one-third of the world’s population. Toxoplasma gondii rhoptry protein 18 (TgROP18) is a key virulence factor determining the parasite’s acute virulence and is secreted into host cells during infection. We previously identified the interaction of TgROP18 and host cell immune-related receptor protein IL20RB, and observed the activation of STAT3 in human keratinocytes (HaCaT) cells infected by the rop16 knockout RH strain, though TgROP16 is regarded as being responsible for host STAT3 activation during T. gondii invasion. Therefore, we hypothesize TgROP18 can activate host STAT3 through binding to IL20RB.

Methods

CRISPR-CAS9 technology was used to generate the ROP16 and ROP18 double knockout RH strain, RH-∆rop16∆rop18. SDS-PAGE and western blot were used to detect STAT3 activation in different HaCaT cells with high endogenous IL20RB expression treated with T. gondii tachyzoites infection, recombinant ROP18, or IL-20. FRET and co-immunoprecipitation (Co-IP) was used to detect the protein-protein interaction.

Results

We observed that TgROP18 was involved in a synergic activation of the host JAK/STAT3 pathway together with TgROP16 in human HaCaT cells infected with T. gondii or treated with recombinant TgROP18 protein, stimulating host proinflammatory immune responses such as expression of TNF-α. The effect of recombinant ROP18 on STAT3 phosphorylation was presented in a dose-dependent manner. Additionally, TgROP18 was identified to target IL20RB on its extracellular domain. When we treated different cell lines with the recombinant ROP18, STAT3 phosphorylation could only be observed in the cells with endogenous IL20RB expression, such as HaCaT cells.

Conclusions

These findings indicate that TgROP18-IL20RB interaction upon T. gondii invasion was involved in STAT3 activation, which is associated with host cell defense.

Strain-specific disruption of interferon-stimulated N-myc and STAT interactor (NMI) function by Toxoplasma gondii type I ROP18 in human cells

Toxoplasma gondii rhoptry protein TgROP18 is a polymorphic virulence effector that targets immunity-related GTPases (IRGs) in rodents. Given that IRGs are uniquely diversified in rodents and not in other T. gondii intermediate hosts, the role of TgROP18 in manipulating non-rodent cells is unclear. Here we show that in human cells TgROP18I interacts with the interferon-gamma-inducible protein N-myc and STAT interactor (NMI) and that this is a property that is unique to the type I TgROP18 allele. Specifically, when expressed ectopically in mammalian cells only TgROP18I co-immunoprecipitates with NMI in IFN-γ-treated cells, while TgROP18II does not. In parasites expressing TgROP18I or TgROP18II, NMI only co-immunoprecipitates with TgROP18I and this is associated with allele-specific immunolocalization of NMI on the parasitophorous vacuolar membrane (PVM). We also found that TgROP18I reduces NMI association with IFN-γ-activated sequences (GAS) in the IRF1 gene promoter. Finally, we determined that polymorphisms in the C-terminal kinase domain of TgROP18I are required for allele-specific effects on NMI. Together, these data further define new host pathway targeted by TgROP18I and provide the first function driven by allelic differences in the highly polymorphic ROP18 locus.

Production and characterization of monoclonal antibodies against Toxoplasma gondii ROP18 with strain-specific reactivity

The rhoptry kinase 18 of Toxoplasma gondii (TgROP18) has been identified as a key virulence factor that allows the parasite to escape from host immune defences and promotes its proliferation in host cells. Although much research is focused on the interaction between host cells and TgROP18, the development of monoclonal antibodies (mAbs) against TgROP18 has not been reported till date. To produce mAbs targeting TgROP18, two hybridomas secreting mAbs against TgROP18, designated as A1 and T2, were generated using cell fusion technology. The subtypes of the A1 and T2 mAbs were identified as IgG3 λ and IgM κ, and peptide scanning revealed that the core sequences of the antigenic epitopes were 180LRAQRRRSELVFE192 and 351NYFLLMMRAEADM363, respectively. The T2 mAb specifically reacted with both T. gondii type I and Chinese I, but not with T. gondii type II, Plasmodium falciparum or Schistosoma japonicum. Finally, the sequences of heavy chain and light chain complementarity-determining regions of T2 were amplified, cloned and characterized, making the modification of the mAb feasible in the future. The development of mAbs against TgROP18 would aid the investigation of the molecular mechanisms underlying the modulation of host cellular functions by TgROP18, and in the development of strategies to diagnose and treat Toxoplasmosis.

Identification of key modules and prognostic markers in adrenocortical carcinoma by weighted gene co-expression network analysis

Adrenocortical carcinoma (ACC) is a rare and aggressive cancer with a high relapse rate and limited treatment options. Therefore, the identification of potential prognostic markers in patients with ACC may improve early detection, survival rates and may additionally provide novel insights into the early detection of recurrence. In the present study, clinical traits and RNA-seq data of 79 patients with ACC were obtained from The Cancer Genome Atlas (TCGA). Weighted gene co-expression network analysis was carried out and 17 distinct co-expression modules were built to examine the association between the modules and the clinical traits. Of the 17 modules, two co-expression modules, which contained 214 and 168 genes, were significantly correlated with two clinical traits, tumor stage and vital status. Functional enrichment analysis was performed on the selected modules. The results showed that one of the modules was primarily enriched in cell division and the other module was enriched in metabolic pathways, suggesting their involvement in tumor progression. Furthermore, cyclin dependent kinase 1 (CDK1) and ubiquitin C (UBC) were identified as hub genes in both modules. Survival analysis revealed that the high expression of the hub genes significantly correlated with the poor survival rate of patients, suggesting that CDK1 and UBC have vital roles in the progression of ACC. In the present study, a co-expression gene module of ACC was provided and the prognostic genes that may serve as new diagnostic markers in the future were defined.

Toxoplasma gondii ROP18 inhibits human glioblastoma cell apoptosis through a mitochondrial pathway by targeting host cell P2X1

Background

Apoptosis plays a critical role in the embryonic development, homeostasis of immune system and host defense against intracellular microbial pathogens. Infection by the obligate intracellular pathogen Toxoplasma gondii can both inhibit and induce host cell apoptosis; however, the parasitic factors involved remain unclear. The T. gondii virulence factor ROP18 (TgROP18) has been reported to regulate host cell apoptosis; nevertheless, results for this regulation have been rarely reported or have provided contradictory findings. Human purinergic receptor 1 (P2X1) is an ATP-gated ion channel that responds to ATP stimulation and functions in cell apoptosis mediation. The precise roles of TgROP18 in T. gondii pathogenesis, and the relationship between TgROP18 and host P2X1 in host cell apoptosis are yet to be revealed.

Methods

Apoptosis rates were determined by flow cytometry (FCM) and TUNEL assay. The interaction between TgROP18 and the host P2X1 was measured by fluorescence resonance energy transfer (FRET) and co-immunoprecipitation (co-IP) assay. Calcium influx and mitochondrial membrane depolarization were determined by FCM after JC-1 staining. The translocation of cytochrome C (Cyt C), Bax and Bcl2 proteins, expression of the apoptotic proteins PARP and caspase activation were detected by western blotting.

Results

The apoptosis rates of glial or immune cells (human SF268, mouse RAW264.7 and human THP-1 cells) infected by any T. gondii strain (RH-type I, ME49-type II and VEG-type III) were significantly inhibited compared with their uninfected controls. TgROP18 inhibited ATP-induced apoptosis of SF268 with P2X1 expression, but had no effect on RAW264.7 or THP-1 cells without detectable P2X1 expression. It was further identified that TgROP18 interacted with P2X1, and overexpression of ROP18 in COS7 cells significantly inhibited cell apoptosis mediated by P2X1. Moreover, TgROP18 also inhibited P2X1-mediated Ca²⁺ influx, translocation of cytochrome C from the mitochondria to the cytosol, and ATP-triggered caspase activation.

Conclusions

Toxoplasma gondii infection inhibits ATP-induced host cell apoptosis, regardless of strain virulence and host cell lines. TgROP18 targets the purinergic receptor P2X1 of the SF268 human neural cells and inhibits ATP-induced apoptosis through the mitochondrial pathway, suggesting a sensor role for the host proapoptotic protein P2X1 in this process.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3529-1) contains supplementary material, which is available to authorized users.

Strategies Developed by Toxoplasma gondii to Survive in the Host

One of the most successful intracellular parasites, Toxoplasma gondii has developed several strategies to avoid destruction by the host. These include approaches such as rapid and efficient cell invasion to avoid phagocytic engulfment, negative regulation of the canonical CD40-CD40L-mediated autophagy pathway, impairment of the noncanonical IFN-γ-dependent autophagy pathway, and modulation of host cell survival and death to obtain lifelong parasite survival. Different virulent strains have even evolved different ways to cope with and evade destruction by the host. This review aims to illustrate every aspect of the game between the host and Toxoplasma during the process of infection. A better understanding of all aspects of the battle between Toxoplasma and its hosts will be useful for the development of better strategies and drugs to control the parasite.

Prediction of Toxoplasma gondii virulence factor ROP18 competitive inhibitors by virtual screening

Background

Rhoptry protein 18 (ROP18) is a key virulence factor of Toxoplasma gondii. The host’s immune responses mediated by immune-related GTPases (IRGs) could be blocked by ROP18’s kinase activity. ROP18 also interacts with various substrates, such as activating transcription factor 6 beta (ATF6β) and affects multiple physiological functions within host cells, thereby inducing intense virulence. In this study, competitive inhibitors targeted to ROP18 were subjected to virtual screening based on the principle of structure-based drug design (SBDD).

Methods

The preparation of the ROP18 structure was conducted using the “Structure Prepare” function of Molecular Operating Environment (MOE) software. The ATP-binding pocket was selected as the starting point for virtual screening. Construction of the pharmacophore model used Extended Hückel Theory (EHT) half-quantitative measurement and construction, as well as the characteristics of Type I kinase inhibitors. The pharmacophore model of ROP18 was imported into the Specs database for small molecule similarity screening using EHT pharmacophore measurement. Hit compounds were selected using the functions of London dG and generalized-born volume integral/weighted surface area (GBVI/WSA) scoring. The top 100 hits were analyzed by molecular docking and structure activity relationships (SAR) analysis.

Results

The final pharmacophore comprised three typical characteristics: three hydrogen bond acceptors/donors, two ring aromatic features occupying the hydrophobic core, and one cation group feature targeted to the terminus of ATP. A total of 1314 hit compounds analogous to ROP18 pharmacophore were passed through the Specs. After two rounds of docking, 25 out of 100 hits were identified as belonging to two main scaffold types: phthalimide ring structure, thiazole ring and styrene structure. Additionally, the screen also identified 13 inhibitors with distinct scaffold types. The docking models and SAR analysis demonstrated that these hits could engage in multiple hydrogen bonds, salt bridges halogen bonds, and hydrophobic interactions with ROP18, and para-position halo substituents on the benzene ring may enhance their affinity scoring.

Conclusions

A pharmacophore against the ROP18 ATP-binding pocket was successfully constructed, and 25 representative inhibitors from 15 scaffold clusters were screened using the Specs database. Our results provide useful scaffold types for the chemical inhibition of ROP18 or alternative conformations to develop new anti-toxoplasmosis drug leads.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3341-y) contains supplementary material, which is available to authorized users.

Protein targets of thiazolidinone derivatives in Toxoplasma gondii and insights into their binding to ROP18

Background

Thiazolidinone derivatives show inhibitory activity (IC₅₀) against the Toxoplasma gondii parasite, as well as high selectivity with high therapeutic index. To disclose the target proteins of the thiazolidinone core in this parasite, we explored in silico the active sites of different T. gondii proteins and estimated the binding-free energy of reported thiazolidinone molecules with inhibitory effect on invasion and replication of the parasite inside host cells. This enabled us to describe some of the most suitable structural characteristics to design a compound derived from the thiazolidinone core.

Results

The best binding affinity was observed in the active site of kinase proteins, we selected the active site of the T. gondii ROP18 kinase, because it is an important factor for the virulence and survival of the parasite. We present the possible effect of a derivative of thiazolidinone core in the active site of T. gondii ROP18 and described some characteristics of substituent groups that could improve the affinity and specificity of compounds derived from the thiazolidinone core against T. gondii.

Conclusions

The results of our study suggest that compounds derived from the thiazolidinone core have a preference for protein kinases of T. gondii, being promising compounds for the development of new drugs with potential anti-toxoplasmosis activity. Our findings highlight the importance of use computational studies for the understanding of the action mechanism of compounds with biological activity.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5223-7) contains supplementary material, which is available to authorized users.