Snake venom glutaminyl cyclases: Purification, cloning, kinetic study, recombinant expression, and comparison with the human enzyme

Proteome analysis of Daboia russelii venom, a medically important snake from the Indian sub-continent

Daboia russelii is a category-I medically important snake throughout the Indian sub-continent contributing to majority of snakebite incidences in this part of the world. As such, extensive studies on its venom composition and search of efficient and appropriate interventions for its treatment become crucial. In this study, the proteome of Daboia russelii venom from Tanore, Rajshahi, Bangladesh was profiled using a combination of chromatographic and mass spectrometric techniques. A total of 37 different proteins belonging to 11 different snake venom protein families were detected. Proteomics analysis revealed the presence of major phospholipase A2 toxins. Daboiatoxin (both A and B subunits), the main lethal PLA2 toxin in the venom of Daboia siamensis (Myanmar viper) which is neurotoxic, myotoxic and cytotoxic was detected. Presence of Daboxin P, which is a major protein in the venom of Indian Daboia russelii with strong anticoagulant activity, was also observed. Inconsistent distribution of such lethal toxins in the venom of same species calls for more investigations of snake venoms from lesser explored regions and formulation of better alternatives to the current antivenom therapy for efficient treatment.

A current perspective on snake venom composition and constituent protein families

Snake venoms are heterogeneous mixtures of proteins and peptides used for prey subjugation. With modern proteomics there has been a rapid expansion in our knowledge of snake venom composition, resulting in the venom proteomes of 30% of vipers and 17% of elapids being characterised. From the reasonably complete proteomic coverage of front-fanged snake venom composition (179 species-68 species of elapids and 111 species of vipers), the venoms of vipers and elapids contained 42 different protein families, although 18 were only reported in < 5% of snake species. Based on the mean abundance and occurrence of the 42 protein families, they can be classified into 4 dominant, 6 secondary, 14 minor, and 18 rare protein families. The dominant, secondary and minor categories account for 96% on average of a snake's venom composition. The four dominant protein families are: phospholipase A₂ (PLA₂), snake venom metalloprotease (SVMP), three-finger toxins (3FTx), and snake venom serine protease (SVSP). The six secondary protein families are: L-amino acid oxidase (LAAO), cysteine-rich secretory protein (CRiSP), C-type lectins (CTL), disintegrins (DIS), kunitz peptides (KUN), and natriuretic peptides (NP). Venom variation occurs at all taxonomic levels, including within populations. The reasons for venom variation are complex, as variation is not always associated with geographical variation in diet. The four dominant protein families appear to be the most important toxin families in human envenomation, being responsible for coagulopathy, neurotoxicity, myotoxicity and cytotoxicity. Proteomic techniques can be used to investigate the toxicological profile of a snake venom and hence identify key protein families for antivenom immunorecognition.

Two-Dimensional Blue Native/SDS Polyacrylamide Gel Electrophoresis for Analysis of Brazilian Bothrops Snake Venoms

Viperidae snakes are the most important agents of snakebites in Brazil. The protein composition of snake venoms has been frequently analyzed by means of electrophoretic techniques, but the interaction of proteins in venoms has barely been addressed. An electrophoretic technique that has gained prominence to study this type of interaction is blue native polyacrylamide gel electrophoresis (BN-PAGE), which allows for the high-resolution separation of proteins in their native form. These protein complexes can be further discriminated by a second-dimension gel electrophoresis (SDS-PAGE) from lanes cut from BN-PAGE. Once there is no study on the use of bidimensional BN/SDS-PAGE with snake venoms, this study initially standardized the BN/SDS-PAGE technique in order to evaluate protein interactions in Bothrops atrox, Bothrops erythromelas, and Bothrops jararaca snake venoms. Results of BN/SDS-PAGE showed that native protein complexes were present, and that snake venom metalloproteinases and venom serine proteinases maintained their enzymatic activity after BN/SDS-PAGE. C-type lectin-like proteins were identified by Western blotting. Therefore, bidimensional BN/SDS-PAGE proved to be an easy, practical, and efficient method for separating functional venom proteins according to their assemblage in complexes, as well as to analyze their biological activities in further details.

Old World Vipers-A Review about Snake Venom Proteomics of Viperinae and Their Variations

Fine-tuned by millions of years of evolution, snake venoms have frightened but also fascinated humanity and nowadays they constitute potential resources for drug development, therapeutics and antivenoms. The continuous progress of mass spectrometry techniques and latest advances in proteomics workflows enabled toxinologists to decipher venoms by modern omics technologies, so-called ‘venomics’. A tremendous upsurge reporting on snake venom proteomes could be observed. Within this review we focus on the highly venomous and widely distributed subfamily of Viperinae (Serpentes: Viperidae). A detailed public literature database search was performed (2003–2020) and we extensively reviewed all compositional venom studies of the so-called Old-World Vipers. In total, 54 studies resulted in 89 venom proteomes. The Viperinae venoms are dominated by four major, four secondary, six minor and several rare toxin families and peptides, respectively. The multitude of different venomics approaches complicates the comparison of venom composition datasets and therefore we differentiated between non-quantitative and three groups of quantitative workflows. The resulting direct comparisons within these groups show remarkable differences on the intra- and interspecies level across genera with a focus on regional differences. In summary, the present compilation is the first comprehensive up-to-date database on Viperinae venom proteomes and differentiating between analytical methods and workflows.

Glutaminyl Cyclase, Diseases, and Development of Glutaminyl Cyclase Inhibitors

Pyroglutamate (pE) modification, catalyzed mainly by glutaminyl cyclase (QC), is prevalent throughout nature and is particularly important in mammals including humans for the maturation of hormones, peptides, and proteins. In humans, the upregulation of QC is involved in multiple diseases and conditions including Alzheimer's disease, Huntington's disease, melanomas, thyroid carcinomas, accelerated atherosclerosis, septic arthritics, etc. This upregulation catalyzes the generation of modified mediators such as pE-amyloid beta (Aß) and pE-chemokine ligand 2 (CCL2) peptides. Not surprisingly, QC has emerged as a reasonable target for the development of therapeutics to combat these diseases and conditions. In this manuscript the deleterious effects of upregulated QC resulting in disease manifestation are reviewed, along with progress on the development of QC inhibitors.

A Unique Carboxylic-Acid Hydrogen-Bond Network (CAHBN) Confers Glutaminyl Cyclase Activity on M28 Family Enzymes

Proteins with sequence or structure similar to those of di-Zn exopeptidases are usually classified as the M28-family enzymes, including the mammalian-type glutaminyl cyclases (QCs). QC catalyzes protein N-terminal pyroglutamate formation, a posttranslational modification important under many physiological and pathological conditions, and is a drug target for treating neurodegenerative diseases, cancers and inflammatory disorders. Without functional characterization, mammalian QCs and their orthologs remain indistinguishable at the sequence and structure levels from other M28-family proteins, leading to few reported QCs. Here, we show that a low-barrier carboxylic-acid hydrogen-bond network (CAHBN) is required for QC activity and discriminates QCs from M28-family peptidases. We demonstrate that the CAHBN-containing M28 peptidases deposited in the PDB are indeed QCs. Our analyses identify several thousands of QCs from the three domains of life, and we enzymatically and structurally characterize several. For the first time, the interplay between a CAHBN and the binuclear metal-binding center of mammalian QCs is made clear. We found that the presence or absence of CAHBN is a key discriminator for the formation of either the mono-Zn QCs or the di-Zn exopeptidases. Our study helps explain the possible roles of QCs in life.

Integrating Top-Down and Bottom-Up Mass Spectrometric Strategies for Proteomic Profiling of Iranian Saw-Scaled Viper, Echis carinatus sochureki, Venom

Saw-scaled or carpet vipers (genus Echis) are considered to cause a higher global snakebite mortality than any other snake. Echis carinatus sochureki (ECS) is a widely distributed snake species, also found across the thirteen provinces of Iran, where it is assumed to be responsible for the most snakebite envenomings. Here, we collected the Iranian specimens of ECS from three different geographically distinct populations, investigated food habits, and performed toxicity assessment and venom proteome profiling to better understand saw-scaled viper life. Our results show that the prey items most commonly found in all populations were arthropods, with scorpions from the family Buthidae particularly well represented. LD₅₀ (median lethal dose) values of the crude venom demonstrate highly comparable venom toxicities in mammals. Consistent with this finding, venom characterization via top-down and bottom-up proteomics, applied to both crude venoms and size-exclusion chromatographic fractions, revealed highly comparable venom compositions among the different populations. By combining all proteomics data, we identified 22 protein families from 102 liquid chromatography and tandem mass spectrometry (LC-MS/MS) raw files, including the most abundant snake venom metalloproteinases (SVMPs, 29-34%); phospholipase A2 (PLA2s, 26-31%); snake venom serine proteinases (SVSPs, 11-12%); l-amino acid oxidases (LAOs, 8-11%), C-type lectins/lectin-like (CTLs, 7-9%) protein families, and many newly detected ones, e.g., renin-like aspartic proteases (RLAPs), fibroblast growth factors (FGFs), peptidyl-prolyl cis-trans isomerases (PPIs), and venom vasodilator peptides (VVPs). Furthermore, we identified and characterized methylated, acetylated, and oxidized proteoforms relating to the PLA2 and disintegrin toxin families and the site of their modifications. It thus seems that post-translational modifications (PTMs) of toxins, particularly target lysine residues, may play an essential role in the structural and functional properties of venom proteins and might be able to influence the therapeutic response of antivenoms, to be investigated in future studies.

Delineating the venom toxin arsenal of Malabar pit viper (Trimeresurus malabaricus) from the Western Ghats of India and evaluating its immunological cross-reactivity and in vitro cytotoxicity

The venom protein components of Malabar pit viper (Trimeresurus malabaricus) were identified by combining SDS-PAGE and ion-exchange chromatography pre-fractionation techniques with LC-MS/MS incorporating Novor and PEAKS-assisted de novo sequencing strategies. Total 97 proteins that belong to 16 protein families such as L-amino acid oxidase, metalloprotease, serine protease, phospholipase A2, 5'-nucleotidase, C-type lectins/snaclecs and disintegrin were recognized from the venom of a single exemplar species. Of the 97 proteins, eighteen were identified through de novo approaches. Immunological cross-reactivity assessed through ELISA and western blot indicate that the Indian antivenoms binds less effectively to Malabar pit viper venom components compared to that of Russell's viper venom. The in vitro cell viability assays suggest that compared to the normal cells, MPV venom induces concentration dependent cell death in various cancer cells. Moreover, crude venom resulted in chromatin condensation and apoptotic bodies implying the induction of apoptosis. Taken together, the present study enabled in dissecting the venom proteome of Trimeresurus malabaricus and revealed the immuno-cross-reactivity profiles of commercially available Indian polyvalent antivenoms that, in turn, is expected to provide valuable insights on the need in improving antivenom preparations against its bite.

Comprehensive Study of the Proteome and Transcriptome of the Venom of the Most Venomous European Viper: Discovery of a New Subclass of Ancestral Snake Venom Metalloproteinase Precursor-Derived Proteins

The nose-horned viper, its nominotypical subspecies Vipera ammodytes ammodytes (Vaa), in particular, is, medically, one of the most relevant snakes in Europe. The local and systemic clinical manifestations of poisoning by the venom of this snake are the result of the pathophysiological effects inflicted by enzymatic and nonenzymatic venom components acting, most prominently, on the blood, cardiovascular, and nerve systems. This venom is a very complex mixture of pharmacologically active proteins and peptides. To help improve the current antivenom therapy toward higher specificity and efficiency and to assist drug discovery, we have constructed, by combining transcriptomic and proteomic analyses, the most comprehensive library yet of the Vaa venom proteins and peptides. Sequence analysis of the venom gland cDNA library has revealed the presence of messages encoding 12 types of polypeptide precursors. The most abundant are those for metalloproteinase inhibitors (MPis), bradykinin-potentiating peptides (BPPs), and natriuretic peptides (NPs) (all three on a single precursor), snake C-type lectin-like proteins (snaclecs), serine proteases (SVSPs), P-II and P-III metalloproteinases (SVMPs), secreted phospholipases A₂ (sPLA₂s), and disintegrins (Dis). These constitute >88% of the venom transcriptome. At the protein level, 57 venom proteins belonging to 16 different protein families have been identified and, with SVSPs, sPLA₂s, snaclecs, and SVMPs, comprise ∼80% of all venom proteins. Peptides detected in the venom include NPs, BPPs, and inhibitors of SVSPs and SVMPs. Of particular interest, a transcript coding for a protein similar to P-III SVMPs but lacking the MP domain was also found at the protein level in the venom. The existence of such proteins, also supported by finding similar venom gland transcripts in related snake species, has been demonstrated for the first time, justifying the proposal of a new P-IIIe subclass of ancestral SVMP precursor-derived proteins.

A Single Mutation Unlocks Cascading Exaptations in the Origin of a Potent Pitviper Neurotoxin

Evolutionary innovations and complex phenotypes seemingly require an improbable amount of genetic change to evolve. Rattlesnakes display two dramatically different venom phenotypes. Type I venoms are hemorrhagic with low systemic toxicity and high expression of tissue-destroying snake venom metalloproteinases. Type II venoms are highly neurotoxic and lack snake venom metalloproteinase expression and associated hemorrhagic activity. This dichotomy hinges on Mojave toxin (MTx), a phospholipase A2 (PLA2) based β-neurotoxin expressed in Type II venoms. MTx is comprised of a nontoxic acidic subunit that undergoes extensive proteolytic processing and allosterically regulates activity of a neurotoxic basic subunit. Evolution of the acidic subunit presents an evolutionary challenge because the need for high expression of a nontoxic venom component and the proteolytic machinery required for processing suggests genetic changes of seemingly little immediate benefit to fitness. We showed that MTx evolved through a cascading series of exaptations unlocked by a single nucleotide change. The evolution of one new cleavage site in the acidic subunit unmasked buried cleavage sites already present in ancestral PLA2s, enabling proteolytic processing. Snake venom serine proteases, already present in the venom to disrupt prey hemostasis, possess the requisite specificities for MTx acidic subunit proteolysis. The dimerization interface between MTx subunits evolved by exploiting a latent, but masked, hydrophobic interaction between ancestral PLA2s. The evolution of MTx through exaptation of existing functional and structural features suggests complex phenotypes that depend on evolutionary innovations can arise from minimal genetic change enabled by prior evolution.

Proteomic Characterization of Two Medically Important Malaysian Snake Venoms, Calloselasma rhodostoma (Malayan Pit Viper) and Ophiophagus hannah (King Cobra)

Calloselasma rhodostoma (CR) and Ophiophagus hannah (OH) are two medically important snakes found in Malaysia. While some studies have described the biological properties of these venoms, feeding and environmental conditions also influence the concentration and distribution of snake venom toxins, resulting in variations in venom composition. Therefore, a combined proteomic approach using shotgun and gel filtration chromatography, analyzed by tandem mass spectrometry, was used to examine the composition of venoms from these Malaysian snakes. The analysis revealed 114 proteins (15 toxin families) and 176 proteins (20 toxin families) in Malaysian Calloselasma rhodostoma and Ophiophagus hannah species, respectively. Flavin monoamine oxidase, phospholipase A₂, phosphodiesterase, snake venom metalloproteinase, and serine protease toxin families were identified in both venoms. Aminopeptidase, glutaminyl-peptide cyclotransferase along with ankyrin repeats were identified for the first time in CR venom, and insulin, c-type lectins/snaclecs, hepatocyte growth factor, and macrophage colony-stimulating factor together with tumor necrosis factor were identified in OH venom for the first time. Our combined proteomic approach has identified a comprehensive arsenal of toxins in CR and OH venoms. These data may be utilized for improved antivenom production, understanding pathological effects of envenoming, and the discovery of biologically active peptides with medical and/or biotechnological value.

Venom Ontogeny in the Mexican Lance-Headed Rattlesnake (Crotalus polystictus)

As trophic adaptations, rattlesnake venoms can vary in composition depending on several intrinsic and extrinsic factors. Ontogenetic changes in venom composition have been documented for numerous species, but little is known of the potential age-related changes in many rattlesnake species found in México. In the current study, venom samples collected from adult and neonate Crotalus polystictus from Estado de México were subjected to enzymatic and electrophoretic analyses, toxicity assays (LD₅₀), and MALDI-TOF mass spectrometry, and a pooled sample of adult venom was analyzed by shotgun proteomics. Electrophoretic profiles of adult males and females were quite similar, and only minor sex-based variation was noted. However, distinct differences were observed between venoms from adult females and their neonate offspring. Several prominent bands, including P-I and P-III snake venom metalloproteinases (SVMPs) and disintegrins (confirmed by MS/MS) were present in adult venoms and absent/greatly reduced in neonate venoms. Age-dependent differences in SVMP, kallikrein-like, phospholipase A₂ (PLA₂), and L-amino acid oxidase (LAAO) activity levels were confirmed by enzymatic activity assays, and like many other rattlesnake species, venoms from adult snakes have higher SVMP activity than neonate venoms. Conversely, PLA₂ activity was approximately 2.5 × greater in venoms from neonates, likely contributing to the increased toxicity (neonate venom LD₅₀ = 4.5 μg/g) towards non-Swiss albino mice when compared to adult venoms (LD₅₀ = 5.5 μg/g). Thrombin-like (TLE) and phosphodiesterase activities did not vary significantly with age. A significant effect of sex (between adult male and adult female venoms) was also observed for SVMP, TLE, and LAAO activities. Analysis of pooled adult venom by LC-MS/MS identified 14 toxin protein families, dominated by bradykinin-inhibitory peptides, SVMPs (P-I, P-II and P-III), disintegrins, PLA₂s, C-type-lectins, CRiSPs, serine proteinases, and LAAOs (96% of total venom proteins). Neonate and adult C. polystictus in this population consume almost exclusively mammals, suggesting that age-based differences in composition are related to physical differences in prey (e.g., surface-to-volume ratio differences) rather than taxonomic differences between prey. Venoms from adult C. polystictus fit a Type I pattern (high SVMP activity, lower toxicity), which is characteristic of many larger-bodied rattlesnakes of North America.

Agkihpin, a Distinct SVTLE from the Venom of Gloydius halys Pallas: Purification, Characterization and Structure-Activity Determination

Blood clots produced by snake-venom thrombin-like enzymes (SVTLEs) are cleared rapidly, which makes SVTLEs attractive as potential candidates for antithrombotic therapy. We isolated a SVTLE, agkihpin, from the venom of Gloydius halys Pallas. Agkihpin was confirmed to a single-chain TLE with molecular mass of 25.5 kD, pI of 7.43, optimal pH of 8.0 (hydrolyzing TAME), linked carbohydrate absent, and weak fibrinogen clotting activity. It was also found that (i) G. halys might be the latest species in SVTLEs phylogenetic tree; (ii) different level of conservation was shown among the SVTLEs from the Viperidae snakes. Some of those site may account for different activities exhibited by those SVTLEs, especially position 181, at which a fibrinogenolytic activity increase was found when a basic and larger amino acid substituted by a neutral and smaller one; (iii) an extra α-helix constructed with a 'Pro + acidic amino acid + aromatic amino acid' pattern was found in the SVTLEs from Gloydius and Agkistrodon snakes, although it does not necessarily imply an effect on the fibrinogenolytic activity of the SVTLEs. This study provided some new insight into the activity of SVTLE.

Proteomic analysis of venom variability and ontogeny across the arboreal palm-pitvipers (genus Bothriechis)

Bothriechis is a genus of eleven currently recognized slender and arboreal venomous snakes, commonly called palm-pitvipers that range from southern Mexico to northern South America. Despite dietary studies suggesting that palm-pitvipers are generalists with an ontogenetic shift toward endothermic prey, venom proteomic analyses have revealed remarkable divergence between the venoms of the Costa Rican species, B. lateralis, B. schlegelii, B. supraciliaris, and B. nigroviridis. To achieve a more complete picture of the venomic landscape across Bothriechis, the venom proteomes of biodiversity of the northern Middle American highland palm-pitvipers, B. thalassinus, B. aurifer, and B. bicolor from Guatemala, B. marchi from Honduras, and neonate Costa Rican B. lateralis and B. schlegelii, were investigated. B. thalassinus and B. aurifer venoms are comprised by similar toxin arsenals dominated by SVMPs (33-39% of the venom proteome), CTLs (11-16%), BPP-like molecules (10-13%), and CRISPs (5-10%), and are characterized by the absence of PLA₂ proteins. Conversely, the predominant (35%) components of B. bicolor are D49-PLA₂ molecules. The venom proteome of B. marchi is similar to B. aurifer and B. thalassinus in that it is rich in SVMPs and BPPs, but also contains appreciable amounts (14.3%) of PLA₂s. The major toxin family found in the venoms of both neonate B. lateralis and B. schlegelii, is serine proteinase (SVSP), comprising about 20% of their toxin arsenals. The venom of neonate B. schlegelii is the only palm-pitviper venom where relative high amounts of Kunitz-type (6.3%) and γPLA₂ (5.2%) inhibitors have been identified. Despite notable differences between their proteomes, neonate venoms are more similar to each other than to adults of their respective species. However, the ontogenetic changes taking place in the venom of B. lateralis strongly differ from those that occur in the venom of B. schlegelii. Thus, the ontogenetic change in B. lateralis produces a SVMP-rich venom, whereas in B. schlegelii the age-dependent compositional shift generates a PLA₂-rich venom. Overall, genus-wide venomics illustrate the high evolvability of palm-pitviper venoms. The integration of the pattern of venom variation across Bothriechis into a phylogenetic and biogeographic framework may lay the foundation for assessing, in future studies, the evolutionary path that led to the present-day variability of the venoms of palm-pitvipers.

SIGNIFICANCE

Bothriechis represents a monophyletic basal genus of eleven arboreal palm-pitvipers that range from southern Mexico to northern South America. Despite palm-pitvipers' putative status as diet generalists, previous proteomic analyses have revealed remarkable divergence between the venoms of Costa Rican species, B. lateralis, B. schlegelii, B. supraciliaris, and B. nigroviridis. Our current proteomic study of Guatemalan species, B. thalassinus, B. aurifer, and B. bicolor, Honduran B. marchi, and neonate B. lateralis and B. schlegelii from Costa Rica was undertaken to deepen our understanding of the evolutionary pattern of venom proteome diversity across Bothriechis. Ancestral characters are often, but not always, preserved in an organism's development. Venoms of neonate B. lateralis and B. schlegelii are more similar to each other than to adults of their respective species, suggesting that the high evolvability of palm-pitviper venoms may represent an inherent feature of Bothriechis common ancestor. Our genus-wide data identified four nodes of venom phenotype differentiation across the phylogeny of Bothriechis. Integrated into a phylogenetic and biogeographic framework, the pattern of venom variation across Bothriechis may lay the groundwork to establish whether divergence was driven by selection for efficient resource exploitation in arboreal 'islands', thereby contributing to the ecological speciation of the genus.

Proteomic Characterization and Comparison of Malaysian Tropidolaemus wagleri and Cryptelytrops purpureomaculatus Venom Using Shotgun-Proteomics

Tropidolaemus wagleri and Cryptelytrops purpureomaculatus are venomous pit viper species commonly found in Malaysia. Tandem mass spectrometry analysis of the crude venoms has detected different proteins in T. wagleri and C. purpureomaculatus. They were classified into 13 venom protein families consisting of enzymatic and nonenzymatic proteins. Enzymatic families detected in T. wagleri and C. purpureomaculatus venom were snake venom metalloproteinase, phospholipase A₂, l-amino acid oxidase, serine proteases, 5′-nucleotidase, phosphodiesterase, and phospholipase B. In addition, glutaminyl cyclotransferase was detected in C. purpureomaculatus. C-type lectin-like proteins were common nonenzymatic components in both species. Waglerin was present and unique to T. wagleri—it was not in C. purpureomaculatus venom. In contrast, cysteine-rich secretory protein, bradykinin-potentiating peptide, and C-type natriuretic peptide were present in C. purpureomaculatus venom. Composition of the venom proteome of T. wagleri and C. purpureomaculatus provides useful information to guide production of effective antivenom and identification of proteins with potential therapeutic applications.