Master Equation: Biological Applications and Thermodynamic Description

It is well known that many realistic mathematical models of biological systems, such as cell growth, cellular development and differentiation, gene expression, gene regulatory networks, enzyme cascades, synaptic plasticity, aging and population growth need to include stochasticity. These systems are not isolated, but rather subject to intrinsic and extrinsic fluctuations, which leads to a quasi equilibrium state (homeostasis). The natural framework is provided by Markov processes and the Master equation (ME) describes the temporal evolution of the probability of each state, specified by the number of units of each species. The ME is a relevant tool for modeling realistic biological systems and allow also to explore the behavior of open systems. These systems may exhibit not only the classical thermodynamic equilibrium states but also the nonequilibrium steady states (NESS). This thesis deals with biological problems that can be treat with the Master equation and also with its thermodynamic consequences. It is organized into six chapters with four new scientific works, which are grouped in two parts: (1) Biological applications of the Master equation: deals with the stochastic properties of a toggle switch, involving a protein compound and a miRNA cluster, known to control the eukaryotic cell cycle and possibly involved in oncogenesis and with the propose of a one parameter family of master equations for the evolution of a population having the logistic equation as mean field limit. (2) Nonequilibrium thermodynamics in terms of the Master equation: where we study the dynamical role of chemical fluxes that characterize the NESS of a chemical network and we propose a one parameter parametrization of BCM learning, that was originally proposed to describe plasticity processes, to study the differences between systems in DB and NESS.

Improved methods for predicting individual leaf area and leaf senescence in maize (Zea mays)

The ability to predict leaf area and leaf area index is crucial in crop simulation models that predict crop growth and yield. Previous studies have shown existing methods of predicting leaf area to be inadequate when applied to a broad range of cultivars with different numbers of leaves. The objectives of the study were to (i) develop generalised methods of modelling individual and total plant leaf area, and leaf senescence, that do not require constants that are specific to environments and/or genotypes, (ii) re-examine the base, optimum, and maximum temperatures for calculation of thermal time for leaf senescence, and (iii) assess the method of calculation of individual leaf area from leaf length and leaf width in experimental work. Five cultivars of maize differing widely in maturity and adaptation were planted in October 1994 in south-eastern Queensland, and grown under non-limiting conditions of water and plant nutrient supplies. Additional data for maize plants with low total leaf number (12-17) grown at Katumani Research Centre, Kenya, were included to extend the range in the total leaf number per plant. The equation for the modified (slightly skewed) bell curve could be generalised for modelling individual leaf area, as all coefficients in it were related to total leaf number. Use of coefficients for individual genotypes can be avoided, and individual and total plant leaf area can be calculated from total leaf number. A single, logistic equation, relying on maximum plant leaf area and thermal time from emergence, was developed to predict leaf senescence. The base, optimum, and maximum temperatures for calculation of thermal time for leaf senescence were 8, 34, and 40 degrees C, and apply for the whole crop-cycle when used in modelling of leaf senescence. Thus, the modelling of leaf production and senescence is simplified, improved, and generalised. Consequently, the modelling of leaf area index (LAI) and variables that rely on LAI will be improved. For experimental purposes, we found that the calculation of leaf area from leaf length and leaf width remains appropriate, though the relationship differed slightly from previously published equations.

Investigation of the electrical and electrochemical properties of nanocomposites from V2O5, polypyrrole, and polyaniline

In this work, an investigation of the electrical and electrochemical properties responsible for the energy storage capability of nanocomposites has been carried out. We demonstrate that, in the case of the V2O5 xerogel and the nanocomposites polypyrrole (Ppy)/V2O5 and polyaniline (PANI)/V2O5, the quadratic logistic equation (QLE) can be used to fit the inverse of the resistance values as a function of the injected charge in non-steady-state conditions. This contributes to a phenomenological understanding of the lithium ion and electron transport. The departure of the experimental curve from the fitting observed for the V2O5 xerogel can be attributed to the trapping sites formed during the lithium electroinsertion, which was observed by electrochemical impedance spectroscopy. The amount of trapping sites was obtained on the basis of the QLE. Similar values used to fit the inverse of the resistance were also used to fit the absorbance changes, which is also associated with the small polaron hopping from the V(IV) to the V(V) sites. On the other hand, there was good agreement between the experimental and the theoretical data when the profile of the inverse of the resistance as a function of the amount of inserted lithium ions of the nanocomposites Ppy/V2O5 and PANI/ V2O5 was concerned. We suggest that the presence of the conducting polymers is responsible for the different electrical profile of the V2O5 xerogel compared with those of the nanocomposites. In the latter case, interactions between the lithium ions and oxygen atoms from V2O5 are shielded, thus decreasing the trapping effect of lithium ions in the V2O5 sites. The different values of the lithium ion diffusion coefficient into these intercalation materials are in agreement with this hypothesis.

Self-assembled films from WO(3): Electrochromism and lithium ion diffusion

WO(3)/chitosan and WO(3)/chitosan/poly(ethylene oxide) (PEO) films were prepared by the layer-by-layer method. The presence of chitosan enabled PEO to be carried into the self-assembled structure, contributing to an increase in the Li(+) diffusion rate. On the basis of the galvanostatic intermittent titration technique (GITT) and the quadratic logistic equation (QLE), a spectroelectrochemical method was used for determination of the ""optical"" diffusion coefficient (D(op)), enabling analysis of the Li(+) diffusion rate and, consequently, the coloration front rate in these host matrices. The D(op) values within the WO(3)/chitosan/PEO film were significantly higher than those within the WO(3)/chitosan film, mainly for higher values of injected charge. The presence of PEO also ensured larger accessibility to the electroactive sites, in accordance with the method employed here. Hence, this spectroelectrochemical method allowed us to separate the contribution of the diffusion process from the number of accessible electroactive sites in the materials, thereby aiding a better understanding of the useful electrochemical and electrochromic properties of these films for use in electrochromic devices. (C) 2010 Elsevier B.V. All rights reserved.

Development and survival of the diamondback moth (Lepidoptera : Plutellidae) at constant and alternating temperatures

Survival and development time from egg to adult emergence of the diamondback moth, Plutella xylostella (L.), were determined at 19 constant and 14 alternating temperature regimes from 4 to 40degreesC. Plutella xylostella developed successfully front egg to adult emergence at constant temperatures from 8 to 32degreesC. At temperatures from 4 to 6degreesC or from 34 to 40degreesC, partial or complete development of individual stages or instars was possible, with third and fourth instars having the widest temperature limits. The insect developed successfully from egg to adult emergence under alternating regimes including temperatures as low as 4degreesC or as high as 38degreesC. The degree-day model, the logistic equation, and the Wang model were used to describe the relationships between temperature and development rate at both constant and alternating temperatures. The degree-day model described the relationships well from 10 to 30degreesC. The logistic equation and the Wang model fit the data well at temperatures 32degreesC. Under alternating regimes, all three models gave good simulations of development in the mid-temperature range, but only the logistic equation gave close simulations in the low temperature range, and none gave close or consistent simulations in the high temperature range. The distribution of development time was described satisfactorily by a Weibull function. These rate and time distribution functions provide tools for simulating population development of P. xylostella over a wide range of temperature conditions.

Aspectos relacionados à utilização da equação logística quadrática em processos eletroquímicos

The concepts of dissipation and feedback are contained in the behavior of many natural dynamical systems. They have been used to predict the evolution of populations leading to the formulation of the quadratic logistic equation (QLE). More recently, the QLE has been used to provide a better understanding of physicochemical systems with promising results. Many physical, chemical and biological dynamic phenomena can be understood on the basis of the QLE and this work describes the main aspects of this equation and some recent applications, with emphasis on electrochemical systems. Also, it is illustrated the concept of potential energy as a convenient way of describing the stability of the fixed points of the QLE.

Cinética e caracterização de ramnolipídeos produzidos por Pseudomonas aeruginosa MSIC02 utilizando glicerol como fonte de carbono

Glycerol, a co-product of biodiesel production, was used as a carbon source for the kinetics studies and production of biosurfactants by P. aeruginosa MSIC02. The highest fermentative parameters (Y PX = 3.04 g g-1; Y PS = 0.189 g g-1, P B = 31.94 mg L-1 h-1 and P X = 10.5 mg L-1 h-1) were obtained at concentrations of 0.4% (w/v) NaNO3 and 2% (w/v) glycerol. The rhamnolipid exhibited 80% of emulsification on kerosene, surface tension of 32.5 mN m-1, CMC = 28.2 mg L-1, C20 (concentration of surfactant in the bulk phase that produces a reduction of 20 dyn/cm in the surface tension of the solvent) = 0.99 mg L-1, Γm (surface concentration excess) = 2.4 x 10-26 mol Å-2 and S (surface area) = 70.4 Ų molecule-1 with solutions containing 10% NaCl. A mathematical model based on logistic equation was considered to representing the process. Model parameters were estimated by non-linear regression method. This approach was able to give a good description of the process.

Use of a predictive model for food insecurity estimates in Brazil

In 2004 the National Household Survey (Pesquisa Nacional par Amostras de Domicilios - PNAD) estimated the prevalence of food and nutrition insecurity in Brazil. However, PNAD data cannot be disaggregated at the municipal level. The objective of this study was to build a statistical model to predict severe food insecurity for Brazilian municipalities based on the PNAD dataset. Exclusion criteria were: incomplete food security data (19.30%); informants younger than 18 years old (0.07%); collective households (0.05%); households headed by indigenous persons (0.19%). The modeling was carried out in three stages, beginning with the selection of variables related to food insecurity using univariate logistic regression. The variables chosen to construct the municipal estimates were selected from those included in PNAD as well as the 2000 Census. Multivariate logistic regression was then initiated, removing the non-significant variables with odds ratios adjusted by multiple logistic regression. The Wald Test was applied to check the significance of the coefficients in the logistic equation. The final model included the variables: per capita income; years of schooling; race and gender of the household head; urban or rural residence; access to public water supply; presence of children; total number of household inhabitants and state of residence. The adequacy of the model was tested using the Hosmer-Lemeshow test (p=0.561) and ROC curve (area=0.823). Tests indicated that the model has strong predictive power and can be used to determine household food insecurity in Brazilian municipalities, suggesting that similar predictive models may be useful tools in other Latin American countries.

Cascades of Hopf bifurcations from boundary delay

We consider a 1-dimensional reaction-diffusion equation with nonlinear boundary conditions of logistic type with delay. We deal with non-negative solutions and analyze the stability behavior of its unique positive equilibrium solution, which is given by the constant function u equivalent to 1. We show that if the delay is small, this equilibrium solution is asymptotically stable, similar as in the case without delay. We also show that, as the delay goes to infinity, this equilibrium becomes unstable and undergoes a cascade of Hopf bifurcations. The structure of this cascade will depend on the parameters appearing in the equation. This equation shows some dynamical behavior that differs from the case where the nonlinearity with delay is in the interior of the domain. (C) 2009 Elsevier Inc. All rights reserved.

Allometry of the secondary sexual characters of the shore crab Pachygrapsus transversus (Gibbes, 1850) (Brachyura, Grapsidae)

The allometric growth of secondary sexual characters in Pachygrapsus transversus is investigated from the 2(nd) crab stage onward. Clear sexual dimorphism is restricted to abdominal morphology, but ANCOVA analyses showed that chelae become larger in males and the carapace becomes wider in females. Size at the puberty moult in both sexes was estimated using Somerton's computer techniques. Mature II analyses applied to bi-log gonopod length vs, carapace length relationships indicated a puberty moult at 5.0 mm in males.In females, Mature I analyses detected the overlapping growth phase lines in bi-log carapace length vs. abdomen width scatterplots. Fitting the logistic equation provided an estimate of 50% maturity at 5.5 mm. The regression lines separate young and resting individuals from the potentially reproductive females, but they do not separate young from adult crabs. Year-round monthly samples showed that the proportion of small adult-like females is higher during the breeding season. After breeding, females may moult to a young-like morphotype, as observed in controlled laboratory conditions. Moulting to a resting condition splits smaller mature females into different growth phase lines. Therefore, estimates of female size at sexual maturity by means of abdomen allometric growth analyses are inadequate in this species.

Das transformadas integrais ao cálculo fracionário aplicado à equação logística

Pós-graduação em Biometria - IBB

A prelude to the fractional calculus applied to tumor dynamic

In order to refine the solution given by the classical logistic equation and extend its range of applications in the study of tumor dynamics, we propose and solve a generalization of this equation, using the so-called Fractional Calculus, i.e., we replace the ordinary derivative of order 1, in one version of the usual equation, by a non-integer derivative of order 0 < α < 1, and recover the classical solution as a particular case. Finally, we analyze the applicability of this model to describe the growth of cancer tumors.

Cifrado de imágenes digitales basado en teoría del caos: mapas logísticos

Las nuevas tendencias de compartir archivos multimedia a través de redes abiertas, demanda el uso de mejores técnicas de encriptación que garanticen la integridad, disponibilidad y confidencialidad, manteniendo y/o mejorando la eficiencia del proceso de cifrado sobre estos archivos. Hoy en día es frecuente la transferencia de imágenes a través de medios tecnológicos, siendo necesario la actualización de las técnicas de encriptación existentes y mejor aún, la búsqueda de nuevas alternativas. Actualmente los algoritmos criptográficos clásicos son altamente conocidos en medio de la sociedad informática lo que provoca mayor vulnerabilidad, sin contar los altos tiempos de procesamiento al momento de ser utilizados, elevando la probabilidad de ser descifrados y minimizando la disponibilidad inmediata de los recursos. Para disminuir estas probabilidades, el uso de la teoría de caos surge como una buena opción para ser aplicada en un algoritmo que tome partida del comportamiento caótico de los sistemas dinámicos, y aproveche las propiedades de los mapas logísticos para elevar el nivel de robustez en el cifrado. Es por eso que este trabajo propone la creación de un sistema criptográfico basado sobre una arquitectura dividida en dos etapas de confusión y difusión. Cada una de ellas utiliza una ecuación logística para generar números pseudoaleatorios que permitan desordenar la posición del píxel y cambiar su intensidad en la escala de grises. Este proceso iterativo es determinado por la cantidad total de píxeles de una imagen. Finalmente, toda la lógica de cifrado es ejecutada sobre la tecnología CUDA que permite el procesamiento en paralelo. Como aporte sustancial, se propone una nueva técnica de encriptación vanguardista de alta sensibilidad ante ruidos externos manteniendo no solo la confidencialidad de la imagen, sino también la disponibilidad y la eficiencia en los tiempos de proceso.---ABSTRACT---New trends to share multimedia files over open networks, demand the best use of encryption techniques to ensure the integrity, availability and confidentiality, keeping and/or improving the efficiency of the encryption process on these files. Today it is common to transfer pictures through technological networks, thus, it is necessary to update existing techniques encryption, and even better, the searching of new alternatives. Nowadays, classic cryptographic algorithms are highly known in the midst of the information society which not only causes greater vulnerability, but high processing times when this algorithms are used. It raise the probability of being deciphered and minimizes the immediate availability of resources. To reduce these odds, the use of chaos theory emerged as a good option to be applied on an algorithm that takes advantage of chaotic behavior of dynamic systems, and take logistic maps’ properties to raise the level of robustness in the encryption. That is why this paper proposes the creation of a cryptographic system based on an architecture divided into two stages: confusion and diffusion. Each stage uses a logistic equation to generate pseudorandom numbers that allow mess pixel position and change their intensity in grayscale. This iterative process is determined by the total number of pixels of an image. Finally, the entire encryption logic is executed on the CUDA technology that enables parallel processing. As a substantial contribution, it propose a new encryption technique with high sensitivity on external noise not only keeping the confidentiality of the image, but also the availability and efficiency in processing times.

Logistic regression analysis of potential prognostic factors for pulmonary thromboembolism

Objective: To identify potential prognostic factors for pulmonary thromboembolism (PTE), establishing a mathematical model to predict the risk for fatal PTE and nonfatal PTE.Method: the reports on 4,813 consecutive autopsies performed from 1979 to 1998 in a Brazilian tertiary referral medical school were reviewed for a retrospective study. From the medical records and autopsy reports of the 512 patients found with macroscopically and/or microscopically,documented PTE, data on demographics, underlying diseases, and probable PTE site of origin were gathered and studied by multiple logistic regression. Thereafter, the jackknife method, a statistical cross-validation technique that uses the original study patients to validate a clinical prediction rule, was performed.Results: the autopsy rate was 50.2%, and PTE prevalence was 10.6%. In 212 cases, PTE was the main cause of death (fatal PTE). The independent variables selected by the regression significance criteria that were more likely to be associated with fatal PTE were age (odds ratio [OR], 1.02; 95% confidence interval [CI], 1.00 to 1.03), trauma (OR, 8.5; 95% CI, 2.20 to 32.81), right-sided cardiac thrombi (OR, 1.96; 95% CI, 1.02 to 3.77), pelvic vein thrombi (OR, 3.46; 95% CI, 1.19 to 10.05); those most likely to be associated with nonfatal PTE were systemic arterial hypertension (OR, 0.51; 95% CI, 0.33 to 0.80), pneumonia (OR, 0.46; 95% CI, 0.30 to 0.71), and sepsis (OR, 0.16; 95% CI, 0.06 to 0.40). The results obtained from the application of the equation in the 512 cases studied using logistic regression analysis suggest the range in which logit p > 0.336 favors the occurrence of fatal PTE, logit p < - 1.142 favors nonfatal PTE, and logit P with intermediate values is not conclusive. The cross-validation prediction misclassification rate was 25.6%, meaning that the prediction equation correctly classified the majority of the cases (74.4%).Conclusions: Although the usefulness of this method in everyday medical practice needs to be confirmed by a prospective study, for the time being our results suggest that concerning prevention, diagnosis, and treatment of PTE, strict attention should be given to those patients presenting the variables that are significant in the logistic regression model.

New Equation for Prediction of Reverse Remodeling after Cardiac Resynchronization Therapy

Objectives: To integrate data from two-dimensional echocardiography (2D ECHO), three-dimensional echocardiography (3D ECHO), and tissue Doppler imaging (TDI) for prediction of left ventricular (LV) reverse remodeling (LVRR) after cardiac resynchronization therapy (CRT). It was also compared the evaluation of cardiac dyssynchrony by TDI and 3D ECHO. Methods: Twenty-four consecutive patients with heart failure, sinus rhythm, QRS = 120 msec, functional class III or IV and LV ejection fraction (LVEF) = 0.35 underwent CRT. 2D ECHO, 3D ECHO with systolic dyssynchrony index (SDI) analysis, and TDI were performed before, 3 and 6 months after CRT. Cardiac dyssynchrony analyses by TDI and SDI were compared with the Pearson's correlation test. Before CRT, a univariate analysis of baseline characteristics was performed for the construction of a logistic regression model to identify the best predictors of LVRR. Results: After 3 months of CRT, there was a moderate correlation between TDI and SDI (r = 0.52). At other time points, there was no strong correlation. Nine of twenty-four (38%) patients presented with LVRR 6 months after CRT. After logistic regression analysis, SDI (SDI > 11%) was the only independent factor in the prediction of LVRR 6 months of CRT (sensitivity = 0.89 and specificity = 0.73). After construction of receiver operator characteristic (ROC) curves, an equation was established to predict LVRR: LVRR =-0.4LVDD (mm) + 0.5LVEF (%) + 1.1SDI (%), with responders presenting values >0 (sensitivity = 0.67 and specificity = 0.87). Conclusions: In this study, there was no strong correlation between TDI and SDI. An equation is proposed for the prediction of LVRR after CRT. Although larger trials are needed to validate these findings, this equation may be useful to candidates for CRT. (Echocardiography 2012;29:678-687)