A generalized Markov sampler

A recent development of the Markov chain Monte Carlo (MCMC) technique is the emergence of MCMC samplers that allow transitions between different models. Such samplers make possible a range of computational tasks involving models, including model selection, model evaluation, model averaging and hypothesis testing. An example of this type of sampler is the reversible jump MCMC sampler, which is a generalization of the Metropolis-Hastings algorithm. Here, we present a new MCMC sampler of this type. The new sampler is a generalization of the Gibbs sampler, but somewhat surprisingly, it also turns out to encompass as particular cases all of the well-known MCMC samplers, including those of Metropolis, Barker, and Hastings. Moreover, the new sampler generalizes the reversible jump MCMC. It therefore appears to be a very general framework for MCMC sampling. This paper describes the new sampler and illustrates its use in three applications in Computational Biology, specifically determination of consensus sequences, phylogenetic inference and delineation of isochores via multiple change-point analysis.

Reduction of β-Endorphin-Containing Immune Cells in Inflamed Paw Tissue Corresponds with a Reduction in Immune-Derived Antinociception: Reversible by Donor Activated Lymphocytes

The functional integrity of the immune system is essential for peripheral antinociception. Previous studies have demonstrated that immune cells elicit potent antinociception in inflamed tissues and that corticotropin-releasing factor-induced antinociception is significantly inhibited in animals that have undergone cyclosporin A (CsA)-induced immunosuppression. In this study, we examined the effect of a single bolus of CsA on inflammatory nociception. CsA-treated rats had substantially increased nociception compared with nonimmunosuppressed rats, consistent with a reduction in circulating and infiltrating lymphocytes. Furthermore, CsA-treated rats had inhibition of corticotropin-releasing factor-induced immune-derived antinociception, which was dose-dependently reversed by IV injection of concanavalin A-activated donor lymphocytes (1.0-7.0 X 10(6) cells/0.1 mL). In conclusion, our findings provided further evidence that opioid-containing immune cells are essential for peripheral analgesia. It is evident from these findings that control of inflammatory pain relies heavily on a functioning immune system.

Bayesian estimation of g-and-k distributions using MCMC

In this paper we investigate a Bayesian procedure for the estimation of a flexible generalised distribution, notably the MacGillivray adaptation of the g-and-κ distribution. This distribution, described through its inverse cdf or quantile function, generalises the standard normal through extra parameters which together describe skewness and kurtosis. The standard quantile-based methods for estimating the parameters of generalised distributions are often arbitrary and do not rely on computation of the likelihood. MCMC, however, provides a simulation-based alternative for obtaining the maximum likelihood estimates of parameters of these distributions or for deriving posterior estimates of the parameters through a Bayesian framework. In this paper we adopt the latter approach, The proposed methodology is illustrated through an application in which the parameter of interest is slightly skewed.

Accessing chain length dependent termination rate coefficients of methyl methacrylate (MMA) via the reversible addition fragmentation chain transfer (RAFT) process

The RAFT-CLD-T methodology is demonstrated to be not only applicable to 1-substituted monomers such as styrene and acrylates, but also to 1,1-disubstituted monomers such as MMA. The chain length of the terminating macromolecules is controlled by CPDB in MMA bulk free radical polymerization at 80 degrees C. The evolution of the chain length dependent termination rate coefficient, k(t)(i,i), was constructed in a step-wise fashion, since the MMA/CPDB system displays hybrid behavior (between conventional and living free radical polymerization) resulting in initial high molecular weight polymers formed at low RAFT agent concentrations. The obtained CLD of k(t) in MMA polymerizations is compatible with the composite model for chain length dependent termination. For the initial chain-length regime, up to a degree of polymerization of 100, k(t) decreases with alpha (in the expression k(t)(i,i) = k(t)(0) . i(-alpha)) being close to 0.65 at 80 degrees C. At chain lengths exceeding 100, the decrease is less pronounced (affording an alpha of 0.15 at 80 degrees C). However, the data are best represented by a continuously decreasing nonlinear functionality implying a chain length dependent alpha.

Modeling the molecular weight distribution of block copolymer formation in a reversible addition-fragmentation chain transfer mediated living radical polymerization

Block copolymers have become an integral part of the preparation of complex architectures through self-assembly. The use of reversible addition-fragmentation chain transfer (RAFT) allows blocks ranging from functional to nonfunctional polymers to be made with predictable molecular weight distributions. This article models block formation by varying many of the kinetic parameters. The simulations provide insight into the overall polydispersities (PDIs) that will be obtained when the chain-transfer constants in the main equilibrium steps are varied from 100 to 0.5. When the first dormant block [polymer-S-C(Z)=S] has a PDI of 1 and the second propagating radical has a low reactivity to the RAFT moiety, the overall PDI will be greater than 1 and dependent on the weight fraction of each block. When the first block has a PDI of 2 and the second propagating radical has a low reactivity to the RAFT moiety, the PDI will decrease to around 1.5 because of random coupling of two broad distributions. It is also shown how we can in principle use only one RAFT agent to obtain block copolymers with any desired molecular weight distribution. We can accomplish this by maintaining the monomer concentration at a constant level in the reactor over the course of the reaction. (c) 2005 Wiley Periodicals, Inc.

Controlled Radical Polymerization of Styrene and Methyl Acrylate in the Presence of Reversible Addition-Fragmentation Chain Transfer Agents, Phenylethyl Phenyl Dithioacetate and Phenyldithioacetic Acid

The use of phenyldithioacetic acid (PDA) in homopolymerizations of styrene or methyl acrylate produced only a small fraction of chains with dithioester end groups. The polymerizations using 1-phenylentyl phenyldithioacetate (PEPDTA) and PDA in the same reaction showed that PDA had little or no influence on the rate or molecular weight distribution even when a 1:1 ratio is used. The mechanistic pathway for the polymerizations in the presence of PDA seemed to be different for each monomer. Styrene favors addition of styrene to PDA via a Markovnikov type addition to form a reactive RAFT agent. The polymer was shown by double detection SEC to contain dithioester end groups over the whole distribution. This polymer was then used in a chain extension experiment and the M-n was close to theory. A unique feature of this work was that PDA could be used to form a RAFT agent in situ by heating a mixture of styrene and PDA for 24 h at 70 degrees C and then polymerizing in the presence of AIBN to give a linear increase in Mn and low values of PDI (< 1.14). In the case of the polymerization of MA with PDA, the mechanism was proposed to be via degradative chain transfer. (c) 2005 Wiley Periodicals, Inc.

Design strategies for controlling the molecular weight and rate using reversible addition-fragmentation chain transfer mediated living radical polymerization

Living radical polymerization has allowed complex polymer architectures to be synthesized in bulk, solution, and water. The most versatile of these techniques is reversible addition-fragmentation chain transfer (RAFT), which allows a wide range of functional and nonfunctional polymers to be made with predictable molecular weight distributions (MWDs), ranging from very narrow to quite broad. The great complexity of the RAFT mechanism and how the kinetic parameters affect the rate of polymerization and MWD are not obvious. Therefore, the aim of this article is to provide useful insights into the important kinetic parameters that control the rate of polymerization and the evolution of the MWD with conversion. We discuss how a change in the chain-transfer constant can affect the evolution of the MWD. It is shown how we can, in principle, use only one RAFT agent to obtain a poly-mer with any MWD. Retardation and inhibition are discussed in terms of (1) the leaving R group reactivity and (2) the intermediate radical termination model versus the slow fragmentation model. (c) 2005 Wiley Periodicals, Inc.

Pitch jump probability measures for the analysis of snoring sounds in apnea

Obstructive sleep apnea (OSA) is a highly prevalent disease in which upper airways are collapsed during sleep, leading to serious consequences. The gold standard of diagnosis, called polysomnography (PSG), requires a full-night hospital stay connected to over ten channels of measurements requiring physical contact with sensors. PSG is inconvenient, expensive and unsuited for community screening. Snoring is the earliest symptom of OSA, but its potential in clinical diagnosis is not fully recognized yet. Diagnostic systems intent on using snore-related sounds (SRS) face the tough problem of how to define a snore. In this paper, we present a working definition of a snore, and propose algorithms to segment SRS into classes of pure breathing, silence and voiced/unvoiced snores. We propose a novel feature termed the 'intra-snore-pitch-jump' (ISPJ) to diagnose OSA. Working on clinical data, we show that ISPJ delivers OSA detection sensitivities of 86-100% while holding specificity at 50-80%. These numbers indicate that snore sounds and the ISPJ have the potential to be good candidates for a take-home device for OSA screening. Snore sounds have the significant advantage in that they can be conveniently acquired with low-cost non-contact equipment. The segmentation results presented in this paper have been derived using data from eight patients as the training set and another eight patients as the testing set. ISPJ-based OSA detection results have been derived using training data from 16 subjects and testing data from 29 subjects.

Resolution and characterisation of multiple isoforms of bovine kappa-casein by 2-DE following a reversible cysteine-tagging enrichment strategy

Visualisation of multiple isoforms of kappa-casein on 2-D gels is restricted by the abundant alpha- and beta-caseins that not only limit gel loading but also migrate to similar regions as the more acidic kappa-casein isoforms. To overcome this problem, we took advantage of the absence of cysteine residues in alpha(S1)- and beta-casein by devising an affinity enrichment procedure based on reversible biotinylation of cysteine residues. Affinity capture of cysteine-containing proteins on avidin allowed the removal of the vast majority of alpha(S1)- and beta-casein, and on subsequent 2-D gel analysis 16 gel spots were identified as kappa-casein by PMF. Further analysis of the C-terminal tryptic peptide along with structural predictions based on mobility on the 2-D gel allowed us to assign identities to each spot in terms of genetic variant (A or B), phosphorylation status (1, 2 or 3) and glycosylation status (from 0 to 6). Eight isoforms of the A and B variants with the same PTMs were observed. When the casein fraction of milk from a single cow, homozygous for the B variant of kappa-casein, was used as the starting material, 17 isoforms from 13 gel spots were characterised. Analysis of isoforms of low abundance proved challenging due to the low amount of material that could be extracted from the gels as well as the lability of the PTMs during MS analysis. However, we were able to identify a previously unrecognised site, T-166, that could be phosphorylated or glycosylated. Despite many decades of analysis of milk proteins, the reasons for this high level of heterogeneity are still not clear.

Reversible active switching of the mechanical properties of a peptide film at a fluid-fluid interface

Designer peptides have recently been developed as building blocks for novel self-assembled materials with stimuli-responsive properties. To date, such materials have been based on self-assembly in bulk aqueous solution or at solid-fluid interfaces. We have designed a 21-residue peptide, AM1, as a stimuli-responsive surfactant that switches molecular architectures at a fluid-fluid interface in response to changes in bulk aqueous solution composition. In the presence of divalent zinc at neutral pH, the peptide forms a mechanically strong 'film state'. In the absence of metal ions or at acid pH, the peptide adsorbs to form a mobile 'detergent state'. The two interfacial states can be actively and reversibly switched. Switching between the two states by a change in pH or the addition of a chelating agent leads to rapid emulsion coalescence or foam collapse. This work introduces a new class of surfactants that offer an environmentally friendly approach to control the stability of interfaces in foams, emulsions and fluid-fluid interfaces more generally.

Synthesis of monocyclic and linear polystyrene using the reversible coupling/cleavage of thiol/disulfide groups

By carefully controlling the concentration of alpha,omega-thiol polystyrene in solution, we achieved formation of unique monocyclic polystyrene chains (i.e., polymer chains with only one disulfide linkage). The presence of cyclic polystyrene was confirmed by its lower than expected molecular weight due to a lower hydrodynamic volume and loss of thiol groups as detected by using Ellman's reagent. The alpha,omega-thiol polystyrene was synthesized by polymerizing styrene in the presence of a difunctional RAFT agent and subsequent conversion of the dithioester end groups to thiols via the addition of hexylamine. Oxidation gave either monocyclic polymer chains (i.e., with only one disulfide linkage) or linear multiblock polymers with many disulfide linkages depending on the concentration of polymer used with greater chance of cyclization in more dilute solutions. At high polymer concentrations, linear multiblock polymers were formed. To control the MWD of these linear multiblocks, monofunctional X-PSTY (X = PhCH2C(S)-S-) was added. It was found that the greatest ratio of X-PSTY to X-PSTY-X resulted in a low M-n and PDI. We have shown that we can control both the structure and MWD using this chemistry, but more importantly such disulfide linkages can be readily reduced back to the starting polystyrene with thiol end groups, which has potential use for a recyclable polymer material.

Submerged and unsubmerged natural hydraulic jumps in a bedrock step-pool mountain channel (Fig. 1 Classical hydraulic jump with partially developed inflow conditions)

Fig. 1. Classical hydraulic jump with partially developed inflow conditions. F1 = 13.6, V1 = 4.7 m/s, B = 0.25 m, h = 0.020 mm, d1 = 0.012 mm, Q = 14 L/s. Photo courtesy of Dr. Hubert Chanson. published in: Geomorphology Volume 82, Issues 1-2, 6 December 2006, Pages 146-159 The Hydrology and Geomorphology of Bedrock Rivers doi:10.1016/j.geomorph.2005.09.024 Submerged and unsubmerged natural hydraulic jumps in a bedrock step-pool mountain channel Brett L. Vallé and Gregory B. Pasternacka