Sequence and phylogenetic analysis of viper venom serine proteases

Snakebites are a major neglected tropical disease responsible for as many as 95000 deaths every year worldwide. Viper venom serine proteases disrupt haemostasis of prey and victims by affecting various stages of the blood coagulation system. A better understanding of their sequence, structure, function and phylogenetic relationships will improve the knowledge on the pathological conditions and aid in the development of novel therapeutics for treating snakebites. A large dataset for all available viper venom serine proteases was developed and analysed to study various features of these enzymes. Despite the large number of venom serine protease sequences available, only a small proportion of these have been functionally characterised. Although, they share some of the common features such as a C-terminal extension, GWG motif and disulphide linkages, they vary widely between each other in features such as isoelectric points, potential N-glycosylation sites and functional characteristics. Some of the serine proteases contain substitutions for one or more of the critical residues in catalytic triad or primary specificity pockets. Phylogenetic analysis clustered all the sequences in three major groups. The sequences with substitutions in catalytic triad or specificity pocket clustered together in separate groups. Our study provides the most complete information on viper venom serine proteases to date and improves the current knowledge on the sequence, structure, function and phylogenetic relationships of these enzymes. This collective analysis of venom serine proteases will help in understanding the complexity of envenomation and potential therapeutic avenues.

Likelihood-free estimation of model evidence

Statistical methods of inference typically require the likelihood function to be computable in a reasonable amount of time. The class of “likelihood-free” methods termed Approximate Bayesian Computation (ABC) is able to eliminate this requirement, replacing the evaluation of the likelihood with simulation from it. Likelihood-free methods have gained in efficiency and popularity in the past few years, following their integration with Markov Chain Monte Carlo (MCMC) and Sequential Monte Carlo (SMC) in order to better explore the parameter space. They have been applied primarily to estimating the parameters of a given model, but can also be used to compare models. Here we present novel likelihood-free approaches to model comparison, based upon the independent estimation of the evidence of each model under study. Key advantages of these approaches over previous techniques are that they allow the exploitation of MCMC or SMC algorithms for exploring the parameter space, and that they do not require a sampler able to mix between models. We validate the proposed methods using a simple exponential family problem before providing a realistic problem from human population genetics: the comparison of different demographic models based upon genetic data from the Y chromosome.

Bayesian parameter estimation for latent Markov random fields and social networks

Undirected graphical models are widely used in statistics, physics and machine vision. However Bayesian parameter estimation for undirected models is extremely challenging, since evaluation of the posterior typically involves the calculation of an intractable normalising constant. This problem has received much attention, but very little of this has focussed on the important practical case where the data consists of noisy or incomplete observations of the underlying hidden structure. This paper specifically addresses this problem, comparing two alternative methodologies. In the first of these approaches particle Markov chain Monte Carlo (Andrieu et al., 2010) is used to efficiently explore the parameter space, combined with the exchange algorithm (Murray et al., 2006) for avoiding the calculation of the intractable normalising constant (a proof showing that this combination targets the correct distribution in found in a supplementary appendix online). This approach is compared with approximate Bayesian computation (Pritchard et al., 1999). Applications to estimating the parameters of Ising models and exponential random graphs from noisy data are presented. Each algorithm used in the paper targets an approximation to the true posterior due to the use of MCMC to simulate from the latent graphical model, in lieu of being able to do this exactly in general. The supplementary appendix also describes the nature of the resulting approximation.

Phylogenetic relationships withinPelargonium sect. Peristera (Geraniaceae) inferred from nrDNA and cpDNA sequence comparisons

Phylogenetic analysis of nrDNA ITS and trnL (UAA) 5′ exon-trnF (GAA) chloroplast DNA sequences from 17 species ofPelargonium sect.Peristera, together with nine putative outgroups, suggests paraphyly for the section and a close relationship between the highly disjunct South African and Australian species of sect.Peristera. Representatives fromPelargonium sectt.Reniformia, Ligularia s. l. andIsopetalum (the St. Helena endemicP. cotyledonis) appear to be nested within thePeristera clade. The close relationship between the South African and AustralianPeristera is interpreted as being caused by long-range dispersal to Australia, probably as recent as the late Pliocene.

Synthesis, crystal structures, magnetic properties and catecholase activity of double phenoxido-bridged penta-coordinated dinuclear nickel(II) complexes derived from reduced Schiff-base ligands: mechanistic inference of catecholase activity

Three double phenoxido-bridged dinuclear nickel(II) complexes, namely [Ni-2(L-1)(2)(NCS)(2)] (1), [Ni-2(L-2)(2)(NCS)(2)] (2), and [Ni-2(L-3)(2)(NCS)(2)] (3) have been synthesized using the reduced tridentate Schiff-base ligands 2-[1-(3-methylamino-propylamino)-ethyl]-phenol (HL1), 2-[1-(2-dimethylamino-ethylamino)-ethyl]-phenol (HL2), and 2-[1-(3-dimethylarnino-propylamino)-ethyl]-phenol (HL3), respectively. The coordination compounds have been characterized by X-ray structural analyses, magnetic-susceptibility measurements, and various spectroscopic methods. In all complexes, the nickel(II) ions are penta-coordinated in a square-pyramidal environment, which is severely distorted in the case of 1 (Addison parameter tau = 0.47) and 3 (tau = 0.29), while it is almost perfect for 2 (tau = 0.03). This arrangement leads to relatively strong antiferromagnetic interactions between the Ni(II) (S = 1) metal centers as mediated by double phenoxido bridges (with J values of -23.32 (1), -35.45 (2), and -34.02 (3) cm(3) K mol(-1), in the convention H = -2JS(1)S(2)). The catalytic activity of these Ni compounds has been investigated for the aerial oxidation of 3,5-di-tert-butylcatechol. Kinetic data analysis following Michaelis-Menten treatment reveals that the catecholase activity of the complexes is influenced by the flexibility of the ligand and also by the geometry around the metal ion. Electrospray ionization mass spectroscopy (ESI-MS) studies (in the positive mode) have been performed for all the coordination compounds in the presence of 3,5-DTBC to characterize potential complex-substrate intermediates. The mass-spectrometry data, corroborated by electron paramagnetic resonance (EPR) measurements, suggest that the metal centers are involved in the catecholase activity exhibited by the complexes.

Excise taxes and commodity promotion bayesian retrieval of the optimum

This article shows how the solution to the promotion problem—the problem of locating the optimal level of advertising in a downstream market—can be derived simply, empirically, and robustly through the application of some simple calculus and Bayesian econometrics. We derive the complete distribution of the level of promotion that maximizes producer surplus and generate recommendations about patterns as well as levels of expenditure that increase net returns. The theory and methods are applied to quarterly series (1978:2S1988:4) on red meats promotion by the Australian Meat and Live-Stock Corporation. A slightly different pattern of expenditure would have profited lamb producers

Bayesian estimation of the double hurdle model in the presence of fixed costs

We present a model of market participation in which the presence of non-negligible fixed costs leads to random censoring of the traditional double-hurdle model. Fixed costs arise when household resources must be devoted a priori to the decision to participate in the market. These costs, usually of time, are manifested in non-negligible minimum-efficient supplies and supply correspondence that requires modification of the traditional Tobit regression. The costs also complicate econometric estimation of household behavior. These complications are overcome by application of the Gibbs sampler. The algorithm thus derived provides robust estimates of the fixed-costs, double-hurdle model. The model and procedures are demonstrated in an application to milk market participation in the Ethiopian highlands.

Bayesian ranking and slection of fishing boat efficiencies

The steadily accumulating literature on technical efficiency in fisheries attests to the importance of efficiency as an indicator of fleet condition and as an object of management concern. In this paper, we extend previous work by presenting a Bayesian hierarchical approach that yields both efficiency estimates and, as a byproduct of the estimation algorithm, probabilistic rankings of the relative technical efficiencies of fishing boats. The estimation algorithm is based on recent advances in Markov Chain Monte Carlo (MCMC) methods— Gibbs sampling, in particular—which have not been widely used in fisheries economics. We apply the method to a sample of 10,865 boat trips in the US Pacific hake (or whiting) fishery during 1987–2003. We uncover systematic differences between efficiency rankings based on sample mean efficiency estimates and those that exploit the full posterior distributions of boat efficiencies to estimate the probability that a given boat has the highest true mean efficiency.

Temporal genetic variation of the red fox, Vulpes vulpes, across western Europe and the British Isles

Quaternary climatic fluctuations have had profound effects on the phylogeographic structure of many species. Classically, species were thought to have become isolated in peninsular refugia, but there is limited evidence that large, non-polar species survived outside traditional refugial areas. We examined the phylogeographic structure of the red fox (Vulpes vulpes), a species that shows high ecological adaptability in the western Palaearctic region. We compared mitochondrial DNA sequences (cytochrome b and control region) from 399 modern and 31 ancient individuals from across Europe. Our objective was to test whether red foxes colonised the British Isles from mainland Europe in the late Pleistocene, or whether there is evidence that they persisted in the region through the Last Glacial Maximum. We found red foxes to show a high degree of phylogeographic structuring across Europe and, consistent with palaeontological and ancient DNA evidence, confirmed via phylogenetic indicators that red foxes were persistent in areas outside peninsular refugia during the last ice age. Bayesian analyses and tests of neutrality indicated population expansion. We conclude that there is evidence that red foxes from the British Isles derived from central European populations that became isolated after the closure of the landbridge with Europe.