Identifying regulational alterations in gene regulatory networks by state space representation of vector autoregressive models and variational annealing

Background: In the analysis of effects by cell treatment such as drug dosing, identifying changes on gene network structures between normal and treated cells is a key task. A possible way for identifying the changes is to compare structures of networks estimated from data on normal and treated cells separately. However, this approach usually fails to estimate accurate gene networks due to the limited length of time series data and measurement noise. Thus, approaches that identify changes on regulations by using time series data on both conditions in an efficient manner are demanded. Methods: We propose a new statistical approach that is based on the state space representation of the vector autoregressive model and estimates gene networks on two different conditions in order to identify changes on regulations between the conditions. In the mathematical model of our approach, hidden binary variables are newly introduced to indicate the presence of regulations on each condition. The use of the hidden binary variables enables an efficient data usage; data on both conditions are used for commonly existing regulations, while for condition specific regulations corresponding data are only applied. Also, the similarity of networks on two conditions is automatically considered from the design of the potential function for the hidden binary variables. For the estimation of the hidden binary variables, we derive a new variational annealing method that searches the configuration of the binary variables maximizing the marginal likelihood. Results: For the performance evaluation, we use time series data from two topologically similar synthetic networks, and confirm that our proposed approach estimates commonly existing regulations as well as changes on regulations with higher coverage and precision than other existing approaches in almost all the experimental settings. For a real data application, our proposed approach is applied to time series data from normal Human lung cells and Human lung cells treated by stimulating EGF-receptors and dosing an anticancer drug termed Gefitinib. In the treated lung cells, a cancer cell condition is simulated by the stimulation of EGF-receptors, but the effect would be counteracted due to the selective inhibition of EGF-receptors by Gefitinib. However, gene expression profiles are actually different between the conditions, and the genes related to the identified changes are considered as possible off-targets of Gefitinib. Conclusions: From the synthetically generated time series data, our proposed approach can identify changes on regulations more accurately than existing methods. By applying the proposed approach to the time series data on normal and treated Human lung cells, candidates of off-target genes of Gefitinib are found. According to the published clinical information, one of the genes can be related to a factor of interstitial pneumonia, which is known as a side effect of Gefitinib.

Constitutive models for biodegradable thermoplastic ropes for ligament repair

The concept behind a biodegradable ligament device is to temporarily replace the biomechanical functions of the ruptured ligament, while it progressively regenerates its capacities. However, there is a lack of methods to predict the mechanical behaviour evolution of the biodegradable devices during degradation, which is an important aspect of the project. In this work, a hyper elastic constitutive model will be used to predict the mechanical behaviour of a biodegradable rope made of aliphatic polyesters. A numerical approach using ABAQUS is presented, where the material parameters of the model proposal are automatically updated in correspondence to the degradation time, by means of a script in PYTHON. In this method we also use a User Material subroutine (UMAT) to apply a failure criterion base on the strength that decreases according to a first order differential equation. The parameterization of the material model proposal for different degradation times were achieved by fitting the theoretical curves with the experimental data of tensile tests on fibres. To model all the rope behaviour we had considered one step of homogenisation considering the fibres architectures in an elementary volume. (C) 2012 Elsevier Ltd. All rights reserved.

Use of recent geoid models to estimate mean dynamic topography and geostrophic currents in South Atlantic and Brazil Malvinas confluence

The use of geoid models to estimate the Mean Dynamic Topography was stimulated with the launching of the GRACE satellite system, since its models present unprecedented precision and space-time resolution. In the present study, besides the DNSC08 mean sea level model, the following geoid models were used with the objective of computing the MDTs: EGM96, EIGEN-5C and EGM2008. In the method adopted, geostrophic currents for the South Atlantic were computed based on the MDTs. In this study it was found that the degree and order of the geoid models affect the determination of TDM and currents directly. The presence of noise in the MDT requires the use of efficient filtering techniques, such as the filter based on Singular Spectrum Analysis, which presents significant advantages in relation to conventional filters. Geostrophic currents resulting from geoid models were compared with the HYCOM hydrodynamic numerical model. In conclusion, results show that MDTs and respective geostrophic currents calculated with EIGEN-5C and EGM2008 models are similar to the results of the numerical model, especially regarding the main large scale features such as boundary currents and the retroflection at the Brazil-Malvinas Confluence.

Cation Distribution and Local Order in Mixed Sodium Metaphosphate Glasses

Two structural properties in mixed alkali metal phosphate glasses that seem to be crucial to the development of the mixed ion effect in dc conductivity were systematically analyzed in Na mixed metaphosphates: the local order around the mobile species, and their distribution and mixing in the glass network. The set of glasses considered here, Na1-xMxPO3 with M = Li, Ag, K, Rb, and Cs and 0 <= x <= 1, encompass a broad degree of size mismatch between the mixed cation species. A comprehensive solid-state nuclear magnetic resonance study was carried out using P-31 MAS, Na-23 triple quantum MAS, Rb-87 QCPMG, P-31-Na-23 REDOR, Na-23-Li-7 and Li-7-Li-6 SEDOR, and Na-23 spin echo decay. It was observed that the arrangement of P atoms around Na in the mixed glasses was indistinguishable from that observed in the NaPO3 glass. However, systematic distortions in the local structure of the 0 environments around Na were observed, related to the presence of the second cation. The average Na-O distances show an expansion/compression When Na+ ions are replaced by cations with respectively smaller/bigger radii. The behavior of the nuclear electric quadrupole coupling. constants indicates that this expansion reduces the local symmetry, while the compression produces the opposite effect These effects become marginally small when the site mismatch between the cations is small, as in Na-Ag mixed glasses. The present study confirms the intimate mixing of cation species at the atomic scale, but clear deviations from random mixing were detected in systems with larger alkali metal ions (Cs-Na, K-Na, Rb-Na). In contrast, no deviations from the statistical ion mixture were found in the systems Ag-Na and Li-Na, where mixed cations are either of radii comparable to (Ag+) or smaller than (Li+) Na+. The set of results supports two fundamental structural features of the models proposed to explain the mixed ion effect: the. structural specificity of the sites occupied by each cation species and their mixing at the atomic scale.

Analysis of elastic functionally graded materials under large displacements via high-order tetrahedral elements

The purpose of this study is to present a position based tetrahedral finite element method of any order to accurately predict the mechanical behavior of solids constituted by functionally graded elastic materials and subjected to large displacements. The application of high-order elements makes it possible to overcome the volumetric and shear locking that appears in usual homogeneous isotropic situations or even in non-homogeneous cases developing small or large displacements. The use of parallel processing to improve the computational efficiency, allows employing high-order elements instead of low-order ones with reduced integration techniques or strain enhancements. The Green-Lagrange strain is adopted and the constitutive relation is the functionally graded Saint Venant-Kirchhoff law. The equilibrium is achieved by the minimum total potential energy principle. Examples of large displacement problems are presented and results confirm the locking free behavior of high-order elements for non-homogeneous materials. (C) 2011 Elsevier B.V. All rights reserved.

Plasma corticosterone levels in mouse models of pain

Background: Pain markedly activates the hypothalamic-pituitary-adrenal (HPA) axis and increases plasma corticosterone release interfering significantly with nociceptive behaviour as well as the mechanism of action of analgesic drugs. Aims/Methods: In the present study, we monitored the time course of circulating corticosterone in two mouse strains (C57Bl/6 and Balb/C) under different pain models. In addition, the stress response was investigated following animal handling, intrathecal (i.t.) manipulation and habituation to environmental conditions commonly used in nociceptive experimental assays. We also examined the influence of within-cage order of testing on plasma corticosterone. Results: Subcutaneous injection of capsaicin precipitated a prompt stress response whereas carrageenan and complete Freund's adjuvant induced an increased corticosterone release around the third hour post-injection. However, carrageenan induced a longer increased corticosterone in C57Bl/6 mice. In partial sciatic nerve ligation, neuropathic pain model corticosterone increased only in the first days whereas mechanical hypersensitivity remained much longer. Animal handling also represents an important stressor whereas the i.t. injection per se does not exacerbate the handling-induced stress response. Moreover, the order of testing animals from the same cage does not interfere with plasma corticosterone levels in the intrathecal procedure. Animal habituation to the testing apparatus also does not reduce the immediate corticosterone increase as compared with non-habituated mice. Conclusion: Our data indicate that HPA axis activation in acute and chronic pain models is time dependent and may be dissociated from evoked hyperalgesia. Therefore, HPA-axis activation represents an important variable to be considered when designing experimental assays of persistent pain as well as for interpretation of data.

Exact cosmological solutions of models with an interacting dark sector

In this work we extend the first order formalism for cosmological models that present an interaction between a fermionic and a scalar field. Cosmological exact solutions describing universes filled with interacting dark energy and dark matter have been obtained. Viable cosmological solutions with an early period of decelerated expansion followed by late acceleration have been found, notably one which presents a dark matter component dominating in the past and a dark energy component dominating in the future. In another one, the dark energy alone is the responsible for both periods, similar to a Chaplygin gas case. Exclusively accelerating solutions have also been obtained.

Weibull and generalised exponential overdispersion models with an application to ozone air pollution

We consider the problem of estimating the mean and variance of the time between occurrences of an event of interest (inter-occurrences times) where some forms of dependence between two consecutive time intervals are allowed. Two basic density functions are taken into account. They are the Weibull and the generalised exponential density functions. In order to capture the dependence between two consecutive inter-occurrences times, we assume that either the shape and/or the scale parameters of the two density functions are given by auto-regressive models. The expressions for the mean and variance of the inter-occurrences times are presented. The models are applied to the ozone data from two regions of Mexico City. The estimation of the parameters is performed using a Bayesian point of view via Markov chain Monte Carlo (MCMC) methods.

Validation of the MAPGEO2010 and comparison with recent geopotential models

This paper presents the offorts to calculate the geoid model for Brazil. It is limited by 6 degrees N and 35 degrees S in latitude and 30 degrees W and 75 degrees W in longitude. The terrestrial gravity data for the continent have been updated by means of the most recent surveys in Brazil and in the neighbour countries. The complete Bouguer and Helmert gravity anomalies have been derived through the Canadian package SHGEO. The short wavelength component was estimated via FFT. The geopotential model EGM2008 was used as a reference field restricted to degree and order 150. The model was validated over 844 GPS observations on Bench Marks of the spirit leveling network. The height anomalies plus a topography dependent correction term derived from EGM2008 (degree 2190 and order 2159), GO_CONS_GCF_2_DIR_R2 (degree and order 240), GOCO02S (degree and order 250), EIGEN 51C (degree and order 359) and EIGEN 6C (degree and order 1420), geoidal height derived from MAPGEO2004 (old official geoid model in Brazil) have also been compared to the GPS points on Bench Marks.

Holographic field theory models of dark energy in interaction with dark matter

We discuss two Lagrangian interacting dark energy models in the context of the holographic principle. The potentials of the interacting fields are constructed. The models are compared with CMB distance information, baryonic acoustic oscillations, lookback time and the Constitution supernovae sample. For both models, the results are consistent with a nonvanishing interaction in the dark sector of the Universe and the sign of coupling is consistent with dark energy decaying into dark matter, alleviating the coincidence problem-with more than 3 standard deviations of confidence for one of them. However, this is because the noninteracting holographic dark energy model is a bad fit to the combination of data sets used in this work as compared to the cosmological constant with cold dark matter model, so that one needs to introduce the interaction in order to improve this model.

Low fatness, reduced fat intake and adequate plasmatic concentrations of LDL-cholesterol are associated with high bone mineral density in women: a cross-sectional study with control group

Background: Several parameters are associated with high bone mineral density (BMD), such as overweight, black background, intense physical activity (PA), greater calcium intake and some medications. The objectives are to evaluate the prevalence and the main aspects associated with high BMD in healthy women. Methods: After reviewing the database of approximately 21,500 BMD scans performed in the metropolitan area of Sao Paulo, Brazil, from June 2005 to October 2010, high BMD (over 1400 g/cm(2) at lumbar spine and/or above 1200 g/cm2 at femoral neck) was found in 421 exams. Exclusion criteria were age below 30 or above 60 years, black ethnicity, pregnant or obese women, disease and/or medications known to interfere with bone metabolism. A total of 40 women with high BMD were included and matched with 40 healthy women with normal BMD, paired to weight, age, skin color and menopausal status. Medical history, food intake and PA were assessed through validated questionnaires. Body composition was evaluated through a GE-Lunar DPX MD + bone densitometer. Radiography of the thoracic and lumbar spine was carried out to exclude degenerative alterations or fractures. Biochemical parameters included both lipid and hormonal profiles, along with mineral and bone metabolism. Statistical analysis included parametric and nonparametric tests and linear regression models. P < 0.05 was considered significant. Results: The mean age was 50.9 (8.3) years. There was no significant difference between groups in relation to PA, smoking, intake of calcium and vitamin D, as well as laboratory tests, except serum C-telopeptide of type I collagen (s-CTX), which was lower in the high BMD group (p = 0.04). In the final model of multivariate regression, a lower fat intake and body fatness as well a better profile of LDL-cholesterol predicted almost 35% of high BMD in women. (adjusted R2 = 0.347; p < 0.001). In addition, greater amounts of lean mass and higher IGF-1 serum concentrations played a protective role, regardless age and weight. Conclusion: Our results demonstrate the potential deleterious effect of lipid metabolism-related components, including fat intake and body fatness and worse lipid profile, on bone mass and metabolism in healthy women.

On Wald Residuals in Generalized Linear Models

A rigorous asymptotic theory for Wald residuals in generalized linear models is not yet available. The authors provide matrix formulae of order O(n(-1)), where n is the sample size, for the first two moments of these residuals. The formulae can be applied to many regression models widely used in practice. The authors suggest adjusted Wald residuals to these models with approximately zero mean and unit variance. The expressions were used to analyze a real dataset. Some simulation results indicate that the adjusted Wald residuals are better approximated by the standard normal distribution than the Wald residuals.

Validation and comparison of two soil-vegetation-atmosphere transfer models for tropical Africa

This study aims to compare and validate two soil-vegetation-atmosphere-transfer (SVAT) schemes: TERRA-ML and the Community Land Model (CLM). Both SVAT schemes are run in standalone mode (decoupled from an atmospheric model) and forced with meteorological in-situ measurements obtained at several tropical African sites. Model performance is quantified by comparing simulated sensible and latent heat fluxes with eddy-covariance measurements. Our analysis indicates that the Community Land Model corresponds more closely to the micrometeorological observations, reflecting the advantages of the higher model complexity and physical realism. Deficiencies in TERRA-ML are addressed and its performance is improved: (1) adjusting input data (root depth) to region-specific values (tropical evergreen forest) resolves dry-season underestimation of evapotranspiration; (2) adjusting the leaf area index and albedo (depending on hard-coded model constants) resolves overestimations of both latent and sensible heat fluxes; and (3) an unrealistic flux partitioning caused by overestimated superficial water contents is reduced by adjusting the hydraulic conductivity parameterization. CLM is by default more versatile in its global application on different vegetation types and climates. On the other hand, with its lower degree of complexity, TERRA-ML is much less computationally demanding, which leads to faster calculation times in a coupled climate simulation.

Percutaneous sciatic nerve block with tramadol induces analgesia and motor blockade in two animal pain models

Local anesthetic efficacy of tramadol has been reported following intradermal application. Our aim was to investigate the effect of perineural tramadol as the sole analgesic in two pain models. Male Wistar rats (280-380 g; N = 5/group) were used in these experiments. A neurostimulation-guided sciatic nerve block was performed and 2% lidocaine or tramadol (1.25 and 5 mg) was perineurally injected in two different animal pain models. In the flinching behavior test, the number of flinches was evaluated and in the plantar incision model, mechanical and heat thresholds were measured. Motor effects of lidocaine and tramadol were quantified and a motor block score elaborated. Tramadol, 1.25 mg, completely blocked the first and reduced the second phase of the flinching behavior test. In the plantar incision model, tramadol (1.25 mg) increased both paw withdrawal latency in response to radiant heat (8.3 +/- 1.1, 12.7 +/- 1.8, 8.4 +/- 0.8, and 11.1 +/- 3.3 s) and mechanical threshold in response to von Frey filaments (459 +/- 82.8, 447.5 +/- 91.7, 320.1 +/- 120, 126.43 +/- 92.8 mN) at 5, 15, 30, and 60 min, respectively. Sham block or contralateral sciatic nerve block did not differ from perineural saline injection throughout the study in either model. The effect of tramadol was not antagonized by intraperitoneal naloxone. High dose tramadol (5 mg) blocked motor function as well as 2% lidocaine. In conclusion, tramadol blocks nociception and motor function in vivo similar to local anesthetics.