The full Bayesian significance test for mixture models: results in gene expression clustering

Gene clustering is a useful exploratory technique to group together genes with similar expression levels under distinct cell cycle phases or distinct conditions. It helps the biologist to identify potentially meaningful relationships between genes. In this study, we propose a clustering method based on multivariate normal mixture models, where the number of clusters is predicted via sequential hypothesis tests: at each step, the method considers a mixture model of m components (m = 2 in the first step) and tests if in fact it should be m - 1. If the hypothesis is rejected, m is increased and a new test is carried out. The method continues (increasing m) until the hypothesis is accepted. The theoretical core of the method is the full Bayesian significance test, an intuitive Bayesian approach, which needs no model complexity penalization nor positive probabilities for sharp hypotheses. Numerical experiments were based on a cDNA microarray dataset consisting of expression levels of 205 genes belonging to four functional categories, for 10 distinct strains of Saccharomyces cerevisiae. To analyze the method's sensitivity to data dimension, we performed principal components analysis on the original dataset and predicted the number of classes using 2 to 10 principal components. Compared to Mclust (model-based clustering), our method shows more consistent results.

Possibility for a full optical determination of photodynamic therapy outcome

The efficacy of photodynamic therapy (PDT) depends on a variety of parameters: concentration of the photosensitizer at the time of treatment, light wavelength, fluence, fluence rate, availability of oxygen within the illuminated volume, and light distribution in the tissue. Dosimetry in PDT requires the congregation of adequate amounts of light, drug, and tissue oxygen. The adequate dosimetry should be able to predict the extension of the tissue damage. Photosensitizer photobleaching rate depends on the availability of molecular oxygen in the tissue. Based on photosensitizers photobleaching models, high photobleaching has to be associated with high production of singlet oxygen and therefore with higher photodynamic action, resulting in a greater depth of necrosis. The purpose of this work is to show a possible correlation between depth of necrosis and the in vivo photosensitizer (in this case, Photogem (R)) photodegradation during PDT. Such correlation allows possibilities for the development of a real time evaluation of the photodynamic action during PDT application. Experiments were performed in a range of fluence (0-450 J/cm(2)) at a constant fluence rate of 250 mW/cm(2) and applying different illumination times (0-1800 s) to achieve the desired fluence. A quantity was defined (psi) as the product of fluorescence ratio (related to the photosensitizer degradation at the surface) and the observed depth of necrosis. The correlation between depth of necrosis and surface fluorescence signal is expressed in psi and could allow, in principle, a noninvasive monitoring of PDT effects during treatment. High degree of correlation is observed and a simple mathematical model to justify the results is presented.

The Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS) - Part 2: Model sensitivity to the biomass burning

We describe an estimation technique for biomass burning emissions in South America based on a combination of remote-sensing fire products and field observations, the Brazilian Biomass Burning Emission Model (3BEM). For each fire pixel detected by remote sensing, the mass of the emitted tracer is calculated based on field observations of fire properties related to the type of vegetation burning. The burnt area is estimated from the instantaneous fire size retrieved by remote sensing, when available, or from statistical properties of the burn scars. The sources are then spatially and temporally distributed and assimilated daily by the Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS) in order to perform the prognosis of related tracer concentrations. Three other biomass burning inventories, including GFEDv2 and EDGAR, are simultaneously used to compare the emission strength in terms of the resultant tracer distribution. We also assess the effect of using the daily time resolution of fire emissions by including runs with monthly-averaged emissions. We evaluate the performance of the model using the different emission estimation techniques by comparing the model results with direct measurements of carbon monoxide both near-surface and airborne, as well as remote sensing derived products. The model results obtained using the 3BEM methodology of estimation introduced in this paper show relatively good agreement with the direct measurements and MOPITT data product, suggesting the reliability of the model at local to regional scales.

Elastic Maier-Saupe-Zwanzig model and some properties of nematic elastomers

We introduce a simple mean-field lattice model to describe the behavior of nematic elastomers. This model combines the Maier-Saupe-Zwanzig approach to liquid crystals and an extension to lattice systems of the Warner-Terentjev theory of elasticity, with the addition of quenched random fields. We use standard techniques of statistical mechanics to obtain analytic solutions for the full range of parameters. Among other results, we show the existence of a stress-strain coexistence curve below a freezing temperature, analogous to the P-V diagram of a simple fluid, with the disorder strength playing the role of temperature. Below a critical value of disorder, the tie lines in this diagram resemble the experimental stress-strain plateau and may be interpreted as signatures of the characteristic polydomain-monodomain transition. Also, in the monodomain case, we show that random fields may soften the first-order transition between nematic and isotropic phases, provided the samples are formed in the nematic state.

Slavnov-Taylor identities for the 2+1 dimensional noncommutative CPN-1 model

In the context of the 1/N expansion, the validity of the Slavnov-Taylor identity relating three- and two-point functions for the 2 + 1-dimensional noncommutative CP(N-1) model is investigated, up to subleading 1/N order, in the Landau gauge.

Dynamical generation of mass in the D = (2+1) Wess-Zumino model

In this work we study the dynamical generation of mass in the massless N = 1 Wess-Zumino model in a three-dimensional spacetime. Using the tadpole method to compute the effective potential, we observe that supersymmetry is dynamically broken together with the discrete symmetry A(x) -> A(x). We show that this model, different from nonsupersymmetric scalar models, exhibits a consistent perturbative dynamical generation of mass after two-loop corrections to the effective potential.

Crystallization and preliminary structural analysis of the giant haemoglobin from Glossoscolex paulistus at 3.2 angstrom

Glossoscolex paulistus is a free-living earthworm encountered in south-east Brazil. Its oxygen transport requirements are undertaken by a giant extracellular haemoglobin, or erythrocruorin (HbGp), which has an approximate molecular mass of 3.6 MDa and, by analogy with its homologue from Lumbricus terrestris (HbLt), is believed to be composed of a total of 180 polypeptide chains. In the present work the full 3.6 MDa particle in its cyanomet state was purified and crystallized using sodium citrate or PEG8000 as precipitant. The crystals contain one-quarter of the full particle in the asymmetric unit of the I222 cell and have parameters of a = 270.8 angstrom, b = 320.3 angstrom and c = 332.4 angstrom. Diffraction data were collected to 3.15 angstrom using synchrotron radiation on beamline X29A at the Brookhaven National Laboratory and represent the highest resolution data described to date for similar erythrocruorins. The structure was solved by molecular replacement using a search model corresponding to one-twelfth of its homologue from HbLt. This revealed that HbGp belongs to the type I class of erythrocruorins and provided an interpretable initial electron density map in which many features including the haem groups and disulfide bonds could be identified.

Integrating structural and input design of a 2-DOF high-speed parallel manipulator: A flexible model-based approach

This paper discusses the integrated design of parallel manipulators, which exhibit varying dynamics. This characteristic affects the machine stability and performance. The design methodology consists of four main steps: (i) the system modeling using flexible multibody technique, (ii) the synthesis of reduced-order models suitable for control design, (iii) the systematic flexible model-based input signal design, and (iv) the evaluation of some possible machine designs. The novelty in this methodology is to take structural flexibilities into consideration during the input signal design; therefore, enhancing the standard design process which mainly considers rigid bodies dynamics. The potential of the proposed strategy is exploited for the design evaluation of a two degree-of-freedom high-speed parallel manipulator. The results are experimentally validated. (C) 2010 Elsevier Ltd. All rights reserved.

AnSBBR Applied to the Treatment of Metalworking Fluid Wastewater: Effect of Organic and Shock Load

An investigation was performed regarding the application of a mechanically stirred anaerobic sequencing batch biofilm reactor containing immobilized biomass on inert polyurethane foam (AnSBBR) to the treatment of soluble metalworking fluids to remove organic matter and produce methane. The effect of increasing organic matter and reactor fill time, as well as shock load, on reactor stability and efficiency have been analyzed. The 5-L AnSBBR was operated at 30 A degrees C in 8-h cycles, agitation of 400 rpm, and treated 2.0 L effluent per cycle. Organic matter was increased by increasing the influent concentration (500, 1,000, 2,000, and 3,000 mg chemical oxygen demand (COD)/L). Fill times investigated were in the batch mode (fill time 10 min) and fed-batch followed by batch (fill time 4 h). In the batch mode, organic matter removal efficiencies were 87%, 86%, and 80% for influent concentrations of 500, 1,000, and 2,000 mgCOD/L (1.50, 3.12, and 6.08 gCOD/L.d), respectively. At 3,000 mgCOD/L (9.38 gCOD/L.d), operational stability could not be achieved. The reactor managed to maintain stability when a shock load twice as high the feed concentration was applied, evidencing the robustness of the reactor to potential concentration variations in the wastewater being treated. Increasing the fill time to 4 h did not improve removal efficiency, which was 72% for 2,000 mgCOD/L. Thus, gradual feeding did not improve organic matter removal. The concentration of methane formed at 6.08 gCOD/L was 5.20 mmolCH(4), which corresponded to 78% of the biogas composition. The behavior of the reactor during batch and fed-batch feeding could be explained by a kinetic model that considers organic matter consumption, production, and consumption of total volatile acids and methane production.

Real time optimization (RTO) with model predictive control (MPC)

This paper studies a simplified methodology to integrate the real time optimization (RTO) of a continuous system into the model predictive controller in the one layer strategy. The gradient of the economic objective function is included in the cost function of the controller. Optimal conditions of the process at steady state are searched through the use of a rigorous non-linear process model, while the trajectory to be followed is predicted with the use of a linear dynamic model, obtained through a plant step test. The main advantage of the proposed strategy is that the resulting control/optimization problem can still be solved with a quadratic programming routine at each sampling step. Simulation results show that the approach proposed may be comparable to the strategy that solves the full economic optimization problem inside the MPC controller where the resulting control problem becomes a non-linear programming problem with a much higher computer load. (C) 2010 Elsevier Ltd. All rights reserved.

Thermodynamic modelling of an Al(2)O(3)-MnO system using the ionic model

The thermodynamic assessment of an Al(2)O(3)-MnO pseudo-binary system has been carried out with the use of an ionic model. The use of the electro-neutrality principles in addition to the constitutive relations, between site fractions of the species on each sub-lattice, the thermodynamics descriptions of each solid phase has been determined to make possible the solubility description. Based on the thermodynamics descriptions of each phase in addition to thermo-chemical data obtained from the literature, the Gibbs energy functions were optimized for each phase of the Al(2)O(3)-MnO system with the support of PARROT(R) module from ThemoCalc(R) package. A thermodynamic database was obtained, in agreement with the thermo-chemical data extracted from the literature, to describe the Al(2)O(3)-MnO system including the solubility description of solid phases. (C) 2009 Elsevier Ltd. All rights reserved.

Stability and in vitro penetration study of rutin incorporated in a cosmetic emulsion through an alternative model biomembrane

Rutin is employed as antioxidant and to prevent the capillary fragility and, when incorporated in cosmetic emulsions, it must target the action site. In vitro cutaneous penetration studies through human skin is the ideal situation, however, there are difficulties to obtain and to maintain this tissue viability. Among the membrane models, shed snake skin presents itself as pure stratum corneum, providing barrier function similar to human and it is obtained without the animal sacrifice. The objectives of this research were the development and stability evaluation of a cosmetic emulsion containing rutin and propylene glycol (penetration enhancer) and the evaluation or rutin in vitro cutaneous penetration and retention from the emulsion, employing an alternative model biomembrane. Emulsion was developed with rutin and propylene glycol, both at 5.0% w/w. Active substance presented on the formulation was quantified by a validated spectrophotometric method at 361.0 nm. Rutin Rutin cutaneous penetration and retention was performed in vertical diffusion cells with shed snake skin of Crotalus durissus, as alternative model biomembrane, and distilled water and ethanol 99.5% (1:1), as receptor fluid. The experiment was conducted for six hours, at 37.0 +/- 0.5 degrees C with constant stirring of 300 rpm. Spectrophotometry at 410.0 nm, previously validated, determined the active substance after cutaneous penetration/ retention. Emulsion did not promote rutin cutaneous penetration through C. durissus skin, retaining 0.931 +/- 0.0391 mu g rutin/mg shed snake skin. The referred formulation was chemically stable for 30 days after stored at 25.0 +/- 2.0 degrees C, 5.0 +/- 0.5 degrees C and 45.0 +/- 0.5 degrees C. In conclusion, it has not been verified the active cutaneous penetration through the model biomembrane, but only its retention on the Crotalus durissus stratum corneum, condition considered stable for 30 days.

Magnetic Investigation of CoFe(2)O(4) Nanoparticles Supported in Biocompatible Polymeric Microsphere

Magnetic investigation of spinel ferrite nanoparticles dispersed in biocompatible polymeric microspheres is reported in this study. X-ray diffraction data analysis confirms the presence of nanosized CoFe(2)O(4) particles (mean size of similar to 8 nm). This finding is corroborated by transmission electron microscopy micrographs. Magnetization isotherms suggest a spin disorder likely occurring at the nanoparticle`s surface. The saturation magnetization value is used to estimate particle concentration of 1.6 x 10(18) cm(-3) dispersed in the polymeric template. A T(1/2) dependence of the coercive field is determined in the low-temperature region (T < 30 K). The model of non-interacting mono-domains is used to estimate an effective magnetic anisotropy of K(eff) = 0.6 x 10(5) J/m(3). The K(eff) value we found is lower than the value reported for spherically-shaped CoFe(2)O(4) nanoparticles, though consistent with the low coercive field observed in the investigated sample.

A Comparison Between Pulse Pressure Variation and Right End Diastolic Volume Index as Guides to Resuscitation in a Model of Hemorrhagic Shock in Pigs

Background: Different hemodynamic parameters including static indicators of cardiac preload as right ventricular end-diastolic volume index (RVEDVI) and dynamic parameters as pulse pressure variation (PPV) have been used in the decision-making process regarding volume expansion in critically ill patients. The objective of this study was to compare fluid resuscitation guided by either PPV or RVEDVI after experimentally induced hemorrhagic shock. Methods: Twenty-six anesthetized and mechanically ventilated pigs were allocated into control (group I), PPV (group II), or RVEDVI (group III) group. Hemorrhagic shock was induced by blood withdrawal to target mean arterial pressure of 40 mm Hg, maintained for 60 minutes. Parameters were measured at baseline, time of shock, 60 minutes after shock, immediately after resuscitation with hydroxyethyl starch 6% (130/0.4), 1 hour and 2 hours thereafter. The endpoint of fluid resuscitation was determined as the baseline values of PPV and RVEDVI. Statistical analysis of data was based on analysis of variance for repeated measures followed by the Bonferroni test (p < 0.05). Results: Volume and time to resuscitation were higher in group III than in group II (group III = 1,305 +/- 331 mL and group II = 965 +/- 245 mL, p < 0.05; and group III = 24.8 +/- 4.7 minutes and group II = 8.8 +/- 1.3 minutes, p < 0.05, respectively). All static and dynamic parameters and biomarkers of tissue oxygenation were affected by hemorrhagic shock and nearly all parameters were restored after resuscitation in both groups. Conclusion: In the proposed model of hemorrhagic shock, resuscitation to the established endpoints was achieved within a smaller amount of time and with less volume when guided by PPV than when guided by pulmonary artery catheter-derived RVEDVI.

Derivation and external validation of a simple prediction model for the diagnosis of type 2 Diabetes Mellitus in the Brazilian urban population

A risk score model was developed based in a population of 1,224 individuals from the general population without known diabetes aging 35 years or more from an urban Brazilian population sample in order to select individuals who should be screened in subsequent testing and improve the efficacy of public health assurance. External validation was performed in a second, independent, population from a different city ascertained through a similar epidemiological protocol. The risk score was developed by multiple logistic regression and model performance and cutoff values were derived from a receiver operating characteristic curve. Model`s capacity of predicting fasting blood glucose levels was tested analyzing data from a 5-year follow-up protocol conducted in the general population. Items independently and significantly associated with diabetes were age, BMI and known hypertension. Sensitivity, specificity and proportion of further testing necessary for the best cutoff value were 75.9, 66.9 and 37.2%, respectively. External validation confirmed the model`s adequacy (AUC equal to 0.72). Finally, model score was also capable of predicting fasting blood glucose progression in non-diabetic individuals in a 5-year follow-up period. In conclusion, this simple diabetes risk score was able to identify individuals with an increased likelihood of having diabetes and it can be used to stratify subpopulations in which performing of subsequent tests is necessary and probably cost-effective.