Revealing molecular determinants governing mambalgin-3 pharmacology at acid-sensing ion channel 1 variants.

Acid-sensing ion channels (ASICs) are trimeric proton-gated cation channels that play a role in neurotransmission and pain sensation. The snake venom-derived peptides, mambalgins, exhibit potent analgesic effects in rodents by inhibiting central ASIC1a and peripheral ASIC1b. Despite their distinct species- and subtype-dependent pharmacology, previous structure-function studies have focussed on the mambalgin interaction with ASIC1a. Currently, the specific channel residues responsible for this pharmacological profile, and the mambalgin pharmacophore at ASIC1b remain unknown. Here we identify non-conserved residues at the ASIC1 subunit interface that drive differences in the mambalgin pharmacology from rat ASIC1a to ASIC1b, some of which likely do not make peptide binding interactions. Additionally, an amino acid variation below the core binding site explains potency differences between rat and human ASIC1. Two regions within the palm domain, which contribute to subtype-dependent effects for mambalgins, play key roles in ASIC gating, consistent with subtype-specific differences in the peptides mechanism. Lastly, there is a shared primary mambalgin pharmacophore for ASIC1a and ASIC1b activity, with certain peripheral peptide residues showing variant-specific significance for potency. Through our broad mutagenesis studies across various species and subtype variants, we gain a more comprehensive understanding of the pharmacophore and the intricate molecular interactions that underlie ligand specificity. These insights pave the way for the development of more potent and targeted peptide analogues required to advance our understating of human ASIC1 function and its role in disease.

Red-on-Yellow Queen: Bio-Layer Interferometry Reveals Functional Diversity Within Micrurus Venoms and Toxin Resistance in Prey Species.

Snakes in the family Elapidae largely produce venoms rich in three-finger toxins (3FTx) that bind to the α 1 subunit of nicotinic acetylcholine receptors (nAChRs), impeding ion channel activity. These neurotoxins immobilize the prey by disrupting muscle contraction. Coral snakes of the genus Micrurus are specialist predators who produce many 3FTx, making them an interesting system for examining the coevolution of these toxins and their targets in prey animals. We used a bio-layer interferometry technique to measure the binding interaction between 15 Micrurus venoms and 12 taxon-specific mimotopes designed to resemble the orthosteric binding region of the muscular nAChR subunit. We found that Micrurus venoms vary greatly in their potency on this assay and that this variation follows phylogenetic patterns rather than previously reported patterns of venom composition. The long-tailed Micrurus tend to have greater binding to nAChR orthosteric sites than their short-tailed relatives and we conclude this is the likely ancestral state. The repeated loss of this activity may be due to the evolution of 3FTx that bind to other regions of the nAChR. We also observed variations in the potency of the venoms depending on the taxon of the target mimotope. Rather than a pattern of prey-specificity, we found that mimotopes modeled after snake nAChRs are less susceptible to Micrurus venoms and that this resistance is partly due to a characteristic tryptophan → serine mutation within the orthosteric site in all snake mimotopes. This resistance may be part of a Red Queen arms race between coral snakes and their prey.

Unveiling Novel Kunitz- and Waprin-Type Toxins in the Micrurus mipartitus Coral Snake Venom Gland: An In Silico Transcriptome Analysis.

Kunitz-type peptide expression has been described in the venom of snakes of the Viperidae, Elapidae and Colubridae families. This work aimed to identify these peptides in the venom gland transcriptome of the coral snake Micrurus mipartitus. Transcriptomic analysis revealed a high diversity of venom-associated Kunitz serine protease inhibitor proteins (KSPIs). A total of eight copies of KSPIs were predicted and grouped into four distinctive types, including short KSPI, long KSPI, Kunitz-Waprin (Ku-WAP) proteins, and a multi-domain Kunitz-type protein. From these, one short KSPI showed high identity with Micrurus tener and Austrelaps superbus. The long KSPI group exhibited similarity within the Micrurus genus and showed homology with various elapid snakes and even with the colubrid Pantherophis guttatus. A third group suggested the presence of Kunitz domains in addition to a whey-acidic-protein-type four-disulfide core domain. Finally, the fourth group corresponded to a transcript copy with a putative 511 amino acid protein, formerly annotated as KSPI, which UniProt classified as SPINT1. In conclusion, this study showed the diversity of Kunitz-type proteins expressed in the venom gland transcriptome of M. mipartitus.

The Cloning and Characterization of a Three-Finger Toxin Homolog (NXH8) from the Coralsnake Micrurus corallinus That Interacts with Skeletal Muscle Nicotinic Acetylcholine Receptors.

Coralsnakes (Micrurus spp.) are the only elapids found throughout the Americas. They are recognized for their highly neurotoxic venom, which is comprised of a wide variety of toxins, including the stable, low-mass toxins known as three-finger toxins (3FTx). Due to difficulties in venom extraction and availability, research on coralsnake venoms is still very limited when compared to that of other Elapidae snakes like cobras, kraits, and mambas. In this study, two previously described 3FTx from the venom of M. corallinus, NXH1 (3SOC1_MICCO), and NXH8 (3NO48_MICCO) were characterized. Using in silico, in vitro, and ex vivo experiments, the biological activities of these toxins were predicted and evaluated. The results showed that only NXH8 was capable of binding to skeletal muscle cells and modulating the activity of nAChRs in nerve-diaphragm preparations. These effects were antagonized by anti-rNXH8 or antielapidic sera. Sequence analysis revealed that the NXH1 toxin possesses eight cysteine residues and four disulfide bonds, while the NXH8 toxin has a primary structure similar to that of non-conventional 3FTx, with an additional disulfide bond on the first loop. These findings add more information related to the structural diversity present within the 3FTx class, while expanding our understanding of the mechanisms of the toxicity of this coralsnake venom and opening new perspectives for developing more effective therapeutic interventions.

High-Voltage Toxin'Roll: Electrostatic Charge Repulsion as a Dynamic Venom Resistance Trait in Pythonid Snakes.

The evolutionary interplay between predator and prey has significantly shaped the development of snake venom, a critical adaptation for subduing prey. This arms race has spurred the diversification of the components of venom and the corresponding emergence of resistance mechanisms in the prey and predators of venomous snakes. Our study investigates the molecular basis of venom resistance in pythons, focusing on electrostatic charge repulsion as a defense against α-neurotoxins binding to the alpha-1 subunit of the postsynaptic nicotinic acetylcholine receptor. Through phylogenetic and bioactivity analyses of orthosteric site sequences from various python species, we explore the prevalence and evolution of amino acid substitutions that confer resistance by electrostatic repulsion, which initially evolved in response to predatory pressure by Naja (cobra) species (which occurs across Africa and Asia). The small African species Python regius retains the two resistance-conferring lysines (positions 189 and 191) of the ancestral Python genus, conferring resistance to sympatric Naja venoms. This differed from the giant African species Python sebae, which has secondarily lost one of these lysines, potentially due to its rapid growth out of the prey size range of sympatric Naja species. In contrast, the two Asian species Python brongersmai (small) and Python bivittatus (giant) share an identical orthosteric site, which exhibits the highest degree of resistance, attributed to three lysine residues in the orthosteric sites. One of these lysines (at orthosteric position 195) evolved in the last common ancestor of these two species, which may reflect an adaptive response to increased predation pressures from the sympatric α-neurotoxic snake-eating genus Ophiophagus (King Cobras) in Asia. All these terrestrial Python species, however, were less neurotoxin-susceptible than pythons in other genera which have evolved under different predatory pressure as: the Asian species Malayopython reticulatus which is arboreal as neonates and juveniles before rapidly reaching sizes as terrestrial adults too large for sympatric Ophiophagus species to consider as prey; and the terrestrial Australian species Aspidites melanocephalus which occupies a niche, devoid of selection pressure from α-neurotoxic predatory snakes. Our findings underline the importance of positive selection in the evolution of venom resistance and suggest a complex evolutionary history involving both conserved traits and secondary evolution. This study enhances our understanding of the molecular adaptations that enable pythons to survive in environments laden with venomous threats and offers insights into the ongoing co-evolution between venomous snakes and their prey.

Structural insight into hormone recognition by the natriuretic peptide receptor-A.

Atrial natriuretic peptide (ANP) plays a central role in the regulation of blood pressure and volume. ANP activities are mediated by natriuretic peptide receptor-A (NPR-A), a single-pass transmembrane receptor harboring intrinsic guanylate cyclase activity. This study investigated the mechanism underlying NPR-A-dependent hormone recognition through the determination of the crystal structures of the NPR-A extracellular hormone-binding domain complexed with full-length ANP, truncated mutants of ANP, and dendroaspis natriuretic peptide (DNP) isolated from the venom of the green Mamba snake, Dendroaspis angusticeps. The bound peptides possessed pseudo-two-fold symmetry, despite the lack of two-fold symmetry in the primary sequences, which enabled the tight coupling of the peptide to the receptor, and evidently contributes to guanylyl cyclase activity. The binding of DNP to the NPR-A was essentially identical to that of ANP; however, the affinity of DNP for NPR-A was higher than that of ANP owing to the additional interactions between distinctive sequences in the DNP and NPR-A. Consequently, our findings provide valuable insights that can be applied to the development of novel agonists for the treatment of various human diseases.

Distinct regulatory networks control toxin gene expression in elapid and viperid snakes.

BACKGROUND: Venom systems are ideal models to study genetic regulatory mechanisms that underpin evolutionary novelty. Snake venom glands are thought to share a common origin, but there are major distinctions between venom toxins from the medically significant snake families Elapidae and Viperidae, and toxin gene regulatory investigations in elapid snakes have been limited. Here, we used high-throughput RNA-sequencing to profile gene expression and microRNAs between active (milked) and resting (unmilked) venom glands in an elapid (Eastern Brown Snake, Pseudonaja textilis), in addition to comparative genomics, to identify cis- and trans-acting regulation of venom production in an elapid in comparison to viperids (Crotalus viridis and C. tigris). RESULTS: Although there is conservation in high-level mechanistic pathways regulating venom production (unfolded protein response, Notch signaling and cholesterol homeostasis), there are differences in the regulation of histone methylation enzymes, transcription factors, and microRNAs in venom glands from these two snake families. Histone methyltransferases and transcription factor (TF) specificity protein 1 (Sp1) were highly upregulated in the milked elapid venom gland in comparison to the viperids, whereas nuclear factor I (NFI) TFs were upregulated after viperid venom milking. Sp1 and NFI cis-regulatory elements were common to toxin gene promoter regions, but many unique elements were also present between elapid and viperid toxins. The presence of Sp1 binding sites across multiple elapid toxin gene promoter regions that have been experimentally determined to regulate expression, in addition to upregulation of Sp1 after venom milking, suggests this transcription factor is involved in elapid toxin expression. microRNA profiles were distinctive between milked and unmilked venom glands for both snake families, and microRNAs were predicted to target a diversity of toxin transcripts in the elapid P. textilis venom gland, but only snake venom metalloproteinase transcripts in the viperid C. viridis venom gland. These results suggest differences in toxin gene posttranscriptional regulation between the elapid P. textilis and viperid C. viridis. CONCLUSIONS: Our comparative transcriptomic and genomic analyses between toxin genes and isoforms in elapid and viperid snakes suggests independent toxin regulation between these two snake families, demonstrating multiple different regulatory mechanisms underpin a venomous phenotype.

Genetic insights: mapping sex-specific loci in Siamese cobra (Naja kaouthia) sheds light on the putative sex determining region.

The location of female-specific/linked loci identified in Siamese cobra (Naja kaouthia) previously has been determined through in silico chromosome mapping of the Indian cobra genome (N. naja) as a reference genome. In the present study, we used in silico chromosome mapping to identify sex-specific and linked loci in Siamese cobra. Many sex-specific and sex-linked loci were successfully mapped on the Z sex chromosome, with 227 of the 475 specific loci frequently mapped in a region covering 57 Mb and positioned at 38,992,675-95,561,177 bp of the Indian cobra genome (N. naja). This suggested the existence of a putative sex-determining region (SDR), with one specific locus (PA100000600) homologous to the TOPBP1 gene. The involvement of TOPBP1 gene may lead to abnormal synaptonemal complexes and meiotic chromosomal defects, resulting in male infertility. These findings offer valuable insights into the genetic basis and functional aspects of sex-specific traits in the Siamese cobra, which will contribute to our understanding of snake genetics and evolutionary biology.

miRNAs derived from cobra venom exosomes contribute to the cobra envenomation.

Currently, there is an increasing amount of evidence indicating that exosomes and the miRNAs they contain are crucial players in various biological processes. However, the role of exosomes and miRNAs in snake venom during the envenomation process remains largely unknown. In this study, fresh venom from Naja atra of different ages (2-month-old, 1-year-old, and 5-year-old) was collected, and exosomes were isolated through ultracentrifugation. The study found that exosomes with inactivated proteins and enzymes can still cause symptoms similar to cobra envenomation, indicating that substances other than proteins and enzymes in exosomes may also play an essential role in cobra envenomation. Furthermore, the expression profiles of isolated exosome miRNAs were analyzed. The study showed that a large number of miRNAs were co-expressed and abundant in cobra venom exosomes (CV-exosomes) of different ages, including miR-2904, which had high expression abundance and specific sequences. The specific miR-2094 derived from CV-exosomes (CV-exo-miR-2904) was overexpressed both in vitro and in vivo. As a result, CV-exo-miR-2904 induced symptoms similar to cobra envenomation in mice and caused liver damage, demonstrating that it plays a crucial role in cobra envenomation. These results reveal that CV-exosomes and the miRNAs they contain play a significant regulatory role in cobra envenomation. Our findings provide new insights for the treatment of cobra bites and the development of snake venom-based medicines.

Domain loss enabled evolution of novel functions in the snake three-finger toxin gene superfamily.

Three-finger toxins (3FTXs) are a functionally diverse family of toxins, apparently unique to venoms of caenophidian snakes. Although the ancestral function of 3FTXs is antagonism of nicotinic acetylcholine receptors, redundancy conferred by the accumulation of duplicate genes has facilitated extensive neofunctionalization, such that derived members of the family interact with a range of targets. 3FTXs are members of the LY6/UPAR family, but their non-toxin ancestor remains unknown. Combining traditional phylogenetic approaches, manual synteny analysis, and machine learning techniques (including AlphaFold2 and ProtT5), we have reconstructed a detailed evolutionary history of 3FTXs. We identify their immediate ancestor as a non-secretory LY6, unique to squamate reptiles, and propose that changes in molecular ecology resulting from loss of a membrane-anchoring domain and changes in gene expression, paved the way for the evolution of one of the most important families of snake toxins.

Comparative Venom Multiomics Reveal the Molecular Mechanisms Driving Adaptation to Diverse Predator-Prey Ecosystems in Closely Related Sea Snakes.

Predator-prey arms races are ideal models for studying the natural selection and adaptive evolution that drive the formation of biological diversity. For venomous snakes, venom is a key bridge linking snakes with their prey, but whether and how venom evolves under the selection of diet remains unclear. Here, we focused on two closely related sea snakes, Hydrophis cyanocinctus and Hydrophis curtus, which show significant differences in prey preferences. Data-independent acquisition (DIA)-based proteomic analysis revealed different degrees of homogeneity in the venom composition of the two snakes, which was consistent with the differential phylogenetic diversity of their prey. By investigating the sequences and structures of three-finger toxins (3FTx), a predominant toxin family in elapid venom, we identified significant differences between the two sea snakes in the binding activity of 3FTx to receptors from different prey populations, which could explain the trophic specialization of H. cyanocinctus. Furthermore, we performed integrated multiomic profiling of the transcriptomes, microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and proteomes of the venom glands; constructed venom-related mRNA-miRNA-lncRNA networks; and identified a series of noncoding RNAs involved in the regulation of toxin gene expression in the two species. These findings are highly informative for elucidating the molecular basis and regulatory mechanisms that account for discrepant venom evolution in response to divergent diets in closely related snakes, providing valuable evidence for the study of coselection and coevolution in predator-prey ecosystems.

Highly Evolvable: Investigating Interspecific and Intraspecific Venom Variation in Taipans (Oxyuranus spp.) and Brown Snakes (Pseudonaja spp.).

Snake venoms are complex mixtures of toxins that differ on interspecific (between species) and intraspecific (within species) levels. Whether venom variation within a group of closely related species is explained by the presence, absence and/or relative abundances of venom toxins remains largely unknown. Taipans (Oxyuranus spp.) and brown snakes (Pseudonaja spp.) represent medically relevant species of snakes across the Australasian region and provide an excellent model clade for studying interspecific and intraspecific venom variation. Using liquid chromatography with ultraviolet and mass spectrometry detection, we analyzed a total of 31 venoms covering all species of this monophyletic clade, including widespread localities. Our results reveal major interspecific and intraspecific venom variation in Oxyuranus and Pseudonaja species, partially corresponding with their geographical regions and phylogenetic relationships. This extensive venom variability is generated by a combination of the absence/presence and differential abundance of venom toxins. Our study highlights that venom systems can be highly dynamical on the interspecific and intraspecific levels and underscores that the rapid toxin evolvability potentially causes major impacts on neglected tropical snakebites.

Polymerase chain reaction-based snake origin tracing in commercial venom crystals by targeting the mitochondrial D-loop.

Snake venom is a valuable raw material for numerous therapeutic formulations because of its life-saving pharmacological potential. However, due to their high price, fake "snake venoms" have captured a significant portion of the global market, and there is currently no reliable reported DNA-based method available for quickly distinguishing between fakes and originals. Therefore, in this study, a set of newly designed snake-specific universal primers targeting mitochondrial D-loop fragments were employed to detect snake origins in commercial venom crystals by only simplex polymerase chain reaction analysis. Under the optimal thermal cycling conditions, only the 145-149 bp snake-specific mitochondrial D-loop fragments from pure and mixed backgrounds were amplified by the newly designed primers. Specificity was achieved by confirming no DNA amplification occurred in the DNA admixture of ten different chordates, and universality by individual DNA amplification of nine different snakes. The primers that efficiently amplified the minimum mitochondrial DNA contained in a total of 10-2 ng in a 10.0 µl reaction were also successfully able to detect the snake origin in commercial cobra venom crystals. These findings suggest that the newly designed primers can be used to differentiate the original and fake commercial snake venom crystals in order to achieve the highest standards of snake venom-based medications through amplifying the snake-specific mitochondrial D-loop fragments.

Virus-like particles displaying conserved toxin epitopes stimulate polyspecific, murine antibody responses capable of snake venom recognition.

Antivenom is currently the first-choice treatment for snakebite envenoming. However, only a low proportion of antivenom immunoglobulins are specific to venom toxins, resulting in poor dose efficacy and potency. We sought to investigate whether linear venom epitopes displayed on virus like particles can stimulate an antibody response capable of recognising venom toxins from diverse medically important species. Bioinformatically-designed epitopes, corresponding to predicted conserved regions of group I phospholipase A2 and three finger toxins, were engineered for display on the surface of hepatitis B core antigen virus like particles and used to immunise female CD1 mice over a 14 weeks. Antibody responses to all venom epitope virus like particles were detectable by ELISA by the end of the immunisation period, although total antibody and epitope specific antibody titres were variable against the different epitope immunogens. Immunoblots using pooled sera demonstrated recognition of various venom components in a diverse panel of six elapid venoms, representing three continents and four genera. Insufficient antibody yields precluded a thorough assessment of the neutralising ability of the generated antibodies, however we were able to test polyclonal anti-PLA2 IgG from three animals against the PLA2 activity of Naja nigricollis venom, all of which showed no neutralising ability. This study demonstrates proof-of-principle that virus like particles engineered to display conserved toxin linear epitopes can elicit specific antibody responses in mice which are able to recognise a geographically broad range of elapid venoms.

Dynamic genetic differentiation drives the widespread structural and functional convergent evolution of snake venom proteinaceous toxins.

BACKGROUND: The explosive radiation and diversification of the advanced snakes (superfamily Colubroidea) was associated with changes in all aspects of the shared venom system. Morphological changes included the partitioning of the mixed ancestral glands into two discrete glands devoted for production of venom or mucous respectively, as well as changes in the location, size and structural elements of the venom-delivering teeth. Evidence also exists for homology among venom gland toxins expressed across the advanced snakes. However, despite the evolutionary novelty of snake venoms, in-depth toxin molecular evolutionary history reconstructions have been mostly limited to those types present in only two front-fanged snake families, Elapidae and Viperidae. To have a broader understanding of toxins shared among extant snakes, here we first sequenced the transcriptomes of eight taxonomically diverse rear-fanged species and four key viperid species and analysed major toxin types shared across the advanced snakes. RESULTS: Transcriptomes were constructed for the following families and species: Colubridae - Helicops leopardinus, Heterodon nasicus, Rhabdophis subminiatus; Homalopsidae - Homalopsis buccata; Lamprophiidae - Malpolon monspessulanus, Psammophis schokari, Psammophis subtaeniatus, Rhamphiophis oxyrhynchus; and Viperidae - Bitis atropos, Pseudocerastes urarachnoides, Tropidolaeumus subannulatus, Vipera transcaucasiana. These sequences were combined with those from available databases of other species in order to facilitate a robust reconstruction of the molecular evolutionary history of the key toxin classes present in the venom of the last common ancestor of the advanced snakes, and thus present across the full diversity of colubroid snake venoms. In addition to differential rates of evolution in toxin classes between the snake lineages, these analyses revealed multiple instances of previously unknown instances of structural and functional convergences. Structural convergences included: the evolution of new cysteines to form heteromeric complexes, such as within kunitz peptides (the beta-bungarotoxin trait evolving on at least two occasions) and within SVMP enzymes (the P-IIId trait evolving on at least three occasions); and the C-terminal tail evolving on two separate occasions within the C-type natriuretic peptides, to create structural and functional analogues of the ANP/BNP tailed condition. Also shown was that the de novo evolution of new post-translationally liberated toxin families within the natriuretic peptide gene propeptide region occurred on at least five occasions, with novel functions ranging from induction of hypotension to post-synaptic neurotoxicity. Functional convergences included the following: multiple occasions of SVMP neofunctionalised in procoagulant venoms into activators of the clotting factors prothrombin and Factor X; multiple instances in procoagulant venoms where kunitz peptides were neofunctionalised into inhibitors of the clot destroying enzyme plasmin, thereby prolonging the half-life of the clots formed by the clotting activating enzymatic toxins; and multiple occasions of kunitz peptides neofunctionalised into neurotoxins acting on presynaptic targets, including twice just within Bungarus venoms. CONCLUSIONS: We found novel convergences in both structural and functional evolution of snake toxins. These results provide a detailed roadmap for future work to elucidate predator-prey evolutionary arms races, ascertain differential clinical pathologies, as well as documenting rich biodiscovery resources for lead compounds in the drug design and discovery pipeline.

Cytotoxin antibody-based colourimetric sensor for field-level differential detection of elapid among big four snake venom.

Development of a rapid, on-site detection tool for snakebite is highly sought after, owing to its clinically and forensically relevant medicolegal significance. Polyvalent antivenom therapy in the management of such envenomation cases is finite due to its poor venom neutralization capabilities as well as diagnostic ramifications manifested as untoward immunological reactions. For precise molecular diagnosis of elapid venoms of the big four snakes, we have developed a lateral flow kit using a monoclonal antibody (AB1; IgG1 - κ chain; Kd: 31 nM) generated against recombinant cytotoxin-7 (rCTX-7; 7.7 kDa) protein of the elapid venom. The monoclonal antibody specifically detected the venoms of Naja naja (p < 0.0001) and Bungarus caeruleus (p<0.0001), without showing any immunoreactivity against the viperidae snakes in big four venomous snakes. The kit developed attained the limit of quantitation of 170 pg/µL and 2.1 ng/µL in spiked buffer samples and 28.7 ng/µL and 110 ng/µL in spiked serum samples for detection of N. naja and B. caeruleus venoms, respectively. This kit holds enormous potential in identification of elapid venom of the big four snakes for effective prognosis of an envenomation; as per the existing medical guidelines.

Venom of the Annulated Sea Snake Hydrophis cyanocinctus: A Biochemically Simple but Genetically Complex Weapon.

Given that the venom system in sea snakes has a role in enhancing their secondary adaption to the marine environment, it follows that elucidating the diversity and function of venom toxins will help to understand the adaptive radiation of sea snakes. We performed proteomic and de novo NGS analyses to explore the diversity of venom toxins in the annulated sea snake (Hydrophis cyanocinctus) and estimated the adaptive molecular evolution of the toxin-coding unigenes and the toxicity of the major components. We found three-finger toxins (3-FTxs), phospholipase A2 (PLA2) and cysteine-rich secretory protein (CRISP) in the venom proteome and 59 toxin-coding unigenes belonging to 24 protein families in the venom-gland transcriptome; 3-FTx and PLA2 were the most abundant families. Nearly half of the toxin-coding unigenes had undergone positive selection. The short- (i.p. 0.09 µg/g) and long-chain neurotoxin (i.p. 0.14 µg/g) presented fairly high toxicity, whereas both basic and acidic PLA2s expressed low toxicity. The toxicity of H. cyanocinctus venom was largely determined by the 3-FTxs. Our data show the venom is used by H. cyanocinctus as a biochemically simple but genetically complex weapon and venom evolution in H. cyanocinctus is presumably driven by natural selection to deal with fast-moving prey and enemies in the marine environment.

Elucidating the Venom Diversity in Sri Lankan Spectacled Cobra (Naja naja) through De Novo Venom Gland Transcriptomics, Venom Proteomics and Toxicity Neutralization.

Inadequate effectiveness of Indian antivenoms in treating envenomation caused by the Spectacled Cobra/Indian Cobra (Naja naja) in Sri Lanka has been attributed to geographical variations in the venom composition. This study investigated the de novo venom-gland transcriptomics and venom proteomics of the Sri Lankan N. naja (NN-SL) to elucidate its toxin gene diversity and venom variability. The neutralization efficacy of a commonly used Indian antivenom product in Sri Lanka was examined against the lethality induced by NN-SL venom in mice. The transcriptomic study revealed high expression of 22 toxin genes families in NN-SL, constituting 46.55% of total transcript abundance. Three-finger toxins (3FTX) were the most diversely and abundantly expressed (87.54% of toxin gene expression), consistent with the dominance of 3FTX in the venom proteome (72.19% of total venom proteins). The 3FTX were predominantly S-type cytotoxins/cardiotoxins (CTX) and α-neurotoxins of long-chain or short-chain subtypes (α-NTX). CTX and α-NTX are implicated in local tissue necrosis and fatal neuromuscular paralysis, respectively, in envenomation caused by NN-SL. Intra-species variations in the toxin gene sequences and expression levels were apparent between NN-SL and other geographical specimens of N. naja, suggesting potential antigenic diversity that impacts antivenom effectiveness. This was demonstrated by limited potency (0.74 mg venom/ml antivenom) of the Indian polyvalent antivenom (VPAV) in neutralizing the NN-SL venom. A pan-regional antivenom with improved efficacy to treat N. naja envenomation is needed.

Two Reference-Quality Sea Snake Genomes Reveal Their Divergent Evolution of Adaptive Traits and Venom Systems.

True sea snakes (Hydrophiini) are among the last and most successful clades of vertebrates that show secondary marine adaptation, exhibiting diverse phenotypic traits and lethal venom systems. To better understand their evolution, we generated the first chromosome-level genomes of two representative Hydrophiini snakes, Hydrophis cyanocinctus and H. curtus. Through comparative genomics we identified a great expansion of the underwater olfaction-related V2R gene family, consisting of more than 1,000 copies in both snakes. A series of chromosome rearrangements and genomic structural variations were recognized, including large inversions longer than 30 megabase (Mb) on sex chromosomes which potentially affect key functional genes associated with differentiated phenotypes between the two species. By integrating multiomics we found a significant loss of the major weapon for elapid predation, three-finger toxin genes, which displayed a dosage effect in H. curtus. These genetic changes may imply mechanisms that drove the divergent evolution of adaptive traits including prey preferences between the two closely related snakes. Our reference-quality sea snake genomes also enrich the repositories for addressing important issues on the evolution of marine tetrapods, and provide a resource for discovering marine-derived biological products.

The omega-loop of cobra cytotoxins tolerates multiple amino acid substitutions.

Cobra cytotoxins (CTs), the three-fingered proteins, feature high amino acid sequence homology in the beta-strands and variations in the loop regions. We selected a pair of cytotoxins from Naja kaouthia crude venom to clarify the sequence-structure relationships. Using chromatography and mass spectroscopy, we separated and identified the mixture of cytotoxins 2 and 3, differentiated by the only Val 41/Ala 41 substitution. Here, using natural abundance 13C, 15N NMR-spectroscopy we performed chemical shift assignments of the signals of the both toxins in aqueous solution in the major and minor forms. Combining NOE and chemical shift data, the toxins' spatial structure was determined. Finally, we proved that the tip of the "finger"-2, or the loop-2 of cytotoxins adopts the shape of an omega-loop with a tightly-bound water molecule in its cavity. Comparison with other NMR and X-ray structures of cytotoxins possessing different amino acid sequences reveals spatial similarity in this family of proteins, including the loop-2 region, previously considered to be flexible.