Algorithms and tools for in silico design of cell factories

PhD thesis in Bioengineering

Testing Nelder-Mead based repulsion algorithms for multiple roots of nonlinear systems via a two-level factorial design of experiments

This paper addresses the challenging task of computing multiple roots of a system of nonlinear equations. A repulsion algorithm that invokes the Nelder-Mead (N-M) local search method and uses a penalty-type merit function based on the error function, known as 'erf', is presented. In the N-M algorithm context, different strategies are proposed to enhance the quality of the solutions and improve the overall efficiency. The main goal of this paper is to use a two-level factorial design of experiments to analyze the statistical significance of the observed differences in selected performance criteria produced when testing different strategies in the N-M based repulsion algorithm. The main goal of this paper is to use a two-level factorial design of experiments to analyze the statistical significance of the observed differences in selected performance criteria produced when testing different strategies in the N-M based repulsion algorithm.

Spatial-temporal modellization of the NO2 concentration data through geostatistical tools

The nitrogen dioxide is a primary pollutant, regarded for the estimation of the air quality index, whose excessive presence may cause significant environmental and health problems. In the current work, we suggest characterizing the evolution of NO2 levels, by using geostatisti- cal approaches that deal with both the space and time coordinates. To develop our proposal, a first exploratory analysis was carried out on daily values of the target variable, daily measured in Portugal from 2004 to 2012, which led to identify three influential covariates (type of site, environment and month of measurement). In a second step, appropriate geostatistical tools were applied to model the trend and the space-time variability, thus enabling us to use the kriging techniques for prediction, without requiring data from a dense monitoring network. This method- ology has valuable applications, as it can provide accurate assessment of the nitrogen dioxide concentrations at sites where either data have been lost or there is no monitoring station nearby.

Bootstrap performance profiles in stochastic algorithms assessment

Optimization with stochastic algorithms has become a relevant research field. Due to its stochastic nature, its assessment is not straightforward and involves integrating accuracy and precision. Performance profiles for the mean do not show the trade-off between accuracy and precision, and parametric stochastic profiles require strong distributional assumptions and are limited to the mean performance for a large number of runs. In this work, bootstrap performance profiles are used to compare stochastic algorithms for different statistics. This technique allows the estimation of the sampling distribution of almost any statistic even with small samples. Multiple comparison profiles are presented for more than two algorithms. The advantages and drawbacks of each assessment methodology are discussed.

Flood pulse and spatial dynamics of mercury in sediments in Puruzinho lake, Brazilian Amazon

The study was conducted in Puruzinho lake (Humaitá, AM) considering seasonal periods of rainy and dry in way to elucidate the flood pulse importance in the deposition, remobilization and distributions of mercury and organic matter in bottom sediments in the Madeira River Basin (Brazilian Amazon). Bottom sediments and soils samples were analyzed for total mercury and organic matter. Mercury concentrations obtained in bottom sediment were 32.20-146.40 ng g-1 and organic matter values were 3.5 - 18.0%. The main region for accumulation of mercury and organic matter was in the central and deepest lake area In the rainy season there was a greater distribution of Hg and organic matter, mainly controlled by means of income of the Madeira river water during flooding, while the predominant process in the dry season was the remobilization of total Hg due to the resuspension of bottom sediments.

Seasonal and spatial variations in limnological conditions of a floodplain lake (lake Catalão) connected to both the Solimões and Negro Rivers, Central Amazonia

Lakes play an important role in biogeochemical, ecological and hydrological processes in the river-floodplain system. The aim of this study was to evaluate the dynamics of the limnological conditions of Catalão Lake, an Amazon floodplain lake. Thus, some of the main limnological environment variables (O2, temperature, pH, nutrient, electrical conductivity) of the Catalão Lake were analyzed under temporal and spacial scales. The study was conducted between November/2004 and August/2005. Sampling excursion were carried out every three months; one excursion for each of the four different hydrological periods (low water, rising water, high water and falling water). Sampling points were chosen so that it could be obtained a gradient of the distance from Negro River. Limnological profiles in Catalão Lake showed generally acidic to slightly alcaline water, with low levels of dissolved oxygen and low concentrations of soluble reactive phosphorous. The Negro River seems to exert the main influence during the rising water period, while the Solimões River is the principal controlling river during peak water. The Principal Component Analysis (PCA) grouped the seasonal collections by hydrological period, showing the formation of a north-south spatial gradient within the lake in relation to the limnological variables. Multivariate dispersion analysis based on distance-to-centroid method demonstrated an increase in similarity over the course of the hydrological cycle, as the lake was inundated in response to the flood pulse of the main river channels. However, the largest spatial homogeneity in the lake was observed in the epilimnion layer, during the falling water period. The daily analysis of variation indicated an oligomitic pattern during the years in which the lake was permanently connected to the Negro River. Although Catalão Lake receives large quantities of both black water from the Negro River and sediment-filled water from the Solimões River, the physical and chemical characteristics of the lake are more similar to those of the Solimões (várzea lake) than the Negro (blackwater lake).

Spatial-temporal variation of dissolved inorganic material in the Amazon basin

The Amazon River basin is important in the contribution of dissolved material to oceans (4% worldwide). The aim of this work was to study the spatial and the temporal variability of dissolved inorganic materials in the main rivers of the Amazon basin. Data from 2003 to 2011 from six gauging stations of the ORE-HYBAM localized in Solimões, Purus, Madeira and Amazon rivers were used for this study. The concentrations of Ca2+, Na+, K+, Mg2+, Cl-, SO4 -2, HCO3 - and SiO2 were analyzed. At the stations of Solimões and Amazon rivers, the concentrations of Ca2+, Mg2+, HCO3 - and SO4 -2 had heterogeneous distribution over the years and did not show seasonality. At the stations of Madeira river, the concentration of these ions had seasonality inversely proportional to water discharge (dilution-concentration effect). Similar behavior was observed for the concentrations of Cl- and Na+ at the stations of the Solimões, Amazon and Madeira rivers, indicating almost constant release of Cl- and Na+ fluxes during the hydrological cycle. K+ and SiO2 showed almost constant concentrations throughout the years and all the stations, indicating that their flows depend on the river discharge variation. Therefore, the temporal variability of the dissolved inorganic material fluxes in the Solimões and Amazon rivers depends on the hydro-climatic factor and on the heterogeneity of the sources. In the Madeira and Purus rivers there is less influence of these factors, indicating that dissolved load fluxes are mainly associated to silicates weathering. As the Solimões basin contributes approximately 84% of the total flux of dissolved materials in the basin and is mainly under the influence of a hydro-climatic factor, we conclude that the temporal variability of this factor controls the temporal variability of the dissolved material fluxes of the Amazon basin.

Spatial and temporal evolution of family-farming land use in the Tapajós region of the Brazilian Amazon

Pressures on the Brazilian Amazon forest have been accentuated by agricultural activities practiced by families encouraged to settle in this region in the 1970s by the colonization program of the government. The aims of this study were to analyze the temporal and spatial evolution of land cover and land use (LCLU) in the lower Tapajós region, in the state of Pará. We contrast 11 watersheds that are generally representative of the colonization dynamics in the region. For this purpose, Landsat satellite images from three different years, 1986, 2001, and 2009, were analyzed with Geographic Information Systems. Individual images were subject to an unsupervised classification using the Maximum Likelihood Classification algorithm available on GRASS. The classes retained for the representation of LCLU in this study were: (1) slightly altered old-growth forest, (2) succession forest, (3) crop land and pasture, and (4) bare soil. The analysis and observation of general trends in eleven watersheds shows that LCLU is changing very rapidly. The average deforestation of old-growth forest in all the watersheds was estimated at more than 30% for the period of 1986 to 2009. The local-scale analysis of watersheds reveals the complexity of LCLU, notably in relation to large changes in the temporal and spatial evolution of watersheds. Proximity to the sprawling city of Itaituba is related to the highest rate of deforestation in two watersheds. The opening of roads such as the Transamazonian highway is associated to the second highest rate of deforestation in three watersheds.

Evaluation of geostatistical techniques to estimate the spatial distribution of aboveground biomass in the Amazon rainforest using high-resolution remote sensing data

ABSTRACT The spatial distribution of forest biomass in the Amazon is heterogeneous with a temporal and spatial variation, especially in relation to the different vegetation types of this biome. Biomass estimated in this region varies significantly depending on the applied approach and the data set used for modeling it. In this context, this study aimed to evaluate three different geostatistical techniques to estimate the spatial distribution of aboveground biomass (AGB). The selected techniques were: 1) ordinary least-squares regression (OLS), 2) geographically weighted regression (GWR) and, 3) geographically weighted regression - kriging (GWR-K). These techniques were applied to the same field dataset, using the same environmental variables derived from cartographic information and high-resolution remote sensing data (RapidEye). This study was developed in the Amazon rainforest from Sucumbíos - Ecuador. The results of this study showed that the GWR-K, a hybrid technique, provided statistically satisfactory estimates with the lowest prediction error compared to the other two techniques. Furthermore, we observed that 75% of the AGB was explained by the combination of remote sensing data and environmental variables, where the forest types are the most important variable for estimating AGB. It should be noted that while the use of high-resolution images significantly improves the estimation of the spatial distribution of AGB, the processing of this information requires high computational demand.

Light responsive multilayer surfaces with controlled spatial extinction capability

Multilayer systems obtained using the Layer-by-Layer (LbL) technology have been proposed for a variety of biomedical applications in tissue engineering and regenerative medicine. LbL assembly is a simple and highly versatile method to modify surfaces and fabricate robust and highly-ordered nanostructured coatings over almost any type of substrates and with a wide range of substances. The incorporation of polyoxometalate (POM) inorganic salts as constituents of the layers presents a possibility of promoting light-stimuli responses in LbL substrates. We propose the design of a biocompatible photo-responsive multilayer system based on a Preyssler-type POM ([NaP5W30O110]14â ) and a natural origin polymer, chitosan, using the LbL methodology. The photo-reduction properties of the POM allow the spatially controlled disruption of the assembled layers due to the weakening of the electrostatic interactions between the layers. This system has found applicability in detaching devices, such as the cell sheet technology, which may solve the drawbacks actually found in other cell treatment proposals.

Deteção remota na avaliação das alterações da cobertura do solo no Parque Nacional da Peneda-Gerês

Dissertação de mestrado integrado em Engenharia Civil

Algorithms and computational tools for metabolic flux analysis

PhD thesis in Biomedical Engineering

A survey of distributed data aggregation algorithms

Distributed data aggregation is an important task, allowing the de- centralized determination of meaningful global properties, that can then be used to direct the execution of other applications. The resulting val- ues result from the distributed computation of functions like count, sum and average. Some application examples can found to determine the network size, total storage capacity, average load, majorities and many others. In the last decade, many di erent approaches have been pro- posed, with di erent trade-o s in terms of accuracy, reliability, message and time complexity. Due to the considerable amount and variety of ag- gregation algorithms, it can be di cult and time consuming to determine which techniques will be more appropriate to use in speci c settings, jus- tifying the existence of a survey to aid in this task. This work reviews the state of the art on distributed data aggregation algorithms, providing three main contributions. First, it formally de nes the concept of aggrega- tion, characterizing the di erent types of aggregation functions. Second, it succinctly describes the main aggregation techniques, organizing them in a taxonomy. Finally, it provides some guidelines toward the selection and use of the most relevant techniques, summarizing their principal characteristics.

Hypertension and clustering of cardiovascular risk factors in a community in Southeast Brazil: the Bambuí Health and Ageing Study

OBJECTIVE - A population-based prospective study was analysed to: a) determine the prevalence of hypertension; b) investigate the clustering of other cardiovascular risk factors and c) verify whether older differed from younger adults in the pattern of clustering. METHODS - The data comprised a representative sample of the population of Bambuí, Brazil. Multiple logistic regression was used to investigate the independent association between hypertension and selected factors. RESULTS - A total of 820 younger adults (82.5%) and 1494 older adults (85.9%) participated in this study. The overall prevalence of hypertension was 24.8% (SE=1.4 %), being higher in women (26.9±1.5%) than in men (22.0± 1.7%) (p=0.033). Hypertension was positively and significantly associated with physical inactivity, overweight, hypercholesterolemia hyperglycemia and hypertriglyceridemia. The coexistence of hypertension with 4 or more of these risk factors occurred 6 times more than expected by chance, after adjusting for age and sex (OR=6.3; 95%CI: 3.4-11.9). The pattern of risk factor clustering in hypertensive individuals differed with age. CONCLUSION - Our results reinforce the need to increase detection and treatment of hypertension and to approach patients' global risk profiles.

Desarrollo de Algoritmos Numéricos para el Estudio de la Propagación no Lineal de Ondas en Aplicaciones Aeroespaciales y Astrofísicas. Develop to Numerical Algorithms to Study non-Linear Waves Propagation in Aerospace and Astrophysics Applications.

En este proyecto se desarrollarán algoritmos numéricos para sistemas no lineales hiperbólicos-parabólicos de ecuaciones diferenciales en derivadas parciales. Dichos sistemas tienen aplicación en propagación de ondas en ámbitos aeroespaciales y astrofísicos.Objetivos generales: 1)Desarrollo y mejora de algoritmos numéricos con la finalidad de incrementar la calidad en la simulación de propagación e interacción de ondas gasdinámicas y magnetogasdinámicas no lineales. 2)Desarrollo de códigos computacionales con la finalidad de simular flujos gasdinámicos de elevada entalpía incluyendo cambios químicos, efectos dispersivos y difusivos.3)Desarrollo de códigos computacionales con la finalidad de simular flujos magnetogasdinámicos ideales y reales.4)Aplicación de los nuevos algoritmos y códigos computacionales a la solución del flujo aerotermodinámico alrededor de cuerpos que ingresan en la atmósfera terrestre. 5)Aplicación de los nuevos algoritmos y códigos computacionales a la simulación del comportamiento dinámico no lineal de arcos magnéticos en la corona solar. 6)Desarrollo de nuevos modelos para describir el comportamiento no lineal de arcos magnéticos en la corona solar.Este proyecto presenta como objetivo principal la introducción de mejoras en algoritmos numéricos para simular la propagación e interacción de ondas no lineales en dos medios gaseosos: aquellos que no poseen carga eléctrica libre (flujos gasdinámicos) y aquellos que tienen carga eléctrica libre (flujos magnetogasdinámicos). Al mismo tiempo se desarrollarán códigos computacionales que implementen las mejoras de las técnicas numéricas.Los algoritmos numéricos se aplicarán con la finalidad de incrementar el conocimiento en tópicos de interés en la ingeniería aeroespacial como es el cálculo del flujo de calor y fuerzas aerotermodinámicas que soportan objetos que ingresan a la atmósfera terrestre y en temas de astrofísica como la propagación e interacción de ondas, tanto para la transferencia de energía como para la generación de inestabilidades en arcos magnéticos de la corona solar. Estos dos temas poseen en común las técnicas y algoritmos numéricos con los que serán tratados. Las ecuaciones gasdinámicas y magnetogasdinámicas ideales conforman sistemas hiperbólicos de ecuaciones diferenciales y pueden ser solucionados utilizando "Riemann solvers" junto con el método de volúmenes finitos (Toro 1999; Udrea 1999; LeVeque 1992 y 2005). La inclusión de efectos difusivos genera que los sistemas de ecuaciones resulten hiperbólicos-parabólicos. La contribución parabólica puede ser considerada como términos fuentes y tratada adicionalmente tanto en forma explícita como implícita (Udrea 1999; LeVeque 2005).Para analizar el flujo alrededor de cuerpos que ingresan en la atmósfera se utilizarán las ecuaciones de Navier-Stokes químicamente activas, mientras la temperatura no supere los 6000K. Para mayores temperaturas es necesario considerar efectos de ionización (Anderson, 1989). Tanto los efectos difusivos como los cambios químicos serán considerados como términos fuentes en las ecuaciones de Euler. Para tratar la propagación de ondas, transferencia de energía e inestabilidades en arcos magnéticos de la corona solar se utilizarán las ecuaciones de la magnetogasdinámica ideal y real. En este caso será también conveniente implementar términos fuente para el tratamiento de fenómenos de transporte como el flujo de calor y el de radiación. Los códigos utilizarán la técnica de volúmenes finitos, junto con esquemas "Total Variation Disminishing - TVD" sobre mallas estructuradas y no estructuradas.