A model for venom-induced consumptive coagulopathy in snake bite

The kidney histopathological spectrum of patients with kidney injury following snakebite envenomation in India: scoping review of five decades

INTRODUCTION

Snakebite is a public health problem leading to about 58,000 deaths every year in India. Kidney injury subsequent to snakebite envenomation is common with a reported prevalence of up to 32%. The current study aims to elucidate the spectrum of kidney histopathology in acute kidney injury (AKI) cases associated with snake bites.

METHODS

We searched seven electronic database studies to identify studies describing the histopathological findings in the kidney with snakebite envenomation. Two reviewers independently conducted titles and abstract screening as well as full-text evaluation for the final inclusion decision. Data were extracted as per the standardized form. We conducted narrative synthesis. Studies done exclusively on autopsy findings, in vitro studies, and case reports were excluded.

RESULTS

We retrieved 1464 studies and finally included 28 studies which met the eligibility criteria in the analysis. Most studies were single-centre and the majority were cross-sectional. Overall we included a total of 534 renal biopsies. Russell's viper bite was the most common cause related to AKI. Acute tubular necrosis was the most common finding followed by acute interstitial nephritis, acute cortical necrosis (ACN), and thrombotic microangiopathy (TMA). Vasculitis changes in vessels were rarely reported. Lesions such as ACN and TMA were associated with poor outcomes.

CONCLUSION

This analysis supports the notion that renal biopsies are important to guide prognosis and increase our knowledge about post-snake bite AKI pathophysiology.

Neutralization, by a polyspecific antivenom, of the coagulopathy induced by the venom of Bothrops asper: Assessment by standard coagulation tests and rotational thromboelastometry in a murine model

Venom-induced consumption coagulopathy and thrombocytopenia are common and potentially severe manifestations of viperid snakebite envenoming since they contribute to local and systemic hemorrhage. Therefore, the assessment of the efficacy of antivenoms to neutralize coagulopathic and thrombocytopenic toxins should be part of the preclinical evaluation of these drugs. To evaluate the efficacy of the polyvalent (Crotalinae) antivenom produced in Costa Rica, in this study we have used a mouse model of coagulopathy and thrombocytopenia induced by the venom of Bothrops asper, based on the bolus intravenous (i.v.) injection of venom. When venom and antivenom were incubated before injection, or when antivenom was administered i.v. immediately after venom injection, venom-induced hemostatic alterations were largely abrogated. We also studied the recovery rate of clotting parameters in conditions where antivenom was administered when mice were coagulopathic. Some parameters recovered more rapidly in antivenom-treated mice than in control envenomed animals, but others showed a spontaneous recovery without antivenom. This is due to a rapid clearance of plasma venom levels in these experimental conditions. This implies that models based on the bolus i.v. injection of venom have limitations for assessing the effect of antivenom in the recovery of clotting alterations once coagulopathy has developed. It is suggested that alternative models should be developed based on a slower systemic absorption of venom. Overall, our findings provide a protocol for the preclinical evaluation of antivenoms and demonstrate that the polyvalent antivenom is effective in neutralizing the toxins of B. asper venom responsible for coagulopathy and thrombocytopenia.

A Review of Quantitative Systems Pharmacology Models of the Coagulation Cascade: Opportunities for Improved Usability

Despite the numerous therapeutic options to treat bleeding or thrombosis, a comprehensive quantitative mechanistic understanding of the effects of these and potential novel therapies is lacking. Recently, the quality of quantitative systems pharmacology (QSP) models of the coagulation cascade has improved, simulating the interactions between proteases, cofactors, regulators, fibrin, and therapeutic responses under different clinical scenarios. We aim to review the literature on QSP models to assess the unique capabilities and reusability of these models. We systematically searched the literature and BioModels database reviewing systems biology (SB) and QSP models. The purpose and scope of most of these models are redundant with only two SB models serving as the basis for QSP models. Primarily three QSP models have a comprehensive scope and are systematically linked between SB and more recent QSP models. The biological scope of recent QSP models has expanded to enable simulations of previously unexplainable clotting events and the drug effects for treating bleeding or thrombosis. Overall, the field of coagulation appears to suffer from unclear connections between models and irreproducible code as previously reported. The reusability of future QSP models can improve by adopting model equations from validated QSP models, clearly documenting the purpose and modifications, and sharing reproducible code. The capabilities of future QSP models can improve from more rigorous validation by capturing a broader range of responses to therapies from individual patient measurements and integrating blood flow and platelet dynamics to closely represent in vivo bleeding or thrombosis risk.

Snakebite Associated Thrombotic Microangiopathy and Recommendations for Clinical Practice

Snakebite is a significant and under-resourced global public health issue. Snake venoms cause a variety of potentially fatal clinical toxin syndromes, including venom-induced consumption coagulopathy (VICC) which is associated with major haemorrhage. A subset of patients with VICC develop a thrombotic microangiopathy (TMA). This article reviews recent evidence regarding snakebite-associated TMA and its epidemiology, diagnosis, outcomes, and effectiveness of interventions including antivenom and therapeutic plasma-exchange. Snakebite-associated TMA presents with microangiopathic haemolytic anaemia (evidenced by schistocytes on the blood film), thrombocytopenia in almost all cases, and a spectrum of acute kidney injury (AKI). A proportion of patients require dialysis, most survive and achieve dialysis free survival. There is no evidence that antivenom prevents TMA specifically, but early antivenom remains the mainstay of treatment for snake envenoming. There is no evidence for therapeutic plasma-exchange being effective. We propose diagnostic criteria for snakebite-associated TMA as anaemia with >1.0% schistocytes on blood film examination, together with absolute thrombocytopenia (<150 × 10⁹/L) or a relative decrease in platelet count of >25% from baseline. Patients are at risk of long-term chronic kidney disease and long term follow up is recommended.

Efficacy of specific antivenom and fresh frozen plasma in reversing Daboia palaestinae venom-induced hemostatic derangement

OBJECTIVES

To evaluate the in vitro efficacy of specific Daboia (Vipera) palaestinae (Dp) antivenom or fresh frozen plasma (FFP) against Dp venom-induced hemostatic changes DESIGN: In vitro study.

SETTING

Laboratory of a university referral hospital.

ANIMALS

Five healthy dogs.

INTERVENTIONS

Rotational thromboelastometry (including recombinant tissue factor or kaolin activators [EXTEM and INTEM, respectively]) and conventional hemostatic tests were evaluated in citrated whole blood samples that underwent 4 treatments: (1) no additives (control); (2) 15 μg of Dp venom; (3) 15 μg of Dp venom and 10 μL of specific Dp antivenom; (4) 15 μg of Dp venom and 0.3 mL of FFP. A linear mixed-effects regression model was used to compare results between each treatment and the control.

MEASUREMENTS AND MAIN RESULTS

Dp-venom engendered statistically significant (P < 0.05) EXTEM changes in 8/17 variables, all indicative of hypercoagulability, which were negated by antivenom but not with FFP. Similarly, Dp-venom induced hypercoagulable, hyperfibrinolytic changes in 12 of 17 INTEM variables, of which only 5 of 12 were negated by antivenom but not with FFP. Fibrinogen concentration was decreased, and the activated partial thromboplastin time (aPTT) was shortened (P < 0.05 for both) in all treatments compared to the control.

CONCLUSIONS

This study demonstrated the ephemeral procoagulant phase of Dp envenomation for the first time. Many venom-induced thromboelastometric changes were reversed by specific antivenom but not with FFP. Prospective clinical studies are warranted to investigate whether the present findings translate to clinical efficacy.

Review article: Let us talk about snakebite management: A discussion on many levels

We want to discuss antivenom use in snakebite clinical practice guidelines. Coronial reviews in Victoria of two cases of snakebite envenomation, one described in detail below, prompted us to submit this paper for a wider audience and debate. Venom and antivenom levels were measured in the case detailed below, but not in the other. The coroner received conflicting and varied advice from experts regarding the dose of antivenom. The Victorian Department of Health and Human Services and the Australasian College for Emergency Medicine were instructed to review snakebite management guidelines, particularly with respect to antivenom dosage. The discussion that took place among medical experts led to considerable media attention. We discuss the potential fallout when there is no consensus among medical experts.

Understanding and reducing complex systems pharmacology models based on a novel input-response index

A growing understanding of complex processes in biology has led to large-scale mechanistic models of pharmacologically relevant processes. These models are increasingly used to study the response of the system to a given input or stimulus, e.g., after drug administration. Understanding the input-response relationship, however, is often a challenging task due to the complexity of the interactions between its constituents as well as the size of the models. An approach that quantifies the importance of the different constituents for a given input-output relationship and allows to reduce the dynamics to its essential features is therefore highly desirable. In this article, we present a novel state- and time-dependent quantity called the input-response index that quantifies the importance of state variables for a given input-response relationship at a particular time. It is based on the concept of time-bounded controllability and observability, and defined with respect to a reference dynamics. In application to the brown snake venom-fibrinogen (Fg) network, the input-response indices give insight into the coordinated action of specific coagulation factors and about those factors that contribute only little to the response. We demonstrate how the indices can be used to reduce large-scale models in a two-step procedure: (i) elimination of states whose dynamics have only minor impact on the input-response relationship, and (ii) proper lumping of the remaining (lower order) model. In application to the brown snake venom-fibrinogen network, this resulted in a reduction from 62 to 8 state variables in the first step, and a further reduction to 5 state variables in the second step. We further illustrate that the sequence, in which a recursive algorithm eliminates and/or lumps state variables, has an impact on the final reduced model. The input-response indices are particularly suited to determine an informed sequence, since they are based on the dynamics of the original system. In summary, the novel measure of importance provides a powerful tool for analysing the complex dynamics of large-scale systems and a means for very efficient model order reduction of nonlinear systems.

Catch a tiger snake by its tail: Differential toxicity, co-factor dependence and antivenom efficacy in a procoagulant clade of Australian venomous snakes

A paradigm of venom research is adaptive evolution of toxins as part of a predator-prey chemical arms race. This study examined differential co-factor dependence, variations relative to dietary preference, and the impact upon relative neutralisation by antivenom of the procoagulant toxins in the venoms of a clade of Australian snakes. All genera were characterised by venoms rich in factor Xa which act upon endogenous prothrombin. Examination of toxin sequences revealed an extraordinary level of conservation, which indicates that adaptive evolution is not a feature of this toxin type. Consistent with this, the venoms did not display differences on the plasma of different taxa. Examination of the prothrombin target revealed endogenous blood proteins are under extreme negative selection pressure for diversification, this in turn puts a strong negative selection pressure upon the toxins as sequence diversification could result in a drift away from the target. Thus this study reveals that adaptive evolution is not a consistent feature in toxin evolution in cases where the target is under negative selection pressure for diversification. Consistent with this high level of toxin conservation, the antivenom showed extremely high-levels of cross-reactivity. There was however a strong statistical correlation between relative degree of phospholipid-dependence and clotting time, with the least dependent venoms producing faster clotting times than the other venoms even in the presence of phospholipid. The results of this study are not only of interest to evolutionary and ecological disciplines, but also have implications for clinical toxinology.

A quantitative systems pharmacology model of blood coagulation network describes in vivo biomarker changes in non-bleeding subjects

Essentials Baseline coagulation activity can be detected in non-bleeding state by in vivo biomarker levels. A detailed mathematical model of coagulation was developed to describe the non-bleeding state. Optimized model described in vivo biomarkers with recombinant activated factor VII treatment. Sensitivity analysis predicted prothrombin fragment 1 + 2 and D-dimer are regulated differently.

SUMMARY

Background Prothrombin fragment 1 + 2 (F1 + 2 ), thrombin-antithrombin III complex (TAT) and D-dimer can be detected in plasma from non-bleeding hemostatically normal subjects or hemophilic patients. They are often used as safety or pharmacodynamic biomarkers for hemostatis-modulating therapies in the clinic, and provide insights into in vivo coagulation activity. Objectives To develop a quantitative systems pharmacology (QSP) model of the blood coagulation network to describe in vivo biomarkers, including F1 + 2 , TAT, and D-dimer, under non-bleeding conditions. Methods The QSP model included intrinsic and extrinsic coagulation pathways, platelet activation state-dependent kinetics, and a two-compartment pharmacokinetics model for recombinant activated factor VII (rFVIIa). Literature data on F1 + 2 and D-dimer at baseline and changes with rFVIIa treatment were used for parameter optimization. Multiparametric sensitivity analysis (MPSA) was used to understand key proteins that regulate F1 + 2 , TAT and D-dimer levels. Results The model was able to describe tissue factor (TF)-dependent baseline levels of F1 + 2 , TAT and D-dimer in a non-bleeding state, and their increases in hemostatically normal subjects and hemophilic patients treated with different doses of rFVIIa. The amount of TF required is predicted to be very low in a non-bleeding state. The model also predicts that these biomarker levels will be similar in hemostatically normal subjects and hemophilic patients. MPSA revealed that F1 + 2 and TAT levels are highly correlated, and that D-dimer is more sensitive to the perturbation of coagulation protein concentrations. Conclusions A QSP model for non-bleeding baseline coagulation activity was established with data from clinically relevant in vivo biomarkers at baseline and changes in response to rFVIIa treatment. This model will provide future mechanistic insights into this system.

Quantitative Systems Pharmacology Model to Predict the Effects of Commonly Used Anticoagulants on the Human Coagulation Network

Warfarin is the anticoagulant of choice for venous thromboembolism (VTE) treatment, although its suppression of the endogenous clot-dissolution complex APC:PS may ultimately lead to longer time-to-clot dissolution profiles, resulting in increased risk of re-thrombosis. This detrimental effect might not occur during VTE treatment using other anticoagulants, such as rivaroxaban or enoxaparin, given their different mechanisms of action within the coagulation network. A quantitative systems pharmacology model was developed describing the coagulation network to monitor clotting factor levels under warfarin, enoxaparin, and rivaroxaban treatment. The model allowed for estimation of all factor rate constants and production rates. Predictions of individual coagulation factor time courses under steady-state warfarin, enoxaparin, and rivaroxaban treatment reflected the suppression of protein C and protein S under warfarin compared to rivaroxaban and enoxaparin. The model may be used as a tool during clinical practice to predict effects of anticoagulants on individual clotting factor time courses and optimize antithrombotic therapy.

An In Silico Knockout Model for Gastrointestinal Absorption Using a Systems Pharmacology Approach - Development and Application for Ketones

Gastrointestinal absorption and disposition of ketones is complex. Recent work describing the pharmacokinetics (PK) of d-β-hydroxybutyrate (BHB) following oral ingestion of a ketone monoester ((R)-3-hydroxybutyl (R)-3-hydroxybutyrate) found multiple input sites, nonlinear disposition and feedback on endogenous production. In the current work, a human systems pharmacology model for gastrointestinal absorption and subsequent disposition of small molecules (monocarboxylic acids with molecular weight < 200 Da) was developed with an application to a ketone monoester. The systems model was developed by collating the information from the literature and knowledge gained from empirical population modelling of the clinical data. In silico knockout variants of this systems model were used to explore the mechanism of gastrointestinal absorption of ketones. The knockouts included active absorption across different regions in the gut and also a passive diffusion knockout, giving 10 gut knockouts in total. Exploration of knockout variants has suggested that there are at least three distinct regions in the gut that contribute to absorption of ketones. Passive diffusion predominates in the proximal gut and active processes contribute to the absorption of ketones in the distal gut. Low doses are predominantly absorbed from the proximal gut by passive diffusion whereas high doses are absorbed across all sites in the gut. This work has provided mechanistic insight into the absorption process of ketones, in the form of unique in silico knockouts that have potential for application with other therapeutics. Future studies on absorption process of ketones are suggested to substantiate findings in this study.

A Systems Pharmacology Model for Predicting Effects of Factor Xa Inhibitors in Healthy Subjects: Assessment of Pharmacokinetics and Binding Kinetics

Factor Xa (FXa) emerged as a promising target for effective anticoagulation and several FXa inhibitors are now available for the prevention of venous thromboembolism. However, in previously reported pharmacokinetic/pharmacodynamic (PK/PD) models, the complex coagulation processes and detailed information of drug action are usually unclear, which makes it difficult to predict clinical outcome at the drug discovery stage. In this study, a large‐scale systems pharmacology model was developed based on several published models and clinical data. It takes into account all pathways of the coagulation network, and captures drug‐specific features: plasma pharmacokinetics and drug‐target binding kinetics (BKs). We aimed to predict the anticoagulation effects of FXa inhibitors in healthy subjects, and to use this model to compare the effects of compounds with different binding properties. Our model predicts the clotting time and anti‐FXa effects and could thus serve as a predictive tool for the anticoagulant potential of a new compound.

The simultaneous occurrence of both hypercoagulability and hypofibrinolysis in blood and serum during systemic inflammation, and the roles of iron and fibrin(ogen)

Although the two phenomena are usually studied separately, we summarise a considerable body of literature to the effect that a great many diseases involve (or are accompanied by) both an increased tendency for blood to clot (hypercoagulability) and the resistance of the clots so formed (hypofibrinolysis) to the typical, 'healthy' or physiological lysis. We concentrate here on the terminal stages of fibrin formation from fibrinogen, as catalysed by thrombin. Hypercoagulability goes hand in hand with inflammation, and is strongly influenced by the fibrinogen concentration (and vice versa); this can be mediated via interleukin-6. Poorly liganded iron is a significant feature of inflammatory diseases, and hypofibrinolysis may change as a result of changes in the structure and morphology of the clot, which may be mimicked in vitro, and may be caused in vivo, by the presence of unliganded iron interacting with fibrin(ogen) during clot formation. Many of these phenomena are probably caused by electrostatic changes in the iron-fibrinogen system, though hydroxyl radical (OH˙) formation can also contribute under both acute and (more especially) chronic conditions. Many substances are known to affect the nature of fibrin polymerised from fibrinogen, such that this might be seen as a kind of bellwether for human or plasma health. Overall, our analysis demonstrates the commonalities underpinning a variety of pathologies as seen in both hypercoagulability and hypofibrinolysis, and offers opportunities for both diagnostics and therapies.

Numerical validation of a synthetic cell-based model of blood coagulation

In Fasano et al. (2012) a new reduced mathematical model for blood coagulation was proposed, incorporating biochemical and mechanical actions of blood flow and including platelets activity. The model was characterized by a considerable simplification of the differential system associated to the biochemical network and it incorporated the role of blood slip at the vessel wall as an extra source of activated platelets. The purpose of this work is to check the validity of the reduced mathematical model, using as a benchmark the model presented in Anand et al. (2008), and to investigate the importance of the blood slip velocity in the blood coagulation process.

Prothrombin activator-like toxin appears to mediate cardiovascular collapse following envenoming by Pseudonaja textilis

Brown snake (Pseudonaja spp.)-induced early cardiovascular collapse is a life-threatening medical emergency in Australia. We have previously shown that this effect can be mimicked in animals and is mediated via the release of endogenous mediators. In the present study, we aimed to purify and characterize the component in Pseudonaja textilis venom which induces cardiovascular collapse following envenoming. The component (fraction 3) was isolated using a combination of techniques including hydroxyapatite and reverse phase chromatography. Fraction 3 (10 or 20 μg/kg, i.v.) produced a rapid decrease in mean arterial pressure (MAP) followed by cardiovascular collapse. Fraction 3-induced early collapse was abolished by prior administration of smaller priming doses of fraction 3 (i.e. 2 and 5 μg/kg, i.v.) or heparin (300 units/kg, i.v.). P. textilis whole venom (1 and 3 μg/ml), but not fraction 3 (1 or 3 μg/ml), induced endothelium-dependent relaxation in isolated rat mesenteric arteries. SDS-PAGE gel indicated the presence of 9-10 protein bands of fraction 3. Using proteomic based analysis some protein bands of fraction 3 were identified as subunits of venom prothrombin activator, pseutarin C of P. textilis venom. Our results conclude that prothrombin activator-like toxin is likely to be a contributor to the rapid collapse induced by P. textilis venom.

Current treatment for venom-induced consumption coagulopathy resulting from snakebite

Venomous snakebite is considered the single most important cause of human injury from venomous animals worldwide. Coagulopathy is one of the commonest important systemic clinical syndromes and can be complicated by serious and life-threatening haemorrhage. Venom-induced consumption coagulopathy (VICC) is the commonest coagulopathy resulting from snakebite and occurs in envenoming by Viperid snakes, certain elapids, including Australian elapids, and a few Colubrid (rear fang) snakes. Procoagulant toxins activate the clotting pathway, causing a broad range of factor deficiencies depending on the particular procoagulant toxin in the snake venom. Diagnosis and monitoring of coagulopathy is problematic, particularly in resource-poor countries where further research is required to develop more reliable, cheap clotting tests. MEDLINE and EMBASE up to September 2013 were searched to identify clinical studies of snake envenoming with VICC. The UniPort database was searched for coagulant snake toxins. Despite preclinical studies demonstrating antivenom binding toxins (efficacy), there was less evidence to support clinical effectiveness of antivenom for VICC. There were no placebo-controlled trials of antivenom for VICC. There were 25 randomised comparative trials of antivenom for VICC, which compared two different antivenoms (ten studies), three different antivenoms (four), two or three different doses or repeat doses of antivenom (five), heparin treatment and antivenom (five), and intravenous immunoglobulin treatment and antivenom (one). There were 13 studies that compared two groups in which there was no randomisation, including studies with historical controls. There have been numerous observational studies of antivenom in VICC but with no comparison group. Most of the controlled trials were small, did not use the same method for assessing coagulopathy, varied the dose of antivenom, and did not provide complete details of the study design (primary outcomes, randomisation, and allocation concealment). Non-randomised trials including comparison groups without antivenom showed that antivenom was effective for some snakes (e.g., Echis), but not others (e.g., Australasian elapids). Antivenom is the major treatment for VICC, but there is currently little high-quality evidence to support effectiveness. Antivenom is not risk free, and adverse reactions can be quite common and potentially severe. Studies of heparin did not demonstrate it improved outcomes in VICC. Fresh frozen plasma appeared to speed the recovery of coagulopathy and should be considered in bleeding patients.

Scale reduction of a systems coagulation model with an application to modeling pharmacokinetic-pharmacodynamic data

Bridging systems biology and pharmacokinetics–pharmacodynamics has resulted in models that are highly complex and complicated. They usually contain large numbers of states and parameters and describe multiple input–output relationships. Based on any given data set relating to a specific input–output process, it is possible that some states of the system are either less important or have no influence at all. In this study, we explore a simplification of a systems pharmacology model of the coagulation network for use in describing the time course of fibrinogen recovery after a brown snake bite. The technique of proper lumping is used to simplify the 62-state systems model to a 5-state model that describes the brown snake venom–fibrinogen relationship while maintaining an appropriate mechanistic relationship. The simplified 5-state model explains the observed decline and recovery in fibrinogen concentrations well. The techniques used in this study can be applied to other multiscale models.

Antivenom therapy of carpet viper (Echis ocellatus) envenoming: effectiveness and strategies for delivery in West Africa

In West Africa, response to specific, geographically appropriate, antivenom is often dramatic following carpet viper (Echis ocellatus) envenoming with rapid restoration of blood coagulability and resolution of spontaneous haemorrhage. Envenoming from Australasian snakes causing similar coagulopathies may respond less dramatically and the effectiveness of antivenom is being questioned. Here we have reviewed and re-analysed all published preclinical and clinical studies on envenoming and antivenom therapy conducted in West Africa to determine the effectiveness of antivenom. 22 studies provided relevant information: 12 observational studies, 4 RCTs and 6 preclinical studies. Four comparative studies confirmed statistically significant protection against mortality ranging from 57 to 87% using specific antivenoms compared to non-specific or no antivenoms. Meta-analysis estimated combined Odds Ratio (95% CI) of 0.25 (0.14-0.45) of dying among those treated with specific antivenom or 75% (95% CI: 55-86%) protection against death. Mortality more than doubled during times when stocks of reliable antivenoms ran out, with Relative Risk (95% CI)] of 2.33 (1.26-4.06). Serum kinetics of venom antigen/antivenom levels also confirmed that decline of venom antigen levels coincided with resolution of coagulopathy while decline of antivenom levels was associated with venom antigen reappearance and recurrence of coagulopathy. Preclinical and antivenomics analysis confirmed efficacy of regionally appropriate antivenoms against E. ocellatus and related species' venoms in Sub-Saharan Africa but not against Asian Echis carinatus venom. Antivenoms raised against E. carinatus were ineffective in human studies. In West Africa, specific antivenom is effective in managing carpet viper envenoming. A centralized hub-and-spoke strategy is suggested for broadening antivenom access to endemic rural areas together with instituting quality assurance, standardization and manpower training. Benefits, risks, cost-effectiveness and feasibility of the approach should be formally assessed.

Effect of Australian elapid venoms on blood coagulation: Australian Snakebite Project (ASP-17)

Snake venoms contain toxins that activate the coagulation network and cause venom-induced consumption coagulopathy. A previously developed mathematical model of the coagulation network was refined and used to describe and predict the time course of changes in the coagulation factors following envenomation by Brown snake (Pseudonaja spp.), Tiger snake (Notechis scutatus), Rough-scaled snake (Tropidechis carinatus) and Hoplocephalus spp. (Stephens banded, Pale headed and Broad headed). Simulations of the time course of the change in coagulation factors were compared to data obtained from a large prospective study of Australian snake bites - the Australian Snakebite Project. The model predictions were also compared against data for partial and complete VICC obtained from the same study. The model simulations were used to understand the differences in consumption and recovery of clotting factors in partial versus complete VICC as well as among bites from different snake types. The model suggested that the venoms were absorbed almost instantaneously and provided a reasonable prediction of the observed concentration of clotting factors over time in patients bitten by Australian elapid snakes. The model predictions suggested a higher consumption of factors (fibrinogen, II and IX in particular) in patients with complete VICC compared to those with partial VICC. The model also predicted that snakes with "Xa-like" venoms may produce a less severe VICC than snakes with "Xa:Va-like" venoms.

The belonging of gpMuc, a glycoprotein from Mucuna pruriens seeds, to the Kunitz-type trypsin inhibitor family explains its direct anti-snake venom activity

In Nigeria, Mucuna pruriens seeds are locally prescribed as an oral prophylactic for snake bite and it is claimed that when two seeds are swallowed they protect the individual for a year against snake bites. In order to understand the Mucuna pruriens antisnake properties, the proteins from the acqueous extract of seeds were purified by three chromatographic steps: ConA affinity chromatography, tandem anionic-cationic exchange and gel filtration, obtaining a fraction conventionally called gpMucB. This purified fraction was analysed by SDS-PAGE obtaining 3 bands with apparent masses ranging from 20 to 24 kDa, and by MALDI-TOF which showed two main peaks of 21 and 23 kDa and another small peak of 19 kDa. On the other hand, gel filtration analysis of the native protein indicated a molecular mass of about 70 kDa suggesting that in its native form, gpMucB is most likely an oligomeric multiform protein. Infrared spectroscopy of gpMucB indicated that the protein is particularly thermostable both at neutral and acidic pHs and that it is an all beta protein. All data suggest that gpMucB belongs to the Kunitz-type trypsin inhibitor family explaining the direct anti-snake venom activity of Mucuna pruriens seeds.

Performance and robustness of the Monte Carlo importance sampling algorithm using parallelized S-ADAPT for basic and complex mechanistic models

The Monte Carlo Parametric Expectation Maximization (MC-PEM) algorithm can approximate the true log-likelihood as precisely as needed and is efficiently parallelizable. Our objectives were to evaluate an importance sampling version of the MC-PEM algorithm for mechanistic models and to qualify the default estimation settings in SADAPT-TRAN. We assessed bias, imprecision and robustness of this algorithm in S-ADAPT for mechanistic models with up to 45 simultaneously estimated structural parameters, 14 differential equations, and 10 dependent variables (one drug concentration and nine pharmacodynamic effects). Simpler models comprising 15 parameters were estimated using three of the ten dependent variables. We set initial estimates to 0.1 or 10 times the true value and evaluated 30 bootstrap replicates with frequent or sparse sampling. Datasets comprised three dose levels with 16 subjects each. For simultaneous estimation of the full model, the ratio of estimated to true values for structural model parameters (median [5-95% percentile] over 45 parameters) was 1.01 [0.94-1.13] for means and 0.99 [0.68-1.39] for between-subject variances for frequent sampling and 1.02 [0.81-1.47] for means and 1.02 [0.47-2.56] for variances for sparse sampling. Imprecision was ≤25% for 43 of 45 means for frequent sampling. Bias and imprecision was well comparable for the full and simpler models. Parallelized estimation was 23-fold (6.9-fold) faster using 48 threads (eight threads) relative to one thread. The MC-PEM algorithm was robust and provided unbiased and adequately precise means and variances during simultaneous estimation of complex, mechanistic models in a 45 dimensional parameter space with rich or sparse data using poor initial estimates.

Development of a new pre- and post-processing tool (SADAPT-TRAN) for nonlinear mixed-effects modeling in S-ADAPT

Mechanistic modeling greatly benefits from automated pre- and post-processing of model code and modeling results. While S-ADAPT provides many state-of-the-art parametric population estimation methods, its pre- and post-processing capabilities are limited. Our objective was to develop a fully automated, open-source pre- and post-processor for nonlinear mixed-effects modeling in S-ADAPT. We developed a new translator tool (SADAPT-TRAN) based on Perl scripts. These scripts (a) automatically translate the core model components into robust Fortran code, (b) perform extensive mutual error checks across all input files and the raw dataset, (c) extend the options of the Monte Carlo Parametric Expectation Maximization (MC-PEM) algorithm, and (d) improve the numerical robustness of the model code. The post-processing scripts automatically summarize the results of one or multiple models as tables and, by generating problem specific R scripts, provide an extended series of standard and covariate-stratified diagnostic plots. The SADAPT-TRAN package substantially improved the efficiency to specify, debug, and evaluate models and enhanced the flexibility of using the MC-PEM algorithm for parallelized estimation in S-ADAPT. The parameter variability model can take any combination of normally, log-normally, or logistically distributed parameters and the SADAPT-TRAN package can automatically generate the Fortran code required to specify between occasion variability. Extended estimation features are available to avoid local minima, estimate means with negligible variances, and estimate variances for fixed means. The SADAPT-TRAN package significantly facilitated model development in S-ADAPT, reduced model specification errors, and provided useful error messages for beginner and advanced users. This benefit was greatest for complex mechanistic models.

In vitro effects of Echis carinatus venom on the human plasma proteome

Echis carinatus venom (EV) is a complex mixture of toxins that contribute to its lethality. EV proteolytic activity was analyzed by zymography, chromogenic assays, and SDS-PAGE. To understand the molecular mechanism of the envenomation, we investigated the in vitro effect of EV on human plasma proteins. We looked for EV protein substrates and their proteolytic fragments. We analyzed EV proteolytic activity on standard proteins such as prothrombin or fibrinogen. To set up the optimal EV:plasma protein ratio conditions, plasma was incubated with EV (treated plasma), depleted of abundant proteins, and subjected to SDS-PAGE. Samples from control and treated plasma were also analyzed by 2-DE/MALDI-TOF MS, leading to the identification of four classes of plasma proteins cleaved by EV: proteases, protease inhibitors, binding proteins, and transporters. EV mainly proteolyzes entire proteins but can also act on physiological fragments. In summary, the physiological effects of EV proteases involve other important processes in addition to blood coagulation; complement activation and hemoglobin metabolism are also affected. In particular, the cleavage of protease inhibitors appears to be the mechanism through which the venom neutralizes the body's defenses.

Factor deficiencies in venom-induced consumption coagulopathy resulting from Australian elapid envenomation: Australian Snakebite Project (ASP-10)

BACKGROUND

Limited information exists on the dynamics of hemostasis in patients with venom-induced consumption coagulopathy (VICC) from snake envenomation.

OBJECTIVE

The aim of the present study was to investigate specific factor deficiencies and their time course in Australasian elapid envenomation.

METHODS

We measured coagulation parameters and factor concentrations in patients recruited to the Australian Snakebite Project, an observational cohort study. There were 112 patients with complete VICC, defined as an international normalized ratio (INR) > 3, and 18 with partial VICC. Serial citrated plasma samples were collected from 0.5 to 60 h post-bite. INR, activated partial thromboplastin time (aPTT), coagulation factors (F)I, II, V, VII, VIII, IX, X, von Willebrand factor antigen (VWF:Ag) and D-dimer concentrations were measured.

RESULTS

Complete VICC was characterized by near/total depletion of fibrinogen, FV and FVIII, with an INR and aPTT that exceeded the upper limits of detection, within 2 h of snakebite. Prothrombin levels never fell below 60% of normal, suggesting that the toxins were rapidly eliminated or inactivated and re-synthesis of clotting factors occurred irrespective of antivenom. Partial VICC caused limited depletion of fibrinogen and FV, and almost complete consumption of FVIII. Onset of VICC was more rapid with brown snake (Pseudonaja spp.) venom, which contains a group C prothrombin activator toxin, compared with the tiger snake group, which contains a group D prothrombin activator toxin and requires human FVa formation. Resolution of VICC occurred within 24-36 h irrespective of snake type.

CONCLUSIONS

These results suggest that Australasian elapid prothrombin activators have a potent but short duration of action. Antivenom is unlikely to be administered in time to prevent VICC.

Snake bite

Snake bite is a common and frequently devastating environmental and occupational disease, especially in rural areas of tropical developing countries. Its public health importance has been largely ignored by medical science. Snake venoms are rich in protein and peptide toxins that have specificity for a wide range of tissue receptors, making them clinically challenging and scientifically fascinating, especially for drug design. Although the full burden of human suffering attributable to snake bite remains obscure, hundreds of thousands of people are known to be envenomed and tens of thousands are killed or maimed by snakes every year. Preventive efforts should be aimed towards education of affected communities to use proper footwear and to reduce the risk of contact with snakes to a minimum through understanding of snakes' behaviour. To treat envenoming, the production and clinical use of antivenom must be improved. Increased collaboration between clinicians, epidemiologists, and laboratory toxinologists should enhance the understanding and treatment of envenoming.

Antivenom efficacy or effectiveness: the Australian experience

Despite widespread use of antivenoms, many questions remain about their effectiveness in the clinical setting. The almost universal acceptance of their value is based mainly on in vitro studies, animal studies and human observational studies. Numerous examples exist where they demonstrate clear benefit, such as consumption coagulopathy in viper envenoming, prevention of neurotoxicity in Australasian elapid bites, systemic effects in scorpion and funnel-web spider envenoming. There are also concerns about the quality and efficacy of some antivenoms. However, it is important not to confuse the efficacy of antivenom, defined as its ability to bind and neutralise venom-mediated effects under ideal conditions, and the effectiveness of antivenom, defined as its ability to reverse or prevent envenoming in human cases. There are numerous potential reasons for antivenom failure in human envenoming, of which antivenom inefficacy is only one. Other important reasons include venom-mediated effects being irreversible, antivenom being unable to reach the site of toxin-mediated injury, or the rapidity of onset of venom-mediated effects. A number of recent studies in Australia bring into question the effectiveness of some antivenoms, including snake antivenom for coagulopathy, redback spider and box jellyfish antivenoms. Despite brown snake antivenom being able to neutralise venom induced clotting in vitro, use of the antivenom in human envenoming does not appear to change the time course of coagulopathy. However, it is important that apparent antivenom ineffectiveness in specific cases is correctly interpreted and does not lead to a universal belief that antivenom is ineffective. It should rather encourage further studies to investigate the underlying pathophysiology of envenoming, the pharmacokinetics of venoms and antivenoms, and ultimately the effectiveness of antivenom based on snake type, clinical effects and timing of administration.

Failure of antivenom to improve recovery in Australian snakebite coagulopathy

BACKGROUND

Venom-induced consumption coagulopathy (VICC) is an important feature of snake envenoming.

AIM

To investigate the effect of antivenom and fresh frozen plasma (FFP) on recovery of VICC in Australian elapid snake envenoming.

DESIGN

Prospective cohort study.

METHODS

Patients with VICC were included from the Australian Snakebite Project (ASP). Time to recovery of VICC (defined as time until INR <2) was investigated using a time to event analysis in WinBUGS. The model considered the effects of age, sex, snake type, time of antivenom after bite, antivenom dose and use of FFP within 4 h.

RESULTS

The study included 167 cases of VICC, median age being 41 [interquartile range (IQR): 28-53) years, and 130 (78%) were males. Antivenom was administered at a median of 3.6 (IQR: 2.2-5.6) h after the bite at a median dose of four vials (IQR: 2-6 vials). Thirteen patients received FFP within 4 h. Recovery of VICC occurred after a median of 14.4 (IQR: 11.5-17.5) h, and only the use of FFP within 4 h influenced the time to recovery. Neither antivenom dose nor time of antivenom administration had an effect on recovery of VICC. In patients administered with FFP, 12% [credible interval (CrI): 6-21%] and 81% (CrI: 61-94%) had recovered at 6 and 12 h, respectively, vs 2.5% (CrI: 1.5-4%) and 28% (CrI: 22-34%) not receiving FFP.

DISCUSSION

Antivenom did not appear to be effective for the coagulopathy in snake envenoming in Australia. FFP appeared to shorten the time of VICC recovery.