Effect of vehicle pretreatment on the flux, retention, and diffusion of topically applied penetrants in vitro

Purpose. The flux of a topically applied drug depends on the activity in the skin and the interaction between the vehicle and skin. Permeation of vehicle into the skin can alter the activity of drug and the properties of the skin barrier. The aim of this in vitro study was to separate and quantify these effects. Methods. The flux of four radiolabeled permeants (water, phenol, diflunisal, and diazepam) with log K-oct/water values from 1.4 to 4.3 was measured over 4 h through heat-separated human epidermis pretreated for 30 min with vehicles having Hildebrand solubility parameters from 7.9 to 23.4 (cal/cm(3))(1/2). Results. Enhancement was greatest after pretreatment with the more lipophilic vehicles. A synergistic enhancement was observed using binary mixtures. The flux of diazepam was not enhanced to the same extent as the other permeants, possibly because its partitioning into the epidermis is close to optimal (log K-oct 2.96). Conclusion. An analysis of the permeant remaining in the epidermis revealed that the enhancement can be the result of either increased partitioning of permeant into the epidermis or an increasing diffusivity of permeants through the epidermis.

Cristalização de Biomateriais Vitrocerâmicos e Mineralização em Meio Fisiológico Simulado

We propose new theoretical models, which generalize the classical Avrami-Nakamura models. These models are suitable to describe the kinetics of nucleation and growth in transient regime, and/or with overlapping of nucleation and growth. Simulations and predictions were performed for lithium disilicate based on data reported in the literature. One re-examined the limitations of the models currently used to interpret DTA or DSC results, and to extract the relevant kinetic parameters. Glasses and glass-ceramics with molar formulation 0.45SiO2? (0.45-x)MgO?xK2O?0.1(3CaO.P2O5) (0?x?0.090) were prepared, crystallized and studied as potential materials for biomedical applications. Substitution of K+ for Mg2+ were used to prevent devritification on cooling, to adjust the kinetics of crystallization and to modify the in vitro behaviour of resulting biomaterials. The crystallization of the glass frits was studied by DTA, XRD and SEM. Exothermic peaks were detected corresponding to bulk crystallization of whitlockite-type phosphate, Ca9MgK(PO4)7, at approximately 900ºC, and surface crystallization of a predominant forsterite phase (Mg2SiO4) at higher temperatures. XRD also revealed the presence of diopside (CaMgSi2O6 in some samples. The predominant microstructure of the phosphate phase is of the plate-type, seemingly crystallizing by a 2-dimensional growth mechanism. Impedance spectroscopy revealed significant changes in electrical behaviour, associated to crystallization of the phosphate phase. This showed that electrical measurements can be used to study the kinetics of crystallization for cases when DTA or DSC experiments reveal limitations, and to extract estimates of relevant parameters from the dependence of crystallization peak temperature, and its width at half height. In vitro studies of glasses and glass-ceramics in acelular SBF media showed bioactivity and the development of apatite layers The morphology, composition and adhesion of the apatite layer could be changed by substitution of Mg2+ by K+. Apatite layers were deposited on the surface of glass-ceramics of the nominal compositions with x=0 and 0.09, in contact with SBF at 37ºC. The adhesion of the apatite layer was quantified by the scratch test technique, having been related with SBF?s immersion time, with composition and structure of the glass phase, and with the morphology of the crystalline phase of the glass-ceramics. The structure of three glasses (x=0, 0.045 and 0.090) were investigated by MAS-NMR ( 29Si and 31P), showing that the fraction of Q3 structural units increases with the contents of Mg, and that the structure of these glasses includes orthophosphate groups (PO43-) preferentially connected to Ca2+ ions. Mg2+ ions show preference towards the silicate network. Substitution of Mg2+ by K+ allowed one to change the bioactivity. FTIR data revealed octacalcium phosphate precipitation (Ca8H2(PO4)6.5H2O) in the glass without K, while the morphology of the layer acquires the shape of partially superimposed hemispheres, spread over the surface. The glasses with K present a layer of acicular hidroxyapatite, whose crystallinity and needles thickness tend to increase along with K content.

Geração de simuladores de filas para diferentes contextos com estudo de casos para redes de computadores

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

A coinfection model for HIV and HCV

We study a mathematical model for the human immunodeficiency virus (HIV) and hepatites C virus (HCV) coinfection. The model predicts four distinct equilibria: the disease free, the HIV endemic, the HCV endemic, and the full endemic equilibria. The local and global stability of the disease free equilibrium was calculated for the full model and the HIV and HCV submodels. We present numerical simulations of the full model where the distinct equilibria can be observed. We show simulations of the qualitative changes of the dynamical behavior of the full model for variation of relevant parameters. From the results of the model, we infer possible measures that could be implemented in order to reduce the number of infected individuals.

Estrutura a partir do Movimento em Sistemas Autónomos

Mestrado em Engenharia Electrotécnica e de Computadores - Ramo de Sistemas Autónomos

Estimation of infection rate in epidemic models with multiple populations

Dissertação para obtenção do Grau de Mestre em Matemática e Aplicações Especialização em Actuariado, Estatística e Investigação Operacional

Predisposition assessment modelling

Dissertation submitted for obtaining the degree of Master in Environmental Engineering

Analysis of output properties of the Peroxiredoxin/Thioredoxin/Thioredoxin reductase system

Dissertação de mestrado em Bioinformática

Modelling Amazonian forest eddy covariance data: a comparison of big leaf versus sun/shade models for the C-14 tower at Manaus I. Canopy photosynthesis

In this study, we concentrate on modelling gross primary productivity using two simple approaches to simulate canopy photosynthesis: "big leaf" and "sun/shade" models. Two approaches for calibration are used: scaling up of canopy photosynthetic parameters from the leaf to the canopy level and fitting canopy biochemistry to eddy covariance fluxes. Validation of the models is achieved by using eddy covariance data from the LBA site C14. Comparing the performance of both models we conclude that numerically (in terms of goodness of fit) and qualitatively, (in terms of residual response to different environmental variables) sun/shade does a better job. Compared to the sun/shade model, the big leaf model shows a lower goodness of fit and fails to respond to variations in the diffuse fraction, also having skewed responses to temperature and VPD. The separate treatment of sun and shade leaves in combination with the separation of the incoming light into direct beam and diffuse make sun/shade a strong modelling tool that catches more of the observed variability in canopy fluxes as measured by eddy covariance. In conclusion, the sun/shade approach is a relatively simple and effective tool for modelling photosynthetic carbon uptake that could be easily included in many terrestrial carbon models.

Desenvolupament d'interfícies HW per a FPGA des de Matlab

Aquest projecte consisteix en el desenvolupament d’estructures hardware digitals, sintetitzables sobre FPGA i realitzades des d’un entorn gràfic de disseny a nivell de sistema (alt nivell). S'ha escollit el Simulink (entorn gràfic que treballa sobre el software matemàtic Matlab de Mathworks) com a entorn de disseny, i que gràcies a la interfície proporcionada per Altera (DSPBuilder) és capaç de generar codi VHDL sintetitzable. Concretament ens centrarem en la gestió d’un sistema capturador d'imatges de comptadors del cabal d'aigua, en el qual volem fer la caracterització del comptador. Aquest capturador consta bàsicament d'un sensor d'imatge i una FPGA. En aquesta caracterització el que es pretén es ajustar els diferents paràmetres del sistema per fer que la lectura sigui òptima per a cada model de comptador que existeixen al mercat, com ara l'exposició del sensor, el guany d'un color, la realització d'un filtrat de la imatge, etc.

Understanding the low-temperature limitations to forest growth through calibration of a forest dynamics model with tree-ring data

The sensitivity of altitudinal and latitudinal tree-line ecotones to climate change, particularly that of temperature, has received much attention. To improve our understanding of the factors affecting tree-line position, we used the spatially explicit dynamic forest model TreeMig. Although well-suited because of its landscape dynamics functions, TreeMig features a parabolic temperature growth response curve, which has recently been questioned. and the species parameters are not specifically calibrated for cold temperatures. Our main goals were to improve the theoretical basis of the temperature growth response curve in the model and develop a method for deriving that curve's parameters from tree-ring data. We replaced the parabola with an asymptotic curve, calibrated for the main species at the subalpine (Swiss Alps: Pinus cembra, Larix decidua, Picea abies) and boreal (Fennoscandia: Pinus sylvestris, Betula pubescens, P. abies) tree-lines. After fitting new parameters, the growth curve matched observed tree-ring widths better. For the subalpine species, the minimum degree-day sum allowing, growth (kDDMin) was lowered by around 100 degree-days; in the case of Larix, the maximum potential ring-width was increased to 5.19 mm. At the boreal tree-line, the kDDMin for P. sylvestris was lowered by 210 degree-days and its maximum ring-width increased to 2.943 mm; for Betula (new in the model) kDDMin was set to 325 degree-days and the maximum ring-width to 2.51 mm; the values from the only boreal sample site for Picea were similar to the subalpine ones, so the same parameters were used. However, adjusting the growth response alone did not improve the model's output concerning species' distributions and their relative importance at tree-line. Minimum winter temperature (MinWiT, mean of the coldest winter month), which controls seedling establishment in TreeMig, proved more important for determining distribution. Picea, P. sylvestris and Betula did not previously have minimum winter temperature limits, so these values were set to the 95th percentile of each species' coldest MinWiT site (respectively -7, -11, -13). In a case study for the Alps, the original and newly calibrated versions of TreeMig were compared with biomass data from the National Forest Inventor), (NFI). Both models gave similar, reasonably realistic results. In conclusion, this method of deriving temperature responses from tree-rings works well. However, regeneration and its underlying factors seem more important for controlling species' distributions than previously thought. More research on regeneration ecology, especially at the upper limit of forests. is needed to improve predictions of tree-line responses to climate change further.

Analisi espacial del canvi climàtic al continent europeu. Incerteses i exposicions principals

El canvi climàtic del segle XXI és una realitat, hi ha moltes evidències científiques que indiquen que l’escalfament del sistema climàtic és inequívoc. Malgrat això, també hi ha moltes incerteses respecte els impactes que pot comportar aquest canvi climàtic global. L’objectiu d’aquest projecte és estudiar la possible evolució futura de tres variables climàtiques, que són el rang de la temperatura diürna a prop de la superfície (DTR), la temperatura mitjana a prop de la superfície (MT) i la precipitació mensual (PL_mes) i valorar l’exposició que poden experimentar diferents cobertes del sòl i diferents regions biogeogràfiques del continent europeu davant d’aquests possibles patrons de canvi. Per això s’han utilitzat Models Climàtics Globals que fan projeccions de variables climàtiques que permeten preveure el possible clima futur. Mitjançant l’aplicatiu informàtic Tetyn s’han extret els paràmetres climàtics dels conjunts de dades del Tyndall Centre for Climate Change Research, del futur (TYN SC) i del passat (CRU TS). Les variables obtingudes s’han processat amb eines de sistemes d’informació geogràfica (SIG) per obtenir els patrons de canvi de les variables a cada coberta del sòl. Els resultats obtinguts mostren que hi ha una gran variabilitat, que augmenta amb el temps, entre els diferents models climàtics i escenaris considerats, que posa de manifest la incertesa associada a la modelització climàtica, a la generació d’escenaris d’emissions i a la naturalesa dinàmica i no determinista del sistema climàtic. Però en general, mostren que les glaceres seran una de les cobertes més exposades al canvi climàtic, i la mediterrània, una de les regions més vulnerables

A gamma camera count rate saturation correction method for whole-body planar imaging.

Whole-body (WB) planar imaging has long been one of the staple methods of dosimetry, and its quantification has been formalized by the MIRD Committee in pamphlet no 16. One of the issues not specifically addressed in the formalism occurs when the count rates reaching the detector are sufficiently high to result in camera count saturation. Camera dead-time effects have been extensively studied, but all of the developed correction methods assume static acquisitions. However, during WB planar (sweep) imaging, a variable amount of imaged activity exists in the detector's field of view as a function of time and therefore the camera saturation is time dependent. A new time-dependent algorithm was developed to correct for dead-time effects during WB planar acquisitions that accounts for relative motion between detector heads and imaged object. Static camera dead-time parameters were acquired by imaging decaying activity in a phantom and obtaining a saturation curve. Using these parameters, an iterative algorithm akin to Newton's method was developed, which takes into account the variable count rate seen by the detector as a function of time. The algorithm was tested on simulated data as well as on a whole-body scan of high activity Samarium-153 in an ellipsoid phantom. A complete set of parameters from unsaturated phantom data necessary for count rate to activity conversion was also obtained, including build-up and attenuation coefficients, in order to convert corrected count rate values to activity. The algorithm proved successful in accounting for motion- and time-dependent saturation effects in both the simulated and measured data and converged to any desired degree of precision. The clearance half-life calculated from the ellipsoid phantom data was calculated to be 45.1 h after dead-time correction and 51.4 h with no correction; the physical decay half-life of Samarium-153 is 46.3 h. Accurate WB planar dosimetry of high activities relies on successfully compensating for camera saturation which takes into account the variable activity in the field of view, i.e. time-dependent dead-time effects. The algorithm presented here accomplishes this task.

Study and implementation of some quantitative trading models

Quantitative or algorithmic trading is the automatization of investments decisions obeying a fixed or dynamic sets of rules to determine trading orders. It has increasingly made its way up to 70% of the trading volume of one of the biggest financial markets such as the New York Stock Exchange (NYSE). However, there is not a signi cant amount of academic literature devoted to it due to the private nature of investment banks and hedge funds. This projects aims to review the literature and discuss the models available in a subject that publications are scarce and infrequently. We review the basic and fundamental mathematical concepts needed for modeling financial markets such as: stochastic processes, stochastic integration and basic models for prices and spreads dynamics necessary for building quantitative strategies. We also contrast these models with real market data with minutely sampling frequency from the Dow Jones Industrial Average (DJIA). Quantitative strategies try to exploit two types of behavior: trend following or mean reversion. The former is grouped in the so-called technical models and the later in the so-called pairs trading. Technical models have been discarded by financial theoreticians but we show that they can be properly cast into a well defined scientific predictor if the signal generated by them pass the test of being a Markov time. That is, we can tell if the signal has occurred or not by examining the information up to the current time; or more technically, if the event is F_t-measurable. On the other hand the concept of pairs trading or market neutral strategy is fairly simple. However it can be cast in a variety of mathematical models ranging from a method based on a simple euclidean distance, in a co-integration framework or involving stochastic differential equations such as the well-known Ornstein-Uhlenbeck mean reversal ODE and its variations. A model for forecasting any economic or financial magnitude could be properly defined with scientific rigor but it could also lack of any economical value and be considered useless from a practical point of view. This is why this project could not be complete without a backtesting of the mentioned strategies. Conducting a useful and realistic backtesting is by no means a trivial exercise since the \laws" that govern financial markets are constantly evolving in time. This is the reason because we make emphasis in the calibration process of the strategies' parameters to adapt the given market conditions. We find out that the parameters from technical models are more volatile than their counterpart form market neutral strategies and calibration must be done in a high-frequency sampling manner to constantly track the currently market situation. As a whole, the goal of this project is to provide an overview of a quantitative approach to investment reviewing basic strategies and illustrating them by means of a back-testing with real financial market data. The sources of the data used in this project are Bloomberg for intraday time series and Yahoo! for daily prices. All numeric computations and graphics used and shown in this project were implemented in MATLAB^R scratch from scratch as a part of this thesis. No other mathematical or statistical software was used.

S-system parameter estimation for noisy metabolic profiles using newton-flow analysis.

Biochemical systems are commonly modelled by systems of ordinary differential equations (ODEs). A particular class of such models called S-systems have recently gained popularity in biochemical system modelling. The parameters of an S-system are usually estimated from time-course profiles. However, finding these estimates is a difficult computational problem. Moreover, although several methods have been recently proposed to solve this problem for ideal profiles, relatively little progress has been reported for noisy profiles. We describe a special feature of a Newton-flow optimisation problem associated with S-system parameter estimation. This enables us to significantly reduce the search space, and also lends itself to parameter estimation for noisy data. We illustrate the applicability of our method by applying it to noisy time-course data synthetically produced from previously published 4- and 30-dimensional S-systems. In addition, we propose an extension of our method that allows the detection of network topologies for small S-systems. We introduce a new method for estimating S-system parameters from time-course profiles. We show that the performance of this method compares favorably with competing methods for ideal profiles, and that it also allows the determination of parameters for noisy profiles.