A semi-analytical particle filter for identification of nonlinear oscillators

Particle filters find important applications in the problems of state and parameter estimations of dynamical systems of engineering interest. Since a typical filtering algorithm involves Monte Carlo simulations of the process equations, sample variance of the estimator is inversely proportional to the number of particles. The sample variance may be reduced if one uses a Rao-Blackwell marginalization of states and performs analytical computations as much as possible. In this work, we propose a semi-analytical particle filter, requiring no Rao-Blackwell marginalization, for state and parameter estimations of nonlinear dynamical systems with additively Gaussian process/observation noises. Through local linearizations of the nonlinear drift fields in the process/observation equations via explicit Ito-Taylor expansions, the given nonlinear system is transformed into an ensemble of locally linearized systems. Using the most recent observation, conditionally Gaussian posterior density functions of the linearized systems are analytically obtained through the Kalman filter. This information is further exploited within the particle filter algorithm for obtaining samples from the optimal posterior density of the states. The potential of the method in state/parameter estimations is demonstrated through numerical illustrations for a few nonlinear oscillators. The proposed filter is found to yield estimates with reduced sample variance and improved accuracy vis-a-vis results from a form of sequential importance sampling filter.

Identification of alpha- and beta-Hydroxy Acid Containing Cyclodepsipeptides in Natural Peptide Mixtures Using Negative Ion Mass Spectrometry

Natural peptide libraries often contain cyclodepsipeptides containing alpha or beta hydroxy residues. Extracts of fungal hyphae of Isaria yield a microheterogenous cyclodepsipeptide mixture in which two classes of molecules can be identified by mass spectral fragmentation of negative ions. In the case of isaridins, which contain an alpha-hydroxy residue and a beta-amino acid residue, a characteristic product ion corresponding to a neutral loss of 72 Da is obtained. hi addition, neutral loss of water followed by a 72 Da loss is also observed. Two distinct modes of fragmentation rationalize the observed product ion distribution. The neutral loss of 72 Da has also been obtained for a roseotoxin component, which is also an alpha-hydroxy residue containing cyclodepsipeptide. In the case of isariins, which contain a beta-hydroxy acid residue, ring opening and subsequent loss of the terminal residue as an unsaturated ketene fragment, rationalizes the observed product ion formation. Fragmentation of negative ions provide characteristic neutral losses, which are diagnostic of the presence of alpha-hydroxy or beta-hydroxy residues.

On the Structure of Contractible Edges in k-connected Partial k-trees

Contraction of an edge e merges its end points into a new single vertex, and each neighbor of one of the end points of e is a neighbor of the new vertex. An edge in a k-connected graph is contractible if its contraction does not result in a graph with lesser connectivity; otherwise the edge is called non-contractible. In this paper, we present results on the structure of contractible edges in k-trees and k-connected partial k-trees. Firstly, we show that an edge e in a k-tree is contractible if and only if e belongs to exactly one (k + 1) clique. We use this characterization to show that the graph formed by contractible edges is a 2-connected graph. We also show that there are at least |V(G)| + k - 2 contractible edges in a k-tree. Secondly, we show that if an edge e in a partial k-tree is contractible then e is contractible in any k-tree which contains the partial k-tree as an edge subgraph. We also construct a class of contraction critical 2k-connected partial 2k-trees.

Identification and Characterization of a Novel Deoxyhypusine Synthase in Leishmania donovani

Deoxyhypusine synthase, an NAD(+)-dependent enzyme, catalyzes the first step in the post-translational synthesis of an unusual amino acid, hypusine (N-epsilon-(4-amino-2-hydroxybutyl)lysine), in the eukaryotic initiation factor 5A precursor protein. Two putative deoxyhypusine synthase (DHS) sequences have been identified in the Leishmania donovani genome, which are present on chromosomes 20: DHSL20 (DHS-like gene from chromosome 20) and DHS34 (DHS from chromosome 34). Although both sequences exhibit an overall conservation of key residues, DHSL20 protein lacks a critical lysine residue, and the recombinant protein showed no DHS activity in vitro. However, DHS34 contains the critical lysine residue, and the recombinant DHS34 effectively catalyzed deoxyhypusine synthesis. Furthermore, in vivo labeling confirmed that hypusination of eukaryotic initiation factor 5A occurs in intact Leishmania parasites. Interestingly, the DHS34 is much longer, with 601 amino acids, compared with the human DHS enzyme (369 amino acids) and contains several unique insertions. To study the physiological role of DHS34 in Leishmania, gene deletion mutations were attempted via targeted gene replacement. However, chromosomal null mutants of DHS34 could only be obtained in the presence of a DHS34-containing episome. The present data provide evidence that DHS34 is essential for L. donovani and that structural differences in the human and leishmanial DHS enzyme may be exploited for designing selective inhibitors against the parasite.

Azospirillum lipoferum tRNAs: Fractionation and identification of tRNA species and the nucleotide sequence of tRNA(Asn) (QUU)

Transfer RNAs of Azospirillum lipoferum were separated by two- dimensional gel electrophoresis and identified by aminoacylation. Thirty-six tRNA spots were resolved by this technique and twenty-six tRNA species have been identified. There are five tRNAs for Leu, four for Val, three for Pro, two each for Arg, Ile, Lys and Tyr, and one each for Ala, Asp, His, Phe, Ser and Thr. The tRNA(Asn) (QUU) was purified and its nucleotide sequence was determined. The A. lipoferum tRNA(Asn) (QUU) is 92% similar to B. subtilis tRNA(Asn) gene and two hypermodified nucleosides, queuosine (Q) and N-(9-beta-D Ribofuranosylpurine-6-YL) carbamoyl)-threonine (t(6)A) are present in this tRNA.

Identification and characterization of thionucleosides in the total tRNA of cucumber cotyledons

Total tRNA isolated from cucumber cotyledons grown in the presence of radioactive sulfur was analyzed for the occurrence of thionucleosides. The analysis revealed the presence of at least five thionucleosides which were identified as 5-methylaminomethyl-2-thiouridine (mnm5s2U), 2-methylthio-N6-isopentenyladenosine (ms2i6A), 2-methylthio-N6-hydroxyisopentenyladenosine (ms2io6A), 5-methyl-2-thiouridine (m5s2U) and N-[(9-beta-ribofuranosyl-2- methylthiopurine-2-yl)-carbamoyl]-threonine (ms2t6A). A comparison of relative amounts of these thionucleosides in the total tRNAs of dark-, and light-grown cotyledons shows that the relative amounts of ms2i6A, ms2io6A and ms2t6A remain unchanged whereas mnm5s2U increases with a concomitant decrease in the relative amounts of m5s2U after light treatment of dark-grown cotyledons.

A Kalman filter based strategy for linear structural system identification based on multiple static and dynamic test data

The problem of identification of stiffness, mass and damping properties of linear structural systems, based on multiple sets of measurement data originating from static and dynamic tests is considered. A strategy, within the framework of Kalman filter based dynamic state estimation, is proposed to tackle this problem. The static tests consists of measurement of response of the structure to slowly moving loads, and to static loads whose magnitude are varied incrementally; the dynamic tests involve measurement of a few elements of the frequency response function (FRF) matrix. These measurements are taken to be contaminated by additive Gaussian noise. An artificial independent variable τ, that simultaneously parameterizes the point of application of the moving load, the magnitude of the incrementally varied static load and the driving frequency in the FRFs, is introduced. The state vector is taken to consist of system parameters to be identified. The fact that these parameters are independent of the variable τ is taken to constitute the set of ‘process’ equations. The measurement equations are derived based on the mechanics of the problem and, quantities, such as displacements and/or strains, are taken to be measured. A recursive algorithm that employs a linearization strategy based on Neumann’s expansion of structural static and dynamic stiffness matrices, and, which provides posterior estimates of the mean and covariance of the unknown system parameters, is developed. The satisfactory performance of the proposed approach is illustrated by considering the problem of the identification of the dynamic properties of an inhomogeneous beam and the axial rigidities of members of a truss structure.

Sensitivity of Eigen Value to Damage and Its Identification

The reduction in natural frequencies,however small, of a civil engineering structure, is the first and the easiest method of estimating its impending damage. As a first level screening for health-monitoring, information on the frequency reduction of a few fundamentalmodes can be used to estimate the positions and the magnitude of damage in a smeared fashion. The paper presents the Eigen value sensitivity equations, derived from first-order perturbation technique, for typical infra-structural systems like a simply supported bridge girder, modelled as a beam, an endbearing pile, modelled as an axial rod and a simply supported plate as a continuum dynamic system. A discrete structure, like a building frame is solved for damage using Eigen-sensitivity derived by a computationalmodel. Lastly, neural network based damage identification is also demonstrated for a simply supported bridge beam, where the known-pairs of damage-frequency vector is used to train a neural network. The performance of these methods under the influence of measurement error is outlined. It is hoped that the developed method could be integrated in a typical infra-structural management program, such that magnitudes of damage and their positions can be obtained using acquired natural frequencies, synthesized from the excited/ambient vibration signatures.

Identification and thermodynamic characterization of molten globule states of periplasmic binding proteins

Molten globule-like intermediates have been shown to occur during protein folding and are thought to be involved in protein translocation and membrane insertion. However, the determinants of molten globule stability and the extent of specific packing in molten globules is currently unclear. Using far- and near-UV CD and intrinsic and ANS fluorescence, we show that four periplasmic binding proteins (LBP, LIVBP, MBP, and RBP) form molten globules at acidic pH values ranging from 3.0 to 3.4. Only two of these (LBP and LIVBP) have similar sequences, but all four proteins adopt similar three-dimensional structures. We found that each of the four molten globules binds to its corresponding ligand without conversion to the native state. Ligand binding affinity measured by isothermal titration calorimetry for the molten globule state of LIVBP was found to be comparable to that of the corresponding native state, whereas for LBP, MBP, and RBP, the molten globules bound ligand with approximately 5-30-fold lower affinity than the corresponding native states. All four molten globule states exhibited cooperative thermal unfolding assayed by DSC. Estimated values of Delta C-p of unfolding show that these molten globule states contain 28-67% of buried surface area relative to the native states. The data suggest that molten globules of these periplasmic binding proteins retain a considerable degree of long range order. The ability of these sequentially unrelated proteins to form highly ordered molten globules may be related to their large size as well as an intrinsic property of periplasmic binding protein folds.

Distributed space-time codes for cooperative networks with partial CSI

Design criteria and full-diversity Distributed Space Time Codes (DSTCs) for the two phase transmission based cooperative diversity protocol of Jing-Hassibi and the Generalized Nonorthogonal Amplify and Forward (GNAF) protocol are reported, when the relay nodes are assumed to have knowledge of the phase component of the source to relay channel gains. It is shown that this under this partial channel state information (CSI), several well known space time codes for the colocated MIMO (Multiple Input Multiple Output) channel become amenable for use as DSTCs. In particular, the well known complex orthogonal designs, generalized coordinate interleaved orthogonal designs (GCIODs) and unitary weight single symbol decodable (UW-SSD) codes are shown to satisfy the required design constraints for DSTCs. Exploiting the relaxed code design constraints, we propose DSTCs obtained from Clifford Algebras which have low ML decoding complexity.

Low PAPR STBCs from Complex Partial-Orthogonal Designs (CPODs)

Space-time codes from complex orthogonal designs (CODs) with no zero entries offer low Peak to Average power ratio (PAPR) and avoid the problem of turning off antennas. But CODs for 2(a) antennas with a + 1 complex variables, with no zero entries are not known in the literature for a >= 4. In this paper, a method of obtaining no zero entry (NZE) codes, called Complex Partial-Orthogonal Designs (CPODs), for 2(a+1) antennas whenever a certain type of NZE code exists for 2(a) antennas is presented. This is achieved with slight increase in the ML decoding complexity for regular QAM constellations and no increase for other complex constellations. Since NZE CODs have been constructed recently for 8 antennas our method leads to NZE CPODs for 16 antennas. Moreover, starting from certain NZE CPODs for n antennas, a construction procedure is given to obtain NZE CPODs for 2n antennas. The class of CPODs do not offer full-diversity for all complex constellations. For the NZE CPODs presented in the paper, conditions on the signal sets which will guarantee full-diversity are identified. Simulations results show that bit error performance of our codes under average power constraint is same as that of the CODs and superior to CODs under peak power constraint.

Identification and characterization of DdRPB4, a subunit of Dictyostelium discoideum RNA polymerase II

Rpb4, the fourth largest subunit of the eukaryotic RNA polymerase II (RNAPII), is required for growth at extreme temperatures and for an appropriate response to nutrient starvation in yeast. Sequence homologs of Rpb4 are found in most sequenced genomes from yeast to humans. To elucidate the role of this subunit in nutrient starvation, we chose Dictyostelium discoideum, a soil amoeba, which responds to nutrient deprivation by undergoing a complex developmental program. Here we report the identification of homolog of Saccharomyces cerevisiae RPB4 in D. discoideum. Localization and complementation studies suggest that Rpb4 is functionally conserved. DdRPB4 transcript and protein levels are developmentally regulated. Although DdRPB4 could not be deleted, overexpression revealed that the Rpb4 protein is essential for cell survival and is regulated stringently at the post-transcriptional level in D. discoideum. Thus maintaining a critical level of Rpb4 is important for this organism.

Novel melody line identification algorithm for polyphonic MIDI music

The problem of automatic melody line identification in a MIDI file plays an important role towards taking QBH systems to the next level. We present here, a novel algorithm to identify the melody line in a polyphonic MIDI file. A note pruning and track/channel ranking method is used to identify the melody line. We use results from musicology to derive certain simple heuristics for the note pruning stage. This helps in the robustness of the algorithm, by way of discarding "spurious" notes. A ranking based on the melodic information in each track/channel enables us to choose the melody line accurately. Our algorithm makes no assumption about MIDI performer specific parameters, is simple and achieves an accuracy of 97% in identifying the melody line correctly. This algorithm is currently being used by us in a QBH system built in our lab.

Nonlinear structural dynamical system identification using adaptive particle filters

The problem of identifying parameters of nonlinear vibrating systems using spatially incomplete, noisy, time-domain measurements is considered. The problem is formulated within the framework of dynamic state estimation formalisms that employ particle filters. The parameters of the system, which are to be identified, are treated as a set of random variables with finite number of discrete states. The study develops a procedure that combines a bank of self-learning particle filters with a global iteration strategy to estimate the probability distribution of the system parameters to be identified. Individual particle filters are based on the sequential importance sampling filter algorithm that is readily available in the existing literature. The paper develops the requisite recursive formulary for evaluating the evolution of weights associated with system parameter states. The correctness of the formulations developed is demonstrated first by applying the proposed procedure to a few linear vibrating systems for which an alternative solution using adaptive Kalman filter method is possible. Subsequently, illustrative examples on three nonlinear vibrating systems, using synthetic vibration data, are presented to reveal the correct functioning of the method. (c) 2007 Elsevier Ltd. All rights reserved.

Identification of protein interaction regions of VINC/NEAT1/Men epsilon RNA

The virus inducible non-coding RNA (VINC) was detected initially in the brain of mice infected with Japanese encephalitis virus (JEV) and rabies virus. VINC is also known as NEAT1 or Men epsilon RNA. It is localized in the nuclear paraspeckles of several murine as well as human cell lines and is essential for paraspeckle formation. We demonstrate that VINC interacts with the paraspeckle protein, P54nrb through three different protein interaction regions (PIRs) one of which (PIR-1) is localized near the 50 end while the other two (PIR-2, PIR-3) are localized near the 30 region of VINC. Our studies suggest that VINC may interact with P54nrb through a novel mechanism which is different from that reported for protein coding RNAs. (C) 2010 Federation of European Biochemical Societies. Published by Elsevier B. V. All rights reserved.