Systemic toxicity of snake venom metalloproteinases: Multi-omics analyses of kidney and blood plasma disturbances in a mouse model

Snakebite envenomation is classified as a Neglected Tropical Disease. Bothrops jararaca venom induces kidney injury and coagulopathy. HF3, a hemorrhagic metalloproteinase of B. jararaca venom, participates in the envenomation pathogenesis. We evaluated the effects of HF3 in mouse kidney and blood plasma after injection in the thigh muscle, mimicking a snakebite. Transcriptomic analysis showed differential expression of 31 and 137 genes related to kidney pathology after 2 h and 6 h, respectively. However, only subtle changes were observed in kidney proteome, with differential abundance of 15 proteins after 6 h, including kidney injury markers. N-terminomic analysis of kidney proteins showed 420 proteinase-generated peptides compatible with meprin specificity, indicating activation of host proteinases. Plasma analysis revealed differential abundance of 90 and 219 proteins, respectively, after 2 h and 6 h, including coagulation-cascade and complement-system components, and creatine-kinase, whereas a semi-specific search of N-terminal peptides indicated activation of endogenous proteinases. HF3 promoted host reactions, altering the gene expression and the proteolytic profile of kidney tissue, and inducing plasma proteome imbalance driven by changes in abundance and proteolysis. The overall response of the mouse underscores the systemic action of a hemorrhagic toxin that transcends local tissue damage and is related to known venom-induced systemic effects.

A snapshot of Bothrops jararaca snake venom gland subcellular proteome

Snake venom protein synthesis undergoes finely regulated processes in the specialized secretory epithelium within the venom gland. Such processes occur within a defined period in the cell and at specific cellular locations. Thus, the determination of subcellular proteomes allows the characterization of protein groups for which the site may be relevant to their biological roles, thereby allowing the deconvolution of complex biological circuits into functional information. In this regard, we performed subcellular fractionation of proteins from B. jararaca venom gland, focusing on nuclear proteins since this cellular compartment comprises key effectors that shape gene expression. Our results provided a snapshot of B. jararaca's subcellular venom gland proteome and pointed to a 'conserved' proteome core among different life stages (newborn and adult) and between sexes (adult male and female). Overall, the top 15 highly abundant proteins identified in B. jararaca venom glands mirrored the panel of highly expressed genes in human salivary glands. Therefore, the expression profile observed for such a protein set could be considered a conserved core signature of salivary gland secretory epithelium. Moreover, the newborn venom gland displayed a unique expression signature of transcription factors involved in regulating transcription and biosynthetic processes and may mirror biological constraints of the ontogenetic development of B. jararaca, contributing to venom proteome diversity.

Profilings of subproteomes of lectin-binding proteins of nine Bothrops venoms reveal variability driven by different glycan types

Snake venom proteomes have long been investigated to explore a multitude of biologically active components that are used for prey capture and defense, and are involved in the pathological effects observed upon mammalian envenomation. Glycosylation is a major protein post-translational modification in venoms and contributes to the diversification of proteomes. We have shown that Bothrops venoms are markedly defined by their content of glycoproteins, and that most N-glycan structures of eight Bothrops venoms contain sialic acid, while bisected N-acetylglucosamine was identified in Bothrops cotiara venom. To further investigate the mechanisms involved in the generation of different venoms by related snakes, here the glycoproteomes of nine Bothrops venoms (Bothrops atrox, B. cotiara, Bothrops erythromelas, Bothrops fonsecai, B. insularis, Bothrops jararaca, Bothrops jararacussu, Bothrops moojeni and Bothrops neuwiedi) were comparatively analyzed by enrichment with three lectins of different specificities, recognizing bisecting N-acetylglucosamine- and sialic acid-containing glycoproteins, and mass spectrometry. The lectin capture strategy generated venom fractions enriched with several glycoproteins, including metalloprotease, serine protease, and L- amino acid oxidase, in addition to various types of low abundant enzymes. The different contents of lectin-enriched proteins underscore novel aspects of the variability of the glycoprotein subproteomes of Bothrops venoms and point to the role of distinct types of glycan chains in generating different venoms by closely related snake species.

Co-Exposure of Cardiomyocytes to IFN-γ and TNF-α Induces Mitochondrial Dysfunction and Nitro-Oxidative Stress: Implications for the Pathogenesis of Chronic Chagas Disease Cardiomyopathy

Infection by the protozoan Trypanosoma cruzi causes Chagas disease cardiomyopathy (CCC) and can lead to arrhythmia, heart failure and death. Chagas disease affects 8 million people worldwide, and chronic production of the cytokines IFN-γ and TNF-α by T cells together with mitochondrial dysfunction are important players for the poor prognosis of the disease. Mitochondria occupy 40% of the cardiomyocytes volume and produce 95% of cellular ATP that sustain the life-long cycles of heart contraction. As IFN-γ and TNF-α have been described to affect mitochondrial function, we hypothesized that IFN-γ and TNF-α are involved in the myocardial mitochondrial dysfunction observed in CCC patients. In this study, we quantified markers of mitochondrial dysfunction and nitro-oxidative stress in CCC heart tissue and in IFN-γ/TNF-α-stimulated AC-16 human cardiomyocytes. We found that CCC myocardium displayed increased levels of nitro-oxidative stress and reduced mitochondrial DNA as compared with myocardial tissue from patients with dilated cardiomyopathy (DCM). IFN-γ/TNF-α treatment of AC-16 cardiomyocytes induced increased nitro-oxidative stress and decreased the mitochondrial membrane potential (ΔΨm). We found that the STAT1/NF-κB/NOS2 axis is involved in the IFN-γ/TNF-α-induced decrease of ΔΨm in AC-16 cardiomyocytes. Furthermore, treatment with mitochondria-sparing agonists of AMPK, NRF2 and SIRT1 rescues ΔΨm in IFN-γ/TNF-α-stimulated cells. Proteomic and gene expression analyses revealed that IFN-γ/TNF-α-treated cells corroborate mitochondrial dysfunction, transmembrane potential of mitochondria, altered fatty acid metabolism and cardiac necrosis/cell death. Functional assays conducted on Seahorse respirometer showed that cytokine-stimulated cells display decreased glycolytic and mitochondrial ATP production, dependency of fatty acid oxidation as well as increased proton leak and non-mitochondrial oxygen consumption. Together, our results suggest that IFN-γ and TNF-α cause direct damage to cardiomyocytes’ mitochondria by promoting oxidative and nitrosative stress and impairing energy production pathways. We hypothesize that treatment with agonists of AMPK, NRF2 and SIRT1 might be an approach to ameliorate the progression of Chagas disease cardiomyopathy.

Cell Cycle, Telomeres, and Telomerase in Leishmania spp.: What Do We Know So Far?

Leishmaniases belong to the inglorious group of neglected tropical diseases, presenting different degrees of manifestations severity. It is caused by the transmission of more than 20 species of parasites of the Leishmania genus. Nevertheless, the disease remains on the priority list for developing new treatments, since it affects millions in a vast geographical area, especially low-income people. Molecular biology studies are pioneers in parasitic research with the aim of discovering potential targets for drug development. Among them are the telomeres, DNA–protein structures that play an important role in the long term in cell cycle/survival. Telomeres are the physical ends of eukaryotic chromosomes. Due to their multiple interactions with different proteins that confer a likewise complex dynamic, they have emerged as objects of interest in many medical studies, including studies on leishmaniases. This review aims to gather information and elucidate what we know about the phenomena behind Leishmania spp. telomere maintenance and how it impacts the parasite’s cell cycle.

The distinct N-terminomes of Bothrops jararaca newborn and adult venoms

Using approaches of transcriptomics and proteomics we have shown that the phenotype of Bothrops jararaca venom undergoes a significant rearrangement upon neonate to adult transition. Most regulatory processes in biology are intrinsically related to modifications of protein structure, function, and abundance. However, it is unclear to which extent intrinsic proteolysis affects toxins and snake venom phenotypes upon ontogenesis. Here we assessed the natural N-terminome of Bothrops jararaca newborn and adult venoms and explored the degree of N-terminal protein truncation in ontogenetic-based proteome variation. To this end we applied the Terminal Amine Isotopic Labeling of Substrates (TAILS) technology to characterize venom collected in the presence of proteinase inhibitors. We identified natural N-terminal sequences in the newborn (71) and adult (84) venoms, from which only 37 were common to both. However, truncated toxins were found in higher number in the newborn (212) than in the adult (140) venom. Moreover, sequences N-terminally blocked by pyroglutamic acid were identified in the newborn (55) and adult (49) venoms. Most toxin classes identified by their natural N-terminal sequences showed a similar number of unique peptides in the newborn and adult venoms, however, those of serine proteinases and C-type lectins were more abundant in the adult venom. Truncated sequences from at least ten toxin classes were detected, however the catalytic and cysteine-rich domains of metalloproteinases were the most prone to proteolysis, mainly in the newborn venom. Our results underscore the pervasiveness of truncations in most toxin classes and highlight variable post-translational events in newborn and adult venoms.

Venom Profiling of the Insular Species Bothrops alcatraz: Characterization of Proteome, Glycoproteome, and N-Terminome Using Terminal Amine Isotopic Labeling of Substrates

Bothrops alcatraz, a species endemic to Alcatrazes Islands, is regarded as critically endangered due to its small area of occurrence and the declining quality of its habitat. We recently reported the identification of N-glycans attached to toxins of Bothrops species, showing similar compositions in venoms of the B. jararaca complex (B. jararaca, B. insularis, and B. alcatraz). Here, we characterized B. alcatraz venom using electrophoretic, proteomic, and glycoproteomic approaches. Electrophoresis showed that B. alcatraz venom differs from B. jararaca and B. insularis; however, N-glycan removal revealed similarities between them, indicating that the occupation of N-glycosylation sites contributes to interspecies variability in the B. jararaca complex. Metalloproteinase was the major toxin class identified in the B. alcatraz venom proteome followed by serine proteinase and C-type lectin, and overall, the adult B. alcatraz venom resembles that of B. jararaca juvenile specimens. The comparative glycoproteomic analysis of B. alcatraz venom with B. jararaca and B. insularis indicated that there may be differences in the utilization of N-glycosylation motifs among their different toxin classes. Furthermore, we prospected for the first time the N-terminome of a snake venom using the terminal amine isotopic labeling of substrates (TAILS) approach and report the presence of ∼30% of N-termini corresponding to truncated toxin forms and ∼37% N-terminal sequences blocked by pyroglutamic acid in B. alcatraz venom. These findings underscore a low correlation between venom gland transcriptomes and proteomes and support the view that post-translational processes play a major role in shaping venom phenotypes.

Community-based network analyses reveal emerging connectivity patterns of protein-protein interactions in murine melanoma secretome

Protein-protein interaction networks (PPINs) are static representations of protein connections in which topological features such as subgraphs (communities) may contain proteins functionally related, revealing an additional layer of interactome complexity. We created two PPINs from the secretomes of a paired set of murine melanocytes (a normal melanocyte and its transformed phenotype). Community structures, identified by a graph clustering algorithm, resulted in the identification of subgraphs in both networks. Interestingly, the underlying structure of such communities revealed shared and exclusive proteins (core and exclusive nodes, respectively), in addition to proteins that changed their location within each community (rewired nodes). Functional enrichment analysis of core nodes revealed conserved biological functions in both networks whereas exclusive and rewired nodes in the tumoral phenotype network were enriched in cancer-related processes, including TGFβ signaling. We found a remarkable shift in the tumoral interactome, resulting in an emerging pattern which was driven by the presence of exclusive nodes and may represent functional network motifs. Our findings suggest that the rearrangement in the tumoral interactome may be correlated with the malignant transformation of melanocytes associated with substrate adhesion impediment. The interactions found in core and new/rewired nodes might potentially be targeted for therapeutic intervention in melanoma treatment. SIGNIFICANCE: Malignant transformation is a result of synergistic action of multiple molecular factors in which genetic alterations as well as protein expression play paramount roles. During oncogenesis, cellular crosstalk through the secretion of soluble mediators modulates the phenotype of transformed cells which ultimately enables them to successfully disrupt important signaling pathways, including those related to cell growth and proliferation. Therefore, in this work we profiled the secretomes of a paired set of normal and transformed phenotypes of a murine melanocyte. After assembling the two interactomes, clusters of functionally related proteins (network communities) were observed as well as emerging patterns of network rewiring which may represent an interactome signature of transformed cells. In summary, the significance of this study relies on the understanding of the repertoire of 'normal' and 'tumoral' secretomes and, more importantly, the set of interacting proteins (the interactome) in both of these conditions, which may reveal key components that might be potentially targeted for therapeutic intervention.

Heterotypic signaling between dermal fibroblasts and melanoma cells induces phenotypic plasticity and proteome rearrangement in malignant cells

The signaling events triggered by soluble mediators released from both transformed and stromal cells shape the phenotype of tumoral cells and have significant implications in cancer development and progression. In this study we performed an in vitro heterotypic signaling assays by evaluating the proteome diversity of human dermal fibroblasts after stimulation with the conditioned media obtained from malignant melanoma cells. In addition, we also evaluated the changes in the proteome of melanoma cells after stimulation with their own conditioned media as well as with the conditioned medium from melanoma-stimulated fibroblasts. Our results revealed a clear rearrangement in the proteome of stromal and malignant cells upon crosstalk of soluble mediators. The main proteome signature of fibroblasts stimulated with melanoma conditioned medium was related to protein synthesis, which indicates that this process might be an early response of stromal cells. In addition, the conditioned medium derived from 'primed' stromal cells (melanoma-stimulated fibroblasts) was more effective in altering the functional phenotype (cell migration) of malignant cells than the conditioned medium from non-stimulated fibroblasts. Collectively, self- and cross-stimulation may play a key role in shaping the tumor microenvironment and enable tumoral cells to succeed in the process of melanoma progression and metastasis. Although the proteome landscape of cells participating in such a heterotypic signaling represents a snapshot of a highly dynamic state, understanding the diversity of proteins and enriched biological pathways resulting from stimulated cell states may allow for targeting specific cell regulatory motifs involved in melanoma progression and metastasis.

Biochemical characterization of venom from Pseudoboa neuwiedii (Neuwied's false boa; Xenodontinae; Pseudoboini)

In this work, we examined the proteolytic and phospholipase A₂ (PLA₂) activities of venom from the opisthoglyphous colubrid Pseudoboa neuwiedii. Proteolytic activity (3 and 10 μg of venom) was comparable to that of Bothrops neuwiedii venom but less than Bothrops atrox. This activity was inhibited by EDTA and 1,10-phenanthroline but only slightly affected (≤30% inhibition) by PMSF and AEBSF, indicating it was mediated by snake venom metalloproteinases (SVMPs). The pH and temperature optima for proteolytic activity were 8.0 and 37 °C, respectively. The venom had no esterase activity, whereas PLA₂ activity was similar to B. atrox, greater than B. neuwiedii but less than B. jararacussu. SDS-PAGE revealed venom proteins >100 kDa, 45-70 kDa, 21-24 kDa and ~15 kDa, and mass spectrometry of protein bands revealed SVMPs, cysteine-rich secretory proteins (CRISPs) and PLA₂, but no serine proteinases. In gelatin zymography, the most active bands occurred at 65-68 kDa (seen with 0.05-0.25 μg of venom). Caseinolytic activity occurred at 50-66 kDa and was generally weaker than gelatinolytic activity. RP-HPLC of venom yielded 15 peaks, five of which showed gelatinolytic activity; peak 7 was the most active and apparently contained a P-III class SVMP. The venom showed α-fibrinogenase activity, without affecting the β and γ chains; this activity was inhibited by EDTA and 1,10-phenanthroline. The venom did not clot rat citrated plasma but reduced the rate and extent of coagulation after plasma recalcification. In conclusion, P. neuwiedii venom is highly proteolytic and could potentially affect coagulation in vivo by degrading fibrinogen via SVMPs.

Structures of N-Glycans of Bothrops Venoms Revealed as Molecular Signatures that Contribute to Venom Phenotype in Viperid Snakes

The complexity of snake venoms has long been investigated to explore a myriad of biologically active proteins and peptides that are used for immobilizing or killing prey, and are responsible for the pathological effects observed on envenomation. Glycosylation is the main post-translational modification (PTM) of viperid venoms but currently there is little understanding of how protein glycosylation impacts the variation of venom proteomes. We have previously reported that Bothrops venom glycoproteomes contain a core of components that markedly define their composition and parallel their phylogenetic classification. Here we extend those observations to eight Bothrops species evaluating the N-glycomes by LC-MS as assigned cartoon structures and detailing those structures separately as methylated analogs using ion-trap mass spectrometry (MSⁿ). Following ion disassembly through multiple steps provided sequence and linkage isomeric details that characterized 52 unique compositions in Bothrops venoms. These occurred as 60 structures, of which 26 were identified in the venoms of the Jararaca Complex (B. alcatraz, B. insularis, and B. jararaca), 20 in B. erythromelas, B. jararacussu, B. moojeni and B. neuwiedi venoms, and 22 in B. cotiara venom. Further, quantitative analysis of these N-glycans showed variable relative abundances in the venoms. For the first time a comprehensive set of N-glycan structures present in snake venoms are defined. Despite the fact that glycosylation is not template-defined, the N-glycomes of these venoms mirror the phylogeny cladograms of South American bothropoid snakes reported in studies on morphological, molecular data and feeding habits, exhibiting distinct molecular signatures for each venom. Considering the complexity of N-glycan moieties generally found in glycoproteins, characterized by different degrees of branching, isomer structures, and variable abundances, our findings point to these factors as another level of complexity in Bothrops venoms, features that could dramatically contribute to their distinct biological activities.

Trends in the Evolution of Snake Toxins Underscored by an Integrative Omics Approach to Profile the Venom of the Colubrid Phalotris mertensi

Only few studies on snake venoms were dedicated to deeply characterize the toxin secretion of animals from the Colubridae family, despite the fact that they represent the majority of snake diversity. As a consequence, some evolutionary trends observed in venom proteins that underpinned the evolutionary histories of snake toxins were based on data from a minor parcel of the clade. Here, we investigated the proteins of the totally unknown venom from Phalotris mertensi (Dipsadinae subfamily), in order to obtain a detailed profile of its toxins and to appreciate evolutionary tendencies occurring in colubrid venoms. By means of integrated omics and functional approaches, including RNAseq, Sanger sequencing, high-resolution proteomics, recombinant protein production, and enzymatic tests, we verified an active toxic secretion containing up to 21 types of proteins. A high content of Kunitz-type proteins and C-type lectins were observed, although several enzymatic components such as metalloproteinases and an L-amino acid oxidase were also present in the venom. Interestingly, an arguable venom component of other species was demonstrated as a true venom protein and named svLIPA (snake venom acid lipase). This finding indicates the importance of checking the actual protein occurrence across species before rejecting genes suggested to code for toxins, which are relevant for the discussion about the early evolution of reptile venoms. Moreover, trends in the evolution of some toxin classes, such as simplification of metalloproteinases and rearrangements of Kunitz and Wap domains, parallel similar phenomena observed in other venomous snake families and provide a broader picture of toxin evolution.

A transcriptomic view of the proteome variability of newborn and adult Bothrops jararaca snake venoms

Background

Snake bite is a neglected public health problem in communities in rural areas of several countries. Bothrops jararaca causes many snake bites in Brazil and previous studies have demonstrated that the pharmacological activities displayed by its venom undergo a significant ontogenetic shift. Similarly, the venom proteome of B. jararaca exhibits a considerable variation upon neonate to adult transition, which is associated with changes in diet from ectothermic prey in early life to endothermic prey in adulthood. Moreover, it has been shown that the Brazilian commercial antibothropic antivenom, which is produced by immunization with adult venom, is less effective in neutralizing newborn venom effects. On the other hand, venom gland transcripts of newborn snakes are poorly known since all transcriptomic studies have been carried out using mRNA from adult specimens.

Methods/Principal Findings

Here we analyzed venom gland cDNA libraries of newborn and adult B. jararaca in order to evaluate whether the variability demonstrated for its venom proteome and pharmacological activities was correlated with differences in the structure of toxin transcripts. The analysis revealed that the variability in B. jararaca venom gland transcriptomes is quantitative, as illustrated by the very high content of metalloproteinases in the newborn venom glands. Moreover, the variability is also characterized by the structural diversity of SVMP precursors found in newborn and adult transcriptomes. In the adult transcriptome, however, the content of metalloproteinase precursors considerably diminishes and the number of transcripts of serine proteinases, C-type lectins and bradykinin-potentiating peptides increase. Moreover, the comparison of the content of ESTs encoding toxins in adult male and female venom glands showed some gender-related differences.

Conclusions/Significance

We demonstrate a substantial shift in toxin transcripts upon snake development and a marked decrease in the metalloproteinase P-III/P-I class ratio which are correlated with changes in the venom proteome complexity and pharmacological activities.

Bothrops jararaca is one of the most abundant venomous snake species in Brazil. It is primarily a nocturnal and generalist animal, however, it exhibits a notable ontogenetic shift in diet, feeding mainly on arthropods, lizards, and amphibians (ectothermic prey) through its juvenile phase and on small mammals (endothermic animals) during adult life. Due to its broad geographical distribution, this species is responsible for the majority of the accidents by Bothrops genus in Brazil. Studies on envenomation cases with newborn and adult B. jararaca snakes have shown distinct patterns, mainly related to blood coagulation disorders, which seems to be prominent in accidents with newborn specimens. Moreover, it has been demonstrated that the Brazilian commercial antibothropic antivenom, which is produced by immunization with adult venom, is less effective in neutralizing newborn venom effects. In this study we analyzed the venom gland transcriptome of newborn snake specimens and compared the content of toxin transcripts with that of adult specimens. We demonstrate that upon B. jararaca development, its repertoire of mRNAs encoding toxins changes both qualitatively and quantitatively and these alterations are associated with the venom proteome profiles and pharmacological activities displayed by newborn and adult specimens.