Exploring the Utility of ssDNA Aptamers Directed against Snake Venom Toxins as New Therapeutics for Snakebite Envenoming

Snakebite is a neglected tropical disease that causes considerable death and disability in the tropical world. Although snakebite can cause a variety of pathologies in victims, haemotoxic effects are particularly common and are typically characterised by haemorrhage and/or venom-induced consumption coagulopathy. Antivenoms are the mainstay therapy for treating the toxic effects of snakebite, but despite saving thousands of lives annually, these therapies are associated with limited cross-snake species efficacy due to venom variation, which ultimately restricts their therapeutic utility to particular geographical regions. In this study, we sought to explore the potential of ssDNA aptamers as toxin-specific inhibitory alternatives to antibodies. As a proof of principle model, we selected snake venom serine protease toxins, which are responsible for contributing to venom-induced coagulopathy following snakebite envenoming, as our target. Using SELEX technology, we selected ssDNA aptamers against recombinantly expressed versions of the fibrinogenolytic SVSPs ancrod from the venom of C. rhodostoma and batroxobin from B. atrox. From the resulting pool of specific ssDNA aptamers directed against each target, we identified candidates that exhibited low nanomolar binding affinities to their targets. Downstream aptamer-linked immobilised sorbent assay, fibrinogenolysis, and coagulation profiling experiments demonstrated that the candidate aptamers were able to recognise native and recombinant SVSP toxins and inhibit the toxin- and venom-induced prolongation of plasma clotting times and the consumption of fibrinogen, with inhibitory potencies highly comparable to commercial polyvalent antivenoms. Our findings demonstrate that rationally selected toxin-specific aptamers can exhibit broad in vitro cross-reactivity against toxin isoforms found in different snake venoms and are capable of inhibiting toxins in pathologically relevant in vitro and ex vivo models of venom activity. These data highlight the potential utility of ssDNA aptamers as novel toxin-inhibiting therapeutics of value for tackling snakebite envenoming.


Introduction
Snakebite envenoming is a significant public health issue as more than 5.4 million people are bitten annually by venomous snakes, and this mainly affects resource-poor of multiple vials (ranging from 3 to 30 vials depending on region and product), resulting in costs exceeding USD 1000 in parts of sub-Saharan Africa, for example [29]. This cost pushes many tropical snakebite victims further into poverty [30,31].
In response to the above, recent research has focused on applying alternate therapeutic strategies to circumvent the limitations currently associated with conventional antivenoms. Promising approaches explored to date include the use of toxin-specific monoclonal antibodies [32], small-molecule drugs [33], decoy receptor binding proteins [34] and aptamers [35][36][37]. Although arguably the least studied of these approaches in the context of snakebite, single-stranded DNA (ssDNA) aptamers have emerged as promising alternatives to antibodies in various biosensing platforms and as the biological recognition element for therapeutic applications and diagnostic tools for different diseases [38][39][40]. Aptamers are short ssDNA (or ssRNA) molecules that can bind to their target and fold into complex and stable three-dimensional shapes [41,42] with high specificity and affinity through hydrogen bonding, van der Waals forces, hydrophobic, salt bridges and other electrostatic interactions [43][44][45][46]. Consequently, aptamers act like antibodies by binding and inhibiting target antigens and, thus, have been referred to as chemical antibodies due to their synthetic production [47,48]. ssDNA aptamers possess several advantages over antibodies as, in addition to seemingly exhibiting comparable high-affinity recognition properties (i.e., dissociation constants in the low nanomolar (nM)/picomolar (pM) range) [43,44,49], they are highly specific, sensitive, and chemically stable with long shelf lives [50][51][52] and require facile artificial synthesis [53], resulting in high purity, low interbatch variability and cost-effective production [54][55][56]. Further, ssDNA aptamers are poorly immunogenic, which aids their tolerance, and due to their small size, they diffuse readily into tissues [57]. However, this latter characteristic also causes challenges, as the in vivo half-life of ssDNA aptamers is much shorter (~20 min) than antibodies (e.g., potentially a few days to weeks) [58], though prolonged circulation can be improved by chemical modification [57,[59][60][61]. Despite the therapeutic promise of aptamers, to date, the only FDA-approved aptamer remains the macular degeneration treatment Macugen ® (Pegaptanib), which was first licensed in 2005 [62]. Since then, no other aptamers have entered the market, although many candidate aptamers have entered clinical trials over the past decade as promising new treatments for heart disease, cancer and type II diabetes [62].
ssDNA aptamers are typically selected by Systematic Evolution of Ligands by Exponential enrichment (SELEX) technology, which was introduced by Ellington and Tuerk three decades ago [43,44]. This process is based on the isolation of high-affinity ligands from a combinatorial ssDNA library (approximately 10 17 random sequence oligonucleotides) through repeated cycles (usually requiring seven to fifteen rounds) of binding, elution, and amplification to select for oligonucleotides with the desired specificity and/or sensitivity [48,63,64]. Following cyclical SELEX enrichment, the final ssDNA aptamer pool is subjected to sequencing to identify the optimal resulting binding sequences, which are then manufactured by chemical synthesis and downstream-optimized by the introduction of chemical modifications to further improve their pharmacokinetic and/or pharmacodynamic profiles [65]. To date, several ssDNA aptamers have been successfully selected for and applied against a wide range of targets that include proteins, peptides, small organic compounds, carbohydrates, metal ions and biological cofactors [66][67][68][69][70][71][72]. In addition, ssDNA aptamers have been shown to bind to highly toxic native antigens, which, due to toxicity, cannot be achieved in animal-based methods to generate specific antibodies without antigen modification [73]. In the context of animal toxins, little work has been done, though aptamers have been selected against the three-finger toxin α-bungarotoxin from Bungarus multicinctus venom [74] and shown to exhibit cross-recognition and inhibitory activity against cytotoxins from Naja atra venom, which share similar tertiary structures [37]. Another recent study selected ssDNA aptamers against Indian Bungarus caeruleus venom and used them in a diagnostic context to discriminate envenomings by this species from those of other medically important snakes found in the region [75].
In this study, we sought to explore the potential utility of ssDNA aptamers against snake venom SVSP toxins responsible for contributing to venom-induced coagulopathy following snakebite envenoming. Using SELEX technology, we selected ssDNA aptamers against recombinantly expressed versions of the fibrinogenolytic SVSPs ancrod from the venom of C. rhodostoma and batroxobin from Bothrops atrox. Downstream ssDNA aptamerlinked immobilised sorbent assay (ALISA), fibrinogenolysis, and coagulation experiments demonstrated that the rationally selected ssDNA aptamers were able to recognise native and recombinant SVSP toxins and inhibit the toxin-and venom-induced prolongation of plasma clotting times and the consumption of fibrinogen. These data highlight the potential utility of ssDNA aptamers as novel toxin-inhibiting therapeutics and alternative binding scaffolds for future treatment strategies targeting tropical snakebites.

ssDNA Aptamer Selection
Here, we report the identification of ssDNA aptamers that target the recombinantly expressed SVSP toxins ancrod from C. rhodostoma and batroxobin from B. atrox, which we previously expressed in human embryonic kidney mammalian cell lines (HEK293) [76]. To do so, we employed a parallel process consisting of 14 rounds of SELEX selections against each toxin target using an ssDNA library pool of 10 17 random sequences. The selection process was monitored throughout by polymerase chain reaction (PCR) amplification of the eluted double-stranded DNA (dsDNA) and visualised by using 2% agarose gel electrophoresis ( Figure 1A,B) and by measuring the fluorescence intensity of the eluted ssDNA solution by using denaturing polyacrylamide gel electrophoresis (PAGE) ( Figure 1C,D).
Our findings revealed that initial DNA recovery against each target was low, though this increased with progressive rounds of selection cycles until a marked increase was observed at round 11, indicating enrichment of the DNA pool with sequences specific to the toxins coupled to the beads ( Figure 1). Notably, the patterns of fluorescent intensities obtained from the sequential selection rounds were highly comparable between the two toxin targets (Figure 1). Counter-selection (CS) was performed at round 12 to eliminate nonspecific ssDNA interacting with the Sepharose beads instead of the toxins, and this resulted in an anticipated drop in the fluorescence intensity of the recovered DNA ( Figure 1C,D). Thereafter, two additional rounds of SELEX selection were performed to secondarily enrich toxin-specific aptamers, which resulted in an increase in fluorescence intensity at round 13 but no further increase at round 14, suggesting that target binding sites might be saturated ( Figure 1C,D). Consequently, SELEX selection was stopped, and the resulting ssDNA from the last round (round 14) was cloned for downstream aptamer identification.
Cloning was performed by ligation into a cloning vector and transformation into E. coli competent cells, with validation performed on ten random colonies to demonstrate that the inserted DNA (300 bp) was successfully ligated into the pCR2.1-TOPO cloning vector (3956 bp) ( Figure 1E,F). Thereafter, 100 colonies were picked for each toxin target and subjected to PCR, followed by Sanger sequencing. Bioinformatic analyses of the resulting DNA identified various identical sequences, indicating enrichment of the DNA pool, and resulted in thirteen and nine unique aptamer sequences against ancrod and batroxobin, respectively (Table 1). (C,D) The fluorescence intensity of the eluted ssDNA solution from each round of selection was measured at an excitation wavelength of 470 nm and an emission of 515 nm. Error bars represent the standard error (SEM) of triplicate measurements, with data for ancrod shown in (C) and batroxobin in (D). (E,F) Colony PCR was performed by ligation into a cloning vector and transformation into E. coli competent cells, with validation performed on ten random colonies to demonstrate that the inserted DNA (Lane 1-10; size ≈ 300 bp) was successfully ligated into the pCR2.1-TOPO cloning vector (Lane 11; size ≈ 3956 bp). Data for ancrod are shown in (E) and batroxobin in (F). The fluorescence intensity of the eluted ssDNA solution from each round of selection was measured at an excitation wavelength of 470 nm and an emission of 515 nm. Error bars represent the standard error (SEM) of triplicate measurements, with data for ancrod shown in (C) and batroxobin in (D). (E,F) Colony PCR was performed by ligation into a cloning vector and transformation into E. coli competent cells, with validation performed on ten random colonies to demonstrate that the inserted DNA (Lane 1-10; size ≈ 300 bp) was successfully ligated into the pCR2.1-TOPO cloning vector (Lane 11; size ≈ 3956 bp). Data for ancrod are shown in (E) and batroxobin in (F).

Dissociation Constants of Identified Aptamers
The dissociation constants (K d ) for each of the selected ssDNA aptamers were determined via the nonlinear regression of saturation graphs resulting from concentrationbased fluorescence binding assays. Of the 22 aptamers tested, only three ssDNA aptamer sequences directed against each of the targets exhibited nanomolar affinity, with those targeting ancrod ranging from 3.0 to 17.8 nM and those against batroxobin ranging from 4.7 to 24.3 nM (Figure 2). The remaining ssDNA aptamers exhibited high K d readings and were excluded from further study. The candidate ssDNA aptamer sequences with the lowest K d were ancrod-55 (K d = 3.0 nM) and batroxobin-26 (K d = 4.7 nM) ( Figure 2); thus, these two aptamers were selected for downstream recombinant toxin/crude venom binding and inhibition experiments. The equilibrium Kds of the highest affinity ssDNA aptamers were derived from the fluorescence binding curves for each toxin target. The data shown represent the highest affinity candidate aptamers for ancrod (ancrod-7, ancrod-11 and ancrod-55) (A) and batroxobin (batroxobin-5, batroxobin-26 and batroxobin-29) (B). Concentration-dependent data were generated at an excitation wavelength of 470 nm and an emission of 515 nm, and Kds were calculated via nonlinear regression analysis of the resulting saturation curves. Data represent means of duplicate measurements, and error bars represent standard deviations (SD).

Quantifying Aptamer Binding to Recombinant Toxins and Venom by ALISA
An ALISA was next used to quantify the binding levels of the two highest affinity aptamers (ancrod-55 and batroxobin-26) against each of the recombinant toxins and the native venoms from which the recombinant toxins were derived (i.e., C. rhodostoma and B. atrox venoms) in a reciprocal and comparative manner. The results of these binding assays, which were performed at two aptamer concentrations (0.4 and 2 nM), revealed that the two aptamers tested exhibited considerable cross-recognition and binding to the two SVSP toxins (i.e., irrespective of which toxin the aptamer was selected against) and the two homologous venoms ( Figure 3). Indeed, the aptamer batroxobin-26 exhibited highly comparable binding levels to batroxobin, ancrod and the two native snake venoms ( Figure  3B). These binding levels were, however, lower than those obtained with the ancrod-55 aptamer for all comparisons ( Figure 3A). The equilibrium K d s of the highest affinity ssDNA aptamers were derived from the fluorescence binding curves for each toxin target. The data shown represent the highest affinity candidate aptamers for ancrod (ancrod-7, ancrod-11 and ancrod-55) (A) and batroxobin (batroxobin-5, batroxobin-26 and batroxobin-29) (B). Concentration-dependent data were generated at an excitation wavelength of 470 nm and an emission of 515 nm, and K d s were calculated via nonlinear regression analysis of the resulting saturation curves. Data represent means of duplicate measurements, and error bars represent standard deviations (SD).

Quantifying Aptamer Binding to Recombinant Toxins and Venom by ALISA
An ALISA was next used to quantify the binding levels of the two highest affinity aptamers (ancrod-55 and batroxobin-26) against each of the recombinant toxins and the native venoms from which the recombinant toxins were derived (i.e., C. rhodostoma and B. atrox venoms) in a reciprocal and comparative manner. The results of these binding assays, which were performed at two aptamer concentrations (0.4 and 2 nM), revealed that the two aptamers tested exhibited considerable cross-recognition and binding to the two SVSP toxins (i.e., irrespective of which toxin the aptamer was selected against) and the two homologous venoms ( Figure 3). Indeed, the aptamer batroxobin-26 exhibited highly comparable binding levels to batroxobin, ancrod and the two native snake venoms ( Figure 3B). These binding levels were, however, lower than those obtained with the ancrod-55 aptamer for all comparisons ( Figure 3A). The data on the left show the binding levels obtained against the two recombinant toxins that the aptamers were selected against and the corresponding snake venoms that these toxins are derived from. The data on the right show binding levels obtained against a broad panel of 11 distinct snake venoms. All data represent the mean of duplicates measured at 450 nm, and error bars represent the standard deviation (SD) of the duplicate measurements.
Next, ALISA was used to explore the capability of each aptamer to bind to a panel of 11 distinct venoms that exhibited considerable variation in toxin composition and were sourced from taxonomically diverse snake species [77]. The binding patterns obtained from these experiments suggest that the two aptamers, each selected against a single SVSP toxin, recognise distinct snake venoms to a similar extent as the toxins and homologous snake venoms described above. Thus, at the 2 nM concentration tested, the aptamer an-  Figure 3B). While the range in binding levels observed suggests a considerable degree of variation in terms of venom recognition, it is possible that this might simply reflect the relative abundance of SVSP toxins present in these different venoms. . The data on the left show the binding levels obtained against the two recombinant toxins that the aptamers were selected against and the corresponding snake venoms that these toxins are derived from. The data on the right show binding levels obtained against a broad panel of 11 distinct snake venoms. All data represent the mean of duplicates measured at 450 nm, and error bars represent the standard deviation (SD) of the duplicate measurements.
Next, ALISA was used to explore the capability of each aptamer to bind to a panel of 11 distinct venoms that exhibited considerable variation in toxin composition and were sourced from taxonomically diverse snake species [77]. The binding patterns obtained from these experiments suggest that the two aptamers, each selected against a single SVSP toxin, recognise distinct snake venoms to a similar extent as the toxins and homologous snake venoms described above. Thus, at the 2 nM concentration tested, the aptamer ancrod-  Figure 3B). While the range in binding levels observed suggests a considerable degree of variation in terms of venom recognition, it is possible that this might simply reflect the relative abundance of SVSP toxins present in these different venoms.

Aptamers Inhibit Toxin-Induced Fibrinogenolysis and Fibrinogen Depletion
Many SVSP toxins exert a functional activity akin to thrombin by cleaving the αor/and β-chains of fibrinogen, typically resulting in the release of fibrinopeptides rather than fibrin and, thus, ultimately contributing to dysregulation of coagulation via depletion of fibrinogen [78,79]. To assess whether the binding exhibited by the highest affinity ssDNA aptamers was capable of inhibiting the functional activities known from SVSPs, we first used SDS-PAGE to visualise aptamer-mediated protection against fibrinogenolysis. Our findings revealed that, in line with the known functional activities of native ancrod and batroxobin in recent work [76], both recombinantly expressed toxins cleaved the α-chain of fibrinogen ( Figure 4A,B). We then assessed the inhibitory capability of the two candidate ssDNA aptamers to prevent this functional activity. Both ancrod-55 and batroxobin-26 aptamers protected against fibrinogenolysis caused by the corresponding recombinant toxin, resulting in the visualisation of the α-chain of fibrinogen on the gel ( Figure 4A,B). We next sought to explore whether these two aptamers could also protect against native venom, which might contain additional related (or unrelated) proteins that could contribute to global fibrinogenolysis. Despite this possibility, our data showed that the aptamer ancrod-55 inhibited the fibrinogenolytic effect of C. rhodostoma venom on both the αand β-chains, while batroxobin-26 aptamer prevented B. atrox venom-induced cleavage of the α-chain of fibrinogen ( Figure 4C,D). These findings suggest that the binding specificities obtained during aptamer selection translate into inhibitory capabilities.

Aptamers Inhibit Toxin-Induced Fibrinogenolysis and Fibrinogen Depletion
Many SVSP toxins exert a functional activity akin to thrombin by cleaving the αor/and β-chains of fibrinogen, typically resulting in the release of fibrinopeptides rather than fibrin and, thus, ultimately contributing to dysregulation of coagulation via depletion of fibrinogen [78,79]. To assess whether the binding exhibited by the highest affinity ssDNA aptamers was capable of inhibiting the functional activities known from SVSPs, we first used SDS-PAGE to visualise aptamer-mediated protection against fibrinogenolysis. Our findings revealed that, in line with the known functional activities of native ancrod and batroxobin in recent work [76], both recombinantly expressed toxins cleaved the α-chain of fibrinogen ( Figure 4A,B). We then assessed the inhibitory capability of the two candidate ssDNA aptamers to prevent this functional activity. Both ancrod-55 and batroxobin-26 aptamers protected against fibrinogenolysis caused by the corresponding recombinant toxin, resulting in the visualisation of the α-chain of fibrinogen on the gel ( Figure 4A,B). We next sought to explore whether these two aptamers could also protect against native venom, which might contain additional related (or unrelated) proteins that could contribute to global fibrinogenolysis. Despite this possibility, our data showed that the aptamer ancrod-55 inhibited the fibrinogenolytic effect of C. rhodostoma venom on both the α-and β-chains, while batroxobin-26 aptamer prevented B. atrox venom-induced cleavage of the α-chain of fibrinogen ( Figure 4C,D). These findings suggest that the binding specificities obtained during aptamer selection translate into inhibitory capabilities. The clotting time of diluted plasma with a standard thrombin concentration is inversely related to the fibrinogen concentration [80]. To further explore the inhibitory effect of the selected high-affinity aptamers against both SVSP toxins and native snake venom, we measured fibrinogen depletion using native venom spiking experiments with plateletpoor plasma (PPP). Our results revealed that both ancrod and batroxobin, as well as native venom from C. rhodostoma and B. atrox, stimulate fibrinogen consumption in platelet-poor plasma (PPP), resulting in reductions versus the thrombin-stimulated saline control (<1.5 The clotting time of diluted plasma with a standard thrombin concentration is inversely related to the fibrinogen concentration [80]. To further explore the inhibitory effect of the selected high-affinity aptamers against both SVSP toxins and native snake venom, we measured fibrinogen depletion using native venom spiking experiments with platelet-poor plasma (PPP). Our results revealed that both ancrod and batroxobin, as well as native venom from C. rhodostoma and B. atrox, stimulate fibrinogen consumption in platelet-poor plasma (PPP), resulting in reductions versus the thrombin-stimulated saline control (<1.5 vs. 3.42 g/L) ( Figure 5A). In both cases, however, the aptamers ancrod-55 and batroxobin-26 inhibited the depletion of fibrinogen stimulated by both toxins and native venoms to near control levels (2.7-2.9 vs. 3.42 g/L) ( Figure 5A). Crucially, the aptamers showed highly comparable inhibitory potencies with those obtained with gold-standard commercial antivenoms, as evidenced by the resulting equivalent fibrinogen concentrations observed ( Figure 5A). vs. 3.42 g/L) ( Figure 5A). In both cases, however, the aptamers ancrod-55 and batroxobin-26 inhibited the depletion of fibrinogen stimulated by both toxins and native venoms to near control levels (2.7-2.9 vs. 3.42 g/L) ( Figure 5A). Crucially, the aptamers showed highly comparable inhibitory potencies with those obtained with gold-standard commercial antivenoms, as evidenced by the resulting equivalent fibrinogen concentrations observed ( Figure 5A).

Aptamers Reduce Toxin-and Venom-Induced Prolongations of Clotting Times
The prothrombin time (PT) measures the time for citrated plasma to clot and specifically assesses the clotting capability of extrinsic and common coagulation cascades. A prolonged PT can result from an absence or deficiency of one of the clotting factors X, VII, V, II or fibrinogen. Our recombinant toxin/native venom spiking experiments demonstrated that both ancrod and batroxobin and the corresponding venoms from C. rhodostoma and B. atrox resulted in a prolongation of the PT in PPP compared to the saline control (i.e., >18 vs. ~14 s, respectively) ( Figure 5B). Noticeably, when co-incubated with the corresponding toxins/native venoms, the aptamers ancrod-55 and batroxobin-26 caused a substantial reduction in the PT compared with the toxin-only/native-venom-only samples, resulting in the PT approaching the level of control readings (14.5-15.5 vs. ~14.0 s, respectively) and, thus, demonstrating the inhibition of toxin activity ( Figure 5B). These reductions were highly comparable to those obtained with commercially available antivenoms (14-16 s).

Aptamers Reduce Toxin-and Venom-Induced Prolongations of Clotting Times
The prothrombin time (PT) measures the time for citrated plasma to clot and specifically assesses the clotting capability of extrinsic and common coagulation cascades. A prolonged PT can result from an absence or deficiency of one of the clotting factors X, VII, V, II or fibrinogen. Our recombinant toxin/native venom spiking experiments demonstrated that both ancrod and batroxobin and the corresponding venoms from C. rhodostoma and B. atrox resulted in a prolongation of the PT in PPP compared to the saline control (i.e., >18 vs.~14 s, respectively) ( Figure 5B). Noticeably, when co-incubated with the corresponding toxins/native venoms, the aptamers ancrod-55 and batroxobin-26 caused a substantial reduction in the PT compared with the toxin-only/native-venom-only samples, resulting in the PT approaching the level of control readings (14.5-15.5 vs.~14.0 s, respectively) and, thus, demonstrating the inhibition of toxin activity ( Figure 5B). These reductions were highly comparable to those obtained with commercially available antivenoms (14-16 s).
We employed the activated partial thromboplastin time (aPTT) in the same manner as the PT, except, in this case, to measure clotting in the context of the combined effect of the intrinsic and common coagulation pathways (i.e., Factors II, V, VIII, IX, X, XI, XII and fibrinogen). Ancrod, batroxobin, C. rhodostoma venom and B. atrox venom all caused a prolongation of the aPTT in PPP when compared with the negative control ( Figure 5C). The aptamers showed complete inhibition of toxin and native venom activity by reducing the aPTT to control levels, except for the ancrod and ancrod-55 aptamer combination, which showed substantial inhibition approaching control levels (ancrod, 57 s; ancrod + ancrod-55 aptamer, 48 s; normal control, 45 s) ( Figure 5C). However, the aptamers slightly outperformed the levels of neutralisation obtained with commercial antivenoms. These findings are noteworthy when considering the differences in doses of these therapeutics (aptamer, 1 pM; antivenom, 4.55 pM).

Discussion
Though current antivenom treatments are life-saving therapeutics, they possess several limitations that restrict their clinical efficacy and utility for tackling tropical snakebites. Consequently, a variety of novel experimental approaches are currently being employed to develop new therapeutic interventions that circumvent these existing limitations with snakebite treatment. Next-generation antivenoms consisting of monoclonal antibodies (mAbs) selected for their desirable specificities show much promise in this regard, particularly since different antibody formats can be adapted to specific pharmacodynamic and pharmacokinetic needs [81]. However, to date, little effort has been devoted to exploring antibody alternatives despite outstanding concerns relating to the potential cost of mAbbased snakebite treatments. In other fields, ssDNA aptamers have emerged as promising alternatives to antibodies for both therapeutic applications and their use as diagnostic tools [38][39][40], and several ssDNA aptamers have been successfully selected and applied against a wide range of targets, including small molecules [82] and toxins [37,[83][84][85][86]. Despite this prior research, ssDNA aptamers remain largely unexplored for their potential utility as snakebite therapeutics. Thus, in this proof-of-concept study, we sought to select ssDNA aptamers against fibrinogenolytic SVSP toxins found in medically important viper venoms and to evaluate their potential as novel toxin-inhibiting molecules via in vitro cross-reactivity and neutralisation studies.
Following fourteen cycles of SELEX selection, we selected two pools of aptamers directed against the recombinantly expressed SVSP toxins ancrod and batroxobin, which are found in native form in the venoms of the South East Asian pit viper C. rhodostoma and the South American lancehead pit viper B. atrox. Although the level of binding to each target increased steadily over SELEX selection cycles, by round 11, major increases in fluorescence intensity, indicating enrichment in binding to the target, occurred. Thereafter, counterselection was employed to deplete binders directed towards the target carrier (beads), and two further cycles of enrichment were employed, ultimately resulting in highly comparable fluorescence intensity patterns obtained during the selection process for the two different targets. Sequencing 100 targets from the resulting two pools of ssDNA aptamers yielded thirteen and nine unique aptamers directed against ancrod and batroxobin, respectively, and analyses of the resulting dissociation constants enabled the selection of the highest affinity aptamers against each, which ultimately exhibited desirable low nanomolar K d s (i.e., ancrod-55, 3.0 nM; batroxobin-26, 4.7 nM). These dissociation constants compare favourably with those obtained after ten rounds of diagnostic ssDNA aptamer selection against the neurotoxin β-bungarotoxin from the venom of B. multicinctus (>65.9 nm), for example [35].
Although ancrod and batroxobin are members of the same gene family (i.e., they are both snake venom serine proteases), they have different molecular masses (26)(27)(28)(29)(30)(31)(32)(33) and share only 67% amino acid sequence similarity [87]. Nonetheless, ALISA binding experiments demonstrated that the highest affinity aptamers selected and directed against one of these SVSP toxins were capable of cross-recognising and binding reciprocally to the other recombinant toxin under study here. Moreover, these experiments also revealed that both aptamers (i.e., ancrod-55 and batroxobin-26) were capable of binding to native toxins present in snake venoms that corresponded to the source of the recombinant toxins (i.e., C. rhodostoma and B. atrox venom) and, indeed, a broader variety of unrelated and distinct native snake venoms. Broadly speaking, each of the aptamers exhibited comparable moderate binding levels to each of the venoms tested, presumably via interaction with other SVSP toxins present within them, though ancrod-55 generally exhibited higher binding levels to these different venoms than batroxobin-26. However, we cannot rule out that nonspecific interactions between the aptamers and other venom toxins might be contributing to the overall binding levels obtained in these experiments.
Both ancrod and batroxobin are coagulopathic SVSP toxins capable of stimulating fibrinogen consumption and prolonging clotting times [13,15,16]. Indeed, most SVSPs are fibrinogenolytic enzymes that result in the depletion of circulating fibrinogen during envenoming. Using a panel of in vitro bioassays, we demonstrated that recombinant ancrod and batroxobin used for aptamer selection exhibit functional activities consistent with native toxins, including fibrinogenolytic activity via cleavage of the α-chain of fibrinogen, the reduction of fibrinogen plasma concentrations, and the prolongation of clotting times in PT and aPTT assays. Crucially, the coincubation of the ssDNA aptamers with the corresponding recombinant toxins abolished the cleavage of fibrinogen and reduced the consumption of fibrinogen in human plasma samples to control levels, demonstrating that the binding observed translates into a degree of functional inhibition. Similar results were observed when measuring the PT and aPTT in toxin-spiked platelet-poor plasma, with the homologous aptamers reducing the prolongation of clotting times stimulated by each toxin to near control levels. The inhibition also extended to native snake venoms, with the aptamers inhibiting their fibrinogenolytic and delayed clotting time effects in a highly comparable manner to those obtained with the recombinant toxins. These findings suggest that SVSP toxins likely dominate the fibrinogenolytic effect of C. rhodostoma and B. atrox venoms and are effectively inhibited by the ssDNA aptamers ancrod-55 and batroxobin-26, though we cannot rule out some contributory interactions between the aptamers and distinct coagulopathic toxins found in these venoms. Irrespectively, it is worth noting that the toxin inhibitory activities demonstrated here by the selected high-affinity aptamers are potent, with low concentrations (1 pM) inhibiting toxin and venom activities observed across the various assays, which compares very favourably to the doses required by the conventional polyclonal antivenoms to effect comparable neutralisation (4.55 pM).
Future work should explore whether the in vitro inhibitory data obtained here might translate into in vivo protection against venom-induced toxicity in preclinical animal models or whether toxin-family-directed aptamers, such as those identified here, need to be combined with other inhibitory molecules (i.e., other aptamers or different inhibitory molecules such as mAbs or small molecule drugs) targeting distinct toxin families to prevent severe envenoming pathology in vivo [6]. This latter point seems highly likely, given that snake venoms contain multiple isoforms from different toxin families that often work in concert to cause pathological effects [6,7]. There is, therefore, a clear trade-off between the toxin-specificity of an inhibitory molecule and the number of inhibitors (whether they be aptamers or mAbs) required to neutralise whole venom, where increased inhibitory breadth against additional toxin isoforms is likely to correspond to some degree of reduction in affinity and/or inhibitory potency against broader targets. Further, increasing the breadth of toxin recognition seems likely to increase the risk of off-target interactions; thus, it is important that aptamer-based therapeutics are also assessed for their potential to interact with cellular receptors or cognate proteins in vitro and evaluated for their potential toxicity in vivo in relevant preclinical models. Such experimental approaches will be key to understanding the number of distinct, well-tolerated aptamers required to neutralise a breadth of snake venoms.
Nonetheless, the proof-of-concept findings described here reveal that rationally selected ssDNA aptamers hold great value for exploration as novel inhibitory molecules capable of broadly inhibiting snake venom toxins such as SVSPs. However, the therapeutic potential of aptamers does not come without several challenges. One major limitation is that ssDNA aptamers are susceptible to nucleases and this, combined with their small size, results in short half-lives in vivo, which can be as low as two minutes due to DNA degradation and rapid clearance [49]. However, aptamers also possess great flexibility and can be modified to increase their stability and their resulting half-life, with conjugation to a higher molecular mass protein carrier, for example, liposomes, cholesterol or PEG; this has previously been demonstrated to successfully decrease the clearance rate from plasma to acceptable levels, though, it adds to cost of goods [57,59,88]. Other limitations include the initial selection process often being labour-intensive and time-consuming, though the downstream production of the resulting identified ssDNA aptamers offers several advantages over mAb production, for example, due to their size and chemical synthesis. Finally, to date, it is worth re-emphasising that only a single aptamer has passed through an FDA-approval process (Macugen ® (Pegaptanib)) [62], and as such, regulatory challenges are likely to remain for exploiting ssDNA aptamers as new therapeutics.
Despite these limitations, aptamers possess several clear advantages in terms of therapeutic characteristics over the polyclonal antibodies currently used for treating snakebite. For example, aptamers are chemically synthesised and, thus, do not require the use of ethically and financially costly experimental animals, and the resulting product exhibits a relative lack of batch-to-batch variation, unlike biologic antivenom [89]. Further, aptamers exhibit desirable stability and are likely more resistant to harsh conditions than antibodies [90], which may be of particular benefit in cold-chain-unstable parts of the tropics that suffer a high burden of snakebite. Finally, the large-scale cost of manufacturing aptamers seems likely to be relatively inexpensive compared with existing polyclonal antivenoms. Currently, manufacturing costs for aptamers range from USD 140-280 per dose [88], and while the dose efficacy of aptamers in the context of snakebite remains to be established, it seems likely that aptamer-based treatments could reduce the high costs currently borne by tropical snakebite victims, which often exceed USD 1000 in sub-Saharan Africa, for example [31]. Although much work needs to be undertaken before such benefits can be realised, the proof-of-concept data described herein add further weight to the potential utility of aptamers as future snakebite treatments and strongly advocate for their continued exploration as a novel therapeutic modality for inhibiting pathogenic snake venom toxins.

DNA Library and Primer Design
The DNA library and PCR primers were chemically synthesised and purified by Integrated DNA Technologies, Inc. The random ssDNA library used in the first SELEX cycle consisted of 2 × 10 17 sequences, specifically containing a central random region of 40 nucleotides flanked by two primers consisting of 16 nucleotides at the 5 and 3 ends (5 -TCCCTACGGCGCTAAC-N40-GTTGTAGCACGGTGGC-3 ) for amplification of the library sequence. To facilitate the separation of ssDNA from amplified dsDNA PCR products and to quantify DNA during selection, the forward and reverse primers were designed and modified with fluorescein and a HEG linker (the spacer is designed to block the PCR extension step during amplification). The modified PCR forward primer used during the selection cycles was 5 -6-FAM-TCCCTACGGCGCTAAC-3 and the reverse primer was (5 -poly dA20-HEG-spacer-GTTGTAGCACGGTGGC-3 ). An unmodified primer set was used for PCR amplification and cloning when SELEX rounds were completed and enriched.

Target Conjugation
The recombinant toxins used in this study, ancrod from C. rhodostoma (UniProt: P26324) and batroxobin from B. atrox (UniProt: P04971), were previously expressed in HEK293 mammalian cell lines [76]. In this study, we used these toxins to conjugate to NHS (Nhydroxysuccinimide)-activated Sepharose ® 4 Fast Flow beads (Sigma-Aldrich, Dorset, UK) to generate ssDNA aptamers as novel toxin-inhibiting therapeutics. To do so, the manufacturer's protocol was followed. First, a 300 µg/mL solution of each toxin was prepared in coupling buffer (0.1 M NaHCO 3 , 0.5 M NaCl, pH 8.3). Next, 2 mL of NHSactivated Sepharose beads were rinsed multiple times with 1 mM HCl for 15 min to remove the additives and preserve the activity of the reactive groups. Each venom toxin was then added immediately to the washed Sepharose beads independently (1:1 v:v ratio) and mixed by end-over-end rotation overnight at 4 • C. After incubation, the beads were centrifuged at 14× g for 5 min, the supernatant discarded, and the samples washed with 1 mL of 1 M NaCl. This centrifugation step was repeated before 2 mL of 1 M NaCl was added, and the samples were incubated with rotation for 2 h at 4 • C to block the unreacted amine groups on the beads. After blocking, the toxin-conjugated Sepharose was washed three times with 4 mL in an alternating manner, first with 0.1 M acetic acid/sodium acetate, pH 4, containing 0.5 M NaCl 0.1 M, followed by Tris-HCl buffer, pH 8, containing 0.5 M NaCl. Thereafter, the toxin-conjugated Sepharose was stored in 50 mM Tris-HCl buffer (pH 7.5), 0.05% sodium azide at 4 • C until use.

PCR Amplification
Each of the recombinant toxin-conjugated Sepharose mixtures was washed five times by adding 500 µL binding buffer (50 mM Tris-HCl, pH 7.5, 150 mM NaCl, 2 mM MgCl 2 ) to 100 µL of Sepharose using centrifugal tube filters (0.45 µm) (Costar, Washington, DC, USA). In the first SELEX cycle, 50 µL (3 nM) of the DNA library solution in a 450 µL binding buffer was used, whereas 50 pM of the DNA elution from each round was used for subsequent cycles. The DNA and binding buffer mixture were heated to 90 • C on a heat block for five minutes, cooled at 4 • C for ten minutes and then held at room temperature (RT) for five minutes. Next, 300 µL of DNA and binding buffer mixture was added to the washed beads and incubated end-over-end for 1 h at RT. After incubation, the mixtures were washed and centrifuged five times with 500 µL of binding buffer at 82× g for 2 min to ensure the elimination of unbound DNA. The flow-through from the first wash was retained for use as a positive control in downstream 2% agarose gel electrophoresis and PAGE experiments. DNA bound to the toxin-conjugated beads was collected from each sample using elution buffer (7 M urea in binding buffer) and incubated at 90 • C for 10 min on a heat block. The elution step has performed a total of five times via the addition of 300 µL elution buffer until no fluorescence was detected in the eluted aliquots, indicating the complete elution of bound DNA from the beads. The eluted aliquots were then desalted with filtered water and concentrated using an ultrafiltration device with a 3 kDa cut-off membrane (Amicon Ultra-0.5 Centrifugal Filter Unit; Merck, Dorset, UK) to remove any urea and salt residues.
A dsDNA product with different lengths for each strand was generated via PCR amplification by using two parallel 50 µL reactions (one for each toxin target), each containing five units of Taq Plus and polymerase buffer (New England, BioLabs ® , Hitchin, UK), 2 mM MgCl 2 , 200 µM deoxyribonucleotide triphosphate (dNTP) and 0.2 µM of the unmodified forward and reverse primers. The PCR conditions were: 94 • C for 10 min, followed by 17 cycles of 94 • C for 1 min, 54 • C for 1 min and 72 • C for 1 min, and a final extension step of 10 min at 72 • C. To confirm the amplification of the dsDNA product, the samples were then analysed electrophoretically using 2% agarose gel electrophoresis.
An ssDNA product with different lengths for each strand was also generated via PCR amplification using 20 parallel 50 µL reactions using the same conditions as outlined above, except with the use of 0.2 µM of the labelled forward primer. Following amplification, the 20 reactions were pooled and aliquoted into three Eppendorf tubes, each containing 250 µL of pooled ssDNA product, 50 µL 3 M sodium acetate, pH 5.2 and 750 µL absolute ethanol.
The tubes were then mixed and stored at −80 • C for 2 h before centrifugation at 13,000× g for 30 min at 4 • C. The supernatant was discarded, and the remaining pellet was dried at 90 • C for a couple of minutes on a heat block before resuspension in 400 µL of water and formamide (1:1 v/v) and incubation at 90 • C for 5 min.

Separation of ssDNA by Denaturing PAGE
We employed a 10% denaturing PAGE to separate the resulting ssDNA from the double-stranded (dsDNA) amplified PCR products. Fifty microlitres of each dsDNA sample were independently loaded into the gel and run at 200 volts for 55 min using a Mini-PROTEAN electrophoresis system (Bio-Rad, Hertfordshire, UK). Next, the ssDNA band was observed under UV light, dissected, cut and transferred to a 15 mL falcon tube containing 4 mL of Tris-ethylenediaminetetraacetic acid (EDTA) buffer solution (TE buffer; 10 mM Tris-HCl, pH 7.4, 1 mM EDTA). Samples were then subjected to freeze-thaw for 30 min at −80 • C, 10 min at 90 • C, followed by incubation on a rotary shaker overnight at 37 • C [91]. Finally, ssDNA was eluted and concentrated using an ultrafiltration device with a 3 kDa cut-off membrane (Amicon Ultra-0.5 Centrifugal Filter Unit; Merck, Dorset, UK), as described earlier. The UV and fluorescence intensity of each dsDNA and ssDNA were then measured for each eluted fraction using a NanoDrop 2000 C spectrophotometer and a FLUOstar Omega 384 well black assay plate reader (B.M.G. LabTech), respectively. The fluorescence spectra were detected at an excitation wavelength of 470 nm and an emission of 515 nm.

Counter-Selection
Following enrichment, counter-selection (negative selection) was employed using blank Sepharose beads instead of those conjugated to toxins. The objective was to increase enrichment to the toxin target via the elimination of aptamers directed towards the beads. Sepharose beads were subjected to the aforementioned conjugation process with 0.1 M Tris-HCl buffer, pH 8.3, containing 0.5 M NaCl only (i.e., in the absence of toxins), and counter-selection followed the same experimental strategy described above and was used to filter out negative targets.

Cloning and Sequencing ssDNA Aptamers
After ssDNA recovery, the eluted ssDNA from the 14th round of selection was amplified using the unlabelled PCR primer (5 -TCCCTACGGCGCTAAC-3 ), following the same conditions as PCR amplification: 94 • C for 10 min, followed by 17 cycles of 94 • C for 1 min, 54 • C for 1 min and 72 • C for 1 min, and a final extension step of 10 min at 72 • C. Next, ligation was performed using the TOPO TA Cloning Kit for Sequencing and One Shot TOP10 Chemically Competent E. coli (Thermo-Fisher Scientific, Altrincham, UK). Ligation reactions consisted of 2 µL of fresh PCR product, 1 µL of salt solution, 1 µL of TOPO ® vector and 1 µL of nuclease-free water, with the resulting samples mixed gently and incubated for 30 min at RT (20-25 • C). For cloning, 5 µL of the ligated product was added to a vial of One Shot TOP10 Chemically Competent E. coli (Thermo-Fisher Scientific, Altrincham, UK) and mixed gently. The cells were then placed on ice for 30 min, heat-shock-treated for 30 s at 42 • C and transferred back to ice. Next, 250 µL of pre-warmed S.O.C medium was added, and samples were incubated at 37 • C for 1 h with horizontal shaking at 200 rpm. Thereafter, 50 µL of the medium was spread on 25 mL of Lauria Broth (LB) agar (VWR, Leicestershire, UK) plates supplemented with 25 µL ampicillin (20 mg/mL) (Sigma-Aldrich, Dorset, UK), 250 µL 5-bromo-4-chloro-3-indolyl-β-D-galactoside (X-Gal) (20 mg/mL) (Sigma-Aldrich, Dorset, UK) and IPTG 25 µL (100 mM) (Sigma-Aldrich, Dorset, UK) and incubated at 37 • C overnight. Following incubation, ten random positive colonies (white colonies) were picked for colony PCR analysis: 5 µL 10 × PCR buffer, 0.5 µL dNTPs mix (50 mM), 0.5 µL of M13 forward primer (5 -GTAAAACGACGGCCAGT-3 ) and M13 reverse primer (5 -CAGGAAACAGCTATGAC-3 ) (0.1 µg/µL each), 41.5 µL nuclease-free water and 1 µL Taq polymerase (1 unit/µL). Colony PCR was performed using the following cycling pa-rameters: initial denaturation at 94 • C for 2 min, followed by 25 cycles of denaturation at 94 • C for 1 min, annealing at 55 • C for 1 min and extension at 72 • C for 1 min, with a final extension step at 72 • C for 7 min. To confirm insert amplification, 2% agarose gel electrophoresis was performed.
One hundred positive bacterial colonies were sterilely selected and inoculated with 5 mL of LB medium containing 50 µL ampicillin (20 mg/mL) and incubated overnight at 37 • C with shaking at 200 rpm. Next, plasmid DNA was purified using a QIAprep Spin Miniprep Kit (Qiagen, Düsseldorf, Germany). The overnight culture was centrifuged at 6000× g for 10 min at RT. The pellet was then resuspended in 250 µL Buffer P1 (with Rnase A) and transferred to a microcentrifuge tube. Then, 250 µL Buffer P2 was added and mixed thoroughly by inverting the tube gently 4-6 times. After that, 350 µL Buffer N3 was added and mixed immediately and thoroughly by inverting the tube 4-6 times. Next, samples were centrifuged for 10 min at 17,900× g, and the supernatant was transferred to QIAprep 2.0 Spin Columns before brief centrifugation at 2817× g for one minute with the eluate discarded. Columns were then washed with 750 µL Buffer PE via centrifugation at 2817× g for one minute, the eluate discarded, and columns centrifuged for an additional minute to remove any residual buffer. Finally, 50 µL Buffer EB (10 mM Tris-HCl, pH 8.5) was added to each column and incubated at RT for two minutes prior to elution via centrifugation at 2817× g for 1 min. The resulting PCR products were then sequenced commercially via Sanger sequencing (SourceBioScience, Nottingham, UK), and the resulting ssDNA aptamer sequences were analysed and aligned using PRALINE [92].

Dissociation Constants
The binding affinities of the identified ssDNA aptamers against each of the corresponding toxins used for selection were evaluated using a fluorescence-based assay. Both forward and reverse primer-binding sites were truncated and labelled with 6-FAM (fluoresceinlabelled) before synthesis. Different concentrations of fluorescein-labelled ssDNA aptamer (0, 5, 10, 50, 100, 200, 300 and 400 nM) were incubated with a constant amount of the corresponding toxin-conjugated beads (50 µL), as described in the selection protocol (by heating to 90 • C for 10 min, 4 • C for 10 min and then at RT for 5 min). After 1 h of incubation at RT with rotation, the unbound ssDNA aptamers were washed using 3 kDa centrifuge filters and 500 µL binding buffer, and the toxin-bound ssDNA aptamers were eluted five times via the addition of 300 µL elution buffer (7 M urea in binding buffer) and incubation at 90 • C for 10 min on a heat block. The five eluted aliquots were then concentrated using an ultrafiltration device with a 3 kDa cut-off membrane (Amicon Ultra-0.5 Centrifugal Filter Unit; Merck, Dorset, UK) and desalted with filtered water to remove urea and salts. The fluorescence intensity of each ssDNA was then measured for each eluted fraction at an excitation wavelength of 470 nm and an emission of 515 nm. The K d of each ssDNA aptamer was calculated via nonlinear regression analysis of the resulting plotted hyperbolic curves using Prism v8 software (GraphPad).

Binding by ALISA
To quantify reciprocal binding to the recombinant toxins (ancrod and batroxobin) by the highest affinity aptamers (determined by the K d analysis above), we used an ALISA assay. We also assessed binding to native venom samples that corresponded to the source of each toxin (i.e., C. rhodostoma, captive-bred and B. atrox, Brazil) and an additional panel of eleven geographically and taxonomically diverse snake species (Bitis arietans, Nigeria; Bothrops asper, Costa Rica; Bothrops jararaca, Brazil; Crotalus atrox, captive-bred; Daboia russelii, Sri Lanka; Deinakgistrodon acutus, captive-bred; Dispholidus typus, South Africa; Echis carinatus, India; Echis ocellatus, Nigeria; Rhabdophis subminiatus, Hong Kong; Trimeresurus albolabris, captive-bred) to assess the breadth of binding. Many of these venoms were collected from animals maintained under controlled environmental and dietary conditions in the UK-Home-Office-licensed and -inspected herpetarium of the Centre for Snakebite Research and Interventions (CSRI) at the Liverpool School of Tropical Medicine (LSTM). The

Sample Preparation for Clotting Profiling Experiments
Blood samples for clotting profiling experiments were obtained according to ethically approved protocols (LSTM research tissue bank, REC ref. 11/H1002/9) from consenting healthy volunteers who confirmed they had not taken any anticoagulant treatments for at least three months prior to blood collection. Blood samples were collected in tubes containing acid citrate dextrose adenine (ACD-A) as an anticoagulant. To prepare fresh frozen plasma (FFP) a fresh blood sample was centrifuged at 2500× g at 20-25 • C for 10 min, and the supernatant was retained and stored at −80 • C until use. Platelet-poor plasma (PPP) was then generated by centrifuging the FFP samples again in an identical manner, with the resulting platelet-depleted supernatant stored at −80 • C until further use.
In addition to the samples prepared above, we also used commercially sourced samples (Diagnostic Reagents Ltd., Thame, UK) to serve as quality control checks for the various reagents employed. To this end, normal and abnormal test controls were implemented for the experiments described below prior to experimentation. These consisted of (1) Diagen control plasma for fibrinogen: normal (RCPN070) and abnormal (RCPA080) and (2) control plasma for aPTT and PT: normal (IQCN130), abnormal 1 (Mild, IQCM140) and abnormal 2 (Severe, IQCS150). All QC samples were aliquoted into 500 µL, stored at −80 • C, and defrosted in a water bath for 5 min at 37 • C immediately prior to use. To act as positive controls for the aptamers, we used commercial antivenoms, specifically the Thai Red Cross monovalent equine Malayan pit viper antivenom (Lot #: CR00316, expiry date 06/2021) for ancrod and C. rhodostoma experiments and the Instituto Butantan polyvalent equine SORO antibotropico/crotalico antivenom (Lot #: 1012308, expiry date: 2013) for batroxobin and B. atrox experiments. These freeze-dried antivenoms were reconstituted in the pharmaceuticalgrade water supplied by the manufacturers, the protein concentrations quantified using a NanoDrop (Thermo-Fisher Scientific, Dorset, UK), and then stored at 4 • C in the short term until use.

Fibrinogen Consumption via the Clauss Method
The Clauss method is a quantitative, clot-based assay. It measures the ability of thrombin to convert fibrinogen to fibrin clots, followed by manual time measurements of clotting [80]. Here, we applied this method in toxin-and venom-spiking experiments to assess the inhibitory capability of the selected ssDNA aptamers against the depletion of fibrinogen. Twenty microlitres of PPP were spiked with 0.6 ng of ancrod, batroxobin, C. rhodostoma or B. atrox venom, or 0.9% saline solution as the negative control. All samples were also repeated in the presence of 1 pM of the selected ssDNA aptamers or 0.5 µg of the commercial antivenoms following preincubation at 37 • C for 5 min. Samples were then diluted tenfold with 0.02 M imidazole buffer (pH 7.35), transferred to glass test tubes (10 × 75 mm), and warmed at 37 • C for 120 s. Thereafter, 100 µL of thrombin reagent (20 units/mL) (Diagnostic Reagents Ltd., Thame, UK) was added and time measurements commenced. Tubes were gently tilted at regular intervals (returning to the water bath between tilting), and the time for the formation of a clot was recorded. All experiments were performed in duplicate.

Prothrombin Time (PT)
To measure the inhibitory capability of selected aptamers against toxins acting on the extrinsic coagulation pathway, we quantified differences in the PT between toxin and venom samples in the presence and absence of aptamers. Measurements of PT were undertaken by first adding 100 µL of calcium rabbit brain thromboplastin (Diagnostic Reagents Ltd., Thame, UK) to a glass test tube (10 × 75 mm) and incubating at 37 • C for 60-120 s in a water bath. Next, 50 µL of PPP was spiked with 50 µL containing 0.6 ng of ancrod, batroxobin, C. rhodostoma or B. atrox venom, or saline solution as the negative control. All samples were also repeated in the presence of 1 pM of the selected ssDNA aptamers or 0.5 µg of the commercial antivenoms, following preincubation at 37 • C for 5 min. Time measurements were commenced, tubes were gently tilted at regular intervals (returning to the water bath between tilting), and the time for the formation of a clot was recorded. All experiments were performed in duplicate.

Activated Partial Thromboplastin Time (aPTT)
To measure the inhibitory capability of selected aptamers against toxins acting on the intrinsic and common coagulation pathways, we quantified differences in the aPTT between toxin and venom samples in the presence and absence of aptamers. To do so, 50 µL of Micronised Silica/Platelet Substitute Mixture (Diagnostic Reagents Ltd., Thame, UK) was placed in glass test tubes (10 × 75 mm) in a water bath at 37 • C and incubated for 60-120 s. Next, 50 µL of PPP was spiked with 0.6 ng of ancrod, batroxobin, C. rhodostoma venom or B. atrox venom, or saline solution as the negative control. All samples were also repeated in the presence of 1 pM of the selected ssDNA aptamers or 0.5 µg of the commercial antivenoms, following preincubation at 37 • C for 5 min. Thereafter, samples were added to the glass test tubes and gently tilted at regular intervals for precisely five minutes at 37 • C. Finally, 50 µL of 25 mM calcium chloride (preincubated at 37 • C) was added to each tube, and the tube was gently tilted until the resulting clot time was recorded. All experiments were performed in duplicate.

Data Availability Statement:
The data presented in this study are available in this article.