Partial automation of database processing of simulation outputs from L-systems models of plant morphogenesis

Models of plant architecture allow us to explore how genotype environment interactions effect the development of plant phenotypes. Such models generate masses of data organised in complex hierarchies. This paper presents a generic system for creating and automatically populating a relational database from data generated by the widely used L-system approach to modelling plant morphogenesis. Techniques from compiler technology are applied to generate attributes (new fields) in the database, to simplify query development for the recursively-structured branching relationship. Use of biological terminology in an interactive query builder contributes towards making the system biologist-friendly. (C) 2002 Elsevier Science Ireland Ltd. All rights reserved.

An approach to verifying and debugging simulation models governed by ordinary differential equations: Part 2. Residuals analysis and a case study

In Part 1 of this paper a methodology for back-to-back testing of simulation software was described. Residuals with error-dependent geometric properties were generated. A set of potential coding errors was enumerated, along with a corresponding set of feature matrices, which describe the geometric properties imposed on the residuals by each of the errors. In this part of the paper, an algorithm is developed to isolate the coding errors present by analysing the residuals. A set of errors is isolated when the subspace spanned by their combined feature matrices corresponds to that of the residuals. Individual feature matrices are compared to the residuals and classified as 'definite', 'possible' or 'impossible'. The status of 'possible' errors is resolved using a dynamic subset testing algorithm. To demonstrate and validate the testing methodology presented in Part 1 and the isolation algorithm presented in Part 2, a case study is presented using a model for biological wastewater treatment. Both single and simultaneous errors that are deliberately introduced into the simulation code are correctly detected and isolated. Copyright (C) 2003 John Wiley Sons, Ltd.

Cobalt complexes with pyrazole ligands as catalyst precursors for the peroxidative oxidation of cyclohexane: Z-ray absorption spectroscopy studies and biological applications

[CoCl(-Cl)(Hpz(Ph))(3)](2) (1) and [CoCl2(Hpz(Ph))(4)] (2) were obtained by reaction of CoCl2 with HC(pz(Ph))(3) and Hpz(Ph), respectively (Hpz(Ph)=3-phenylpyrazole). The compounds were isolated as air-stable solids and fully characterized by IR and far-IR spectroscopy, MS(ESI+/-), elemental analysis, cyclic voltammetry (CV), controlled potential electrolysis, and single-crystal X-ray diffraction. Electrochemical studies showed that 1 and 2 undergo single-electron irreversible (CoCoIII)-Co-II oxidations and (CoCoI)-Co-II reductions at potentials measured by CV, which also allowed, in the case of dinuclear complex 1, the detection of electronic communication between the Co centers through the chloride bridging ligands. The electrochemical behavior of models of 1 and 2 were also investigated by density functional theory (DFT) methods, which indicated that the vertical oxidation of 1 and 2 (that before structural relaxation) affects mostly the chloride and pyrazolyl ligands, whereas adiabatic oxidation (that after the geometry relaxation) and reduction are mostly metal centered. Compounds 1 and 2 and, for comparative purposes, other related scorpionate and pyrazole cobalt complexes, exhibit catalytic activity for the peroxidative oxidation of cyclohexane to cyclohexanol and cyclohexanone under mild conditions (room temperature, aqueous H2O2). Insitu X-ray absorption spectroscopy studies indicated that the species derived from complexes 1 and 2 during the oxidation of cyclohexane (i.e., Ox-1 and Ox-2, respectively) are analogous and contain a Co-III site. Complex 2 showed low invitro cytotoxicity toward the HCT116 colorectal carcinoma and MCF7 breast adenocarcinoma cell lines.

Synchronization in Von Bertalanffy’s models

Many data have been useful to describe the growth of marine mammals, invertebrates and reptiles, seabirds, sea turtles and fishes, using the logistic, the Gom-pertz and von Bertalanffy's growth models. A generalized family of von Bertalanffy's maps, which is proportional to the right hand side of von Bertalanffy's growth equation, is studied and its dynamical approach is proposed. The system complexity is measured using Lyapunov exponents, which depend on two biological parameters: von Bertalanffy's growth rate constant and the asymptotic weight. Applications of synchronization in real world is of current interest. The behavior of birds ocks, schools of fish and other animals is an important phenomenon characterized by synchronized motion of individuals. In this work, we consider networks having in each node a von Bertalanffy's model and we study the synchronization interval of these networks, as a function of those two biological parameters. Numerical simulation are also presented to support our approaches.

From weak Allee effect to no Allee effect in Richards’ growth models

Population dynamics have been attracting interest since many years. Among the considered models, the Richards’ equations remain one of the most popular to describe biological growth processes. On the other hand, Allee effect is currently a major focus of ecological research, which occurs when positive density dependence dominates at low densities. In this chapter, we propose the dynamical study of classes of functions based on Richards’ models describing the existence or not of Allee effect. We investigate bifurcation structures in generalized Richards’ functions and we look for the conditions in the (β, r) parameter plane for the existence of a weak Allee effect region. We show that the existence of this region is related with the existence of a dovetail structure. When the Allee limit varies, the weak Allee effect region disappears when the dovetail structure also disappears. Consequently, we deduce the transition from the weak Allee effect to no Allee effect to this family of functions. To support our analysis, we present fold and flip bifurcation curves and numerical simulations of several bifurcation diagrams.

Models for pheromone evaluation in Ant Systems for Mobile Ad-hoc networks

On a mobile ad-hoc network environment, where the resources are scarce, the knowledge about the network's link state is essential to optimize the routing procedures. This paper presents a study about different pheromone evaluation models and how they react to possible changes in traffic rate. Observing how the pheromone value on a link changes, it could be possible to identify certain patterns which can indicate the path status. For this study, the behavior of the Ant System evaluation model was compared with a Temporal Active Pheromone model (a biological approach) and a Progressive Pheromone Reduction model with and without a maximum pheromone limit.

Participatory Surveillance and Mathematical Models in Epidemiologic Research: Successes and Challenges

Theoretical epidemiology aims to understand the dynamics of diseases in populations and communities. Biological and behavioral processes are abstracted into mathematical formulations which aim to reproduce epidemiological observations. In this thesis a new system for the self-reporting of syndromic data — Influenzanet — is introduced and assessed. The system is currently being extended to address greater challenges of monitoring the health and well-being of tropical communities.(...)

The effect of key process operational conditions on enhanced biological phosphorus removal from wastewater

Enhanced biological phosphorus removal (EBPR) is the most economic and sustainable option used in wastewater treatment plants (WWTPs) for phosphorus removal. In this process it is important to control the competition between polyphosphate accumulating organisms (PAOs) and glycogen accumulating organisms (GAOs), since EBPR deterioration or failure can be related with the proliferation of GAOs over PAOs. This thesis is focused on the effect of operational conditions (volatile fatty acid (VFA) composition, dissolved oxygen (DO) concentration and organic carbon loading) on PAO and GAO metabolism. The knowledge about the effect of these operational conditions on EBPR metabolism is very important, since they represent key factors that impact WWTPs performance and sustainability. Substrate competition between the anaerobic uptake of acetate and propionate (the main VFAs present in WWTPs) was shown in this work to be a relevant factor affecting PAO metabolism, and a metabolic model was developed that successfully describes this effect. Interestingly, the aerobic metabolism of PAOs was not affected by different VFA compositions, since the aerobic kinetic parameters for phosphorus uptake, polyhydroxyalkanoates (PHAs) degradation and glycogen production were relatively independent of acetate or propionate concentration. This is very relevant for WWTPs, since it will simplify the calibration procedure for metabolic models, facilitating their use for full-scale systems. The DO concentration and aerobic hydraulic retention time (HRT) affected the PAO-GAO competition, where low DO levels or lower aerobic HRT was more favourable for PAOs than GAOs. Indeed, the oxygen affinity coefficient was significantly higher for GAOs than PAOs, showing that PAOs were far superior at scavenging for the often limited oxygen levels in WWTPs. The operation of WWTPs with low aeration is of high importance for full-scale systems, since it decreases the energetic costs and can potentially improve WWTP sustainability. Extended periods of low organic carbon load, which are the most common conditions that exist in full-scale WWTPs, also had an impact on PAO and GAO activity. GAOs exhibited a substantially higher biomass decay rate as compared to PAOs under these conditions, which revealed a higher survival capacity for PAOs, representing an advantage for PAOs in EBPR processes. This superior survival capacity of PAOs under conditions more closely resembling a full-scale environment was linked with their ability to maintain a residual level of PHA reserves for longer than GAOs, providing them with an effective energy source for aerobic maintenance processes. Overall, this work shows that each of these key operational conditions play an important role in the PAO-GAO competition and should be considered in WWTP models in order to improve EBPR processes.

Development of 3D in vitro models for prediction of hepatic metabolism and toxicity

The development of human cell models that recapitulate hepatic functionality allows the study of metabolic pathways involved in toxicity and disease. The increased biological relevance, cost-effectiveness and high-throughput of cell models can contribute to increase the efficiency of drug development in the pharmaceutical industry. Recapitulation of liver functionality in vitro requires the development of advanced culture strategies to mimic in vivo complexity, such as 3D culture, co-cultures or biomaterials. However, complex 3D models are typically associated with poor robustness, limited scalability and compatibility with screening methods. In this work, several strategies were used to develop highly functional and reproducible spheroid-based in vitro models of human hepatocytes and HepaRG cells using stirred culture systems. In chapter 2, the isolation of human hepatocytes from resected liver tissue was implemented and a liver tissue perfusion method was optimized towards the improvement of hepatocyte isolation and aggregation efficiency, resulting in an isolation protocol compatible with 3D culture. In chapter 3, human hepatocytes were co-cultivated with mesenchymal stem cells (MSC) and the phenotype of both cell types was characterized, showing that MSC acquire a supportive stromal function and hepatocytes retain differentiated hepatic functions, stability of drug metabolism enzymes and higher viability in co-cultures. In chapter 4, a 3D alginate microencapsulation strategy for the differentiation of HepaRG cells was evaluated and compared with the standard 2D DMSO-dependent differentiation, yielding higher differentiation efficiency, comparable levels of drug metabolism activity and significantly improved biosynthetic activity. The work developed in this thesis provides novel strategies for 3D culture of human hepatic cell models, which are reproducible, scalable and compatible with screening platforms. The phenotypic and functional characterization of the in vitro systems performed contributes to the state of the art of human hepatic cell models and can be applied to the improvement of pre-clinical drug development efficiency of the process, model disease and ultimately, development of cell-based therapeutic strategies for liver failure.

Biological colonization on majolica glazed tiles: biodeterioration, bioreceptivity and mitigation strategies

The impact of microbial activity on the deterioration of cultural heritage is a well-recognized global problem. Glazed wall tiles constitute an important part of the worldwide cultural heritage. When exposed outdoors, biological colonization and consequently biodeterioration may occur. Few studies have dealt with this issue, as shown in the literature review on biodiversity, biodeterioration and bioreceptivity of architectural ceramic materials. Due to the lack of knowledge on the biodeteriogens affecting these assets, the characterization of microbial communities growing on Portuguese majolica glazed tiles, from Pena National Palace (Sintra, Portugal) and another from Casa da Pesca (Oeiras, Portugal) was carried out by culture and molecular biology techniques. Microbial communities were composed of microalgae, cyanobacteria, bacteria and fungi, including a new fungal species (Devriesia imbrexigena) described for the first time. Laboratory-based colonization experiments were performed to assess the biodeterioration patterns and bioreceptivity of glazed wall tiles produced in laboratory. Microorganisms previously identified on glazed tiles were inoculated on pristine and artificially aged tile models and incubated under laboratory conditions for 12 months. Phototrophic microorganisms were able to grow into glaze fissures and the tested fungus was able to form oxalates over the glaze. The bioreceptivity of artificially aged tiles was higher for phototrophic microorganisms than pristine tile models. A preliminary approach on mitigation strategies based on in situ application of commercial biocides and titanium dioxide (TiO2) nanoparticles on glazed tiles demonstrated that commercial biocides did not provide long term protection. In contrast, TiO2 treatment caused biofilm detachment. In addition, the use of TiO2 thin films on glazed wall tiles as a protective coating to prevent biological colonization was analysed under laboratorial conditions. Finally, conservation notes on tiles exposed to biological colonization were presented.

Exploring the role of O-Mannosylation in the modulation of E-Cadherin function in Gastric Cancer cell line models

Dissertação de mestrado em Bioquímica Aplicada – Biomedicina

Modeling sequencing errors by combining Hidden Markov models.

Among the largest resources for biological sequence data is the large amount of expressed sequence tags (ESTs) available in public and proprietary databases. ESTs provide information on transcripts but for technical reasons they often contain sequencing errors. Therefore, when analyzing EST sequences computationally, such errors must be taken into account. Earlier attempts to model error prone coding regions have shown good performance in detecting and predicting these while correcting sequencing errors using codon usage frequencies. In the research presented here, we improve the detection of translation start and stop sites by integrating a more complex mRNA model with codon usage bias based error correction into one hidden Markov model (HMM), thus generalizing this error correction approach to more complex HMMs. We show that our method maintains the performance in detecting coding sequences.

Caracteritació de models funcionals d'àmfores romanes mitjançant Finite Element Analysis

Projecte de recerca elaborat a partir d’una estada al Laboratory of Archaeometry del National Centre of Scientific Research “Demokritos” d’Atenes, Grècia, entre juny i setembre 2006. Aquest estudi s’emmarca dins d’un context més ampli d’estudi del canvi tecnològic que es documenta en la producció d’àmfores de tipologia romana durant els segles I aC i I dC en els territoris costaners de Catalunya. Una part d’aquest estudi contempla el càlcul de les propietats mecàniques d’aquestes àmfores i la seva avaluació en funció de la tipologia amforal, a partir de l’Anàlisi d’Elements Finits (AEF). L’AEF és una aproximació numèrica que té el seu origen en les ciències d’enginyeria i que ha estat emprada per estimar el comportament mecànic d’un model en termes, per exemple, de deformació i estrès. Així, un objecte, o millor dit el seu model, es dividit en sub-dominis anomenats elements finits, als quals se’ls atribueixen les propietats mecàniques del material en estudi. Aquests elements finits estan connectats formant una xarxa amb constriccions que pot ser definida. En el cas d’aplicar una força determinada a un model, el comportament de l’objecte pot ser estimat mitjançant el conjunt d’equacions lineals que defineixen el rendiment dels elements finits, proporcionant una bona aproximació per a la descripció de la deformació estructural. Així, aquesta simulació per ordinador suposa una important eina per entendre la funcionalitat de ceràmiques arqueològiques. Aquest procediment representa un model quantitatiu per predir el trencament de l’objecte ceràmic quan aquest és sotmès a diferents condicions de pressió. Aquest model ha estat aplicat a diferents tipologies amforals. Els resultats preliminars mostren diferències significatives entre la tipologia pre-romana i les tipologies romanes, així com entre els mateixos dissenys amforals romans, d’importants implicacions arqueològiques.

Two different 3D biogeochemical models for the NW mediterranean sea: validation with Meris data and intercomparison

Report for the scientific sojourn carried out at the University of New South Wales from February to June the 2007. Two different biogeochemical models are coupled to a three dimensional configuration of the Princeton Ocean Model (POM) for the Northwestern Mediterranean Sea (Ahumada and Cruzado, 2007). The first biogeochemical model (BLANES) is the three-dimensional version of the model described by Bahamon and Cruzado (2003) and computes the nitrogen fluxes through six compartments using semi-empirical descriptions of biological processes. The second biogeochemical model (BIOMEC) is the biomechanical NPZD model described in Baird et al. (2004), which uses a combination of physiological and physical descriptions to quantify the rates of planktonic interactions. Physical descriptions include, for example, the diffusion of nutrients to phytoplankton cells and the encounter rate of predators and prey. The link between physical and biogeochemical processes in both models is expressed by the advection-diffusion of the non-conservative tracers. The similarities in the mathematical formulation of the biogeochemical processes in the two models are exploited to determine the parameter set for the biomechanical model that best fits the parameter set used in the first model. Three years of integration have been carried out for each model to reach the so called perpetual year run for biogeochemical conditions. Outputs from both models are averaged monthly and then compared to remote sensing images obtained from sensor MERIS for chlorophyll.

Transcriptional response to cardiac injury in the zebrafish: systematic identification of genes with highly concordant activity across in vivo models.

BACKGROUND: Zebrafish is a clinically-relevant model of heart regeneration. Unlike mammals, it has a remarkable heart repair capacity after injury, and promises novel translational applications. Amputation and cryoinjury models are key research tools for understanding injury response and regeneration in vivo. An understanding of the transcriptional responses following injury is needed to identify key players of heart tissue repair, as well as potential targets for boosting this property in humans. RESULTS: We investigated amputation and cryoinjury in vivo models of heart damage in the zebrafish through unbiased, integrative analyses of independent molecular datasets. To detect genes with potential biological roles, we derived computational prediction models with microarray data from heart amputation experiments. We focused on a top-ranked set of genes highly activated in the early post-injury stage, whose activity was further verified in independent microarray datasets. Next, we performed independent validations of expression responses with qPCR in a cryoinjury model. Across in vivo models, the top candidates showed highly concordant responses at 1 and 3 days post-injury, which highlights the predictive power of our analysis strategies and the possible biological relevance of these genes. Top candidates are significantly involved in cell fate specification and differentiation, and include heart failure markers such as periostin, as well as potential new targets for heart regeneration. For example, ptgis and ca2 were overexpressed, while usp2a, a regulator of the p53 pathway, was down-regulated in our in vivo models. Interestingly, a high activity of ptgis and ca2 has been previously observed in failing hearts from rats and humans. CONCLUSIONS: We identified genes with potential critical roles in the response to cardiac damage in the zebrafish. Their transcriptional activities are reproducible in different in vivo models of cardiac injury.