A Kushner-Stratonovich Monte Carlo filter applied to nonlinear dynamical system identification

A Monte Carlo filter, based on the idea of averaging over characteristics and fashioned after a particle-based time-discretized approximation to the Kushner-Stratonovich (KS) nonlinear filtering equation, is proposed. A key aspect of the new filter is the gain-like additive update, designed to approximate the innovation integral in the KS equation and implemented through an annealing-type iterative procedure, which is aimed at rendering the innovation (observation prediction mismatch) for a given time-step to a zero-mean Brownian increment corresponding to the measurement noise. This may be contrasted with the weight-based multiplicative updates in most particle filters that are known to precipitate the numerical problem of weight collapse within a finite-ensemble setting. A study to estimate the a-priori error bounds in the proposed scheme is undertaken. The numerical evidence, presently gathered from the assessed performance of the proposed and a few other competing filters on a class of nonlinear dynamic system identification and target tracking problems, is suggestive of the remarkably improved convergence and accuracy of the new filter. (C) 2013 Elsevier B.V. All rights reserved.

Identification of heat shock factor binding protein in Plasmodium falciparum

Background: Heat shock factor binding protein (HSBP) was originally discovered in a yeast two-hybrid screen as an interacting partner of heat shock factor (HSF). It appears to be conserved in all eukaryotes studied so far, with yeast being the only exception. Cell biological analysis of HSBP in mammals suggests its role as a negative regulator of heat shock response as it appears to interact with HSF only during the recovery phase following exposure to heat stress. While the identification of HSF in the malaria parasite is still eluding biologists, this study for the first time, reports the presence of a homologue of HSBP in Plasmodium falciparum. Methods: PfHSBP was cloned and purified as his-tag fusion protein. CD (Circular dichroism) spectroscopy was performed to predict the secondary structure. Immunoblots and immunofluorescence approaches were used to study expression and localization of HSBP in P. falciparum. Cellular fractionation was performed to examine subcellular distribution of PfHSBP. Immunoprecipitation was carried out to identify HSBP interacting partner in P. falciparum. Results: PfHSBP is a conserved protein with a high helical content and has a propensity to form homo-oligomers. PfHSBP was cloned, expressed and purified. The in vivo protein expression profile shows maximal expression in trophozoites. The protein was found to exist in oligomeric form as trimer and hexamer. PfHSBP is predominantly localized in the parasite cytosol, however, upon heat shock, it translocates to the nucleus. This study also reports the interaction of PfHSBP with PfHSP70-1 in the cytoplasm of the parasite. Conclusions: This study emphasizes the structural and biochemical conservation of PfHSBP with its mammalian counterpart and highlights its potential role in regulation of heat shock response in the malaria parasite. Analysis of HSBP may be an important step towards identification of the transcription factor regulating the heat shock response in P. falciparum.

Titanium promoted reduction of imines with Grignards, silanes, and zinc: identification of a new mechanism with silanes

Aldimines react with reducing agents, such as Grignards, phenylsilane or zinc in the presence of titanium(IV) isopropoxide to form amines and reductively coupled imines (diamines). Using deuterium labeled reagents, the mechanism of reduction to form amines is described. Reducing agents, such as the Grignard and zinc result in the formation of low valent titanium (LVT), which in turn reduces the imine. On the other hand, phenylsilane reacts by a distinctly different mechanism and where a hydrogen atom from silicon is directly transferred to the titanium coordinated imine. (c) 2014 Elsevier Ltd. All rights reserved.

Identification of miRNAs from French bean (Phaseolus vulgaris) under low nitrate stress

Objective: In this study, we report the role of miRNAs involved under nitrogen starvation from widely grown vegetable crop, French bean. In recent years, a great deal of attention has been paid to the elucidation of miRNAs involved in low nitrate stress. Methods: To identify miRNAs expressed under stress, cDNA libraries were analyzed. Results: We reported the nine potential miRNAs with 67 targets involved in nutrient transporters and other stress specific genes. Among the miRNA sequences obtained 6 sequences belong to miR172 family, one with miR169. RT-PCR analysis of expression of miR172 family was induced upon low nitrate stress while miR169 family was repressed. In addition, Pvu-SN7b and Pvu-miR16 may be new members of miRNA172 and miR169 families, respectively. Conclusion: The targets of Pvu-SN7b were major protein kinases, one among which is the Protein Kinase CK2. CK2 Kinase is found to involve in transcription-directed signaling, gene control and cell-cycle regulation. Other targets of Pvu-SN7b were involved in DNA-dependent transcription regulation, photo-periodism, calcium-mediated signaling. Pvu-miR16 targets Thymidine kinase, the key enzyme of deoxy-nucleotide synthesis. The cleavage of these targets affects cell proliferation there by affecting nodule formation. Pvu-miR8 inhibits translation of its target protein Pre-protein translocase, a membrane-bound protein transporter involved in trans-membrane protein transportation. Together these results denote the response and role of miRNAs to nitrate-limiting conditions in French bean.

Automatic Tonic Identification in Indian Art Music: Approaches and Evaluation

The tonic is a fundamental concept in Indian art music. It is the base pitch, which an artist chooses in order to construct the melodies during a rg(a) rendition, and all accompanying instruments are tuned using the tonic pitch. Consequently, tonic identification is a fundamental task for most computational analyses of Indian art music, such as intonation analysis, melodic motif analysis and rg recognition. In this paper we review existing approaches for tonic identification in Indian art music and evaluate them on six diverse datasets for a thorough comparison and analysis. We study the performance of each method in different contexts such as the presence/absence of additional metadata, the quality of audio data, the duration of audio data, music tradition (Hindustani/Carnatic) and the gender of the singer (male/female). We show that the approaches that combine multi-pitch analysis with machine learning provide the best performance in most cases (90% identification accuracy on average), and are robust across the aforementioned contexts compared to the approaches based on expert knowledge. In addition, we also show that the performance of the latter can be improved when additional metadata is available to further constrain the problem. Finally, we present a detailed error analysis of each method, providing further insights into the advantages and limitations of the methods.

Temporal variability of forest communities: empirical estimates of population change in 4000 tree species

Long-term surveys of entire communities of species are needed to measure fluctuations in natural populations and elucidate the mechanisms driving population dynamics and community assembly. We analysed changes in abundance of over 4000 tree species in 12 forests across the world over periods of 6-28years. Abundance fluctuations in all forests are large and consistent with population dynamics models in which temporal environmental variance plays a central role. At some sites we identify clear environmental drivers, such as fire and drought, that could underlie these patterns, but at other sites there is a need for further research to identify drivers. In addition, cross-site comparisons showed that abundance fluctuations were smaller at species-rich sites, consistent with the idea that stable environmental conditions promote higher diversity. Much community ecology theory emphasises demographic variance and niche stabilisation; we encourage the development of theory in which temporal environmental variance plays a central role.

Role of Vertex Index in Substructure Identification and Activity Prediction: A Study on Antitubercular Activity of a Series of Acid Alkyl Ester Derivatives

Tuberculosis (TB) is a life threatening disease caused due to infection from Mycobacterium tuberculosis (Mtb). That most of the TB strains have become resistant to various existing drugs, development of effective novel drug candidates to combat this disease is a need of the day. In spite of intensive research world-wide, the success rate of discovering a new anti-TB drug is very poor. Therefore, novel drug discovery methods have to be tried. We have used a rule based computational method that utilizes a vertex index, named `distance exponent index (D-x)' (taken x = -4 here) for predicting anti-TB activity of a series of acid alkyl ester derivatives. The method is meant to identify activity related substructures from a series a compounds and predict activity of a compound on that basis. The high degree of successful prediction in the present study suggests that the said method may be useful in discovering effective anti-TB compound. It is also apparent that substructural approaches may be leveraged for wide purposes in computer-aided drug design.

A shortest path tree based algorithm for relay placement in a wireless sensor network and its performance analysis

In this paper, we study a problem of designing a multi-hop wireless network for interconnecting sensors (hereafter called source nodes) to a Base Station (BS), by deploying a minimum number of relay nodes at a subset of given potential locations, while meeting a quality of service (QoS) objective specified as a hop count bound for paths from the sources to the BS. The hop count bound suffices to ensure a certain probability of the data being delivered to the BS within a given maximum delay under a light traffic model. We observe that the problem is NP-Hard. For this problem, we propose a polynomial time approximation algorithm based on iteratively constructing shortest path trees and heuristically pruning away the relay nodes used until the hop count bound is violated. Results show that the algorithm performs efficiently in various randomly generated network scenarios; in over 90% of the tested scenarios, it gave solutions that were either optimal or were worse than optimal by just one relay. We then use random graph techniques to obtain, under a certain stochastic setting, an upper bound on the average case approximation ratio of a class of algorithms (including the proposed algorithm) for this problem as a function of the number of source nodes, and the hop count bound. To the best of our knowledge, the average case analysis is the first of its kind in the relay placement literature. Since the design is based on a light traffic model, we also provide simulation results (using models for the IEEE 802.15.4 physical layer and medium access control) to assess the traffic levels up to which the QoS objectives continue to be met. (C) 2014 Elsevier B.V. All rights reserved.

Identification of a novel AMPK-PEA15 axis in the anoikis-resistant growth of mammary cells

Introduction: Matrix detachment triggers anoikis, a form of apoptosis, in most normal epithelial cells, while acquisition of anoikis resistance is a prime requisite for solid tumor growth. Of note, recent studies have revealed that a small population of normal human mammary epithelial cells (HMECs) survive in suspension and generate multicellular spheroids termed `mammospheres'. Therefore, understanding how normal HMECs overcome anoikis may provide insights into breast cancer initiation and progression. Methods: Primary breast tissue-derived normal HMECs were grown as adherent monolayers or mammospheres. The status of AMP-activated protein kinase (AMPK) and PEA15 signaling was investigated by immunoblotting. Pharmacological agents and an RNA interference (RNAi) approach were employed to gauge their roles in mammosphere formation. Immunoprecipitation and in vitro kinase assays were undertaken to evaluate interactions between AMPK and PEA15. In vitro sphere formation and tumor xenograft assays were performed to understand their roles in tumorigenicity. Results: In this study, we show that mammosphere formation by normal HMECs is accompanied with an increase in AMPK activity. Inhibition or knockdown of AMPK impaired mammosphere formation. Concomitant with AMPK activation, we detected increased Ser(116) phosphorylation of PEA15, which promotes its anti-apoptotic functions. Inhibition or knockdown of AMPK impaired PEA15 Ser(116) phosphorylation and increased apoptosis. Knockdown of PEA15, or overexpression of the nonphosphorylatable S116A mutant of PEA15, also abrogated mammosphere formation. We further demonstrate that AMPK directly interacts with and phosphorylates PEA15 at Ser(116) residue, thus identifying PEA15 as a novel AMPK substrate. Together, these data revealed that AMPK activation facilitates mammosphere formation by inhibition of apoptosis, at least in part, through Ser(116) phosphorylation of PEA15. Since anoikis resistance plays a critical role in solid tumor growth, we investigated the relevance of these findings in the context of breast cancer. Significantly, we show that the AMPK-PEA15 axis plays an important role in the anchorage-independent growth of breast cancer cells both in vitro and in vivo. Conclusions: Our study identifies a novel AMPK-PEA15 signaling axis in the anchorage-independent growth of both normal and cancerous mammary epithelial cells, suggesting that breast cancer cells may employ mechanisms of anoikis resistance already inherent within a subset of normal HMECs. Thus, targeting the AMPK-PEA15 axis might prevent breast cancer dissemination and metastasis.

Multiscale Symmetry Detection in Scalar Fields by Clustering Contours

The complexity in visualizing volumetric data often limits the scope of direct exploration of scalar fields. Isocontour extraction is a popular method for exploring scalar fields because of its simplicity in presenting features in the data. In this paper, we present a novel representation of contours with the aim of studying the similarity relationship between the contours. The representation maps contours to points in a high-dimensional transformation-invariant descriptor space. We leverage the power of this representation to design a clustering based algorithm for detecting symmetric regions in a scalar field. Symmetry detection is a challenging problem because it demands both segmentation of the data and identification of transformation invariant segments. While the former task can be addressed using topological analysis of scalar fields, the latter requires geometry based solutions. Our approach combines the two by utilizing the contour tree for segmenting the data and the descriptor space for determining transformation invariance. We discuss two applications, query driven exploration and asymmetry visualization, that demonstrate the effectiveness of the approach.

And yet it shrinks: A novel method for correcting bias in forest tree growth estimates caused by water-induced fluctuations

Accuracy in tree woody growth estimates is important to global carbon budget estimation and climate-change science. Tree growth in permanent sampling plots (PSPs) is commonly estimated by measuring stem diameter changes, but this method is susceptible to bias resulting from water-induced reversible stem shrinkage. In the absence of bias correction, temporal variability in growth is likely to be overestimated and incorrectly attributed to fluctuations in resource availability, especially in forests with high seasonal and inter-annual variability in water. We propose and test a novel approach for estimating and correcting this bias at the community level. In a 50-ha PSP from a seasonally dry tropical forest in southern India, where tape measurements have been taken every four years from 1988 to 2012, for nine trees we estimated bias due to reversible stem shrinkage as the difference between woody growth measured using tree rings and that estimated from tape. We tested if the bias estimated from these trees could be used as a proxy to correct bias in tape-based growth estimates at the PSP scale. We observed significant shrinkage-related bias in the growth estimates of the nine trees in some censuses. This bias was strongly linearly related to tape-based growth estimates at the level of the PSP, and could be used as a proxy. After bias was corrected, the temporal variance in growth rates of the PSP decreased, while the effect of exceptionally dry or wet periods was retained, indicating that at least a part of the temporal variability arose from reversible shrinkage-related bias. We also suggest that the efficacy of the bias correction could be improved by measuring the proxy on trees that belong to different size classes and census timing, but not necessarily to different species. Our approach allows for reanalysis - and possible reinterpretation of temporal trends in tree growth, above ground biomass change, or carbon fluxes in forests, and their relationships with resource availability in the context of climate change. (C) 2014 Elsevier B.V. All rights reserved.

Iterated stochastic filters with additive updates for dynamic system identification: Annealing-type iterations and the filter bank

A nonlinear stochastic filtering scheme based on a Gaussian sum representation of the filtering density and an annealing-type iterative update, which is additive and uses an artificial diffusion parameter, is proposed. The additive nature of the update relieves the problem of weight collapse often encountered with filters employing weighted particle based empirical approximation to the filtering density. The proposed Monte Carlo filter bank conforms in structure to the parent nonlinear filtering (Kushner-Stratonovich) equation and possesses excellent mixing properties enabling adequate exploration of the phase space of the state vector. The performance of the filter bank, presently assessed against a few carefully chosen numerical examples, provide ample evidence of its remarkable performance in terms of filter convergence and estimation accuracy vis-a-vis most other competing filters especially in higher dimensional dynamic system identification problems including cases that may demand estimating relatively minor variations in the parameter values from their reference states. (C) 2014 Elsevier Ltd. All rights reserved.

On Whitespace Identification Using Randomly Deployed Sensors

This work considers the identification of the available whitespace, i.e., the regions that do not contain any existing transmitter within a given geographical area. To this end, n sensors are deployed at random locations within the area. These sensors detect for the presence of a transmitter within their radio range r(s) using a binary sensing model, and their individual decisions are combined to estimate the available whitespace. The limiting behavior of the recovered whitespace as a function of n and r(s) is analyzed. It is shown that both the fraction of the available whitespace that the nodes fail to recover as well as their radio range optimally scale as log(n)/n as n gets large. The problem of minimizing the sum absolute error in transmitter localization is also analyzed, and the corresponding optimal scaling of the radio range and the necessary minimum transmitter separation is determined.

Bayesian parameter identification in dynamic state space models using modified measurement equations

When Markov chain Monte Carlo (MCMC) samplers are used in problems of system parameter identification, one would face computational difficulties in dealing with large amount of measurement data and (or) low levels of measurement noise. Such exigencies are likely to occur in problems of parameter identification in dynamical systems when amount of vibratory measurement data and number of parameters to be identified could be large. In such cases, the posterior probability density function of the system parameters tends to have regions of narrow supports and a finite length MCMC chain is unlikely to cover pertinent regions. The present study proposes strategies based on modification of measurement equations and subsequent corrections, to alleviate this difficulty. This involves artificial enhancement of measurement noise, assimilation of transformed packets of measurements, and a global iteration strategy to improve the choice of prior models. Illustrative examples cover laboratory studies on a time variant dynamical system and a bending-torsion coupled, geometrically non-linear building frame under earthquake support motions. (C) 2015 Elsevier Ltd. All rights reserved.

Identification of local variations within secondary structures of proteins

Secondary-structure elements (SSEs) play an important role in the folding of proteins. Identification of SSEs in proteins is a common problem in structural biology. A new method, ASSP (Assignment of Secondary Structure in Proteins), using only the path traversed by the C atoms has been developed. The algorithm is based on the premise that the protein structure can be divided into continuous or uniform stretches, which can be defined in terms of helical parameters, and depending on their values the stretches can be classified into different SSEs, namely -helices, 3(10)-helices, -helices, extended -strands and polyproline II (PPII) and other left-handed helices. The methodology was validated using an unbiased clustering of these parameters for a protein data set consisting of 1008 protein chains, which suggested that there are seven well defined clusters associated with different SSEs. Apart from -helices and extended -strands, 3(10)-helices and -helices were also found to occur in substantial numbers. ASSP was able to discriminate non--helical segments from flanking -helices, which were often identified as part of -helices by other algorithms. ASSP can also lead to the identification of novel SSEs. It is believed that ASSP could provide a better understanding of the finer nuances of protein secondary structure and could make an important contribution to the better understanding of comparatively less frequently occurring structural motifs. At the same time, it can contribute to the identification of novel SSEs. A standalone version of the program for the Linux as well as the Windows operating systems is freely downloadable and a web-server version is also available at .