Estimation of Hidden Markov Models parameters using differential evolution

Hidden Markov Models (HMMs) have been successfully applied to different modelling and classification problems from different areas over the recent years. An important step in using HMMs is the initialisation of the parameters of the model as the subsequent learning of HMM’s parameters will be dependent on these values. This initialisation should take into account the knowledge about the addressed problem and also optimisation techniques to estimate the best initial parameters given a cost function, and consequently, to estimate the best log-likelihood. This paper proposes the initialisation of Hidden Markov Models parameters using the optimisation algorithm Differential Evolution with the aim to obtain the best log-likelihood.

Exploration vs exploitation in differential evolution

Differential Evolution (DE) is a tool for efficient optimisation, and it belongs to the class of evolutionary algorithms, which include Evolution Strategies and Genetic Algorithms. DE algorithms work well when the population covers the entire search space, and they have shown to be effective on a large range of classical optimisation problems. However, an undesirable behaviour was detected when all the members of the population are in a basin of attraction of a local optimum (local minimum or local maximum), because in this situation the population cannot escape from it. This paper proposes a modification of the standard mechanisms in DE algorithm in order to change the exploration vs. exploitation balance to improve its behaviour.

Application of auto regressive models of wavelet sub-bands for classifying Terahertz pulse measurements

This paper presents an approach for automatic classification of pulsed Terahertz (THz), or T-ray, signals highlighting their potential in biomedical, pharmaceutical and security applications. T-ray classification systems supply a wealth of information about test samples and make possible the discrimination of heterogeneous layers within an object. In this paper, a novel technique involving the use of Auto Regressive (AR) and Auto Regressive Moving Average (ARMA) models on the wavelet transforms of measured T-ray pulse data is presented. Two example applications are examined - the classi. cation of normal human bone (NHB) osteoblasts against human osteosarcoma (HOS) cells and the identification of six different powder samples. A variety of model types and orders are used to generate descriptive features for subsequent classification. Wavelet-based de-noising with soft threshold shrinkage is applied to the measured T-ray signals prior to modeling. For classi. cation, a simple Mahalanobis distance classi. er is used. After feature extraction, classi. cation accuracy for cancerous and normal cell types is 93%, whereas for powders, it is 98%.

On the spectra of certain integro-differential-delay problems with applications in neurodynamics

We investigate the spectrum of certain integro-differential-delay equations (IDDEs) which arise naturally within spatially distributed, nonlocal, pattern formation problems. Our approach is based on the reformulation of the relevant dispersion relations with the use of the Lambert function. As a particular application of this approach, we consider the case of the Amari delay neural field equation which describes the local activity of a population of neurons taking into consideration the finite propagation speed of the electric signal. We show that if the kernel appearing in this equation is symmetric around some point a= 0 or consists of a sum of such terms, then the relevant dispersion relation yields spectra with an infinite number of branches, as opposed to finite sets of eigenvalues considered in previous works. Also, in earlier works the focus has been on the most rightward part of the spectrum and the possibility of an instability driven pattern formation. Here, we numerically survey the structure of the entire spectra and argue that a detailed knowledge of this structure is important within neurodynamical applications. Indeed, the Amari IDDE acts as a filter with the ability to recognise and respond whenever it is excited in such a way so as to resonate with one of its rightward modes, thereby amplifying such inputs and dampening others. Finally, we discuss how these results can be generalised to the case of systems of IDDEs.

Contribution of the double layer to transient faradaic processes: implications for hydrodynamic modulated voltammetry of nanostructures

The electrochemistry of Pt nanostructured electrodes is investigated using hydrodynamic modulated voltammetry (HMV). Here a liquid crystal templating process is used to produce platinum-modified electrodes with a range of surface areas (roughness factor 42.4-280.8). The electroreduction of molecular oxygen at these nanostructured platinum surfaces is used to demonstrate the ability of HMV to discriminate between faradaic and nonfaradaic electrode reactions. The HMV approach shows that the reduction of molecular oxygen experiences considerable signal loss within the high pseudocapacitive region of the voltammetry. Evidence for the contribution of the double layer to transient mass transfer events is presented. In addition, a model circuit and appropriate theoretical analysis are used to illustrate the transient responses of a time variant faradaic component. This in conjunction with the experimental evidence shows that, far from being a passive component in this system, the double layer can contribute to HMV faradaic reactions under certain conditions.

Differential changes in transforming growth factor β isoform expression during postnatal cardiac growth

Transforming growth factor-β (TGF-β) is synthesised as an inactive precursor protein; this is cleaved to produce the mature peptide and a latency associated protein (LAP), which remains associated with the mature peptide until activation by LAP degradation. Isoform specific antibodies raised against the LAPs for TGF-β2and -β3were used to determine the myocardial levels of LAP (activatable TGF-β) and full length precursor (inactive TGF-β) forms during post-natal development in the rat. TGF-β2was present predominantly as the precursor in 2 day old myocardium. There was an age-dependent shift from precursor protein to LAP between 2 and 28 days. A corresponding increase in the level of mature (activatable) TGF-β2was found. TGF-β3was detected in significant quantities only as LAP. However, a four-fold increase in the expression of TGF-β3LAP was observed between 2 and 28 days. The substantial increases in activatable forms of TGF-β2and -β3that occur in myocardium during the first 28 days of life in the rat support a role for these proteins in post-natal cardiac development.

Differential protein expression and subcellular distribution of TGFβ1, β2 and β3 in cardiomyocytes during pressure overload-induced hypertrophy

The transforming growth factorβ(TGFβ) superfamily plays an important role in the myocardial response to hypertrophy. We have investigated the protein expression of TGFβ1,β2andβ3in left ventricular tissue, and determined their subcellular distribution in myocytes by immunoblotting and immunocytochemistry during the development of left ventricular hypertrophy (LVH), using isoform specific antibodies to TGFβ1,β2andβ3. LVH was produced in rats by aortic constriction (AC) and LV tissue was obtained at days (d)0, 1, 3, 7, 14, 21 and 42 following operation. Compared with age matched sham-operated controls (SH), TGFβ1levels in LV tissue of AC rats increased significantly from d1–d14 (P<0.03) concomitant with the adaptive growth of LV tissue. In contrast, TGFβ3levels decreased in LV tissue of AC rats from d3 post-operation (significant from d14–d42,P<0.03). No significant difference in TGFβ2levels were observed from SH and AC rats after operation. Antibodies to TGFβ1stained intercalated disks, sarcolemmal membranes and cytoplasm, but not nuclei, of cardiomyocytes on LV sections from untreated and SH rats. However, a trans-localisation of TGFβ1to the nuclei of cardiomyocytes was observed in AC hearts. Antibodies to TGFβ3stained T tubules, cytoplasm and the nuclei of cardiomyocytes from untreated and SH rats. However, by d7 post-AC operation, TGFβ3expression was lost rapidly from nuclei of cardiomyocytes followed by a reduction in total TGFβ3immunofluorescence in myocytes. Antibodies to TGFβ2stained sarcolemmal membranes of cardiomyocytes from both SH and AC rats without significant difference between groups. Thus, the differential pattern of protein expression and subcellular distribution of TGFβ1,β2andβ3in myocytes during the development of LVH suggests that these molecules play different roles in the response of cardiomyocytes to LVH.

Differential induction of 12-O-tetradecanoylphorbol 13-acetate sequence gene expression in murine melanocytes and melanoma cells

We previously showed that growth of the nontumorigenic, immortal murine melanocyte line Mel-ab correlates with the depletion of protein kinase C (PKC), whereas quiescence is associated with elevated levels of this enzyme (Brooks G, et al., Cancer Res 51: 3281–3288, 1991). Here we report responses that occur in these cells downstream of PKC activation or downregulation. We examined induction of 12-O-tetradecanoylphorbol-13-acetate (TPA)-inducible sequence (TIS) gene expression in Mel-ab melanocytes and in their transformed counterparts, B16 melanoma cells. Exposure of quiescent Mel-ab cells to the PKC-activating phorbol esters TPA or sapintoxin A at 81 nM for 2 h increased levels of mRNA for six of seven TIS genes examined (twofold to 80-fold increase in steady-state RNA levels for TIS 1, 7, 8, 11, 21, and 28 (c-fos); TIS 10 expression was not affected). No induction of 115 gene expression was observed either in growing Mel-ab cells maintained in 324 nM phorbol 12,13-dibutyrate or in B16 cells previously unexposed to phorbol esters, in which normal PKC levels were endogenously depressed. The cAMP-elevating agents choleratoxin (10 nM) and dibutyryl cyclic AMP (2.5 mM) increased levels of TIS mRNA (with the exception of TIS 10) in both proliferating Mel-ab and B16 cells, suggesting that downregulation of the PKC pathway is specific and not a consequence of a general inhibition of all signalling pathways.

Multiplexed quantitative proteomics using differential metabolic 15N-labeling

There are several advantages of using metabolic labeling in quantitative proteomics. The early pooling of samples compared to post-labeling methods eliminates errors from different sample processing, protein extraction and enzymatic digestion. Metabolic labeling is also highly efficient and relatively inexpensive compared to commercial labeling reagents. However, methods for multiplexed quantitation in the MS-domain (or ‘non-isobaric’ methods), suffer from signal dilution at higher degrees of multiplexing, as the MS/MS signal for peptide identification is lower given the same amount of peptide loaded onto the column or injected into the mass spectrometer. This may partly be overcome by mixing the samples at non-uniform ratios, for instance by increasing the fraction of unlabeled proteins. We have developed an algorithm for arbitrary degrees of nonisobaric multiplexing for relative protein abundance measurements. We have used metabolic labeling with different levels of 15N, but the algorithm is in principle applicable to any isotope or combination of isotopes. Ion trap mass spectrometers are fast and suitable for LC-MS/MS and peptide identification. However, they cannot resolve overlapping isotopic envelopes from different peptides, which makes them less suitable for MS-based quantitation. Fourier-transform ion cyclotron resonance (FTICR) mass spectrometry is less suitable for LC-MS/MS, but provides the resolving power required to resolve overlapping isotopic envelopes. We therefore combined ion trap LC-MS/MS for peptide identification with FTICR LC-MS for quantitation using chromatographic alignment. We applied the method in a heat shock study in a plant model system (A. thaliana) and compared the results with gene expression data from similar experiments in literature.

Differential recognition of Staphylococcus aureus quorum-sensing signals depends on both extracellular loops 1 and 2 of the transmembrane sensor AgrC.

Virulence in Staphylococcus aureus is regulated via agr-dependent quorum sensing in which an autoinducing peptide (AIP) activates AgrC, a histidine protein kinase. AIPs are usually thiolactones containing seven to nine amino acid residues in which the thiol of the central cysteine is linked to the alpha-carboxyl of the C-terminal amino acid residue. The staphylococcal agr locus has diverged such that the AIPs of the four different S. aureus agr groups self-activate but cross-inhibit. Consequently, although the agr system is conserved among the staphylococci, it has undergone significant evolutionary divergence whereby to retain functionality, any changes in the AIP-encoding gene (agrD) that modifies AIP structure must be accompanied by corresponding changes in the AgrC receptor. Since AIP-1 and AIP-4 only differ by a single amino acid, we compared the transmembrane topology of AgrC1 and AgrC4 to identify amino acid residues involved in AIP recognition. As only two of the three predicted extracellular loops exhibited amino acid differences, site-specific mutagenesis was used to exchange the key AgrC1 and AgrC4 amino acid residues in each loop either singly or in combination. A novel lux-based agrP3 reporter gene fusion was constructed to evaluate the response of the mutated AgrC receptors. The data obtained revealed that while differential recognition of AIP-1 and AIP-4 depends primarily on three amino acid residues in loop 2, loop 1 is essential for receptor activation by the cognate AIP. Furthermore, a single mutation in the AgrC1 loop 2 resulted in conversion of (Ala5)AIP-1 from a potent antagonist to an activator, essentially resulting in the forced evolution of a new AIP group. Taken together, our data indicate that loop 2 constitutes the predicted hydrophobic pocket that binds the AIP thiolactone ring while the exocyclic amino acid tail interacts with loop 1 to facilitate receptor activation.