Identifying the preferred subset of enzymatic profiles in nonlinear kinetic metabolic models via multiobjective global optimization and pareto filters

Optimization models in metabolic engineering and systems biology focus typically on optimizing a unique criterion, usually the synthesis rate of a metabolite of interest or the rate of growth. Connectivity and non-linear regulatory effects, however, make it necessary to consider multiple objectives in order to identify useful strategies that balance out different metabolic issues. This is a fundamental aspect, as optimization of maximum yield in a given condition may involve unrealistic values in other key processes. Due to the difficulties associated with detailed non-linear models, analysis using stoichiometric descriptions and linear optimization methods have become rather popular in systems biology. However, despite being useful, these approaches fail in capturing the intrinsic nonlinear nature of the underlying metabolic systems and the regulatory signals involved. Targeting more complex biological systems requires the application of global optimization methods to non-linear representations. In this work we address the multi-objective global optimization of metabolic networks that are described by a special class of models based on the power-law formalism: the generalized mass action (GMA) representation. Our goal is to develop global optimization methods capable of efficiently dealing with several biological criteria simultaneously. In order to overcome the numerical difficulties of dealing with multiple criteria in the optimization, we propose a heuristic approach based on the epsilon constraint method that reduces the computational burden of generating a set of Pareto optimal alternatives, each achieving a unique combination of objectives values. To facilitate the post-optimal analysis of these solutions and narrow down their number prior to being tested in the laboratory, we explore the use of Pareto filters that identify the preferred subset of enzymatic profiles. We demonstrate the usefulness of our approach by means of a case study that optimizes the ethanol production in the fermentation of Saccharomyces cerevisiae.

Electrophysiological measurements of spectral sensitivities: a review

Spectral sensitivities of visual systems are specified as the reciprocals of the intensities of light (quantum fluxes) needed at each wavelength to elicit the same criterion amplitude of responses. This review primarily considers the methods that have been developed for electrophysiological determinations of criterion amplitudes of slow-wave responses from single retinal cells. Traditional flash methods can require tedious dark adaptations and may yield erroneous spectral sensitivity curves which are not seen in such modifications as ramp methods. Linear response methods involve interferometry, while constant response methods involve manual or automatic adjustments of continuous illumination to keep response amplitudes constant during spectral scans. In DC or AC computerized constant response methods, feedback to determine intensities at each wavelength is derived from the response amplitudes themselves. Although all but traditional flash methods have greater or lesser abilities to provide on-line determinations of spectral sensitivities, computerized constant response methods are the most satisfactory due to flexibility, speed and maintenance of a constant adaptation level

Entropy application to improve construction finance decisions

In financial decision-making processes, the adopted weights of the objective functions have significant impacts on the final decision outcome. However, conventional rating and weighting methods exhibit difficulty in deriving appropriate weights for complex decision-making problems with imprecise information. Entropy is a quantitative measure of uncertainty and has been useful in exploring weights of attributes in decision making. A fuzzy and entropy-based mathematical approach is employed to solve the weighting problem of the objective functions in an overall cash-flow model. The multiproject being undertaken by a medium-size construction firm in Hong Kong was used as a real case study to demonstrate the application of entropy. Its application in multiproject cash flow situations is demonstrated. The results indicate that the overall before-tax profit was HK$ 0.11 millions lower after the introduction of appropriate weights. In addition, the best time to invest in new projects arising from positive cash flow was identified to be two working months earlier than the nonweight system.

Microbial−mammalian cometabolites dominate the age-associated urinary metabolic phenotype in Taiwanese and American populations

Understanding the metabolic processes associated with aging is key to developing effective management and treatment strategies for age-related diseases. We investigated the metabolic profiles associated with age in a Taiwanese and an American population. 1H NMR spectral profiles were generated for urine specimens collected from the Taiwanese Social Environment and Biomarkers of Aging Study (SEBAS; n = 857; age 54–91 years) and the Mid-Life in the USA study (MIDUS II; n = 1148; age 35–86 years). Multivariate and univariate linear projection methods revealed some common age-related characteristics in urinary metabolite profiles in the American and Taiwanese populations, as well as some distinctive features. In both cases, two metabolites—4-cresyl sulfate (4CS) and phenylacetylglutamine (PAG)—were positively associated with age. In addition, creatine and β-hydroxy-β-methylbutyrate (HMB) were negatively correlated with age in both populations (p < 4 × 10–6). These age-associated gradients in creatine and HMB reflect decreasing muscle mass with age. The systematic increase in PAG and 4CS was confirmed using ultraperformance liquid chromatography–mass spectrometry (UPLC–MS). Both are products of concerted microbial–mammalian host cometabolism and indicate an age-related association with the balance of host–microbiome metabolism.

Measurements of stratospheric ozone at a mid-latitude observing station Valentia, Ireland (51.94° N, 10.25° W), using ground-based and ozonesonde observations from 1994 to 2009

Sixteen years (1994 – 2009) of ozone profiling by ozonesondes at Valentia Meteorological and Geophysical Observatory, Ireland (51.94° N, 10.23° W) along with a co-located MkIV Brewer spectrophotometer for the period 1993–2009 are analyzed. Simple and multiple linear regression methods are used to infer the recent trend, if any, in stratospheric column ozone over the station. The decadal trend from 1994 to 2010 is also calculated from the monthly mean data of Brewer and column ozone data derived from satellite observations. Both of these show a 1.5 % increase per decade during this period with an uncertainty of about ±0.25 %. Monthly mean data for March show a much stronger trend of ~ 4.8 % increase per decade for both ozonesonde and Brewer data. The ozone profile is divided between three vertical slots of 0–15 km, 15–26 km, and 26 km to the top of the atmosphere and a 11-year running average is calculated. Ozone values for the month of March only are observed to increase at each level with a maximum change of +9.2 ± 3.2 % per decade (between years 1994 and 2009) being observed in the vertical region from 15 to 26 km. In the tropospheric region from 0 to 15 km, the trend is positive but with a poor statistical significance. However, for the top level of above 26 km the trend is significantly positive at about 4 % per decade. The March integrated ozonesonde column ozone during this period is found to increase at a rate of ~6.6 % per decade compared with the Brewer and satellite positive trends of ~5 % per decade.

Detection of the Vandibular Canal and the Mental Foramen in Panoramic Radiographs: Intraexaminer Agreement

The aim of this study was to evaluate the intraexaminer agreement in the detection of the mandibular canal roof (MCR) and mental foramen (MF) in panoramic radiographs. Forty panoramic radiographs of edentulous patients were used. Two calibrated examiners (A and B) read the images 2 times, for both sides independently, under blind conditions. The interval between the readings was 10 days. The intraexaminer agreement in the interpretation of MCR and MF was performed by kappa statistics with linear weighting (x). The intraexaminer agreement for the detection of MCR, in the left side, was good for both examiners (A: kappa = 0.67; B: kappa = 0.71). Related to the right side, it was found to be kappa = 0.47 and kappa = 0.62, respectively to A and B. The intraexaminer agreement for the detection of MF was good for both examiners interpreting the left side (A: kappa = 0.61; B: kappa = 0.63), and in relation to the right side, it was moderate (A: kappa = 0.51) and fair (B: kappa = 0.38). The intraexaminer agreement in the detection of MCR was good and from good to fair in the detection of MF.

Efeitos agudos do tabagismo sobre a modulação autonômica: análise por meio do plot de poincaré

FUNDAMENTO: O tabagismo altera a função autonômica. OBJETIVO: Investigar os efeitos agudos do tabagismo sobre a modulação autonômica e a recuperação dos índices de variabilidade de frequência cardíaca (VFC) pós-fumo, por meio do plot de Poincaré e índices lineares. MÉTODOS: Foram avaliados 25 fumantes jovens, os quais tiveram a frequência cardíaca analisada, batimento a batimento, na posição sentada, após 8 horas de abstinência, por 30 minutos em repouso, 20 minutos durante o fumo e 30 minutos pós-fumo. Análise de variância para medidas repetidas, seguido do teste de Tukey, ou teste de Friedman seguido do teste de Dunn foram aplicados dependendo da normalidade dos dados, com p < 0,05. RESULTADOS: Durante o fumo, houve redução dos índices SD1 (23,4 ± 9,2 vs 13,8 ± 4,8), razão SD1/SD2 (0,31 ± 0,08 vs 0,2 ± 0,04), RMSSD (32,7 ± 13 vs 19,1 ± 6,8), SDNN (47,6 ± 14,8 vs 35,5 ± 8,4), HFnu (32,5 ± 11,6 vs 19 ± 8,1) e do intervalo RR (816,8 ± 89 vs 696,5 ± 76,3) em relação ao repouso, enquanto que aumentos do índice LFnu (67,5 ± 11,6 vs 81 ± 8,1) e da razão LF/HF (2,6 ± 1,7 vs 5,4 ± 3,1) foram observados. A análise visual do plot mostrou menor dispersão dos intervalos RR durante o fumo. Com exceção da razão SD1/SD2, os demais índices apresentaram recuperação dos valores, 30 minutos após o tabagismo. CONCLUSÃO: O tabagismo produziu agudamente modificações no controle autonômico, caracterizadas por ativação simpática e retirada vagal, com recuperação 30 minutos após o fumo.

Atitudes em relação à Bioestatística de discentes e docentes da Faculdade de Ciências Farmacêuticas de Araraquara-UNESP

In this study, the attitudes of students and teachers at the Faculty of Pharmaceutical Sciences of Araraquara (UNESP University) towards Biostatistics were assessed. The Survey of Attitudes Toward Statistics (SATS) scale was used as the measuring instrument. The reproducibility of the scale was estimated by Cohen's Kappa (κ) coefficient with linear weighting and its internal consistency by Cronbach's alpha coefficient (α). The individuals were first placed in two groups, according to their positive or negative attitude toward Statistics; then, the association of their attitude with the variables of interest was tested by the chi-squared (χ 2) test at a significance level of 5%. The sample consisted of 272 undergraduate students, 83 graduate students and 24 teachers, predominantly female (78.2%). Among the students, 67.5% participated in the scientific research initial training program. Reproducibility and internal consistency of the scale were adequate (κ=0.7093; α=0.9334). Most of the subjects (74.4%) had a positive attitude toward Statistics. Significant association was found between attitude and functional activity (p=0.0204), the course taken (p=0.0316) and effort (p=0.0002). Thus, it was concluded that the great majority of the participants had a positive attitude towards Biostatistics and that undergraduate students and those who reported good performance in Biostatistics showed a significantly higher proportion of positive attitude than the other students.

Métodos numéricos para o retoque digital

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Inversão gravimétrica rápida do relevo do embasamento aplicando o funcional da variação total

Apresentamos um novo método de inversão linear bidimensional de dados gravimétricos produzidos por bacias sedimentares com relevo do embasamento descontínuo. O método desenvolvido utiliza um modelo interpretativo formado por um conjunto de fitas horizontais bidimensionais justapostas cujas espessuras são os parâmetros a serem estimados. O contraste de densidade entre o embasamento e os sedimentos é presumido constante e conhecido. As estimativas das espessuras foram estabilizadas com o funcional da Variação Total (VT) que permite soluções apresentando descontinuidades locais no relevo do embasamento. As estimativas do relevo são obtidas através da resolução de um sistema de equações lineares, resolvido na norma L₁. Como métodos lineares subestimam as estimativas de profundidade do embasamento de bacias maiores que cerca de 500 m, amplificamos as estimativas de profundidade através da modificação da matriz associada ao modelo interpretativo de fitas. As estimativas obtidas através deste procedimento são em geral ligeiramente superestimadas. Desse modo, elas são corrigidas através de uma correção definida pela expressão da placa Bouguer. Testes em dados sintéticos e reais produziram resultados comparáveis aos produzidos pelo método não linear, mas exigiram menor tempo computacional. A razão R entre os tempos exigidos pelo método não linear e o método proposto cresce com o número de observações e parâmetros. Por exemplo, para 60 observações e 60 parâmetros, R é igual a 4, enquanto para 2500 observações e 2500 parâmetros R cresce para 16,8. O método proposto e o método de inversão não linear foram aplicados também em dados reais do Steptoe Valley, Nevada, Estados Unidos, e da ponte do POEMA, no Campus do Guamá em Belém, produzindo soluções similares às obtidas com o método não linear exigindo menor tempo computacional.

Influência da miliamperagem sobre a qualidade do diagnóstico das fraturas radiculares por meio de tomografia computadorizada de feixe cônico: estudo in vitro

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Influência da miliamperagem sobre a qualidade do diagnóstico das fraturas radiculares por meio de tomografia computadorizada de feixe cônico: estudo in vitro

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Finite element models of coseismic deformation due to the 2009 L'Aquila (Italy) and 2008 Wenchuan(China) earthquakes

The topic of my Ph.D. thesis is the finite element modeling of coseismic deformation imaged by DInSAR and GPS data. I developed a method to calculate synthetic Green functions with finite element models (FEMs) and then use linear inversion methods to determine the slip distribution on the fault plane. The method is applied to the 2009 L’Aquila Earthquake (Italy) and to the 2008 Wenchuan earthquake (China). I focus on the influence of rheological features of the earth's crust by implementing seismic tomographic data and the influence of topography by implementing Digital Elevation Models (DEM) layers on the FEMs. Results for the L’Aquila earthquake highlight the non-negligible influence of the medium structure: homogeneous and heterogeneous models show discrepancies up to 20% in the fault slip distribution values. Furthermore, in the heterogeneous models a new area of slip appears above the hypocenter. Regarding the 2008 Wenchuan earthquake, the very steep topographic relief of Longmen Shan Range is implemented in my FE model. A large number of DEM layers corresponding to East China is used to achieve the complete coverage of the FE model. My objective was to explore the influence of the topography on the retrieved coseismic slip distribution. The inversion results reveals significant differences between the flat and topographic model. Thus, the flat models frequently adopted are inappropriate to represent the earth surface topographic features and especially in the case of the 2008 Wenchuan earthquake.

Inversion of low frequency seismo-volcanic events

The inversion of seismo-volcanic events is performed to retrieve the source geometry and to determine volumetric budgets of the source. Such observations have shown to be an important tool for the seismological monitoring of volcanoes. We developed a novel technique for the non-linear constrained inversion of low frequency seismo-volcanic events. Unconstrained linear inversion methods work well when a dense network of broadband seismometers is available. We propose a new constrained inversion technique, which has shown to be efficient also in a reduced network configuration and a low signal-noise ratio. The waveform inversion is performed in the frequency domain, constraining the source mechanism during the event to vary only in its magnitude. The eigenvectors orientation and the eigenvalue ratio are kept constant. This significantly reduces the number of parameters to invert, making the procedure more stable. The method has been tested over a synthetic dataset, reproducing realistic very-long-period (VLP) signals of Stromboli volcano. The information obtained by performing the synthetic tests is used to assess the reliability of the results obtained on a VLP dataset recorded on Stromboli volcano and on a low frequency events recorded at Vesuvius volcano.

Improved Localization of Implanted Subdural Electrode Contacts on MRI Using an Elastic Image Fusion Algorithm in Invasive EEG Recording

BACKGROUND: Accurate projection of implanted subdural electrode contacts in presurgical evaluation of pharmacoresistant epilepsy cases by invasive EEG is highly relevant. Linear fusion of CT and MRI images may display the contacts in the wrong position due to brain shift effects. OBJECTIVE: A retrospective study in five patients with pharmacoresistant epilepsy was performed to evaluate whether an elastic image fusion algorithm can provide a more accurate projection of the electrode contacts on the pre-implantation MRI as compared to linear fusion. METHODS: An automated elastic image fusion algorithm (AEF), a guided elastic image fusion algorithm (GEF), and a standard linear fusion algorithm (LF) were used on preoperative MRI and post-implantation CT scans. Vertical correction of virtual contact positions, total virtual contact shift, corrections of midline shift and brain shifts due to pneumencephalus were measured. RESULTS: Both AEF and GEF worked well with all 5 cases. An average midline shift of 1.7mm (SD 1.25) was corrected to 0.4mm (SD 0.8) after AEF and to 0.0mm (SD 0) after GEF. Median virtual distances between contacts and cortical surface were corrected by a significant amount, from 2.3mm after LF to 0.0mm after AEF and GEF (p<.001). Mean total relative corrections of 3.1 mm (SD 1.85) after AEF and 3.0mm (SD 1.77) after GEF were achieved. The tested version of GEF did not achieve a satisfying virtual correction of pneumencephalus. CONCLUSION: The technique provided a clear improvement in fusion of pre- and post-implantation scans, although the accuracy is difficult to evaluate.