Characterization of subgraph relationships and distribution in complex networks

A network can be analyzed at different topological scales, ranging from single nodes to motifs, communities, up to the complete structure. We propose a novel approach which extends from single nodes to the whole network level by considering non-overlapping subgraphs (i.e. connected components) and their interrelationships and distribution through the network. Though such subgraphs can be completely general, our methodology focuses on the cases in which the nodes of these subgraphs share some special feature, such as being critical for the proper operation of the network. The methodology of subgraph characterization involves two main aspects: (i) the generation of histograms of subgraph sizes and distances between subgraphs and (ii) a merging algorithm, developed to assess the relevance of nodes outside subgraphs by progressively merging subgraphs until the whole network is covered. The latter procedure complements the histograms by taking into account the nodes lying between subgraphs, as well as the relevance of these nodes to the overall subgraph interconnectivity. Experiments were carried out using four types of network models and five instances of real-world networks, in order to illustrate how subgraph characterization can help complementing complex network-based studies.

Bayesian online algorithms for learning in discrete Hidden Markov Models

We propose and analyze two different Bayesian online algorithms for learning in discrete Hidden Markov Models and compare their performance with the already known Baldi-Chauvin Algorithm. Using the Kullback-Leibler divergence as a measure of generalization we draw learning curves in simplified situations for these algorithms and compare their performances.

Upconversion in Nd(3+)-doped glasses: Microscopic theory and spectroscopic measurements

In this work, we report a systematic investigation of upconversion losses and their effects on fluorescence quantum efficiency and fractional thermal loading in Nd(3+)-doped fluoride glasses. The energy transfer upconversion (gamma(up)) parameter, which describes upconversion losses, was experimentally determined using different methods: thermal lens (TL) technique and steady state luminescence (SSL) measurements. Additionally, the upconversion parameter was also obtained from energy transfer models and excited state absorption measurements. The results reveal that the microscopic treatment provided by the energy transfer models is similar to the macroscopic ones achieved from the TL and SSL measurements because similar gamma(up) parameters were obtained. Besides, the achieved results also point out the migration-assisted energy transfer according to diffusion-limited regime rather than hopping regime as responsible for the upconversion losses in Nd-doped glasses. (c) 2008 American Institute of Physics.

Correlation between MMPs and their inhibitors in breast cancer tumor tissue specimens and in cell lines with different metastatic potential

Background: The metastatic disease rather than the primary tumor itself is responsible for death in most solid tumors, including breast cancer. The role of matrix metalloproteinases ( MMPs), tissue inhibitors of MMPs (TIMPs) and Reversion-inducing cysteine-rich protein with Kazal motifs ( RECK) in the metastatic process has previously been established. However, in all published studies only a limited number of MMPs/MMP inhibitors was analyzed in a limited number of cell lines. Here, we propose a more comprehensive approach by analyzing the expression levels of several MMPs (MMP-2, MMP-9 and MMP-14) and MMP inhibitors (TIMP-1, TIMP-2 and RECK) in different models ( five human breast cancer cell lines, 72 primary breast tumors and 30 adjacent normal tissues). Methods: We analyzed the expression levels of MMP-2, MMP-9 and MMP-14 and their inhibitors (TIMP-1, TIMP-2 and RECK) by quantitative RT-PCR (qRT-PCR) in five human breast cancer cell lines presenting increased invasiveness and metastatic potential, 72 primary breast tumors and 30 adjacent normal tissues. Moreover, the role of cell-extracellular matrix elements interactions in the regulation of expression and activity of MMPs and their inhibitors was analyzed by culturing these cell lines on plastic or on artificial ECM (Matrigel). Results: The results demonstrated that MMPs mRNA expression levels displayed a positive and statistically significant correlation with the transcriptional expression levels of their inhibitors both in the cell line models and in the tumor tissue samples. Furthermore, the expression of all MMP inhibitors was modulated by cell-Matrigel contact only in highly invasive and metastatic cell lines. The enzyme/inhibitor balance at the transcriptional level significantly favors the enzyme which is more evident in tumor than in adjacent non-tumor tissue samples. Conclusion: Our results suggest that the expression of MMPs and their inhibitors, at least at the transcriptional level, might be regulated by common factors and signaling pathways. Therefore, the multi-factorial analysis of these molecules could provide new and independent prognostic information contributing to the determination of more adequate therapy strategies for each patient.`

Laser-induced breakdown spectroscopy and chemometrics for classification of toys relying on toxic elements

Quality control of toys for avoiding children exposure to potentially toxic elements is of utmost relevance and it is a common requirement in national and/or international norms for health and safety reasons. Laser-induced breakdown spectroscopy (LIBS) was recently evaluated at authors` laboratory for direct analysis of plastic toys and one of the main difficulties for the determination of Cd. Cr and Pb was the variety of mixtures and types of polymers. As most norms rely on migration (lixiviation) protocols, chemometric classification models from LIBS spectra were tested for sampling toys that present potential risk of Cd, Cr and Pb contamination. The classification models were generated from the emission spectra of 51 polymeric toys and by using Partial Least Squares - Discriminant Analysis (PLS-DA), Soft Independent Modeling of Class Analogy (SIMCA) and K-Nearest Neighbor (KNN). The classification models and validations were carried out with 40 and 11 test samples, respectively. Best results were obtained when KNN was used, with corrected predictions varying from 95% for Cd to 100% for Cr and Pb. (C) 2011 Elsevier B.V. All rights reserved.

THREE-DIMENSIONAL PHOTOIONIZATION STRUCTURE AND DISTANCES OF PLANETARY NEBULAE. IV. NGC 40

Continuing our series of papers on the three-dimensional (3D) structure and accurate distances of planetary nebulae (PNe), we present here the results obtained for PN NGC 40. Using data from different sources and wavelengths, we construct 3D photoionization models and derive the physical quantities of the ionizing source and nebular gas. The procedure, discussed in detail in the previous papers, consists of the use of 3D photoionization codes constrained by observational data to derive the 3D nebular structure, physical and chemical characteristics, and ionizing star parameters of the objects by simultaneously fitting the integrated line intensities, the density map, the temperature map, and the observed morphologies in different emission lines. For this particular case we combined hydrodynamical simulations with the photoionization scheme in order to obtain self-consistent distributions of density and velocity of the nebular material. Combining the velocity field with the emission-line cubes we also obtained the synthetic position-velocity plots that are compared to the observations. Finally, using theoretical evolutionary tracks of intermediate-and low-mass stars, we derive the mass and age of the central star of NGC 40 as (0.567 +/- 0.06) M(circle dot) and (5810 +/- 600) yr, respectively. The distance obtained from the fitting procedure was (1150 +/- 120) pc.

Financial distress, financial constraint and investment decision: Evidence from Brazil

This paper analyses the presence of financial constraint in the investment decisions of 367 Brazilian firms from 1997 to 2004, using a Bayesian econometric model with group-varying parameters. The motivation for this paper is the use of clustering techniques to group firms in a totally endogenous form. In order to classify the firms we used a hybrid clustering method, that is, hierarchical and non-hierarchical clustering techniques jointly. To estimate the parameters a Bayesian approach was considered. Prior distributions were assumed for the parameters, classifying the model in random or fixed effects. Ordinate predictive density criterion was used to select the model providing a better prediction. We tested thirty models and the better prediction considers the presence of 2 groups in the sample, assuming the fixed effect model with a Student t distribution with 20 degrees of freedom for the error. The results indicate robustness in the identification of financial constraint when the firms are classified by the clustering techniques. (C) 2010 Elsevier B.V. All rights reserved.

Fast experimental identification of suspension models for control design

This paper presents both the theoretical and the experimental approaches of the development of a mathematical model to be used in multi-variable control system designs of an active suspension for a sport utility vehicle (SUV), in this case a light pickup truck. A complete seven-degree-of-freedom model is successfully quickly identified, with very satisfactory results in simulations and in real experiments conducted with the pickup truth. The novelty of the proposed methodology is the use of commercial software in the early stages of the identification to speed up the process and to minimize the need for a large number of costly experiments. The paper also presents major contributions to the identification of uncertainties in vehicle suspension models and in the development of identification methods using the sequential quadratic programming, where an innovation regarding the calculation of the objective function is proposed and implemented. Results from simulations of and practical experiments with the real SUV are presented, analysed, and compared, showing the potential of the method.

Adaptation for tandem cold mill models

The ideal conditions for the operation of tandem cold mills are connected to a set of references generated by models and used by dynamic regulators. Aiming at the optimization of the friction and yield stress coefficients an adaptation algorithm is proposed in this paper. Experimental results obtained from an industrial cold rolling mill are presented. (C) 2008 Elsevier B.V. All rights reserved.

Cloning, expression, molecular modelling and docking analysis of glutathione transferase from Saccharum officinarum

Sugarcane yield and quality are affected by a number of biotic and abiotic stresses. In response to such stresses, plants may increase the activities of some enzymes such as glutathione transferase (GST), which are involved in the detoxification of xenobiotics. Thus, a sugarcane GST was modelled and molecular docked using the program LIGIN to investigate the contributions of the active site residues towards the binding of reduced glutathione (GSH) and 1-chloro-2,4-dinitrobenzene (CDNB). As a result, W13 and I119 were identified as key residues for the specificity of sugarcane GSTF1 (SoGSTF1) towards CDNB. To obtain a better understanding of the catalytic specificity of sugarcane GST (SoGSTF1), two mutants were designed, W13L and I119F. Tertiary structure models and the same docking procedure were performed to explain the interactions between sugarcane GSTs with GSH and CDNB. An electron-sharing network for GSH interaction was also proposed. The SoGSTF1 and the mutated gene constructions were cloned and expressed in Escherichia coli and the expressed protein purified. Kinetic analyses revealed different Km values not only for CDNB, but also for GSH. The Km values were 0.2, 1.3 and 0.3 mM for GSH, and 0.9, 1.2 and 0.5 mM for CDNB, for the wild type, W13L mutant and I119F mutant, respectively. The V(max) values were 297.6, 224.5 and 171.8 mu mol min(-1) mg(-1) protein for GSH, and 372.3, 170.6 and 160.4 mu mol min(-1) mg(-1) protein for CDNB.

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.

Within-stand and seasonal variations of specific leaf area in a clonal Eucalyptus plantation in the Republic of Congo

Specific leaf area (SLA; m(leaf)(2) kg(leaf)(-1)) is a key ecophysiological parameter influencing leaf physiology, photosynthesis, and whole plant carbon gain. Both individual tree-based models and other forest process-based models are generally highly sensitive to this parameter, but information on its temporal or within-stand variability is still scarce. In a 2-4-year-old Eucalyptus plantation in Congo, prone to seasonal drought, the within-stand and seasonal variability in SLA were investigated by means of destructive sampling carried out at 2-month intervals, over a 2-year period. Within-crown vertical gradients of SLA were small. Highly significant relationships were found between tree-average SLA (SLA(t)) and tree size (tree height, H(t), or diameter at breast height, DBH): SLA(t) ranged from about 9 m(2) kg(-1) for dominant trees to about 14-15 m(2) kg(-1) for the smallest trees. The decrease in SLA(t) with increasing tree size was accurately predicted from DBH using power functions. Stand-average SLA varied by about 20% during the year, with lowest values at the end of the 5-month dry season, and highest values about 2-3 months after the onset of the wet season. Variability in leaf water status according to tree size and season is discussed as a possible determinant of both the within-stand and seasonal variations in SM. (C) 2009 Elsevier B.V. All rights reserved.

Rational Design and 3D-Pharmacophore Mapping of 5 `-Thiourea-Substituted alpha-Thymidine Analogues as Mycobacterial TMPK Inhibitors

Thymidine monophosphate kinase (TMPK) has emerged as an attractive target for developing inhibitors of Mycobacterium tuberculosis growth. In this study the receptor-independent (RI) 4D-QSAR formalism has been used to develop QSAR models and corresponding 3D-pharmacophores for a set of 5`-thiourea-substituted alpha-thymidine inhibitors. Models were developed for the entire training set and for a subset of the training set consisting of the most potent inhibitors. The optimized (RI) 4D-QSAR models are statistically significant (r(2) = 0.90, q(2) = 0.83 entire set, r(2) = 0.86, q(2) = 0.80 high potency subset) and also possess good predictivity based on test set predictions. The most and least potent inhibitors, in their respective postulated active conformations derived from the models, were docked in the active site of the TMPK crystallographic structure. There is a solid consistency between the 3D-pharmacophore sites defined by the QSAR models and interactions with binding site residues. This model identifies new regions of the inhibitors that contain pharmacophore sites, such as the sugar-pyrimidine ring structure and the region of the 5`-arylthiourea moiety. These new regions of the ligands can be further explored and possibly exploited to identify new, novel, and, perhaps, better antituberculosis inhibitors of TMPKmt. Furthermore, the 3D-pharmacophores defined by these models can be used as a starting point for future receptor-dependent antituberculosis drug design as well as to elucidate candidate sites for substituent addition to optimize ADMET properties of analog inhibitors.

Polytomous item response theory models

Polytomous Item Response Theory Models provides a unified, comprehensive introduction to the range of polytomous models available within item response theory (IRT). It begins by outlining the primary structural distinction between the two major types of polytomous IRT models. This focuses on the two types of response probability that are unique to polytomous models and their associated response functions, which are modeled differently by the different types of IRT model. It describes, both conceptually and mathematically, the major specific polytomous models, including the Nominal Response Model, the Partial Credit Model, the Rating Scale model, and the Graded Response Model. Important variations, such as the Generalized Partial Credit Model are also described as are less common variations, such as the Rating Scale version of the Graded Response Model. Relationships among the models are also investigated and the operation of measurement information is described for each major model. Practical examples of major models using real data are provided, as is a chapter on choosing an appropriate model. Figures are used throughout to illustrate important elements as they are described.

The kinetics of NO and N2O reduction over coal chars in fluidised-bed combustion

This paper presents a comprehensive and critical review of the mechanisms and kinetics of NO and N2O reduction reaction with coal chars under fluidised-bed combustion conditions (FBC). The heterogeneous reactions of NO and N2O with char/carbon surface have been well recognised as the most important processes in reducing both NOx and N2O in situ FBC. Compared to NO-carbon reactions in FBC, the reactions of N2O with chars have been relatively less understood and studied. Beginning with the overall reaction schemes for both NO and N2O reduction, the paper extensively discusses the reaction mechanisms including the effects of active surface sites. Generally, NO- and N2O-carbon reactions follow a series of step reactions. However, questions remain concerning the role of adsorbed phases of NO and N2O, and the behaviour of different surface sites. Important kinetics factors such as the rate expressions, kinetics parameters as well as the effects of surface area and pore structure are discussed in detail. The main factors influencing the reduction of NO and N2O in FBC conditions are the chemical and physical properties of chars, and the operating parameters of FBC such as temperature, presence of CO, O-2 and pressure. It is shown that under similar conditions, N2O is more readily reduced on the char surface than NO. Temperature was found to be a very important parameter in both NO and N2O reduction. It is generally agreed that both NO- and N2O-carbon reactions follow first-order reaction kinetics with respect to the NO and N2O concentrations. The kinetic parameters for NO and N2O reduction largely depend on the pore structure of chars. The correlation between the char surface area and the reactivities of NO/N2O-char reactions is considered to be of great importance to the determination of the reaction kinetics. The rate of NO reduction by chars is strongly enhanced by the presence of CO and O-2, but these species may not have significant effects on the rate of N2O reduction. However, the presence of these gases in FBC presents difficulties in the study of kinetics since CO cannot be easily eliminated from the carbon surface. In N2O reduction reactions, ash in chars is found to have significant catalytic effects, which must be accounted for in the kinetic models and data evaluation. (C) 1997 Elsevier Science Ltd.