Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Bacillus subtilis Type I antitoxin SR6 Promotes Degradation of Toxin yonT mRNA and Is Required to Prevent Toxic yoyJ Overexpression
Previous Article in Journal
PhcrTx2, a New Crab-Paralyzing Peptide Toxin from the Sea Anemone Phymanthus crucifer
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 10
Issue 2
10.3390/toxins10020073
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Pharmacokinetics of Snake Venom

by
Suchaya Sanhajariya
1,2,
Stephen B. Duffull
2 and
Geoffrey K. Isbister
1,*
1
Clinical Toxicology Research Group, University of Newcastle, Newcastle 2298, Australia
2
School of Pharmacy, University of Otago, Dunedin 9016, New Zealand
*
Author to whom correspondence should be addressed.
Toxins 2018, 10(2), 73; https://doi.org/10.3390/toxins10020073
Submission received: 9 December 2017 / Revised: 31 January 2018 / Accepted: 3 February 2018 / Published: 7 February 2018
(This article belongs to the Section Animal Venoms)
Browse Figures
Versions Notes

Abstract

:
Understanding snake venom pharmacokinetics is essential for developing risk assessment strategies and determining the optimal dose and timing of antivenom required to bind all venom in snakebite patients. This review aims to explore the current knowledge of snake venom pharmacokinetics in animals and humans. Literature searches were conducted using EMBASE (1974–present) and Medline (1946–present). For animals, 12 out of 520 initially identified studies met the inclusion criteria. In general, the disposition of snake venom was described by a two-compartment model consisting of a rapid distribution phase and a slow elimination phase, with half-lives of 5 to 48 min and 0.8 to 28 h, respectively, following rapid intravenous injection of the venoms or toxins. When the venoms or toxins were administered intramuscularly or subcutaneously, an initial absorption phase and slow elimination phase were observed. The bioavailability of venoms or toxins ranged from 4 to 81.5% following intramuscular administration and 60% following subcutaneous administration. The volume of distribution and the clearance varied between snake species. For humans, 24 out of 666 initially identified publications contained sufficient information and timed venom concentrations in the absence of antivenom therapy for data extraction. The data were extracted and modelled in NONMEM. A one-compartment model provided the best fit, with an elimination half-life of 9.71 ± 1.29 h. It is intended that the quantitative information provided in this review will provide a useful basis for future studies that address the pharmacokinetics of snakebite in humans.
Keywords:
snakes; venom; pharmacokinetics; elapid; viper; toxins
Key Contribution:
This review investigates the current knowledge of snake venom pharmacokinetics. In addition we conducted an exploratory population pharmacokinetic analysis of timed venom concentration data extracted from the literature, which can be used as a prior for the future development of substantial snake venom pharmacokinetic models.

Graphical Abstract

1. Introduction

In 2009, WHO listed snakebite as a neglected tropical disease, recognising its importance alongside many infectious diseases [1]. This is a particularly important public health issue in tropical and subtropical regions [2,3] mostly affecting those who live in rural areas, including the agricultural workforce in developing countries [3,4,5]. Snake envenomation causes significant morbidity and mortality usually requiring hospitalisation and occasionally causing permanent disabilities and in severe cases death [6]. Delayed access to appropriate medical facilities, lack of antivenom, and limited supportive treatments are considered the main contributors to the high morbidity and mortality [3,4].
Venomous snakes with medical importance are predominantly front-fanged, originating from three families: Atractaspididae, Elapidae, and Viperidae [7]. Snakes from the Viperidae family can be further divided into two subfamilies, Viperinae (true vipers) and Crotalinae (pit vipers). In all families, the venom is produced in a specialised gland and typically delivered to the target organism via modified teeth [8,9]. The major functions of snake venoms are to facilitate prey capture by immobilisation or death, to assist with snake digestion, and to act as a defensive mechanism against predators. Snake venoms comprise a mixture of biologically active proteins and polypeptides (comprising approximately 90–95% of a venom load), along with other non-protein components including carbohydrates, lipids, amines, and inorganic salts (Table 1 and Table 2) [7,10,11]. Proteins and polypeptides can be classified into (1) enzymes [e.g., phospholipase A2 (PLA2), metalloproteases (SVMP), serine proteases (SVSP), and L-amino acid oxidases (LAAO)], and (2) non-enzymatic substances [e.g., three-finger toxins (3FTx), kunitz peptides (KUN), and disintegrins (DIS)]. The composition of snake venoms varies depending on a variety of factors, including snake family, genus and species, geographical location, typical prey type, and age and size of the snake [12]. For instance, 3FTx and PLA2 generally predominate in venoms of elapid snakes, while PLA2, SVMP, and SVSP are the most common in vipers, and 3FTx are essentially absent in vipers [7]. Recently, Tasoulis and Isbister (2017) comprehensively reviewed the snake venom proteomes from various compositional studies and classified proteins into five groups based on their abundance and ubiquity [7]. In this work, it was determined that more than 90% of elapid and viper venoms are comprised of 10 protein families. Four of these were considered dominant protein families (PLA2, SVMP, SVSP, and 3FTx) and six were secondary protein families (Cysteine-rich secretory proteins (CRiSP), LAAO, KUN, DIS, C-type lectins (CTL), and natriuretic peptides (NP)) [7]. Understanding the variations in snake venom composition is important and will affect studies of snake venoms, their pharmacokinetics, and ultimately the development, production, and treatments of venomous snakebites using antivenoms.
Systemic snake envenomation in humans involves the injected venom being absorbed and entering the systemic circulation, which results in a range of clinical effects depending on the particular snake and the venom components [13]. Snakes from the family Elapidae (cobras, kraits, mambas, Australian venomous terrestrial snakes, coral snakes, and sea snakes) and from the family Viperidae (true vipers and pit vipers) are responsible for most severe envenomation in humans [2,6]. The cocktail of biologically active proteins in venoms gives rise to a wide variety of local and systemic effects (Table 1 and Table 2). The local effects range from local pain and swelling to tissue necrosis [13,14]. Local tissue damage is caused by a range of myotoxic and cytotoxic components such as non-catalytic PLA2 myotoxins [15,16], SVMP [15,17], and cytotoxins [15], which can lead to extensive tissue necrosis and potentially require debridement and amputation [5]. The systemic effects are potentially life-threatening and include coagulopathy, neurotoxicity, acute kidney injury, and myotoxicity. Venom-induced consumptive coagulopathy is arguably the most common systemic envenoming syndrome worldwide and is caused by many vipers and Australasian elapids [18]. Neurotoxicity occurs mainly following elapids bite and results from a blockade at the neuromuscular junction, resulting in the paralysis of the facial, bulbar, respiratory, and limb muscles that can lead to respiratory failure and death [2,19].
The measurement of snake venom in blood is rarely available in the clinical setting, but has been undertaken in research studies of snakebite. A number of assay techniques have been developed for the detection of snake venom and venom antigens (toxins) in biological samples, including enzyme-linked immunosorbent assay (ELISA) and radioimmunoassay (RIA). Venom assays have been used in research mainly to identify the snake type involved in envenoming and to evaluate the efficacy of an antivenom in binding the free venom [20]. A venom-specific ELISA was first developed by Theakston and colleagues in 1977 [21] and has good sensitivity and specificity, except in some cases where cross-reactivity may occur from closely related species of snakes [22,23,24]. The limits of detection for a snake venom or toxin are reported between 0.1 mcg/L and 20 mcg/L [25,26,27,28,29,30,31,32,33,34]. The intra- and inter-assay coefficients of variability reported are 1 to 20% and 4 to 20%, respectively [26,27,35]. ELISA, in comparison to RIA, is more practical, easier to perform, more cost-efficient, less time-consuming, and confers no health hazard from handling radioisotopes [21,23,24,36]. The ELISA technique has been used to develop snake venom detection kits, which are used as diagnostic tools for assisting clinicians to make decisions about the appropriate antivenom, but tests swabs from the skin for venom, rather than blood [13,23].
The treatment of snake envenomation principally involves the administration of specific antivenoms for the snake species or type involved, and supportive care. Snake antivenoms are antibodies or antibody fragments which are derived from the plasma of animals (typically horses and sheep) that have been immunised with a snake venom [63]. The intravenously injected antibodies (IgG immunoglobulin) or antibody fragments [F(ab’)2 or Fab] bind and neutralise free venom in the patient’s plasma, which aims to reverse or prevent further toxic effects. However, animal-derived antivenoms also carry a risk of hypersensitivity reactions because they are foreign proteins, which can result in cutaneous and multiorgan reactions that are potentially life-threatening with severe anaphylaxis [5,63]. Antivenoms remain expensive, and their supply is often limited in some regions, so more effective dosing and antivenom use may improve their cost-effectiveness.
The dose and timing of antivenom administration is still largely empirical and often based on in vitro neutralisation studies in animals, in which the toxic effects differ compared to humans [64]. In practice, the dose administered is determined by the treating clinicians on the basis of subjective symptoms or clinical or laboratory results. New approaches are required to improve antivenom dosing and determine the optimal dose and timing of antivenom administration. The knowledge of snake venom pharmacokinetics (the study of the time course of venom distribution in biological systems) will provide important information about the time course of the exposure to venom. This will allow for a better determination of a sufficient dose and timing of antivenom administration. In particular, it remains unclear how long post-bite the administration of antivenom remains an effective therapeutic intervention.
This paper aims to bring together the current knowledge of the pharmacokinetics of snake venom by reviewing: (1) laboratory studies performed in animals to investigate the time course of venom concentrations in plasma; and (2) reports describing the time course of venom concentrations in humans.

2. Results

2.1. Pharmacokinetic Studies in Animals

2.1.1. Literature Search

The initial search generated 493 studies for screening, after duplicates were removed. Of these, 481 studies were excluded on the basis of the inclusion and exclusion criteria (see Figure 1). This left 12 studies available for data extraction and to investigate the characteristics and pharmacokinetic parameters of snake venom in animals.

2.1.2. Demographics of Snake Venom Pharmacokinetic Studies in Animals

Of the 12 studies included, eight studies investigated the pharmacokinetics of whole venom, three investigated the pharmacokinetics of isolated snake toxins, and one looked at both venom and purified toxins. Eight of the 12 studies were performed in rabbits. The pharmacokinetic parameters reported in these studies were estimated from the plasma concentration time course following the administration of snake venoms or toxins via the intravenous (IV), intramuscular (IM), or subcutaneous (SC) routes. A summary of the pharmacokinetic parameters reported in these 12 studies are shown in Table 3, Table 4 and Table 5.

2.1.3. Pharmacokinetic Parameters of Snake Venoms and Toxins Following Intravenous Administration

In most studies, the investigators identified a two-compartment output model as the best fit for their data. Overall, regardless of the type of snake, the concentration initially decreased rapidly with a half-life ranging from, for whole venoms, 15 (Cryptelytrops purpureomaculatus) to 48 min (Naja sumatrana and Hypnale hypnale), for toxins, 5 (habutobin of Trimeresurus flavoviridis) to 42 min (PLA2 of N. sumatrana). In the terminal phase, the concentrations decreased more slowly with longer half-lives ranging from, for whole venoms, 12 (Vipera aspis and Bothrops alternatus) to 28 h (C. purpureomaculatus), for toxins, 0.8 (habutobin of T. flavoviridis) to 12 h (PLA2 of N. sumatrana). It is assumed that the initial phase is the distribution and the terminal phase is the elimination, but this cannot be determined from the data presented. Habutobin, a toxin purified from the venom of T. flavoviridis, had the shortest distribution and elimination half-lives compared to the rest of the snake venoms and toxins. Also of note, N. naja atra cytotoxin had a shorter half-life than the venom and toxins from two other types of cobra studied (N. sumatrana and N. sputatrix).
Hart et al. (2014) [65] and Paniagua et al. (2012) [66] only reported single values for the half-lives from a one-compartment model. It is interpreted that the disposition phase half-life would be a composite of the distribution and elimination phases reported with the two-compartment models. Tan et al. (2014) [67] reported a three-compartment model fit for their data, with three half-lives (24 min, 48 min, and 19 h, respectively). The half-lives were similar to the two-compartment models.
The reported volume of distribution at steady state (Vss) for snake venoms and toxins ranged widely. The Vss of snake venoms ranged from as low as 0.12 L·kg−1 in M. fulvius venom to 1.2 L·kg−1 in V. aspis venom. The Vss was notably higher for the thrombin-like enzyme (TLE) from A. halys ussuriensis Emelianov venom (1.8 L·kg−1).
The systemic clearance (CL) of snake venoms and toxins reported also varied from species to species. The CL of snake venom ranged from 0.007 L·h−1·kg−1 for H. hypnale venom to 0.093 L·h−1·kg−1 for M. fulvius venom, and from 0.048 to as high as 0.324 L·h−1·kg−1 for snake toxins (PLA2 in N. sumatrana and TLE from A. halys ussuriensis Emelianov venom, respectively). It appears that the reported values of CL of snake toxins are generally larger than those reported for whole snake venoms.

2.1.4. Pharmacokinetic Parameters of Snake Venoms and Toxins Following Intramuscular or Subcutaneous Administration

When snake venom or toxins are administered via the IM or SC routes, the additional influence of absorption rate and bioavailability (F) of venom and toxins needs to be considered. Guo et al., 1993 [27], Audebert et al., 1994 [35], Sim et al., 2013 [74], Yap et al., 2013 [73], Tan et al., 2014 [67], Yap et al., 2014 [70] investigated the pharmacokinetics of snake venom following IM administration, while Zhao et al., 2001 [68] and Paniagua et al., 2012 [66] investigated the pharmacokinetics of snake venom following SC administration. All studies demonstrated evidence of an absorption phase of venom and toxins. Not all venom or toxins were fully absorbed into the circulation following IM administration [F ranging from 4 (H. hypnale venom) to 67% (V. aspis venom), and from 40 to 82% for N. sumatrana toxins] and SC administration [F = 60% for M. fulvius venom]. In these studies, the distribution half-life, which is expected to be related to the ongoing absorption that occurs during the initial distribution phase, was not commonly seen. The studies by Zhao et al., 2001 [68] and Guo et al., 1993 [27] are the only two studies that report the half-life of the absorption phase as well as two phases of disposition half-lives. In these two studies, the half-lives following SC administration were reported to be 2.5, 4.8, and 125 h, respectively [68], and those following IM administration were 0.077, 0.37, and 5.9 h, respectively [27].
Guo et al., 1993 [27], Audebert et al., 1994 [35], and Yap et al., 2013 [73] found that the terminal half-life following IM injection was longer than that following IV administration. However, in the studies by Sim et al., 2013 [74], Yap et al., 2013 [73], Yap et al., 2014 [70], and Tan et al., 2014 [67], the terminal half-lives reported following IV and IM administration appeared similar. Both Zhao et al., 2001 [68] and Paniagua et al., 2012 [66] found the terminal half-life to be longer following SC administration than after IV injection. In these cases, it is expected that the prolonged duration of the half-life would be a function of prolonged absorption since the disposition processes should be unaffected by the input processes.

2.2. Pharmacokinetic Studies in Humans

2.2.1. Literature Search

The initial search generated 576 studies for screening, after duplicates were removed. A total of 552 studies were further removed on the basis of the inclusion and exclusion criteria (Figure 2), leaving 26 studies for data extraction and ultimately post-hoc (model-based) analysis.

2.2.2. Data Extraction

Timed venom concentration data were extracted from 24 (Table 6) of the 26 studies included from the literature review. One study was not included for further analysis, as the serum antigen concentrations appeared to have increased over time post-bite for all patients over a period of 24 h, which was not consistent with data reported in the other studies. Another study was excluded from data extraction because of the difficulty in distinguishing individual data points from the overlapping concentration-time profiles of the 37 patients, and, hence, we were unable to digitise the data from the figure. Nine studies reported elapid envenomation, and sixteen studies reported viper envenomation.
From the 24 included studies, we were able to retrieve data for 145 individuals. A total of 218 timed concentration data were obtained from the text of the case reports. For studies in which the results were reported graphically, we successfully recreated similar plots by digitising the data points from the figures using the WebPlotDigitizer software.

2.2.3. Data Analysis

A plot of the extracted data is shown in Figure 3 and Figure 4. Most data were collected within the first few hours since snakebite, and most individuals provided only a single timed venom concentration. This is because most studies primarily reported a clinical presentation of snakebite, and only one sample was taken prior to the administration of the antivenom. Only five studies reported serial venom concentrations over time in the absence of antivenom. The dataset obtained was sparse in terms of observations per individual, and there were large magnitudes of difference between the snake venom concentrations. In some cases, snake venom concentrations could still be detected 24 h after the bite. Figure 4 compares the concentration-time profiles of patients bitten by snakes from the Elapidae family (Figure 4a) and the Viperidae family (Figure 4b). It can be seen that the venom concentrations detected in patients bitten by vipers are typically higher than those in patients bitten by elapids.
The extracted data were modelled using NONMEM. The concentration-time data were best described by a one-compartmental model with zero-order input and first-order elimination. A prediction-corrected visual predictive check is presented in Figure 5. The observed percentiles are mostly within the 95% confidence interval of simulated percentiles, suggesting that the model provides a good description of the data. The half-life calculated from the estimated values of clearance and volume of distribution was 9.71 ± 1.29 h (see Table A1 of Appendix A for parameter estimates from the model). The model also accounted for the relative exposure (i.e., “dose”) of patients, resulting from bites of snakes from the Elapidae and Viperidae family. It was apparent from these data that vipers inject approximately 75% more venom than elapids, but it should be noted the between-subject variability in the exposure was large, corresponding to 275%, meaning that the variability between bites within a family is at least as significant as the differences between families of snakes. It is possible that some of this variability could be further explained by the snake species, but these data were not available.

3. Discussion

3.1. Pharmacokinetic Studies in Animals

To date, there is only a handful of experimental studies on the pharmacokinetics of snake venom in the absence of antivenom. Most of these studies were performed in large animals, commonly rabbits, as they allowed for the analysis of multiple samples over a longer period of time, consistent with the duration of envenomation. Since a snake venom is comprised of a cocktail of protein components, it is common that snake venoms are investigated as whole venom rather than as individual toxins. Some investigators took the additional step of isolating the major toxins of interest in order to study their profiles [68,69,70]. For the purpose of this review, we only assessed the pharmacokinetic parameters that arose from studies carried out using ELISA, although there were some other studies performed using RIA methods [94,95,96].
The pharmacokinetic profiles of snake venoms and isolated toxins were investigated in animals following intravenous injection. This avoids issues related to bioavailability and provides profiles that are dependent only on distribution and clearance. Most studies demonstrated that the concentration-time profiles of snake venoms and toxins from both snake families were best described by a two-compartment model. The concentrations reportedly declined rapidly during the initial phase. This suggests that the venoms (or toxins) were rapidly distributed out of the central compartment to peripheral compartments. The initial half-life values reported for the whole venoms and the isolated toxins ranged from 5 min to 48 min, which is significantly faster than the terminal half-lives which could be as long as 28 h. Given the low CL values and the relatively high Vss values of the venoms, we suspect that the terminal phase half-life from the two-compartment model, may be predominantly contributed to by elimination. This does not rule out further distribution components. The volume of distribution of snake venom in animals was relatively large when compared to their blood volume (approximately 0.054–0.070 L·kg−1 for rats, 0.057–0.065 L·kg−1 for rabbits, and 0.058–0.064 L·kg−1 for sheep), which supports the significant distribution of the venom outside of blood. This means that an antivenom will be much less effective (potentially ineffective) in binding the venom once the venom leaves the circulation. Understanding the distribution of a venom between the central compartment and the peripheral compartments and, importantly, the rate at which a venom moves out of the central compartment is essential to determining the effective timing of antivenom administration.
Pharmacokinetic studies in animals were also performed following intramuscular or subcutaneous injection in order to mimic the real-life nature of snakebites. The initial phase of the concentration-time profile was governed by the absorption and simultaneous distribution of the snake venom (and toxins). Some studies observed multiple absorption phases which may be due to the difference in the absorption rate of venom components with a wide range of molecular weights (Table 1 and Table 2). In addition, the lymphatic absorption from the injection site to the blood compartment may be associated with a delayed increase in venom concentration observed in these studies. The absorption of venom appeared to be relatively slow, on the basis of the prolonged elimination half-lives reported, especially following SC administration. This indicates that the absorption is likely to be the rate-limiting step for snake venom disposition, and that the ‘flip-flop’ phenomenon may govern the pharmacokinetics of snake venom. In contrast, some studies reported a similar elimination half-life following IV and IM injection, which indicates that the slow elimination is likely to be the rate-limiting factor of the terminal phase of snake venom time course.
Not all snake venoms are absorbed from the injection site, with their bioavailability being as low as 4% in some cases and as high as 81.5% in other cases. The unabsorbed components of snake venom retained at the injection site may possibly be responsible for local tissue damage. Some confounding estimates of Vss characterized by large values were observed in some studies (Guo et al., 1993 [27], Zhao et al., 2001 [68] and Paniagua et al., 2012 [66]), compared with the Vss following IV injection.
It is important to keep in mind that snake venom is a mixture of large and small molecular weight molecules and that there are limitations in using ELISA to measure such a mixture, because multiple different antibodies are used to detect multiple different antigens. Therefore, the venom concentrations measured and the pharmacokinetic parameters obtained in these studies are an averaged representation of the venom profile and may not capture the true characteristics or variability of the various snake toxins in the biological system. High-performance liquid chromatography with mass spectrometry would provide a more accurate measure of each toxin type, but such an approach to venom measurement has not been undertaken to date.

3.2. Pharmacokinetic Studies in Humans

In this meta-analysis, we were able to model literature-derived data to explore the pharmacokinetic parameters of snake venom for the first time. Timed venom concentration data from the literature were obtained from pre-antivenom samples described mainly in studies investigating patients who were intended to be administered an antivenom. The data extracted from the literature and modelled were best described by a one-compartment model. It is likely that venom may well display multicompartment pharmacokinetics in humans, but the available data that could be extracted were too sparse to support more than a one-compartment model. In this work, we also identified a bite-dependent exposure related to snakes of different families, indicating that the Viperidae family injects approximately a 75% higher dose. In the model, we noted a large between-subject variability in the relative dose injected (indicated as the between-subject variability in F1 from Table A1 of Appendix A). The variability between subjects in the remaining parameters, e.g., CL and V, was not wholly different from that of therapeutic drugs for which the normal variability in CL is approximately 50% in the case of small drug molecules [97]. In the model, the elimination half-life was calculated to be approximately 10 h. Of note, Audebert et al., 1994 also reported similar findings, showing a mono-exponential decline of snake venom concentration, with elimination half-lives ranging from 5 h to 12 h [75]. However, our results are in contrast to those of most animal studies in which two-compartment model fits are reported and hence the terminal half-life is often prolonged more than what observed here. Larger and richer venom concentration datasets in human envenoming are required to support multicompartment models and explore between-subject and between-snake variability, as well as other potential pharmacokinetic covariates.
When comparing our results to an allometrically scaled human half-life from animal data following IV administration of venom, a significant interspecies variability is observed. The allometric scaling of half-life is calculated as:
t 1 / 2   ( h u m a n ) = t 1 / 2 ( a n i m a l ) × ( W T h u m a n W T a n i m a l ) 0.25
where WT is the weight of the human (assumed 70 kg) or animal (assumed 0.31 kg for rats and 2.3 kg for rabbits) examined. Considering a terminal half-life of 5.3 h in rats (derived as the arithmetic average from the data in Table 3) and 16 h in rabbits (derived as the arithmetic average from the data in Table 3), the allometrically scaled half-life in humans was estimated to be 21 h and 39 h. The terminal phase half-lives calculated from rat and rabbit data were significantly longer compared to those in our study. The discrepancy in these results may be due to the sparse human data only supporting a one-compartment model compared to the two-compartment models for the rats and rabbits, considering that the half-life from the one-compartment model will be shorter than the terminal phase half-life from the two-compartment model. It is possible, therefore, that the terminal half-life in humans could be longer than one day.
Further analysis needs to be performed to account for these differences and to identify outstanding covariates that may inform the terminal phase of exposure.

4. Conclusions

While there is a limited number of studies that investigated the pharmacokinetics of snake venom in animals and reported venom concentrations in humans, this review yields important quantitative information which can be used as a basis for the future development of snake venom pharmacokinetic models. A meta-analysis of snake venom pharmacokinetics in humans has led to the development of a preliminary model, which was able to describe literature-derived data. Further exploration of the relationship between parameters and covariates is necessary to examine the factors that may influence the pharmacokinetic profiles of snake venom and construct a substantial model from data in a larger series of snake envenomation cases. It is surprising that a relatively meaningful analysis was possible, despite the limitations of ELISA, the fact that venom is a mixture of small and large molecules, and the differences in snakes. The development of improved methods to measure venom is required, as well as studies focusing on single snake species or on similar types of snakes.

5. Methods

5.1. Pharmacokinetic Studies in Animals

We reviewed pharmacokinetic studies of the administration of snake venom to animals and summarised relevant data.

5.1.1. Literature Search Strategy

EMBASE (1974–present) and Medline (1946–present) were used to identify relevant articles. Keywords and text terms included “exp envenomation”, “exp snakebite”, “exp snake venom”, “exp pharmacokinetics”, “exp toxicokinetics”. The search was limited to the English language and animals.

5.1.2. Inclusion and Exclusion Criteria

The publications were only included if they identified the venom or toxins administered, provided timed venom and toxin concentrations (in the absence of antivenom) or a compartmental analysis of their data measured by ELISA, and conducted a pharmacokinetic analysis. Studies that measured venom concentration using RIA were not included because of the possibility of an inaccurate determination of the pharmacokinetic parameters as a result of the metabolism or of the degradation of the radio-labelled proteins [98]. Review articles, commentaries, editorial papers, and conference abstracts were excluded.

5.1.3. Data Extraction and Creation of a Summary Table

The studies meeting all the inclusion criteria were grouped according to animal species. The pharmacokinetic parameters extracted from these publications included: half-life, bioavailability, clearance, and volume of distribution. A summary table was created showing the demographics of the animal studies, snake identification or venom or toxin type, type and number of animals, route of administration, dose administered, and pharmacokinetic parameters. When different units of doses or parameters were reported in the studies, they were converted to mcg·kg−1 for the dose administered, hours for the half-lives, L·kg−1 for the volume of distribution at steady state (Vss), and L·h−1·kg−1 for the systemic clearance (CL).

5.2. Pharmacokinetic Studies in Humans

We reviewed all pharmacokinetic studies of snake venom in humans and extracted relevant data for further analysis.

5.2.1. Literature Search Strategy

Both EMBASE (1974–present) and Medline (1946–present) databases were used for the search. Keywords and text words used for the search were “exp envenomation”, “exp snakebite”, “exp snake venom”, “exp venom antiserum”, “exp snake venom antiserum”, “exp antivenin”, “exp enzymelinked immunosorbent assay”, “exp enzyme immunoassay”, “exp immunoenzyme technique”.
The publications were only included for data extraction if they reported the identity of the snake, timed (post-bite) venom concentrations (in the absence of antivenom), and used ELISA to measure the snake venom concentrations. The review was limited to English language and humans. Review articles, commentaries, editorial papers, and conference abstracts were excluded.

5.2.2. Data Extraction and Synthesis

The snake venom concentration-time data post-envenomation were extracted by (1) obtaining the reported concentrations from the written text, or (2) digitising relevant data points from the concentration-time graphs published in the paper. The graphical data were digitised using WebPlotDigitizer v3.12 software (available at: http://arohatgi.info/WebPlotDigitizer/). All data points considered to have occurred after the administration of an antivenom were excluded.

5.2.3. Data Analysis

The extracted data were transferred to a Microsoft Excel spreadsheet and replotted as a concentration-time profile of different snake venoms from different studies. Since the data included repeated measures on many individuals, the data were modelled within a nonlinear mixed-effects modelling framework using NONMEM 7.2 (ICON Development Solutions, Ellicott City, MD, USA), and standard pharmacokinetic models were considered. The best model was determined using the likelihood ratio test, with degrees of freedom equal to the number of parameters different in successive models. The final model was evaluated using a prediction-corrected visual predictive check.

Acknowledgments

Suchaya Sanhajariya is supported by the University of Newcastle International Postgraduate Research Scholarship and the National Health and Medical Research Council Centre for Research Excellence grant 1110343. Geoffrey K. Isbister is supported by a National Health and Medical Research Council Senior Research Fellowship, ID1061041.

Author Contributions

Suchaya Sanhajariya, Stephen B. Duffull, and Geoffrey K. Isbister designed the research; Suchaya Sanhajariyaperformed the literature search and extracted the data; Suchaya Sanhajariya and Stephen B. Duffull analysed data; Suchaya Sanhajariya, Stephen B. Duffull, and Geoffrey K. Isbister wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table
Table A1. Parameter estimates from the pharmacokinetics models.
Table A1. Parameter estimates from the pharmacokinetics models.
ParametersParameter Estimates (RSE%)
Base ModelCovariate Model
CL (L/h)15.2 (11%)13.3 (14%)
V (L)215 (9%)184 (10%)
D1 (h)1 FIX1 FIX
F1 (Viperidae)1 FIX1 FIX
F1 (Elapidae)1 FIX0.569 (43%)
Between subject variability
CL (CV %)52.6% (30%)43.7% (52%)
V (CV %)15.5% (237%)29.8% (101%)
D1 (CV %)44.1% (21%)44.1% (17%)
F1 (CV %)286.5% (7%)275.4% (7%)
Residual error
Proportional error0.047 (25%)0.047 (25%)
Abbreviations: CL, clearance; V, apparent volume of distribution; D1, duration for central compartment; F1, bioavailability for central compartment; RSE, relative standard error; CV%, percentage coefficient of variation.

References

  1. Warrell, D.A. Guidelines for the Management of Snake-Bites; World Health Organization: Geneva, Switzerland, 2010; p. 162. [Google Scholar]
  2. Gutiérrez, J.M.; Theakston, R.D.G.; Warrell, D.A. Confronting the neglected problem of snake bite envenoming: The need for a global partnership. PLoS Med. 2006, 3, e150. [Google Scholar] [CrossRef] [PubMed]
  3. Kasturiratne, A.; Wickremasinghe, A.R.; de Silva, N.; Gunawardena, N.K.; Pathmeswaran, A.; Premaratna, R.; Savioli, L.; Lalloo, D.G.; de Silva, H.J. The global burden of snakebite: A literature analysis and modelling based on regional estimates of envenoming and deaths. PLoS Med. 2008, 5, e218. [Google Scholar] [CrossRef] [PubMed]
  4. Chippaux, J.P. Snake-bites: Appraisal of the global situation. Bull. World Health Organ. 1998, 76, 515–524. [Google Scholar] [PubMed]
  5. Warrell, D.A. Snake bite. Lancet 2010, 375, 77–88. [Google Scholar] [CrossRef]
  6. Warrell, D.A.; Gutiérrez, J.M.; Organization, W.H. Rabies and Envenomings: A Neglected Public Health Issue: Report of a Consultative Meeting; World Health Organization: Geneva, Switzerland, 2007. [Google Scholar]
  7. Tasoulis, T.; Isbister, K.G. A review and database of snake venom proteomes. Toxins 2017, 9, 290. [Google Scholar] [CrossRef] [PubMed]
  8. Fry, B.G.; Vidal, N.; van der Weerd, L.; Kochva, E.; Renjifo, C. Evolution and diversification of the toxicofera reptile venom system. J. Proteom. 2009, 72, 127–136. [Google Scholar] [CrossRef] [PubMed]
  9. Vidal, N. Colubroid systematics: Evidence for an early appearance of the venom apparatus followed by extensive evolutionary tinkering. J. Toxicol. Toxin Rev. 2002, 21, 21–41. [Google Scholar] [CrossRef]
  10. Bieber, A.L. Metal and nonprotein constituents in snake venoms. In Snake Venoms; Lee, C.-Y., Ed.; Springer: Berlin, Heidelberg, 1979; pp. 295–306. [Google Scholar]
  11. Bon, C. Pharmacokinetics of venom toxins and their modification by antivenom therapy. J. Toxicol. Toxin Rev. 2003, 22, 129–138. [Google Scholar] [CrossRef]
  12. Chippaux, J.P.; Williams, V.; White, J. Snake venom variability: Methods of study, results and interpretation. Toxicon 1991, 29, 1279–1303. [Google Scholar] [CrossRef]
  13. Isbister, G.K.; Brown, S.G.A.; Page, C.B.; McCoubrie, D.L.; Greene, S.L.; Buckley, N.A. Snakebite in australia: A practical approach to diagnosis and treatment. Med. J. Aust. 2013, 199, 763–768. [Google Scholar] [CrossRef] [PubMed]
  14. Mehta, S.R.; Sashindran, V.K. Clinical features and management of snake bite. Med. J. Armed Forces India 2002, 58, 247–249. [Google Scholar] [CrossRef]
  15. Gasanov, S.E.; Dagda, R.K.; Rael, E.D. Snake venom cytotoxins, phospholipase A2s, and Zn2+-dependent metalloproteinases: Mechanisms of action and pharmacological relevance. J. Clin. Toxicol. 2014, 4, 1000181. [Google Scholar] [CrossRef] [PubMed]
  16. Kini, R.M. Excitement ahead: Structure, function and mechanism of snake venom phospholipase A2 enzymes. Toxicon 2003, 42, 827–840. [Google Scholar] [CrossRef] [PubMed]
  17. Gutiérrez, J.M.; Rucavado, A. Snake venom metalloproteinases:Their role in the pathogenesis of local tissue damage. Biochimie 2000, 82, 841–850. [Google Scholar] [CrossRef]
  18. Isbister, G.K. Snakebite doesn’t cause disseminated intravascular coagulation: Coagulopathy and thrombotic microangiopathy in snake envenoming. Semin. Thromb. Hemost. 2010, 36, 444–451. [Google Scholar] [CrossRef] [PubMed]
  19. Ranawaka, U.K.; Lalloo, D.G.; de Silva, H.J. Neurotoxicity in snakebite—The limits of our knowledge. PLoS Negl. Trop. Dis. 2013, 7, e2302. [Google Scholar] [CrossRef] [PubMed]
  20. Theakston, R.D.; Laing, G.D. Diagnosis of snakebite and the importance of immunological tests in venom research. Toxins 2014, 6, 1667–1695. [Google Scholar] [CrossRef] [PubMed]
  21. Theakston, R.D.; Lloyd-Jones, M.J.; Reid, H.A. Micro-elisa for detecting and assaying snake venom and venom-antibody. Lancet 1977, 2, 639–641. [Google Scholar] [CrossRef]
  22. Theakston, R.D.G. The application of immunoassay techniques, including enzyme-linked immunosorbent assay (elisa), to snake venom research. Toxicon 1983, 21, 341–352. [Google Scholar] [CrossRef]
  23. Dhananjaya, B.L.; Menon, J.C.; Joseph, J.K.; Raveendran, D.K.; Oommen, O.V. Snake venom detection kit (svdk): Update on current aspects and challenges. In Clinical Toxinology in Asia Pacific and Africa; Gopalakrishnakone, P., Faiz, A., Fernando, R., Gnanathasan, C.A., Habib, A.G., Yang, C.-C., Eds.; Springer: Dordrecht, The Netherlands, 2015; pp. 379–400. [Google Scholar]
  24. Selvanayagam, Z.E.; Gopalakrishnakone, P. Tests for detection of snake venoms, toxins and venom antibodies: Review on recent trends (1987–1997). Toxicon 1999, 37, 565–586. [Google Scholar] [CrossRef]
  25. Allen, G.E.; Brown, S.G.A.; Buckley, N.A.; O’Leary, M.A.; Page, C.B.; Currie, B.J.; White, J.; Isbister, G.K. Clinical effects and antivenom dosing in brown snake (pseudonaja spp.) envenoming—Australian snakebite project (asp-14). PLoS ONE 2012, 7. [Google Scholar] [CrossRef] [PubMed]
  26. Audebert, F.; Grosselet, O.; Sabouraud, A.; Bon, C. Quantitation of venom antigens from european vipers in human serum or urine by elisa. J. Anal. Toxicol. 1993, 17, 236–240. [Google Scholar] [CrossRef] [PubMed]
  27. Guo, M.P.; Wang, Q.C.; Liu, G.F. Pharmacokinetics of cytotoxin from chinese cobra (naja naja atra) venom. Toxicon 1993, 31, 339–343. [Google Scholar] [CrossRef]
  28. Hanvivatvong, O.; Mahasandana, S.; Karnchanachetanee, C. Kinetic study of russell’s viper venom in envenomed patients. Am. J. Trop. Med. Hyg. 1997, 57, 605–609. [Google Scholar] [CrossRef] [PubMed]
  29. Khin Ohn, L.; Aye Aye, M.; Tun, P.; Theingie, N.; Min, N. Russell’s viper venom levels in serum of snake bite victims in burma. Trans. R. Soc. Trop. Med. Hyg. 1984, 78, 165–168. [Google Scholar]
  30. Sano-Martins, I.S.; Tomy, S.C.; Campolina, D.; Dias, M.B.; de Castro, S.C.; de Sousa-e-Silva, M.C.; Amaral, C.F.; Rezende, N.A.; Kamiguti, A.S.; Warrell, D.A.; et al. Coagulopathy following lethal and non-lethal envenoming of humans by the south american rattlesnake (crotalus durissus) in brazil. QJM 2001, 94, 551–559. [Google Scholar] [CrossRef] [PubMed]
  31. Barral-Netto, M.; Schriefer, A.; Barral, A.; Almeida, A.R.P.; Mangabeira, A. Serum levels of bothropic venom in patients without antivenom intervention. Am. J. Trop. Med. Hyg. 1991, 45, 751–754. [Google Scholar] [CrossRef] [PubMed]
  32. Kulawickrama, S.; O’Leary, M.A.; Hodgson, W.C.; Brown, S.G.; Jacoby, T.; Davern, K.; Isbister, G.K. Development of a sensitive enzyme immunoassay for measuring taipan venom in serum. Toxicon 2010, 55, 1510–1518. [Google Scholar] [CrossRef] [PubMed]
  33. O’Leary, M.A.; Isbister, G.K.; Schneider, J.J.; Brown, S.G.A.; Currie, B.J. Enzyme immunoassays in brown snake (pseudonaja spp.) envenoming: Detecting venom, antivenom and venom-antivenom complexes. Toxicon 2006, 48, 4–11. [Google Scholar] [CrossRef] [PubMed]
  34. Sjostrom, L.; Karlson-Stiber, C.; Persson, H.; Al-Abdulla, I.H.; Smith, D.C. Development and clinical application of immunoassays for european adder (vipera berus berus) venom and antivenom. Toxicon 1996, 34, 91–98. [Google Scholar] [CrossRef]
  35. Audebert, F.; Urtizberea, M.; Sabouraud, A.; Scherrmann, J.M.; Bon, C. Pharmacokinetics of vipera aspis venom after experimental envenomation in rabbits. J. Pharmacol. Exp. Ther. 1994, 268, 1512–1517. [Google Scholar] [PubMed]
  36. Minton, S.A. Present tests for detection of snake venom: Clinical applications. Ann. Emerg. Med. 1987, 16, 932–937. [Google Scholar] [CrossRef]
  37. Toyama, M.H.; Soares, A.M.; Wen-Hwa, L.; Polikarpov, I.; Giglio, J.R.; Marangoni, S. Amino acid sequence of piratoxin-ii, a myotoxic lys49 phospholipase A2 homologue from bothrops pirajai venom. Biochimie 2000, 82, 245–250. [Google Scholar] [CrossRef]
  38. Burke, J.E.; Dennis, E.A. Phospholipase A2 biochemistry. Cardiovasc. Drugs Ther. 2009, 23, 49–59. [Google Scholar] [CrossRef] [PubMed]
  39. Arni, R.K.; Ward, R.J. Phospholipase A2—A structural review. Toxicon 1996, 34, 827–841. [Google Scholar] [CrossRef]
  40. Bjarnason, J.B.; Fox, J.W. Hemorrhagic metalloproteinases from snake venoms. Pharmacol. Ther. 1994, 62, 325–372. [Google Scholar] [CrossRef]
  41. Takeda, S.; Takeya, H.; Iwanaga, S. Snake venom metalloproteinases: Structure, function and relevance to the mammalian adam/adamts family proteins. Biochim. Biophys. Acta Proteins Proteom. 2012, 1824, 164–176. [Google Scholar] [CrossRef] [PubMed]
  42. Serrano, S.M.T.; Maroun, R.C. Snake venom serine proteinases: Sequence homology vs. Substrate specificity, a paradox to be solved. Toxicon 2005, 45, 1115–1132. [Google Scholar] [CrossRef] [PubMed]
  43. Matsui, T.; Fujimura, Y.; Titani, K. Snake venom proteases affecting hemostasis and thrombosis. Biochim. Biophys. Acta Protein Struct. Mol. Enzymol. 2000, 1477, 146–156. [Google Scholar] [CrossRef]
  44. Du, X.-Y.; Clemetson, K.J. Snake venom l-amino acid oxidases. Toxicon 2002, 40, 659–665. [Google Scholar] [CrossRef]
  45. Fox, J.W. A brief review of the scientific history of several lesser-known snake venom proteins: L-amino acid oxidases, hyaluronidases and phosphodiesterases. Toxicon 2013, 62, 75–82. [Google Scholar] [CrossRef] [PubMed]
  46. Guo, C.; Liu, S.; Yao, Y.; Zhang, Q.; Sun, M.-Z. Past decade study of snake venom l-amino acid oxidase. Toxicon 2012, 60, 302–311. [Google Scholar] [CrossRef] [PubMed]
  47. Mackessy, S.P. The field of reptile toxinology. In Handbook of Venoms and Toxins of Reptiles; CRC Press: Boca Raton, FL, USA, 2009; pp. 3–23. [Google Scholar]
  48. Iwanaga, S.; Suzuki, T. Enzymes in snake venom. In Snake Venoms; Lee, C.-Y., Ed.; Springer: Berlin/Heidelberg, Germany, 1979; pp. 61–158. [Google Scholar]
  49. Ouyang, C.; Huang, T.-F. Inhibition of platelet aggregation by 5′-nucleotidase purified from trimeresurus gramineus snake venom. Toxicon 1983, 21, 491–501. [Google Scholar] [CrossRef]
  50. Aird, S.D. Ophidian envenomation strategies and the role of purines. Toxicon 2002, 40, 335–393. [Google Scholar] [CrossRef]
  51. Kang, T.S.; Georgieva, D.; Genov, N.; Murakami, M.T.; Sinha, M.; Kumar, R.P.; Kaur, P.; Kumar, S.; Dey, S.; Sharma, S.; et al. Enzymatic toxins from snake venom: Structural characterization and mechanism of catalysis. FEBS J. 2011, 278, 4544–4576. [Google Scholar] [CrossRef] [PubMed]
  52. Frobert, Y.; Créminon, C.; Cousin, X.; Rémy, M.-H.; Chatel, J.-M.; Bon, S.; Bon, C.; Grassi, J. Acetylcholinesterases from elapidae snake venoms: Biochemical, immunological and enzymatic characterization. Biochim. Biophys. Acta Protein Struct. Mol. Enzymol. 1997, 1339, 253–267. [Google Scholar] [CrossRef]
  53. Kini, R.M.; Doley, R. Structure, function and evolution of three-finger toxins: Mini proteins with multiple targets. Toxicon 2010, 56, 855–867. [Google Scholar] [CrossRef] [PubMed]
  54. Nirthanan, S.; Gwee, M.C.E. Three-finger α-neurotoxins and the nicotinic acetylcholine receptor, forty years on. J. Pharmacol. Sci. 2004, 94, 1–17. [Google Scholar] [CrossRef] [PubMed]
  55. Willmott, N.; Gaffney, P.; Masci, P.; Whitaker, A. A novel serine protease inhibitor from the australian brown snake, pseudonaja textilis textilis: Inhibition kinetics. Fibrinolysis 1995, 9, 1–8. [Google Scholar] [CrossRef]
  56. Mukherjee, A.K.; Mackessy, S.P.; Dutta, S. Characterization of a kunitz-type protease inhibitor peptide (rusvikunin) purified from daboia russelii russelii venom. Int. J. Biol. Macromol. 2014, 67, 154–162. [Google Scholar] [CrossRef] [PubMed]
  57. Yamazaki, Y.; Morita, T. Structure and function of snake venom cysteine-rich secretory proteins. Toxicon 2004, 44, 227–231. [Google Scholar] [CrossRef] [PubMed]
  58. Morita, T. Structures and functions of snake venom clps (c-type lectin-like proteins) with anticoagulant-, procoagulant-, and platelet-modulating activities. Toxicon 2005, 45, 1099–1114. [Google Scholar] [CrossRef] [PubMed]
  59. Ogawa, T.; Chijiwa, T.; Oda-Ueda, N.; Ohno, M. Molecular diversity and accelerated evolution of c-type lectin-like proteins from snake venom. Toxicon 2005, 45, 1–14. [Google Scholar] [CrossRef] [PubMed]
  60. Kamiguti, A.S.; Zuzel, M.; Theakston, R.D.G. Snake venom metalloproteinases and disintegrins: Interactions with cells. Braz. J. Med. Biol. Res. 1998, 31, 853–862. [Google Scholar] [CrossRef] [PubMed]
  61. St Pierre, L.; Flight, S.; Masci, P.P.; Hanchard, K.J.; Lewis, R.J.; Alewood, P.F.; de Jersey, J.; Lavin, M.F. Cloning and characterisation of natriuretic peptides from the venom glands of australian elapids. Biochimie 2006, 88, 1923–1931. [Google Scholar] [CrossRef] [PubMed]
  62. Levin, E.R.; Gardner, D.G.; Samson, W.K. Natriuretic peptides. N. Engl. J. Med. 1998, 339, 321–328. [Google Scholar] [PubMed]
  63. Gutierrez, J.M.; Leon, G.; Lomonte, B.; Angulo, Y. Antivenoms for snakebite envenomings. Inflamm. Allergy Drug Targets 2011, 10, 369–380. [Google Scholar] [CrossRef] [PubMed]
  64. Maduwage, K.; Silva, A.; O’Leary, M.A.; Hodgson, W.C.; Isbister, G.K. Efficacy of indian polyvalent snake antivenoms against sri lankan snake venoms: Lethality studies or clinically focussed in vitro studies. Sci. Rep. 2016, 6, 26778. [Google Scholar] [CrossRef] [PubMed]
  65. Hart, A.J.; Hodgson, W.C.; O’Leary, M.; Isbister, G.K. Pharmacokinetics and pharmacodynamics of the myotoxic venom of pseudechis australis (mulga snake) in the anesthetised rat. Clin. Toxicol. 2014, 52, 604–610. [Google Scholar] [CrossRef] [PubMed]
  66. Paniagua, D.; Jimenez, L.; Romero, C.; Vergara, I.; Calderon, A.; Benard, M.; Bernas, M.J.; Rilo, H.; De Roodt, A.; D’Suze, G.; et al. Lymphatic route of transport and pharmacokinetics of micrurus fulvius (coral snake) venom in sheep. Lymphology 2012, 45, 144–153. [Google Scholar] [PubMed]
  67. Tan, C.H.; Sim, S.M.; Gnanathasan, C.A.; Fung, S.Y.; Tan, N.H. Pharmacokinetics of the sri lankan hump-nosed pit viper (hypnale hypnale) venom following intravenous and intramuscular injections of the venom into rabbits. Toxicon 2014, 79, 37–44. [Google Scholar] [CrossRef] [PubMed]
  68. Zhao, H.; Zheng, J.; Jiang, Z. Pharmacokinetics of thrombin-like enzyme from venom of agkistrodon halys ussuriensis emelianov determined by elisa in the rat. Toxicon 2001, 39, 1821–1826. [Google Scholar] [CrossRef]
  69. Nakamura, M.; Kinjoh, K.; Miyagi, C.; Oka, U.; Sunagawa, M.; Yamashita, S.; Kosugi, T. Pharmacokinetics of habutobin in rabbits. Toxicon 1995, 33, 1201–1206. [Google Scholar] [CrossRef]
  70. Yap, M.K.; Tan, N.H.; Sim, S.M.; Fung, S.Y.; Tan, C.H. Pharmacokinetics of naja sumatrana (equatorial spitting cobra) venom and its major toxins in experimentally envenomed rabbits. PLoS Negl. Trop. Dis. 2014, 8, e2890. [Google Scholar] [CrossRef] [PubMed]
  71. Mello, S.M.; Linardi, A.; Renno, A.L.; Tarsitano, C.A.B.; Pereira, E.M.; Hyslop, S. Renal kinetics of bothrops alternatus (urutu) snake venom in rats. Toxicon 2010, 55, 470–480. [Google Scholar] [CrossRef] [PubMed]
  72. Riviere, G.; Choumet, V.; Saliou, B.; Debray, M.; Bon, C. Absorption and elimination of viper venom after antivenom administration. J. Pharmacol. Exp. Ther. 1998, 285, 490–495. [Google Scholar] [PubMed]
  73. Yap, M.K.K.; Tan, N.H.; Sim, S.M.; Fung, S.Y. Toxicokinetics of naja sputatrix (javan spitting cobra) venom following intramuscular and intravenous administrations of the venom into rabbits. Toxicon 2013, 68, 18–23. [Google Scholar] [CrossRef] [PubMed]
  74. Sim, S.M.; Saremi, K.; Tan, N.H.; Fung, S.Y. Pharmacokinetics of cryptelytrops purpureomaculatus (mangrove pit viper) venom following intravenous and intramuscular injections in rabbits. Int. Immunopharmacol. 2013, 17, 997–1001. [Google Scholar] [CrossRef] [PubMed]
  75. Audebert, F.; Sorkine, M.; Robbe-Vincent, A.; Bon, C. Viper bites in france: Clinical and biological evaluation; kinetics of envenomations. Hum. Exp. Toxicol. 1994, 13, 683–688. [Google Scholar] [CrossRef] [PubMed]
  76. Barraviera, B.; Sartori, A.; Pereira Da Silva, M.F.; Kaneno, R.; Peracoli, M.T.S. Use of an elisa assay to evaluate venom, antivenom, igg and igm human antibody levels in serum and cerebrospinal fluid from patients bitten by crotalus durissus terrificus in brazil. J. Venom. Anim. Toxins 1996, 2, 14–27. [Google Scholar] [CrossRef]
  77. Brvar, M.; Kurtovic, T.; Grenc, D.; Lang Balija, M.; Krizaj, I.; Halassy, B. Vipera ammodytes bites treated with antivenom viperatab: A case series with pharmacokinetic evaluation. Clin. Toxicol. 2017, 55, 241–248. [Google Scholar] [CrossRef] [PubMed]
  78. Bucaretchi, F.; Hyslop, S.; Mello, S.M.; Vieira, R.J. Bothrops snakebite on the head: Case report and review of the literature. Ann. Trop. Med. Parasitol. 2007, 101, 733–743. [Google Scholar] [CrossRef] [PubMed]
  79. Bucaretchi, F.; Capitani, E.M.D.; Hyslop, S.; Mello, S.M.; Madureira, P.R.; Zanardi, V.; Ferreira, D.M.; Meirelles, G.V.; Fernandes, L.C.R. Compartment syndrome after bothrops jararaca snakebite: Monitoring, treatment, and outcome. Clin. Toxicol. 2010, 48, 57–60. [Google Scholar] [CrossRef] [PubMed]
  80. Bucaretchi, F.; De Capitani, E.M.; Hyslop, S.; Mello, S.M.; Fernandes, C.B.; Bergo, F.; Nascimento, F.B.P. Compartment syndrome after south american rattlesnake (crotalus durissus terrificus) envenomation. Clin. Toxicol. 2014, 52, 639–641. [Google Scholar] [CrossRef] [PubMed]
  81. Bucher, B.; Canonge, D.; Robbe-Vincent, A.; Choumet, V.; Bon, C.; Lang, J.; Biao, T.; Brunod, M.; Garnier, D.; Guerin, J.F.; et al. Clinical indicators of envenoming and serum levels of venom antigens in patients bitten by bothrops lanceolatus in martinique. Trans. R. Soc. Trop. Med. Hyg. 1997, 91, 186–190. [Google Scholar] [CrossRef]
  82. Hung, D.Z.; Liau, M.Y.; Lin-Shiau, S.Y. The clinical significance of venom detection in patients of cobra snakebite. Toxicon 2003, 41, 409–415. [Google Scholar] [CrossRef]
  83. Hung, D.Z.; Yu, Y.J.; Hsu, C.L.; Lin, T.J. Antivenom treatment and renal dysfunction in russell’s viper snakebite in taiwan: A case series. Trans. R. Soc. Trop. Med. Hyg. 2006, 100, 489–494. [Google Scholar] [CrossRef] [PubMed]
  84. Isbister, G.K.; O’Leary, M.A.; Schneider, J.J.; Brown, S.G.A.; Currie, B.J. Efficacy of antivenom against the procoagulant effect of australian brown snake (pseudonaja sp.) venom: In vivo and in vitro studies. Toxicon 2007, 49, 57–67. [Google Scholar] [CrossRef] [PubMed]
  85. Isbister, G.K.; Halkidis, L.; O’Leary, M.A.; Whitaker, R.; Cullen, P.; Mulcahy, R.; Bonnin, R.; Brown, S.G.A. Human anti-snake venom igg antibodies in a previously bitten snake-handler, but no protection against local envenoming. Toxicon 2010, 55, 646–649. [Google Scholar] [CrossRef] [PubMed]
  86. Isbister, G.K.; Gault, A.; Tasoulis, T.; O’Leary, M.A. A definite bite by the ornamental snake (denisonia maculata) causing mild envenoming. Clin. Toxicol. 2016, 54, 241–244. [Google Scholar] [CrossRef] [PubMed]
  87. Lim, A.Y.L.; Singh, P.N.; Isbister, G.K. Severe rhabdomyolysis from red-bellied black snake (pseudechis porphyriacus) envenoming despite antivenom. Toxicon 2016, 117, 46–48. [Google Scholar] [CrossRef] [PubMed]
  88. McNally, T.; Conway, G.S.; Jackson, L.; Theakston, R.D.G.; Marsh, N.A.; Warrell, D.A.; Young, L.; Mackie, I.J.; Machin, S.J. Accidental envenoming by a gaboon viper (bitis gabonica): The haemostatic disturbances observed and investigation of in vitro haemostatic properties of whole venom. Trans. R. Soc. Trop. Med. Hyg. 1993, 87, 66–70. [Google Scholar] [CrossRef]
  89. Phillips, R.E.; Theakston, R.D.G.; Warrell, D.A.; Galigedara, Y.; Abeysekera, D.T.D.J.; Dissanayaka, P.; Hutton, R.A.; Aloysius, D.J. Paralysis, rhabdomyolysis and haemolysis caused by bites of russell’s viper (vipera russelli pulchella) in sri lanka: Failure of indian (haffkine) antivenom. Q. J. Med. 1988, 68, 691–716. [Google Scholar] [PubMed]
  90. Schneemann, M.; Cathomas, R.; Laidlaw, S.T.; El Nahas, A.M.; Theakston, R.D.G.; Warrell, D.A. Life-threatening envenoming by the saharan horned viper (cerastes cerastes) causing micro-angiopathic haemolysis, coagulopathy and acute renal failure: Clinical cases and review. QJM Mon. J. Assoc. Phys. 2004, 97, 717–727. [Google Scholar] [CrossRef] [PubMed]
  91. Silva, A.; Maduwage, K.; Sedgwick, M.; Pilapitiya, S.; Weerawansa, P.; Dahanayaka, N.J.; Buckley, N.A.; Johnston, C.; Siribaddana, S.; Isbister, G.K. Neuromuscular effects of common krait (bungarus caeruleus) envenoming in sri lanka. PLoS Negl. Trop. Dis. 2016, 10. [Google Scholar] [CrossRef] [PubMed]
  92. Tun, P.; Aye Aye, M.; Warrell, D.A.; Tin, M. King cobra (ophiophagus hannah) bites in myanmar: Venom antigen levels and development of venom antibodies. Toxicon 1995, 33, 379–382. [Google Scholar] [CrossRef]
  93. Warrell, D.A.; Hudson, B.J.; Lalloo, D.G.; Trevett, A.J.; Whitehead, P.; Bamler, P.R.; Ranaivoson, M.; Wiyono, A.; Richie, T.L.; Fryauff, D.J.; et al. The emerging syndrome of envenoming by the new guinea small-eyed snake micropechis ikaheka. QJM Mon. J. Assoc. Phys. 1996, 89, 523–530. [Google Scholar] [CrossRef]
  94. Ismail, M.; Abd-Elsalam, M.A.; Al-Ahaidib, M.S. Pharmacokinetics of 125i-labelled walterinnesia aegyptia venom and its distribution of the venom and its toxin versus slow absorption and distribution of igg, f(ab’)2 and f(ab) of the antivenin. Toxicon 1998, 36, 93–114. [Google Scholar] [CrossRef]
  95. Ismail, M.; Aly, M.H.M.; Abd-Elsalam, M.A.; Morad, A.M. A three-compartment open pharmacokinetic model can explain variable toxicities of cobra venoms and their alpha toxins. Toxicon 1996, 34, 1011–1026. [Google Scholar] [CrossRef]
  96. Maung Maung, T.; Khin Mee, M.; Mi Mi, K.; Thein, T. Kinetics of envenomation with russell’s viper (vipera russelli) venom and of antivenom use in mice. Toxicon 1988, 26, 373–378. [Google Scholar] [CrossRef]
  97. Al-Sallami, H.S.; Cheah, S.L.; Han, S.Y.; Liew, J.; Lim, J.; Ng, M.A.; Solanki, H.; Soo, R.J.; Tan, V.; Duffull, S.B. Between-subject variability: Should high be the new normal? Eur. J. Clin. Pharmacol. 2014, 70, 1403–1404. [Google Scholar] [CrossRef] [PubMed]
  98. Krifi, M.N.; Savin, S.; Debray, M.; Bon, C.; El Ayeb, M.; Choumet, V. Pharmacokinetic studies of scorpion venom before and after antivenom immunotherapy. Toxicon 2005, 45, 187–198. [Google Scholar] [CrossRef] [PubMed]
Toxins 10 00073 g001 550
Figure 1. Schematic diagram of article identification (animal studies).
Figure 1. Schematic diagram of article identification (animal studies).
Toxins 10 00073 g001
Toxins 10 00073 g002 550
Figure 2. Schematic diagram of article identification (human studies).
Figure 2. Schematic diagram of article identification (human studies).
Toxins 10 00073 g002
Toxins 10 00073 g003 550
Figure 3. Snake venom concentrations versus time post-bite in 145 patients, extracted from 24 studies. The closed triangles represent data of patients bitten by snakes of the Elapidae family, and the closed circles represents data of patients bitten by snakes of the Viperidae family. The lines are used to join samples from the same patient.
Figure 3. Snake venom concentrations versus time post-bite in 145 patients, extracted from 24 studies. The closed triangles represent data of patients bitten by snakes of the Elapidae family, and the closed circles represents data of patients bitten by snakes of the Viperidae family. The lines are used to join samples from the same patient.
Toxins 10 00073 g003
Toxins 10 00073 g004a 550Toxins 10 00073 g004b 550
Figure 4. Snake venom concentration versus time data post-bite of patients bitten by snakes of the Elapidae family (a) and the Viperidae family (b).
Figure 4. Snake venom concentration versus time data post-bite of patients bitten by snakes of the Elapidae family (a) and the Viperidae family (b).
Toxins 10 00073 g004aToxins 10 00073 g004b
Toxins 10 00073 g005 550
Figure 5. Simulation-based diagnostic plots for snake venom concentration (mcg/L) versus time post-bite (h) for the pharmacokinetic model. The solid line represents the median observed snake venom concentration, and the dashed lines represent the observed 10% and 90% percentile. The shaded areas are the corresponding 95% confidence intervals of the percentiles predicted by the model.
Figure 5. Simulation-based diagnostic plots for snake venom concentration (mcg/L) versus time post-bite (h) for the pharmacokinetic model. The solid line represents the median observed snake venom concentration, and the dashed lines represent the observed 10% and 90% percentile. The shaded areas are the corresponding 95% confidence intervals of the percentiles predicted by the model.
Toxins 10 00073 g005
Table
Table 1. List of common venom components (enzymatic) found in snakes of the Elapidae and Viperidae families, their size and activities.
Table 1. List of common venom components (enzymatic) found in snakes of the Elapidae and Viperidae families, their size and activities.
Enzymatic ComponentsApproximate Molecular Mass (kDa)Mechanism of ActionExamples of Biological EffectsSnake Families
Phospholipase A2 (PLA2)12–14 [37]Hydrolyses the ester bond at sn-2 position of phospholipids producing free fatty acids and lysophospholipid. Toxic effects can result from this enzymatic action or may be the results of non-enzymatic activity [38].Myotoxicity, oedema formation, anticoagulant effects, hypotension, presynaptic neurotoxicity [16,39]Elapidae (type I PLA2); Viperidae (type II PLA2) [7]
Snake-venom metalloproteases (SVMP)Classified into 3 groups based on domain organisation [40]; P-I: 20–30; P-II: 30–60; P-III: 60–100Proteolytic activities leading to degradation of protein structures e.g., basement membranes of blood vessel and extracellular matrix components [17,41]. Disintegrin-like domain of SVMP may also contribute to the haemorrhagic effects [17].Induce local and systemic bleeding and disrupt haemostasis through its pro-/anticoagulation properties. Extravasation of blood, inflammation and tissue necrosis [17,41]Major protein family in viper venoms, but less abundant in elapid venom [7]
Serine proteases (SVSP) e.g., thrombin-like enzymes26–67 [42]Hydrolyse peptide bonds mainly in pro-enzymes in the coagulation cascade, causing procoagulant, fibrinolytic and/or fibrinogenolytic activities. Some SVSPs have kallikrein-like activity leading to release of bradykinin [42,43].Disruption of haemostasis and hypotension [42]Almost all Viperidae, uncommon in Elapidae except Australian snakes [7]
l-amino acid oxidases (LAAO)110–150 when measured by gel-filtration method under non-denaturing conditions; 50–70 when measured by SDS/PAGE method under reducing and non-reducing conditions [44]Catalyse stereospecific oxidative deamination of l-amino acid, resulting in production of α-keto acid, ammonia and hydrogen peroxide [45].Effects on platelet aggregation, inducing cell apoptosis, and antimicrobial activities [46]Both Elapidae and Viperidae. Most common in Crotalinae [7].
5′-Nucleotidases53–82 [47]Hydrolyse phosphate monoester linked to 5′-position of DNA and RNA [48].Platelet aggregation inhibition [49,50]Both Elapidae and Viperidae [48]
Acetylcholinesterases55–60 [47]Hydrolyse acetylcholine to choline and acetate group [51].Termination of neurotransmission by acetylcholine [51,52]Elapidae except Dendroaspis genus [52]
Hyaluronidases33–110 [45]Hydrolyse hyaluronan into oligosaccharides and N-acetylglucosamine [45].“Spreading factor” alters the structural, rheological, and chemical properties of the extracellular matrix [45]Both Elapidae and Viperidae [48]
Abbreviations: PLA2, phospholipase A2; SVMP, snake venom metalloproteases; SVSP, snake-venom serine proteases; LAAO, l-amino acid oxidases.
Table
Table 2. List of common venom components (non-enzymatic) found in snakes of the Elapidae and Viperidae families, their size and activities.
Table 2. List of common venom components (non-enzymatic) found in snakes of the Elapidae and Viperidae families, their size and activities.
Non-Enzymatic ComponentsApproximate Molecular Mass (kDa)Mechanism of ActionExamples of Biological EffectsSnake Families
Three-finger toxins (3FTx) e.g., α-neurotoxins6–9 [47]Inhibit postsynaptic nicotinic acetylcholine receptors in neuromuscular junction and interfere with neuromuscular transmission [53,54]. Other activities include cardiotoxins, L-type calcium channel blockage, inhibition of platelet aggregation [53].Postsynaptic neurotoxicityElapidae and very rare in Viperidae [7]
Kunitz peptides (KUN)7 [55,56]Inhibit serine protease (e.g., trypsin and plasmin) activities, interfering with blood coagulation and fibrinolysis [55,56]. Other activities include ion channel blockade and inflammation [56].Disruption of haemostasisElapidae and Viperinae (absent in Crotalinae) [7]
Cysteine-rich secretory proteins (CRiSP)20–30 [57]L-type calcium and cyclic nucleotide-gated (CNG) channel blockade [57].Inhibit smooth muscle contraction [57]More common and abundant in Viperidae [7]
C-type lectins (CTL)Composed of two subunits [58]; α (A chain): 14–15; β (B chain): 13–14Bind to, inhibit, or activate specific platelet membrane receptors or blood coagulation factors [59].Anticoagulation, promote or inhibit platelet aggregation [59]More abundant in Viperidae [7]
Disintegrins (DIS)5–10 [60]Bind to glycoprotein IIb/IIIa (αIIbβ3 integrin) expressed on activated platelet to prevent interaction with fibrinogen [60].Inhibit platelet aggregation [60]Viperidae, absent in Elapidae [7]
Natriuretic peptides (NP)3.5–4 [61]Interaction between NPs and guanylyl cyclase receptors leads to an increase in cGMP levels and subsequent signalling cascade [62]. NPs can also affect renin-angiotensin system by inhibiting angiotensin-converting enzyme [61].Vasodilation, diuresis, and natriuresis leading to hypotension, and promote sodium and water excretion [62]Both Elapidae (atrial-type and brain-type) and Viperidae (C-type) [5]. More common and abundant in Viperidae than Elapidae [7]
Abbreviations: 3FTx, three-finger toxins; KUN, Kunitz peptides; CRiSP, Cysteine-rich secretory proteins; CTL, C-type lectins; DIS, disintegrins; NP, natriuretic peptides.
Table
Table 3. Pharmacokinetic parameters of snake venom and toxins following intravenous (IV) injection expressed as mean (±SD).
Table 3. Pharmacokinetic parameters of snake venom and toxins following intravenous (IV) injection expressed as mean (±SD).
Snake SpeciesAnimal ModelNo.Dose (mcg·kg−1)t1/2α (h)t1/2β (h)Vss (L·kg−1)CL (L·h−1·kg−1)Ref
Isolated Toxin
Rats
Agkistrodon halys ussuriensis EmelianovSprague-Dawley rats of either sex (180–200 g)550 (thrombin-like enzyme)0.3 (±0.12 *)3.9 (±1.63 *)1.8 (±1.03 *)0.324 (±0.067 *)[68]
Rabbits
Naja naja atraNew Zealand rabbits of either sex (1.82 ± 0.09 kg)6200 (cytotoxin)0.097 (±0.01)3.5 (±0.2)1.7 a (±0.3)0.185 b[27]
Trimeresurus flavoviridisJapanese white rabbits (3.2–4.4 kg)550 (habutobin)0.074 (±0.021 **)0.84 (±0.13 **)0.031 b,c, 0.021 b,d0.061 b[69]
Naja sumatranaNew Zealand white rabbits (approx. 2 kg)350 (PLA2)0.7 (±0.03)11.7 (±0.8)0.25 b,c, 0.45 b,d0.048 b[70]
50 (neurotoxin)0.5 (±0.1)8.8 (±0.9)0.45 b,c, 0.5 b,d0.082 b
50 (cardiotoxin)0.6 (±0.1)8.6 (±0.1)0.5 b,c, 0.55 b,d0.087 b
100 (cardiotoxin in whole venom)0.5 (±0.01)11.0 (±0.2)0.4 b,c, 0.5 b,d0.060 b
Whole Venom
Rats
Bothrops alternatusMale Wistar rats (200–250 g)68000.38 (±0.03)12.1 (±6.4)0.50 (±0.12)0.033 (±0.011)[71]
Pseudechis australisMale Sprague-Dawley rats (320–420 g)8100-0.27 ***--[65]
Rabbits
Vipera aspisCharles de Bouscat HY rabbits (2.5–3 kg)52500.71 (±0.2 *)12 (±2.24 *)1.2 (±0.089 *)0.084 (±0.013 *)[35]
Vipera aspisNew Zealand rabbits (2.75–3 kg)52500.53 (±0.31 *)14.2 (±2.68 *)0.7 (±0.11 *)0.040 (±0.002 *)[72]
Naja sputatrixNew Zealand rabbits (2 kg)3900.5 (±0.3)15.4 (±2.5)0.8 b0.034 b[73]
Cryptelytrops purpureomaculatusMale New Zealand rabbits (1.7–2.1 kg)32000.25 (±0.01)27.7 (±0.0)0.39 c (±0.01), 1.80 d (±0.11)0.055 (±0.003)[74]
Naja sumatranaNew Zealand white rabbits (approx. 2 kg)31000.8 (±0.3)13.6 (±1.1)0.5 b,c, 0.4 b,d0.046 b[70]
Hypnale hypnaleNew Zealand white rabbits (1.95 ± 0.05 kg)3100.8 (±0.17 *)19.3 (±3.29 *)0.13 b0.007 (±0.001 *)[67]
Sheep
Micrurus fulviusSheep (36–60 kg)41000 mcg-0.42 (±0.11 *)0.12 b0.093 b[66]
* SD calculated from standard errors reported in the original papers using the following equation: S D = S E   ×   √ n ; ** Assumed to be reported as SD; *** SD not reported; a = Apparent volume of distribution; b = Unit by weight (·kg−1) calculated on the basis of the mean of the reported animal weight; c = Apparent volume of central compartment by weight (Vc); d = Apparent volume of peripheral compartment by weight (Vp); Abbreviations: t1/2α, half-life of the distribution phase; t1/2β, half-life of the elimination phase; Vss, apparent volume of distribution at steady state; CL, systemic clearance.
Table
Table 4. Pharmacokinetic parameters of snake venom and toxins following intramuscular (IM) injection expressed as mean (±SD).
Table 4. Pharmacokinetic parameters of snake venom and toxins following intramuscular (IM) injection expressed as mean (±SD).
Snake SpeciesAnimal ModelNo.Dose (mcg·kg−1)t1/2ka (h)F (%)t1/2α (h)t1/2β (h)Vss (L·kg−1)CL (L·h−1·kg−1)Ref
Isolated Toxin
Rats
Naja naja atraNew Zealand rabbits of either sex (1.82 ± 0.09 kg)6500 (cytotoxin)0.077 (±0.018)-0.37 (±0.12)5.9 (±0.9)9 a (±4)0.56 b[27]
Rabbits
Naja sumatranaNew Zealand white rabbits (approx. 2 kg)3100 (PLA2)-68.6 (±0.8)-10.18 (±1.18)-0.048 b[70]
70 (neurotoxin)-81.5 (±0.6)-8.6 (±0.5)-0.082 b
150 (cardiotoxin)-45.6 (±0.1)-8.2 (±0.1)-0.087 b
500 (cardiotoxin in whole venom)-39.5 (±1.1)-11.6 (±0.9)-0.061 b
Whole Venom
Rabbits
Vipera aspisCharles de Bouscat HY rabbits (2.5–3 kg)5300-63 (±17.89 *)-32 (±8.94 *)--[35]
500-67 (±11.18 *)-36 (±8.94 *)--
700-63 (±38.01 *)-29 (±4.47 *)--
Naja sputatrixNew Zealand rabbits (2 kg)3500-41.7 **-18.9 (±4.6)-0.034 b[73]
Cryptelytrops purpureomaculatusMale New Zealand rabbits (1.7–2.1 kg)3500-41.6 (±3.0)-27 (±0.6)-0.055 (±0.004)[74]
Naja sumatranaNew Zealand white rabbits (approx. 2 kg)3500-41.9 (±0.2)-12.5 (±0.9)-0.047 b[70]
Hypnale hypnaleNew Zealand white rabbits (2.03 ± 0.06 kg)31000-4 **-19.3 (±1.21 *)-0.007 (±0.003 *)[67]
* SD calculated from standard errors reported in the original papers using the following equation: S D = S E   ×   √ n ; ** Assumed to be reported as SD; a = Apparent volume of distribution; b = Unit by weight (·kg−1) calculated on the basis of the mean of the reported animal weight; Abbreviations: t1/2ka, half-life of the absorption phase; F, bioavailability; t1/2α, half-life of the distribution phase; t1/2β, half-life of the elimination phase; Vss, apparent volume of distribution at steady state; CL, systemic clearance.
Table
Table 5. Pharmacokinetic parameters of snake venom and toxins following subcutaneous (SC) injection expressed as mean (±SD).
Table 5. Pharmacokinetic parameters of snake venom and toxins following subcutaneous (SC) injection expressed as mean (±SD).
Snake SpeciesAnimal ModelNo.Dose (mcg·kg−1)t1/2ka (h)F (%)t1/2α (h)t1/2β (h)Vss (L·kg−1)CL (L·h−1·kg−1)Ref
Isolated Toxin
Rats
Agkistrodon halys ussuriensis EmelianovSprague-Dawley rats of either sex (180–200 g)6750 (thrombin-like enzyme)2.5 (±0.73 *)-4.8 (±4.16 *)125 (±181.26 *)19 (±49 *)0.294 (±0.103 *)[68]
Whole Venom
Rabbits
Micrurus fulviusSheep (36–60 kg)45000 mcg-60 (±10 *)-4.35 (±1.83 *)0.56 b0.084 b[66]
* SD calculated from standard errors reported in the original papers using the following equation: S D = S E   ×   √ n ; b = Unit by weight (·kg−1) calculated on the basis of the mean of the reported animal weight; Abbreviations: t1/2ka, half-life of the absorption phase; F, bioavailability; t1/2α, half-life of the distribution phase; t1/2β, half-life of the elimination phase; Vss, apparent volume of distribution at steady state; CL, systemic clearance.
Table
Table 6. List of publications obtained from literature search and the demographics of snakebite cases.
Table 6. List of publications obtained from literature search and the demographics of snakebite cases.
Snake SpeciesCountryNumber of Patients with Timed Concentration Post-Bite for Data RetrievalRef
Pseudonaja spp.Australia5[25]
Vipera aspis, Vipera berus, and Vipera ammodytesFrance3[26]
Vipera aspis and Vipera berusFrance6[75]
Crotalus durissus terrificusBrazil9[76]
Vipera ammodytesSlovenia3[77]
Bothrops jararacaBrazil1[78]
Bothrops jararacaBrazil1[79]
Crotalus durissus terrificusBrazil1[80]
Bothrops lanceolatusMartinique1[81]
Daboia russelli siamensisThailand24[28]
Naja atraTaiwan14[82]
Daboia russelli siamensisTaiwan10[83]
Pseudonaja spp.Australia1[84]
Acanthophis spp.Australia1[85]
Denisonia maculataAustralia1[86]
Vipera russelliMyanmar38[29]
Pseudechis porphyriacusAustralia1[87]
Bitis gabonicaUK1[88]
Vipera russelli pulchellaSri Lanka1[89]
Crotalus durissusBrazil11[30]
Cerastes cerastes mutilaSwitzerland & UK2[90]
Bungarus caeruleusSri Lanka8[91]
Ophiophagus HannahMyanmar1[92]
Micropechis ikahekaPapua New Guinea1[93]

Share and Cite

MDPI and ACS Style

Sanhajariya, S.; Duffull, S.B.; Isbister, G.K. Pharmacokinetics of Snake Venom. Toxins 2018, 10, 73. https://doi.org/10.3390/toxins10020073

AMA Style

Sanhajariya S, Duffull SB, Isbister GK. Pharmacokinetics of Snake Venom. Toxins. 2018; 10(2):73. https://doi.org/10.3390/toxins10020073

Chicago/Turabian Style

Sanhajariya, Suchaya, Stephen B. Duffull, and Geoffrey K. Isbister. 2018. "Pharmacokinetics of Snake Venom" Toxins 10, no. 2: 73. https://doi.org/10.3390/toxins10020073

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop