A DISCRETIZATION SCHEME FOR AN ONE-DIMENSIONAL REACTION-DIFFUSION EQUATION WITH DELAY AND ITS DYNAMICS

The goal of this paper is to present an approximation scheme for a reaction-diffusion equation with finite delay, which has been used as a model to study the evolution of a population with density distribution u, in such a way that the resulting finite dimensional ordinary differential system contains the same asymptotic dynamics as the reaction-diffusion equation.

Reaction-diffusion systems coupled at the boundary and the Morse-Smale property

We study an one-dimensional nonlinear reaction-diffusion system coupled on the boundary. Such system comes from modeling problems of temperature distribution on two bars of same length, jointed together, with different diffusion coefficients. We prove the transversality property of unstable and stable manifolds assuming all equilibrium points are hyperbolic. To this end, we write the system as an equation with noncontinuous diffusion coefficient. We then study the nonincreasing property of the number of zeros of a linearized nonautonomous equation as well as the Sturm-Liouville properties of the solutions of a linear elliptic problem. (C) 2008 Elsevier Inc. All rights reserved.

Generic hyperbolicity of stationary solutions for a reaction-diffusion system

In this paper, we study the generic hyperbolicity of equilibria of a reaction-diffusion system with respect to nonlinear terms in the set of C(2)-functions equipped with the Whitney Topology. To accomplish this, we combine Baire`s Lemma and the usual Transversality Theorem. (C) 2010 Elsevier Ltd. All rights reserved.

Reaction-diffusion stochastic lattice model for a predator-prey system

We have the purpose of analyzing the effect of explicit diffusion processes in a predator-prey stochastic lattice model. More precisely we wish to investigate the possible effects due to diffusion upon the thresholds of coexistence of species, i. e., the possible changes in the transition between the active state and the absorbing state devoid of predators. To accomplish this task we have performed time dependent simulations and dynamic mean-field approximations. Our results indicate that the diffusive process can enhance the species coexistence.

Dynamics in dumbbell domains III. Continuity of attractors

In this paper we conclude the analysis started in [J.M. Arrieta, AN Carvalho, G. Lozada-Cruz, Dynamics in dumbbell domains I. Continuity of the set of equilibria, J. Differential Equations 231 (2006) 551-597] and continued in [J.M. Arrieta, AN Carvalho, G. Lozada-Cruz, Dynamics in dumbbell domains II. The limiting problem, J. Differential Equations 247 (1) (2009) 174-202 (this issue)] concerning the behavior of the asymptotic dynamics of a dissipative reaction-diffusion equation in a dumbbell domain as the channel shrinks to a line segment. In [J.M. Arrieta, AN Carvalho. G. Lozada-Cruz, Dynamics in dumbbell domains I. Continuity of the set of equilibria, J. Differential Equations 231 (2006) 551-597], we have established an appropriate functional analytic framework to address this problem and we have shown the continuity of the set of equilibria. In [J.M. Arrieta, AN Carvalho, G. Lozada-Cruz. Dynamics in dumbbell domains II. The limiting problem, J. Differential Equations 247 (1) (2009) 174-202 (this issue)], we have analyzed the behavior of the limiting problem. In this paper we show that the attractors are Upper semicontinuous and, moreover, if all equilibria of the limiting problem are hyperbolic, then they are lower semicontinuous and therefore, continuous. The continuity is obtained in L(p) and H(1) norms. (C) 2008 Elsevier Inc. All rights reserved.

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.

ON THE SUPERCRITICALITY OF THE FIRST HOPF BIFURCATION IN A DELAY BOUNDARY PROBLEM

The goal of this paper is to analyze the character of the first Hopf bifurcation (subcritical versus supercritical) that appears in a one-dimensional reaction-diffusion equation with nonlinear boundary conditions of logistic type with delay. We showed in the previous work [Arrieta et al., 2010] that if the delay is small, the unique non-negative equilibrium solution is asymptotically stable. We also showed that, as the delay increases and crosses certain critical value, this equilibrium becomes unstable and undergoes a Hopf bifurcation. This bifurcation is the first one of a cascade occurring as the delay goes to infinity. The structure of this cascade will depend on the parameters appearing in the equation. In this paper, we show that the first bifurcation that occurs is supercritical, that is, when the parameter is bigger than the delay bifurcation value, stable periodic orbits branch off from the constant equilibrium.

Modeling the skin pattern of fishes

Complicated patterns showing various spatial scales have been obtained in the past by coupling Turing systems in such a way that the scales of the independent systems resonate. This produces superimposed patterns with different length scales. Here we propose a model consisting of two identical reaction-diffusion systems coupled together in such a way that one of them produces a simple Turing pattern of spots or stripes, and the other traveling wave fronts that eventually become stationary. The basic idea is to assume that one of the systems becomes fixed after some time and serves as a source of morphogens for the other system. This mechanism produces patterns very similar to the pigmentation patterns observed in different species of stingrays and other fishes. The biological mechanisms that support the realization of this model are discussed.

Dynamics in dumbbell domains II. The limiting problem

In this work we continue the analysis of the asymptotic dynamics of reaction-diffusion problems in a dumbbell domain started in [J.M. Arrieta, AN Carvalho, G. Lozada-Cruz, Dynamics in dumbbell domains I. Continuity of the set of equilibria, J. Differential Equations 231 (2) (2006) 551-597]. Here we study the limiting problem, that is, an evolution problem in a ""domain"" which consists of an open, bounded and smooth set Omega subset of R(N) with a curve R(0) attached to it. The evolution in both parts of the domain is governed by a parabolic equation. In Omega the evolution is independent of the evolution in R(0) whereas in R(0) the evolution depends on the evolution in Omega through the continuity condition of the solution at the junction points. We analyze in detail the linear elliptic and parabolic problem, the generation of linear and nonlinear semigroups, the existence and structure of attractors. (C) 2009 Elsevier Inc. All rights reserved.

Gene Expression Noise in Spatial Patterning: hunchback Promoter Structure Affects Noise Amplitude and Distribution in Drosophila Segmentation

Positional information in developing embryos is specified by spatial gradients of transcriptional regulators. One of the classic systems for studying this is the activation of the hunchback (hb) gene in early fruit fly (Drosophila) segmentation by the maternally-derived gradient of the Bicoid (Bcd) protein. Gene regulation is subject to intrinsic noise which can produce variable expression. This variability must be constrained in the highly reproducible and coordinated events of development. We identify means by which noise is controlled during gene expression by characterizing the dependence of hb mRNA and protein output noise on hb promoter structure and transcriptional dynamics. We use a stochastic model of the hb promoter in which the number and strength of Bcd and Hb (self-regulatory) binding sites can be varied. Model parameters are fit to data from WT embryos, the self-regulation mutant hb(14F), and lacZ reporter constructs using different portions of the hb promoter. We have corroborated model noise predictions experimentally. The results indicate that WT (self-regulatory) Hb output noise is predominantly dependent on the transcription and translation dynamics of its own expression, rather than on Bcd fluctuations. The constructs and mutant, which lack self-regulation, indicate that the multiple Bcd binding sites in the hb promoter (and their strengths) also play a role in buffering noise. The model is robust to the variation in Bcd binding site number across a number of fly species. This study identifies particular ways in which promoter structure and regulatory dynamics reduce hb output noise. Insofar as many of these are common features of genes (e. g. multiple regulatory sites, cooperativity, self-feedback), the current results contribute to the general understanding of the reproducibility and determinacy of spatial patterning in early development.

Sweet Basil (Ocimum basilicum) Extracts Obtained by Supercritical Fluid Extraction (SFE): Global Yields, Chemical Composition, Antioxidant Activity, and Estimation of the Cost of Manufacturing

In this work, supercritical technology was used to obtain extracts from Ocimum basilicum (sweet basil) with CO(2) and the cosolvent H(2)O at 1, 10, and 20% (w/w). The raw material was obtained from hydroponic cultivation. The extract`s global yield isotherms, chemical compositions, antioxidant activity, and cost of manufacturing were determined. The extraction assays were done for pressures of 10 to 30 MPa at 303 to 323 K. The identification of the compounds present in the extracts was made by GC-MS and ESI-MS. The antioxidant activity of extracts was determined using the coupled reaction of beta-carotene and linolenic acid. At 1% of cosolvent, the largest global yield was obtained at 10 MPa and 303 K (2%, dry basis-d.b.); at 10% of cosolvent the largest global yield was obtained at 10 and 15 MPa (11%, d.b.), and at 20% of cosolvent the largest global yield was detected at 30 MPa and 303 K (24%, d.b.). The main components identified in the extracts were eugenol, germacrene-D, epi-alpha-cadinol, malic acid, tartaric acid, ramnose, caffeic acid, quinic acid, kaempferol, caffeoylquinic acid, and kaempferol 3-O-glucoside. Sweet basil extracts exhibited high antioxidant activity compared to beta-carotene. Three types of SFE extracts from sweet basil were produced, for which the estimated cost of manufacturing (class 5 type) varied from US$ 47.96 to US$ 1,049.58 per kilogram of dry extract.

Lower semicontinuity of attractors for non-autonomous dynamical systems

This paper is concerned with the lower semicontinuity of attractors for semilinear non-autonomous differential equations in Banach spaces. We require the unperturbed attractor to be given as the union of unstable manifolds of time-dependent hyperbolic solutions, generalizing previous results valid only for gradient-like systems in which the hyperbolic solutions are equilibria. The tools employed are a study of the continuity of the local unstable manifolds of the hyperbolic solutions and results on the continuity of the exponential dichotomy of the linearization around each of these solutions.

Semilinear parabolic problems in thin domains with a highly oscillatory boundary

In this paper, we study the behavior of the solutions of nonlinear parabolic problems posed in a domain that degenerates into a line segment (thin domain) which has an oscillating boundary. We combine methods from linear homogenization theory for reticulated structures and from the theory on nonlinear dynamics of dissipative systems to obtain the limit problem for the elliptic and parabolic problems and analyze the convergence properties of the solutions and attractors of the evolutionary equations. (C) 2011 Elsevier Ltd. All rights reserved.

A discontinuous-Galerkin-based immersed boundary method

A numerical method to approximate partial differential equations on meshes that do not conform to the domain boundaries is introduced. The proposed method is conceptually simple and free of user-defined parameters. Starting with a conforming finite element mesh, the key ingredient is to switch those elements intersected by the Dirichlet boundary to a discontinuous-Galerkin approximation and impose the Dirichlet boundary conditions strongly. By virtue of relaxing the continuity constraint at those elements. boundary locking is avoided and optimal-order convergence is achieved. This is shown through numerical experiments in reaction-diffusion problems. Copyright (c) 2008 John Wiley & Sons, Ltd.

Radical acetylation of 2`-deoxyguanosine and L-histidine coupled to the reaction of diacetyl with peroxynitrite in aerated medium

Diacetyl, like other alpha-dicarbonyl compounds, is reportedly cytotoxic and genotoxic. A food and cigarette contaminant, it is related with alcohol hepatotoxicity and lung disease. Peroxynitrite is a potent oxidant formed in vivo by the diffusion-controlled reaction of the superoxide radical anion with nitric oxide, which is able to form adducts with carbon dioxide and carbonyl compounds. Here, we investigate the nucleophilic addition of peroxynitrite to diacetyl forming acetyl radicals, whose reaction with molecular oxygen leads to acetate. Peroxynitrite is shown to react with diacetyl in phosphate buffer (bell-shaped pH profile with maximum at 7.2) at a very high rate constant (k(2) = 1.0 X 10(4) M-1 s(-1)) when compared with monocarbonyl substrates (k(2) < 10(3) M-1 s(-1)). Phosphate ions (100-500 MM) do not affect the rate of spontaneous peroxynitrite decay, but the H2PO4- anion catalyzes the nucleophilic addition of the peroxynitrite anion to diacetyl. The intermediacy of acetyl radicals is suggested by a three-line spectrum (a(N) = a(H) = 0.83 mT) obtained by EPR spin trapping of the reaction mixture with 2-methyl-2-nitrosopropane. The peroxynitrite reaction is accompanied by concentration-dependent oxygen uptake. Stoichiometric amounts of acetate from millimolar amounts of peroxynitrite and diacetyl were obtained under nonlimiting conditions of dissolved oxygen. In the presence of either L-histidine or 2`-deoxyguanosine, the peroxynitrite/diacetyl system afforded the corresponding acetylated molecules identified by HPLC-MS"". These studies provide evidence that the peroxynitrite/diacetyl reaction yields acetyl radicals and raise the hypothesis that protein and DNA nonenzymatic acetylation may occur in cells and be implicated in aging and metabolic disorders in which oxygen and nitrogen reactive species are putatively involved.