Emergence and loss of assortative mating in sympatric speciation

We have studied an agent model which presents the emergence of sexual barriers through the onset of assortative mating, a condition that might lead to sympatric speciation. In the model, individuals are characterized by two traits, each determined by a single locus A or B. Heterozygotes on A are penalized by introducing an adaptive difference from homozygotes. Two niches are available. Each A homozygote is adapted to one of the niches. The second trait, called the marker trait has no bearing on the fitness. The model includes mating preferences, which are inherited from the mother and subject to random variations. A parameter controlling recombination probabilities of the two loci is also introduced. We study the phase diagram by means of simulations, in the space of parameters (adaptive difference, carrying capacity, recombination probability). Three phases are found, characterized by (i) assortative mating, (ii) extinction of one of the A alleles and (iii) Hardy-Weinberg like equilibrium. We also make perturbations of these phases to see how robust they are. Assortative mating can be gained or lost with changes that present hysteresis loops, showing the resulting equilibrium to have partial memory of the initial state and that the process of going from a polymorphic panmictic phase to a phase where assortative mating acts as sexual barrier can be described as a first-order transition. (C) 2009 Published by Elsevier Ltd.

Self-dual hopfions

We construct static and time-dependent exact soliton solutions with nontrivial Hopf topological charge for a field theory in 3 + 1 dimensions with the target space being the two dimensional sphere S(2). The model considered is a reduction of the so-called extended Skyrme-Faddeev theory by the removal of the quadratic term in derivatives of the fields. The solutions are constructed using an ansatz based on the conformal and target space symmetries. The solutions are said self-dual because they solve first order differential equations which together with some conditions on the coupling constants, imply the second order equations of motion. The solutions belong to a sub-sector of the theory with an infinite number of local conserved currents. The equation for the profile function of the ansatz corresponds to the Bogomolny equation for the sine-Gordon model.

Phototransformation of hematoporphyrin in aqueous solution: Anomalous behavior at low oxygen concentration

The photoactivation of a photosensitizer is the initial step in photodynamic therapy (PDT) where photochemical reactions result in the production of reactive oxygen species and eventually cell death. In addition to oxidizing biomolecules, some of these photochemical reactions lead to photosensitizer degradation at a rate dependent on the oxygen concentration among other factors. We investigated photodegradation of Photogem A (R) (28 mu M), a hematoporphyrin derivative, at different oxygen concentrations (9.4 to 625.0 mu M) in aqueous solution. The degradation was monitored by fluorescence spectroscopy. The degradation rate (M/s) increases as the oxygen concentration increases when the molar ratio of oxygen to PhotogemA (R) is greater than 1. At lower oxygen concentrations (< 25 mu M) an inversion of this behavior was observed. The data do not fit a simple kinetic model of first-order dependence on oxygen concentration. This inversion of the degradation rate at low oxygen concentration has not previously been demonstrated and highlights the relationship between photosensitizer and oxygen concentrations in determining the photobleaching mechanism(s). The findings demonstrate that current models for photobleaching are insufficient to explain completely the effects at low oxygen concentration.

The concept of quasi-integrability: a concrete example

We use the deformed sine-Gordon models recently presented by Bazeia et al [1] to take the first steps towards defining the concept of quasi-integrability. We consider one such definition and use it to calculate an infinite number of quasi-conserved quantities through a modification of the usual techniques of integrable field theories. Performing an expansion around the sine-Gordon theory we are able to evaluate the charges and the anomalies of their conservation laws in a perturbative power series in a small parameter which describes the ""closeness"" to the integrable sine-Gordon model. We show that in the case of the two-soliton scattering the charges, up to first order of perturbation, are conserved asymptotically, i.e. their values are the same in the distant past and future, when the solitons are well separated. We indicate that this property may hold or not to higher orders depending on the behavior of the two-soliton solution under a special parity transformation. For closely bound systems, such as breather-like field configurations, the situation however is more complex and perhaps the anomalies have a different structure implying that the concept of quasi-integrability does not apply in the same way as in the scattering of solitons. We back up our results with the data of many numerical simulations which also demonstrate the existence of long lived breather-like and wobble-like states in these models.

Phase diagram and spectral properties of a new exactly integrable spin-1 quantum chain

The spectral properties and phase diagram of the exactly integrable spin-1 quantum chain introduced by Alcaraz and Bariev are presented. The model has a U(1) symmetry and its integrability is associated with an unknown R-matrix whose dependence on the spectral parameters is not of a different form. The associated Bethe ansatz equations that fix the eigenspectra are distinct from those associated with other known integrable spin models. The model has a free parameter t(p). We show that at the special point t(p) = 1, the model acquires an extra U(1) symmetry and reduces to the deformed SU(3) Perk-Schultz model at a special value of its anisotropy q = exp(i2 pi/3) and in the presence of an external magnetic field. Our analysis is carried out either by solving the associated Bethe ansatz equations or by direct diagonalization of the quantum Hamiltonian for small lattice sizes. The phase diagram is calculated by exploring the consequences of conformal invariance on the finite-size corrections of the Hamiltonian eigenspectrum. The model exhibits a critical phase ruled by the c = 1 conformal field theory separated from a massive phase by first-order phase transitions.

Recording from Two Neurons: Second-Order Stimulus Reconstruction from Spike Trains and Population Coding

We study the reconstruction of visual stimuli from spike trains, representing the reconstructed stimulus by a Volterra series up to second order. We illustrate this procedure in a prominent example of spiking neurons, recording simultaneously from the two H1 neurons located in the lobula plate of the fly Chrysomya megacephala. The fly views two types of stimuli, corresponding to rotational and translational displacements. Second-order reconstructions require the manipulation of potentially very large matrices, which obstructs the use of this approach when there are many neurons. We avoid the computation and inversion of these matrices using a convenient set of basis functions to expand our variables in. This requires approximating the spike train four-point functions by combinations of two-point functions similar to relations, which would be true for gaussian stochastic processes. In our test case, this approximation does not reduce the quality of the reconstruction. The overall contribution to stimulus reconstruction of the second-order kernels, measured by the mean squared error, is only about 5% of the first-order contribution. Yet at specific stimulus-dependent instants, the addition of second-order kernels represents up to 100% improvement, but only for rotational stimuli. We present a perturbative scheme to facilitate the application of our method to weakly correlated neurons.

Second-order negative-curvature methods for box-constrained and general constrained optimization

A Nonlinear Programming algorithm that converges to second-order stationary points is introduced in this paper. The main tool is a second-order negative-curvature method for box-constrained minimization of a certain class of functions that do not possess continuous second derivatives. This method is used to define an Augmented Lagrangian algorithm of PHR (Powell-Hestenes-Rockafellar) type. Convergence proofs under weak constraint qualifications are given. Numerical examples showing that the new method converges to second-order stationary points in situations in which first-order methods fail are exhibited.

Density-Profile Processes Describing Biological Signaling Networks: Almost Sure Convergence to Deterministic Trajectories

We introduce jump processes in R(k), called density-profile processes, to model biological signaling networks. Our modeling setup describes the macroscopic evolution of a finite-size spin-flip model with k types of spins with arbitrary number of internal states interacting through a non-reversible stochastic dynamics. We are mostly interested on the multi-dimensional empirical-magnetization vector in the thermodynamic limit, and prove that, within arbitrary finite time-intervals, its path converges almost surely to a deterministic trajectory determined by a first-order (non-linear) differential equation with explicit bounds on the distance between the stochastic and deterministic trajectories. As parameters of the spin-flip dynamics change, the associated dynamical system may go through bifurcations, associated to phase transitions in the statistical mechanical setting. We present a simple example of spin-flip stochastic model, associated to a synthetic biology model known as repressilator, which leads to a dynamical system with Hopf and pitchfork bifurcations. Depending on the parameter values, the magnetization random path can either converge to a unique stable fixed point, converge to one of a pair of stable fixed points, or asymptotically evolve close to a deterministic orbit in Rk. We also discuss a simple signaling pathway related to cancer research, called p53 module.

INSECT CHYMOTRYPSINS: CHLOROMETHYL KETONE INACTIVATION AND SUBSTRATE SPECIFICITY RELATIVE TO POSSIBLE COEVOLUTIONAL ADAPTATION OF INSECTS AND PLANTS

Insect digestive chymotrypsins are present in a large variety of insect orders but their substrate specificity still remains unclear. Ewer insect chymotrypsins from 3 different insect orders (Dictyoptera, Coleoptera and two Lepidoptera) were isolated using affinity chromatography. Enzymes presented molecular masses in the range of 20 to 31 kDa and pH optima in the range of 7.5 to 10.0. Kinetic characterization. using different, colorimetric and fluorescent substrates indicated that insect chymotrypsins differ from, bovine chymotrypsin in their primary specificity toward small substrates (like N-benzoyl-L-Tyr p-nitroanilide) rather than on their preference for large substrates (exemplified by Succynil-Ala-Ala-Pro-Phe P-nitroanilide). Chloromethyl ketones (TPCK, N-alpha-tosyl-L-Phe chloromethyl ketone and Z-GGF-CK, N-carbobenzoxy-Gly-Gly-phe-CK) inactivated all chymotrypsins legated. Inactivation rates follow apparent first-order kinetics with variable second order rates (TPCK, 42 to 130 M(-1)s(-1); Z-GGF-CK, 150 to 450 M(-1)s(-1) that may be remarkably low for S. frugiperda chymotrypsin (TPCK, 6 M(-1)s(-1); Z-GGF-CK, 6.1 M(-1) s(-1)). Homology modelling and sequence alignment showed that. in lepidopteran chymotrypsins, differences in the amino acid residues in the neighborhood of the catalytic His 57 may affect its pKa, value. This is Proposed as the cause of the decrease in His 57 reactivity toward chloromethyl ketones. Such amino acid replacement in the active site is proposed. to be an adaptation to the presence of dietary ketones. (C) 2009 Wiley Periodicals, Inc.

Heterogeneous photocatalytic degradation of acid dyes in aqueous TiO2 suspensions: Kinetics and toxicity

This work assesses the photocatalytic (TiO2/UV) degradation of a simulated acid dye bath (Yellow 3, Red 51, Blue 74, and auxiliary chemicals). Color and phytotoxicity removal were monitored by spectrophotometry and lettuce (Lactuca sativa) seeds as the test organism, respectively. Mineralization was determined by DOC analyses. Photocatalytic, photolytic, and adsorption experiments were performed, showing that adsorption was negligible. After 240 minutes of irradiation, it was achieved 96% and 78% of color removal with photocatalysis and photolysis, respectively. 37% of mineralization occurred with photocatalysis only. The dye bath was rendered completely non-toxic after 60 minutes of photocatalytic treatment; the same result was only achieved with photolysis after 90 minutes. A kinetic model composed of two first-order in series reactions was used. The first photocatalytic decolorization rate constant was k(1) = 0.062 min(-1) and the second k(2) = 0.0043 min(-1), approximately two times greater than the photolytic ones.

Heterogeneous photocatalytic degradation of reactive dyes in aqueous TiO(2) suspensions: Decolorization kinetics

This work assesses the photocatalytic (TiO(2)/UV) degradation of a simulated reactive dye bath (Black 5, Red 239, Yellow 17, and auxiliary chemicals). Color removal was monitored by spectrophotometry. Mineralization was determined by DOC analyses. Photocatalytic, photolytic, and adsorption experiments were performed, showing that adsorption was negligible. After 30 min of irradiation, it was achieved 97% and 40% of color removal with photocatalysis and photolysis, respectively. No mineralization occurred within 30 min. A kinetic model composed of two, first-order in-series reactions was used. The first photocatalytic decolorization rate constant was k(1) = 2.6 min(-1) and the second k(2) = 0.011 min(-1). The fast decolorization of Reactive Black 5 dye is an indication that the number of azo and vinylsulfone groups in the dye molecule maybe a determining factor for the increased photolytic and photocatalytic color removal and degradation rates. (C) 2008 Elsevier B.V. All rights reserved.