Exploring Toxin Genes of Myanmar Russell's Viper, Daboia siamensis, through De Novo Venom Gland Transcriptomics.

The Russell's viper (Daboia siamensis) is a medically important venomous snake in Myanmar. Next-generation sequencing (NGS) shows potential to investigate the venom complexity, giving deeper insights into snakebite pathogenesis and possible drug discoveries. mRNA from venom gland tissue was extracted and sequenced on the Illumina HiSeq platform and de novo assembled by Trinity. The candidate toxin genes were identified via the Venomix pipeline. Protein sequences of identified toxin candidates were compared with the previously described venom proteins using Clustal Omega to assess the positional homology among candidates. Candidate venom transcripts were classified into 23 toxin gene families including 53 unique full-length transcripts. C-type lectins (CTLs) were the most highly expressed, followed by Kunitz-type serine protease inhibitors, disintegrins and Bradykinin potentiating peptide/C-type natriuretic peptide (BPP-CNP) precursors. Phospholipase A2, snake venom serine proteases, metalloproteinases, vascular endothelial growth factors, L-amino acid oxidases and cysteine-rich secretory proteins were under-represented within the transcriptomes. Several isoforms of transcripts which had not been previously reported in this species were discovered and described. Myanmar Russell's viper venom glands displayed unique sex-specific transcriptome profiles which were correlated with clinical manifestation of envenoming. Our results show that NGS is a useful tool to comprehensively examine understudied venomous snakes.

Biogeographic venom variation in Russell's viper (Daboia russelii) and the preclinical inefficacy of antivenom therapy in snakebite hotspots.

BACKGROUND: Snakebite in India results in over 58,000 fatalities and a vast number of morbidities annually. The majority of these clinically severe envenomings are attributed to Russell's viper (Daboia russelii), which has a near pan-India distribution. Unfortunately, despite its medical significance, the influence of biogeography on the composition and potency of venom from disparate D. russelii populations, and the repercussions of venom variation on the neutralisation efficacy of marketed Indian antivenoms, remain elusive. METHODS: Here, we employ an integrative approach comprising proteomic characterisation, biochemical analyses, pharmacological assessment, and venom toxicity profiling to elucidate the influence of varying ecology and environment on the pan-Indian populations of D. russelii. We then conducted in vitro venom recognition experiments and in vivo neutralisation assays to evaluate the efficacy of the commercial Indian antivenoms against the geographically disparate D. russelii populations. FINDINGS: We reveal significant intraspecific variation in the composition, biochemical and pharmacological activities and potencies of D. russelii venoms sourced from five distinct biogeographic zones across India. Contrary to our understanding of the consequences of venom variation on the effectiveness of snakebite therapy, commercial antivenom exhibited surprisingly similar neutralisation potencies against the majority of the investigated populations, with the exception of low preclinical efficacy against the semi-arid population from northern India. However, the ability of Indian antivenoms to counter the severe morbid effects of Daboia envenoming remains to be evaluated. CONCLUSION: The concerning lack of antivenom efficacy against the north Indian population of D. russelii, as well as against two other 'big four' snake species in nearby locations, underscores the pressing need to develop pan-India effective antivenoms with improved efficacy in high snakebite burden locales.

Widespread and Differential Neurotoxicity in Venoms from the Bitis Genus of Viperid Snakes.

Research into the neurotoxic activity of venoms from species within the snake family Viperidae is relatively neglected compared with snakes in the Elapidae family. Previous studies into venoms from the Bitis genus of vipers have identified the presence of presynaptic phospholipase A2 neurotoxins in B. atropos and B. caudalis, as well as a postsynaptic phospholipase A2 in B. arietans. Yet, no studies have investigated how widespread neurotoxicity is across the Bitis genus or if they exhibit prey selectivity of their neurotoxins. Utilising a biolayer interferometry assay, we were able to assess the binding of crude venom from 14 species of Bitis to the neuromuscular α-1 nAChR orthosteric site across a wide range of vertebrate taxa mimotopes. Postsynaptic binding was seen for venoms from B. arietans, B. armata, B. atropos, B. caudalis, B. cornuta, B. peringueyi and B. rubida. To further explore the types of neurotoxins present, venoms from the representatives B. armata, B. caudalis, B. cornuta and B. rubida were additionally tested in the chick biventer cervicis nerve muscle preparation, which showed presynaptic and postsynaptic activity for B. caudalis and only presynaptic neurotoxicity for B. cornuta and B. rubida, with myotoxicity also evident for some species. These results, combined with the biolayer interferometry results, indicate complex neurotoxicity exerted by Bitis species, which varies dramatically by lineage tested upon. Our data also further support the importance of sampling across geographical localities, as significant intraspecific variation of postsynaptic neurotoxicity was reported across the different localities.

Mass spectrometric analysis to unravel the venom proteome composition of Indian snakes: opening new avenues in clinical research.

INTRODUCTION: The 'Big Four' venomous snakes - Daboia russelii, Naja naja, Bungarus caeruleus, and Echis carinatus - are primarily responsible for the majority of snake envenomation in India. Several other lesser-known venomous snake species also inflict severe envenomation in the country. AREAS COVERED: A comprehensive analysis of the venom proteome composition of the 'Big Four' and other medically important venomous snakes of India and the effect of regional variation in venom composition on immunorecognition and/or neutralization by commercial antivenom was undertaken by searching the literature (from 1985 to date) available in large public databases. Further, mass spectrometric identification of poorly immunogenic toxins of snake venom (against which commercial polyvalent antivenom contains a significantly lower proportion of antibodies) and its impact on antivenom therapy against snakebite are discussed. The application of mass spectrometry to identify protein (toxin) complexes as well as drug prototypes from Indian snake venoms and the clinical importance of such studies are also highlighted. EXPERT OPINION: Further detailed clinical and proteomic research is warranted to better understand the effects of regional snake venom composition on the clinical manifestation of envenomation and antivenom therapy and to improve the production of antibodies against poorly immunogenic venom components.

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 A2 (sPLA2s), 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, sPLA2s, 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.

Proteomic analysis reveals geographic variation in venom composition of Russell's Viper in the Indian subcontinent: implications for clinical manifestations post-envenomation and antivenom treatment.

INTRODUCTION: The Russell's Viper (RV) (Daboia russelii), a category I medically important snake, is responsible for a significant level of morbidity and mortality in the Indian sub-continent. Areas covered: The current review highlights the variation in RV venom (RVV) composition from different geographical locales on the Indian sub-continent, as revealed by biochemical and proteomic analyses. A comparison of these RVV proteomes revealed significant differences in the number of toxin isoforms and relative toxin abundances, highlighting the impact of geographic location on RVV composition. Antivenom efficacy studies have shown differential neutralization of toxicity and enzymatic activity of different RVV samples from the Indian sub-continent by commercial polyvalent antivenom (PAV). The proteome analysis has provided deeper insights into the variation of RVV composition leading to differences in antivenom efficacy and severity of clinical manifestations post RV-envenomation across the Indian sub-continent. Expert commentary: Variation in RVV antigenicity due to geographical differences and poor recognition of low molecular mass (<20 kDa) RVV toxins by PAV are serious concerns for effective antivenom treatment against RV envenomation. Improvements in immunization protocols that take into account the poorly immunogenic components and geographic variation in RVV composition, can lead to better hospital management of RV bite patients.

Comparative Profiling of Three Atheris Snake Venoms: A. squamigera, A. nitschei and A. chlorechis.

A proteomic and transcriptomic comparative analysis of the venoms of three Atheris species (A. squamigera, A. nitschei and A. chlorechis) was carried out by size exclusion liquid chromatography, gel electrophoresis, mass spectrometry, and mRNA sequencing. The improved proteomic profiling utilised in this work was combined with transcript studies, advancing our insights into venom composition, protein distribution and inter-species variation among the three bush vipers. Crude venoms of all three samples contained at least 10-20 protein components, ranging in size from ≤ 3 to > 98 kDa. Both approaches yielded converging overall information, pointing to phospholipases, disintegrins, serine proteases and metalloproteases as the major toxin classes, which are likely to explain the local and systemic symptoms observed in envenomation by Atheris genus. Being considered as the main factors involved in the distinct venom-induced pathologies, these identified snake venom proteins are of particular interest in terms of understanding their physiological and biological function as well as for their contribution in potential medical treatments.

Accelerated evolution of toxin genes: Exonization and intronization in snake venom disintegrin/metalloprotease genes.

Toxin genes in animals undergo accelerated evolution compared to non-toxin genes to be effective and competitive in prey capture, as well as to enhance their predator defense. Several mechanisms have been proposed to explain this unusual phenomenon. These include (a) frequent mutations in exons compared to introns and nonsynonymous substitutions in exons; (b) high frequency of point mutations are due to the presence of more unstable triplets in exons compared to introns; (c) Accelerated Segment Switch in Exons to alter Targeting (ASSET); (d) Rapid Accumulation of Variations in Exposed Residues (RAVERs); (e) alteration in intron-exon boundary; (f) deletion of exon; and (g) loss/gain of domains through recombination. By systematic analyses of snake venom disintegrin/metalloprotease genes, I describe a new mechanism in the evolution of these genes through exonization and intronization. In the evolution of RTS/KTS disintegrins, a new exon (10a) is formed in intron 10 of the disintegrin/metalloprotease gene. Unlike more than 90% new exons that are from repetitive elements in introns, exon 10a originated from a non-repetitive element. To incorporate exon 10a, part of the exon 11 is intronized to retain the open reading frame. This is the first case of simultaneous exonization and intronization within a single gene. This new mechanism alters the function of toxins through drastic changes to the molecular surface via insertion of new exons and deletion of exons.

Functional characterization of recombinant snake venom rhodocytin: rhodocytin mutant blocks CLEC-2/podoplanin-dependent platelet aggregation and lung metastasis.

Essentials We generated recombinant rhodocytin that could aggregate platelets via CLEC-2. Recombinant wild-type rhodocytin formed heterooctamer with four α- and ß-subunits. Asp 4 in α-subunit of rhodocytin was required for binding to CLEC-2. Inhibitory mutant of rhodocytin blocked podoplanin-dependent hematogenous metastasis. SUMMARY: Background Rhodocytin, a disulfide-linked heterodimeric C-type lectin from Calloselasma rhodostoma consisting of α-subunits and ß-subunits, induces platelet aggregation through C-type lectin-like receptor 2 (CLEC-2). CLEC-2 is a physiological binding partner of podoplanin (PDPN), which is expressed on some tumor cell types, and is involved in tumor cell-induced platelet aggregation and tumor metastasis. Thus, modified rhodocytin may be a possible source of anti-CLEC-2 drugs for both antiplatelet and antimetastasis therapy. However, its molecular function has not been well characterized, because of the lack of recombinant rhodocytin that induces platelet aggregation. Objective To produce recombinant rhodocytin, in order to verify its function with mutagenesis, and to develop an anti-CLEC-2 drug based on the findings. Methods We used Chinese hamster ovary cells to express recombinant rhodocytin (wild-type [WT] and mutant), which was analyzed for induction/inhibition of platelet aggregation with light transmission aggregometry, the formation of multimers with blue native PAGE, and binding to CLEC-2 with flow cytometry. Finally, we investigated whether mutant rhodocytin could suppress PDPN-induced metastasis in an experimental lung metastasis mouse model. Results Functional WT] rhodocytin (αWTßWT) was obtained by coexpression of both subunits. Asp4 in α-subunits of rhodocytin was required for CLEC-2 binding. αWTßWT formed a heterooctamer similarly to native rhodocytin. Moreover, an inhibitory mutant of rhodocytin (αWTßK53A/R56A), forming a heterotetramer, bound to CLEC-2 without inducing platelet aggregation, and blocked CLEC-2-PDPN interaction-dependent platelet aggregation and experimental lung metastasis. Conclusion These findings provide molecular characterization information on rhodocytin, and suggest that mutant rhodocytin could be used as a therapeutic agent to target CLEC-2.

Population Genomic Analysis of a Pitviper Reveals Microevolutionary Forces Underlying Venom Chemistry.

Venoms are among the most biologically active secretions known, and are commonly believed to evolve under extreme positive selection. Many venom gene families, however, have undergone duplication, and are often deployed in doses vastly exceeding the LD50 for most prey species, which should reduce the strength of positive selection. Here, we contrast these selective regimes using snake venoms, which consist of rapidly evolving protein formulations. Though decades of extensive studies have found that snake venom proteins are subject to strong positive selection, the greater action of drift has been hypothesized, but never tested. Using a combination of de novo genome sequencing, population genomics, transcriptomics, and proteomics, we compare the two modes of evolution in the pitviper, Protobothrops mucrosquamatus. By partitioning selective constraints and adaptive evolution in a McDonald-Kreitman-type framework, we find support for both hypotheses: venom proteins indeed experience both stronger positive selection, and lower selective constraint than other genes in the genome. Furthermore, the strength of selection may be modulated by expression level, with more abundant proteins experiencing weaker selective constraint, leading to the accumulation of more deleterious mutations. These findings show that snake venoms evolve by a combination of adaptive and neutral mechanisms, both of which explain their extraordinarily high rates of molecular evolution. In addition to positive selection, which optimizes efficacy of the venom in the short term, relaxed selective constraints for deleterious mutations can lead to more rapid turnover of individual proteins, and potentially to exploration of a larger venom phenotypic space.

Understanding Russell's viper venom factor V activator's substrate specificity by surface plasmon resonance and in-silico studies.

Blood coagulation factor V (FV) is activated either by Factor X or thrombin, cleaving at three different sites viz., Site I (Arg709-Ser710), site II (Arg1018-Thr1019), and site III (Arg1545-Ser1546). Russell's viper venom factor V activator (RVV-V) is a thrombin-like serine proteinase that activates FV with selective, single cleavage at site III. A long lasting effort is being pending in understanding the 'selective' binding specificity of the RVV-V towards site III. Here, we present the binding kinetic study of RVV-V with two designed peptides corresponding to the regions from site I (Gln699-Asn713) and site II (1008Lys-Pro1022), respectively, that include 15 amino acids. Our investigation for justifying the binding efficacy and kinetics of peptides includes SPR method, protein-peptide docking, molecular dynamics simulation, and principal component analysis (PCA). Surprisingly, the SPR experiment disclosed that the Peptide II showed a lower binding affinity with KD of 2.775 mM while the Peptide I showed none. Docking and simulation of both the peptides with RVV-V engaged either rooted or shallow binding for Peptide II and Peptide I respectively. The peptide binding resulted in global conformational changes in the native fold of RVV-V, whereas the similar studies for thrombin failed to make major changes in the native fold. In support, the PCA analysis for RVV-V showed the dislocation of catalytic triad upon binding both the peptides. Hence, RVV-V, a serine protease, is incompetent in cleaving these two sites. This study suggests a transition in RVV-V from the native rigid to the distorted flexible structure and paves a way to design a new peptide substrate/inhibitor.

Neurotoxicity in Sri Lankan Russell's Viper (Daboia russelii) Envenoming is Primarily due to U1-viperitoxin-Dr1a, a Pre-Synaptic Neurotoxin.

Russell's vipers are snakes of major medical importance in Asia. Russell's viper (Daboia russelii) envenoming in Sri Lanka and South India leads to a unique, mild neuromuscular paralysis, not seen in other parts of the world where the snake is found. This study aimed to identify and pharmacologically characterise the major neurotoxic components of Sri Lankan Russell's viper venom. Venom was fractionated using size exclusion chromatography and reverse-phase high-performance liquid chromatography (RP-HPLC). In vitro neurotoxicities of the venoms, fractions and isolated toxins were measured using chick biventer and rat hemidiaphragm preparations. A phospholipase A2 (PLA2) toxin, U1-viperitoxin-Dr1a (13.6 kDa), which constitutes 19.2 % of the crude venom, was isolated and purified using HPLC. U1-viperitoxin-Dr1a produced concentration-dependent in vitro neurotoxicity abolishing indirect twitches in the chick biventer nerve-muscle preparation, with a t 90 of 55 ± 7 min only at 1 µM. The toxin did not abolish responses to acetylcholine and carbachol indicating pre-synaptic neurotoxicity. Venom, in the absence of U1-viperitoxin-Dr1a, did not induce in vitro neurotoxicity. Indian polyvalent antivenom, at the recommended concentration, only partially prevented the neurotoxic effects of U1-viperitoxin-Dr1a. Liquid chromatography mass spectrometry analysis confirmed that U1-viperitoxin-Dr1a was the basic S-type PLA2 toxin previously identified from this venom (NCBI-GI: 298351762; SwissProt: P86368). The present study demonstrates that neurotoxicity following Sri Lankan Russell's viper envenoming is primarily due to the pre-synaptic neurotoxin U1-viperitoxin-Dr1a. Mild neurotoxicity observed in severely envenomed Sri Lankan Russell's viper bites is most likely due to the low potency of U1-viperitoxin-Dr1a, despite its high relative abundance in the venom.

Strong and widespread action of site-specific positive selection in the snake venom Kunitz/BPTI protein family.

S1 family of serine peptidases is the largest family of peptidases. They are specifically inhibited by the Kunitz/BPTI inhibitors. Kunitz domain is characterized by the compact 3D structure with the most important inhibitory loops for the inhibition of S1 peptidases. In the present study we analysed the action of site-specific positive selection and its impact on the structurally and functionally important parts of the snake venom Kunitz/BPTI family of proteins. By using numerous models we demonstrated the presence of large numbers of site-specific positively selected sites that can reach between 30-50% of the Kunitz domain. The mapping of the positively selected sites on the 3D model of Kunitz/BPTI inhibitors has shown that these sites are located in the inhibitory loops 1 and 2, but also in the Kunitz scaffold. Amino acid replacements have been found exclusively on the surface, and the vast majority of replacements are causing the change of the charge. The consequence of these replacements is the change in the electrostatic potential on the surface of the Kunitz/BPTI proteins that may play an important role in the precise targeting of these inhibitors into the active site of S1 family of serine peptidases.

Insights into the Evolution of a Snake Venom Multi-Gene Family from the Genomic Organization of Echis ocellatus SVMP Genes.

The molecular events underlying the evolution of the Snake Venom Metalloproteinase (SVMP) family from an A Disintegrin And Metalloproteinase (ADAM) ancestor remain poorly understood. Comparative genomics may provide decisive information to reconstruct the evolutionary history of this multi-locus toxin family. Here, we report the genomic organization of Echis ocellatus genes encoding SVMPs from the PII and PI classes. Comparisons between them and between these genes and the genomic structures of Anolis carolinensis ADAM28 and E. ocellatus PIII-SVMP EOC00089 suggest that insertions and deletions of intronic regions played key roles along the evolutionary pathway that shaped the current diversity within the multi-locus SVMP gene family. In particular, our data suggest that emergence of EOC00028-like PI-SVMP from an ancestral PII(e/d)-type SVMP involved splicing site mutations that abolished both the 3' splice AG acceptor site of intron 12* and the 5' splice GT donor site of intron 13*, and resulted in the intronization of exon 13* and the consequent destruction of the structural integrity of the PII-SVMP characteristic disintegrin domain.

Full-Length Venom Protein cDNA Sequences from Venom-Derived mRNA: Exploring Compositional Variation and Adaptive Multigene Evolution.

Envenomation of humans by snakes is a complex and continuously evolving medical emergency, and treatment is made that much more difficult by the diverse biochemical composition of many venoms. Venomous snakes and their venoms also provide models for the study of molecular evolutionary processes leading to adaptation and genotype-phenotype relationships. To compare venom complexity and protein sequences, venom gland transcriptomes are assembled, which usually requires the sacrifice of snakes for tissue. However, toxin transcripts are also present in venoms, offering the possibility of obtaining cDNA sequences directly from venom. This study provides evidence that unknown full-length venom protein transcripts can be obtained from the venoms of multiple species from all major venomous snake families. These unknown venom protein cDNAs are obtained by the use of primers designed from conserved signal peptide sequences within each venom protein superfamily. This technique was used to assemble a partial venom gland transcriptome for the Middle American Rattlesnake (Crotalus simus tzabcan) by amplifying sequences for phospholipases A2, serine proteases, C-lectins, and metalloproteinases from within venom. Phospholipase A2 sequences were also recovered from the venoms of several rattlesnakes and an elapid snake (Pseudechis porphyriacus), and three-finger toxin sequences were recovered from multiple rear-fanged snake species, demonstrating that the three major clades of advanced snakes (Elapidae, Viperidae, Colubridae) have stable mRNA present in their venoms. These cDNA sequences from venom were then used to explore potential activities derived from protein sequence similarities and evolutionary histories within these large multigene superfamilies. Venom-derived sequences can also be used to aid in characterizing venoms that lack proteomic profiles and identify sequence characteristics indicating specific envenomation profiles. This approach, requiring only venom, provides access to cDNA sequences in the absence of living specimens, even from commercial venom sources, to evaluate important regional differences in venom composition and to study snake venom protein evolution.

VTBuilder: a tool for the assembly of multi isoform transcriptomes.

BACKGROUND: Within many research areas, such as transcriptomics, the millions of short DNA fragments (reads) produced by current sequencing platforms need to be assembled into transcript sequences before they can be utilized. Despite recent advances in assembly software, creating such transcripts from read data harboring isoform variation remains challenging. This is because current approaches fail to identify all variants present or they create chimeric transcripts within which relationships between co-evolving sites and other evolutionary factors are disrupted. We present VTBuilder, a tool for constructing non-chimeric transcripts from read data that has been sequenced from sources containing isoform complexity. RESULTS: We validated VTBuilder using reads simulated from 54 Sanger sequenced transcripts (SSTs) expressed in the venom gland of the saw scaled viper, Echis ocellatus. The SSTs were selected to represent genes from major co-expressed toxin groups known to harbor isoform variants. From the simulated reads, VTBuilder constructed 55 transcripts, 50 of which had a greater than 99% sequence similarity to 48 of the SSTs. In contrast, using the popular assembler tool Trinity (r2013-02-25), only 14 transcripts were constructed with a similar level of sequence identity to just 11 SSTs. Furthermore VTBuilder produced transcripts with a similar length distribution to the SSTs while those produced by Trinity were considerably shorter. To demonstrate that our approach can be scaled to real world data we assembled the venom gland transcriptome of the African puff adder Bitis arietans using paired-end reads sequenced on Illumina's MiSeq platform. VTBuilder constructed 1481 transcripts from 5 million reads and, following annotation, all major toxin genes were recovered demonstrating reconstruction of complex underlying sequence and isoform diversity. CONCLUSION: Unlike other approaches, VTBuilder strives to maintain the relationships between co-evolving sites within the constructed transcripts, and thus increases transcript utility for a wide range of research areas ranging from transcriptomics to phylogenetics and including the monitoring of drug resistant parasite populations. Additionally, improving the quality of transcripts assembled from read data will have an impact on future studies that query these data. VTBuilder has been implemented in java and is available, under the GPL GPU V0.3 license, from http:// http://www.lstmed.ac.uk/vtbuilder .

Folding molecular dynamics simulations accurately predict the effect of mutations on the stability and structure of a vammin-derived peptide.

Folding molecular dynamics simulations amounting to a grand total of 4 µs of simulation time were performed on two peptides (with native and mutated sequences) derived from loop 3 of the vammin protein and the results compared with the experimentally known peptide stabilities and structures. The simulations faithfully and accurately reproduce the major experimental findings and show that (a) the native peptide is mostly disordered in solution, (b) the mutant peptide has a well-defined and stable structure, and (c) the structure of the mutant is an irregular ß-hairpin with a non-glycine ß-bulge, in excellent agreement with the peptide's known NMR structure. Additionally, the simulations also predict the presence of a very small ß-hairpin-like population for the native peptide but surprisingly indicate that this population is structurally more similar to the structure of the native peptide as observed in the vammin protein than to the NMR structure of the isolated mutant peptide. We conclude that, at least for the given system, force field, and simulation protocol, folding molecular dynamics simulations appear to be successful in reproducing the experimentally accessible physical reality to a satisfactory level of detail and accuracy.

Recombinant snake venom metalloproteinase inhibitor BJ46A inhibits invasion and metastasis of B16F10 and MHCC97H cells through reductions of matrix metalloproteinases 2 and 9 activities.

Studies have shown that the recombinant BJ46a (rBJ46a) protein can reduce matrix metalloproteinase (MMP) activities and inhibit invasion and metastasis of melanoma cells. Here, we optimized the Pichia pastoris system to evaluate rBJ46a protein as an anticancer agent. The Enzchek gelatinase/collagenase assay showed that rBJ46a inhibited MMP activities (IC50=0.119 mg/ml). Kinetic analyses using a series of double reciprocal Lineweaver-Burk plots (1/V vs. 1/S) showed a competitive mode of inhibition with rBJ46a with inhibitory efficiency against MMPs (Ki=13.6 nmol/l). Matrigel invasion assays showed significant activity of rBJ46a on tumor cells. For lung colonization assays, C57BL/6 mice were inoculated in the lateral tail vein with B16F10 cells and were treated with three i.v. injections of rBJ46a (1, 2, and 4 mg/kg) 24 h before cell inoculation, and 2 and 24 h after cell inoculation. Administration of rBJ46a suppressed lung tumor colony formation significantly. For spontaneous metastasis assays, MHCC97H cells were inoculated subcutaneously into nude mice. After 24 h, rBJ46a was administered by i.p. injections: 1, 2, and 4 mg/kg once daily for 6 days. rBJ46a decreased lung tumor colony formation significantly. Gelatin zymography showed that MMP2/MMP9 enzymatic activities in tumor cells were suppressed by rBJ46a in a dose-dependent manner, and the Km values of rBJ46a against MMP2 and MMP9 activities that were expressed in both B16F10 and MHCC97H cells were 3.6 and 1.4 µmol/l, respectively. Thus, rBJ46a can inhibit the invasion and metastasis of tumor cells by reducing MMP2/MMP9 activities, indicating that rBJ46a may be a novel therapeutic agent for antimetastasis of tumor cells.

Purification, characterization and molecular cloning of chymotrypsin inhibitor peptides from the venom of Burmese Daboia russelii siamensis.

One novel Kunitz BPTI-like peptide designated as BBPTI-1, with chymotrypsin inhibitory activity was identified from the venom of Burmese Daboia russelii siamensis. It was purified by three steps of chromatography including gel filtration, cation exchange and reversed phase. A partial N-terminal sequence of BBPTI-1, HDRPKFCYLPADPGECLAHMRSF was obtained by automated Edman degradation and a Ki value of 4.77nM determined. Cloning of BBPTI-1 including the open reading frame and 3' untranslated region was achieved from cDNA libraries derived from lyophilized venom using a 3' RACE strategy. In addition a cDNA sequence, designated as BBPTI-5, was also obtained. Alignment of cDNA sequences showed that BBPTI-5 exhibited an identical sequence to BBPTI-1 cDNA except for an eight nucleotide deletion in the open reading frame. Gene variations that represented deletions in the BBPTI-5 cDNA resulted in a novel protease inhibitor analog. Amino acid sequence alignment revealed that deduced peptides derived from cloning of their respective precursor cDNAs from libraries showed high similarity and homology with other Kunitz BPTI proteinase inhibitors. BBPTI-1 and BBPTI-5 consist of 60 and 66 amino acid residues respectively, including six conserved cysteine residues. As these peptides have been reported to have influence on the processes of coagulation, fibrinolysis and inflammation, their potential application in biomedical contexts warrants further investigation.

Ammodytoxins efficiently release arachidonic acid and induce apoptosis in a motoneuronal cell line in an enzymatic activity-dependent manner.

Secreted phospholipases A2 (sPLA2s) are phospholipolytic enzymes and receptor ligands whose action affects cell death and survival. We have previously shown that ammodytoxin A (AtxA), a snake venom sPLA2, is rapidly internalized into motoneuronal NSC34 cells, inducing characteristic neurotoxic sPLA2 cell damage and apoptosis. In this study, we have analyzed the role of sPLA2 enzymatic activity, including arachidonic acid (AA) release, in the induction of motoneuronal apoptosis by AtxA and homologous recombinant sPLA2s with different enzymatic properties: an AtxA mutant (V31W) with very high enzymatic activity, enzymatically inactive S49-sPLA2 (ammodytin L, AtnL), its mutant (LW) with restored enzymatic activity, and non-toxic, enzymatically active sPLA2 (AtnI2). Addition of AA, AtxA, AtxA-V31W and AtnL-LW, but not AtnL and AtnI2, to NSC34 cells resulted in caspase-3 activation, DNA fragmentation and disruption of mitochondrial membrane potential, leading to a significant and rapid decrease in motoneuronal cell viability that was not observed in C2C12 myoblasts and HEK293 cells. AtxA, AtxA-V31W and AtnL-LW, but not AtnL and AtnI2, also liberated large amounts of AA specifically from motoneuronal cells, and this ability correlated well with the ability to induce apoptotic changes and decrease cell viability. The enzymatic activity of AtxA and similar sPLA2s is thus necessary, but not sufficient, for inducing motoneuronal apoptosis. This suggests that specific binding to the motoneuronal cell surface, followed by internalization and enzymatic activity-dependent induction of apoptosis, possibly as a consequence of extensive extra- and intracellular AA release, is necessary for Atx-induced motoneuronal cell death.