Reaction-diffusion wave fronts: Multigeneration biological species under climate change

A generalization of reaction-diffusion models to multigeneration biological species is presented. It is based on more complex random walks than those in previous approaches. The new model is developed analytically up to infinite order. Our predictions for the speed agree to experimental data for several butterfly species better than existing models. The predicted dependence for the speed on the number of generations per year allows us to explain the change in speed observed for a specific invasion

Speed of reaction-diffusion fronts in spatially heterogeneous media

The front speed problem for nonuniform reaction rate and diffusion coefficient is studied by using singular perturbation analysis, the geometric approach of Hamilton-Jacobi dynamics, and the local speed approach. Exact and perturbed expressions for the front speed are obtained in the limit of large times. For linear and fractal heterogeneities, the analytic results have been compared with numerical results exhibiting a good agreement. Finally we reach a general expression for the speed of the front in the case of smooth and weak heterogeneities

Speed of reaction-transport processes

We present an approach to determining the speed of wave-front solutions to reaction-transport processes. This method is more accurate than previous ones. This is explicitly shown for several cases of practical interest: (i) the anomalous diffusion reaction, (ii) reaction diffusion in an advective field, and (iii) time-delayed reaction diffusion. There is good agreement with the results of numerical simulations

Speed of wave-front solutions to hyperbolic reaction-diffusion equations

The asymptotic speed problem of front solutions to hyperbolic reaction-diffusion (HRD) equations is studied in detail. We perform linear and variational analyses to obtain bounds for the speed. In contrast to what has been done in previous work, here we derive upper bounds in addition to lower ones in such a way that we can obtain improved bounds. For some functions it is possible to determine the speed without any uncertainty. This is also achieved for some systems of HRD (i.e., time-delayed Lotka-Volterra) equations that take into account the interaction among different species. An analytical analysis is performed for several systems of biological interest, and we find good agreement with the results of numerical simulations as well as with available observations for a system discussed recently

Time-delayed reaction-diffusion fronts

A time-delayed second-order approximation for the front speed in reaction-dispersion systems was obtained by Fort and Méndez [Phys. Rev. Lett. 82, 867 (1999)]. Here we show that taking proper care of the effect of the time delay on the reactive process yields a different evolution equation and, therefore, an alternate equation for the front speed. We apply the new equation to the Neolithic transition. For this application the new equation yields speeds about 10% slower than the previous one

Reaction-diffusion waves of advance in the transition to agricultural economics

In a previous paper [J.Fort and V.Méndez, Phys. Rev. Lett. 82, 867 (1999)], the possible importance of higher-order terms in a human population wave of advance has been studied. However, only a few such terms were considered. Here we develop a theory including all higher-order terms. Results are in good agreement with the experimental evidence involving the expansion of agriculture in Europe

Acute inflammatory reaction associated with endoluminal bypass grafts.

PURPOSE: Nonspecific inflammatory reactions characterized by local tenderness, fever, and flu-like discomfort have been seen in patients undergoing endoluminal graft placement in the abdominal aorta or the femoral arteries. We undertook a study to assess the clinical and laboratory parameters of this inflammation. METHODS: Ten patients with femoropopliteal artery (n = 9) or aortic (n = 1) lesions were treated with EndoPro System 1 stent-grafts made of nitinol alloy and covered with a polyester (Dacron) fabric. Eleven patients implanted with a bare nitinol stent served as the control group. RESULTS: In the stent-graft group, four patients showed clinical signs of acute inflammation manifested by fever and local tenderness. Three of these patients suffered thrombosis of the stent-grafts during the first month of follow-up. Plasma levels of interleukin-1 beta and interleukin-6 in all stent-graft patients were markedly increased 1 day after intervention (7.3 +/- 2.8 versus 90.2 +/- 34.1 pg/mL and 15.6 +/- 5.8 versus 175.5 +/- 66.3 pg/mL, respectively; p < 0.01). This was followed by an increase in fibrinogen (3.0 +/- 0.2 versus 5.0 +/- 0.2 g/L; p < 0.05) and C-reactive protein (14.6 +/- 3.3 versus 77.5 +/- 15.0 mg/L; p < 0.01) at 1 week. No direct correlation between the inflammatory markers and symptoms could be found. In vitro analysis showed that individual components of the stent-graft did not activate human neutrophils, whereas the intact stent-graft itself induced a marked neutrophil activation. CONCLUSIONS: The component of the self-expanding stent-graft responsible for the nonspecific inflammatory reaction was not identified in this study. It is likely that the stent-graft itself or some as yet unrecognized element of the device other than the Dacron fabric or metal alloy may be a potent in vivo inducer of cytokine reaction by neutrophils.

Covalent assembly of a soluble T cell receptor-peptide-major histocompatibility class I complex.

We used stepwise photochemical cross-linking for specifically assembling soluble and covalent complexes made of a T-cell antigen receptor (TCR) and a class I molecule of the major histocompatibility complex (MHC) bound to an antigenic peptide. For that purpose, we have produced in myeloma cells a single-chain Fv construct of a TCR specific for a photoreactive H-2Kd-peptide complex. Photochemical cross-linking of this TCR single-chain Fv with a soluble form of the photoreactive H-2Kd-peptide ligand resulted in the formation of a ternary covalent complex. We have characterized the soluble ternary complex and showed that it reacted with antibodies specific for epitopes located either on the native TCR or on the Kd molecules. By preventing the fast dissociation kinetics observed with most T cell receptors, this approach provides a means of preparing soluble TCR-peptide-MHC complexes on large-scale levels.

Establishment of animal models to analyze the kinetics and distribution of human tumor antigen-specific CD8⁺ T cells.

Many patients develop tumor antigen-specific T cell responses detectable in peripheral blood mononuclear cells (PBMCs) following cancer vaccine. However, measurable tumor regression is observed in a limited number of patients receiving cancer vaccines. There is a need to re-evaluate systemically the immune responses induced by cancer vaccines. Here, we established animal models targeting two human cancer/testis antigens, NY-ESO-1 and MAGE-A4. Cytotoxic T lymphocyte (CTL) epitopes of these antigens were investigated by immunizing BALB/c mice with plasmids encoding the entire sequences of NY-ESO-1 or MAGE-A4. CD8(+) T cells specific for NY-ESO-1 or MAGE-A4 were able to be detected by ELISPOT assays using antigen presenting cells pulsed with overlapping peptides covering the whole protein, indicating the high immunogenicity of these antigens in mice. Truncation of these peptides revealed that NY-ESO-1-specific CD8(+) T cells recognized D(d)-restricted 8mer peptides, NY-ESO-181-88. MAGE-A4-specific CD8(+) T cells recognized D(d)-restricted 9mer peptides, MAGE-A4265-273. MHC/peptide tetramers allowed us to analyze the kinetics and distribution of the antigen-specific immune responses, and we found that stronger antigen-specific CD8(+) T cell responses were required for more effective anti-tumor activity. Taken together, these animal models are valuable for evaluation of immune responses and optimization of the efficacy of cancer vaccines.

Simulation methods with extended stability for stiff biochemical kinetics

Background: With increasing computer power, simulating the dynamics of complex systems in chemistry and biology is becoming increasingly routine. The modelling of individual reactions in (bio)chemical systems involves a large number of random events that can be simulated by the stochastic simulation algorithm (SSA). The key quantity is the step size, or waiting time, τ, whose value inversely depends on the size of the propensities of the different channel reactions and which needs to be re-evaluated after every firing event. Such a discrete event simulation may be extremely expensive, in particular for stiff systems where τ can be very short due to the fast kinetics of some of the channel reactions. Several alternative methods have been put forward to increase the integration step size. The so-called τ-leap approach takes a larger step size by allowing all the reactions to fire, from a Poisson or Binomial distribution, within that step. Although the expected value for the different species in the reactive system is maintained with respect to more precise methods, the variance at steady state can suffer from large errors as τ grows. Results: In this paper we extend Poisson τ-leap methods to a general class of Runge-Kutta (RK) τ-leap methods. We show that with the proper selection of the coefficients, the variance of the extended τ-leap can be well-behaved, leading to significantly larger step sizes.Conclusions: The benefit of adapting the extended method to the use of RK frameworks is clear in terms of speed of calculation, as the number of evaluations of the Poisson distribution is still one set per time step, as in the original τ-leap method. The approach paves the way to explore new multiscale methods to simulate (bio)chemical systems.

Propagation velocity kinetics and repolarization alternans in a free-behaving sheep model of pacing-induced atrial fibrillation.

AIMS: Experimental models have reported conflicting results regarding the role of dispersion of repolarization in promoting atrial fibrillation (AF). Repolarization alternans, a beat-to-beat alternation in action potential duration, enhances dispersion of repolarization when propagation velocity is involved. METHODS AND RESULTS: In this work, original electrophysiological parameters were analysed to study AF susceptibility in a chronic sheep model of pacing-induced AF. Two pacemakers were implanted, each with a single right atrial lead. Right atrial depolarization and repolarization waves were documented at 2-week intervals. A significant and gradual decrease in the propagation velocity at all pacing rates and in the right atrial effective refractory period (ERP) was observed during the weeks of burst pacing before sustained AF developed when compared with baseline conditions. Right atrial repolarization alternans was observed, but because of the development of 2/1 atrioventricular block with far-field ventricular interference, its threshold could not be precisely measured. Non-sustained AF was not observed at baseline, but appeared during the electrical remodelling in association with a decrease in both ERP and propagation velocity. CONCLUSION: We report here on the feasibility of measuring ERP, atrial repolarization alternans, and propagation velocity kinetics and their potential in predicting susceptibility to AF in a free-behaving model of pacing-induced AF using the standard pacemaker technology.

Reverse transcriptase-dependent and -independent phases of infection with mouse mammary tumor virus: implications for superantigen function.

Mouse mammary tumor virus (MMTV) encodes a superantigen (SAg) that promotes stable infection and virus transmission. Upon subcutaneous MMTV injection, infected B cells present SAg to SAg-reactive T cells leading to a strong local immune response in the draining lymph node (LN) that peaks after 6 d. We have used the reverse transcriptase inhibitor 3'-azido-3'-deoxythymidine (AZT) to dissect in more detail the mechanism of SAg-dependent enhancement of MMTV infection in this system. Our data show that no detectable B or T cell response to SAg occurs in AZT pretreated mice. However, if AZT treatment is delayed 1-2 d after MMTV injection, a normal SAg-dependent local immune response is observed on day 6. Quantitation of viral DNA in draining LN of these infected mice indicates that a 4,000-fold increase in the absolute numbers of infected cells occurs between days 2 and 6 despite the presence of AZT. Furthermore MMTV DNA was found preferentially in surface IgG+ B cells of infected mice and was not detectable in SAg-reactive T cells. Collectively our data suggest that MMTV infection occurs preferentially in B cells without SAg involvement and is completed 1-2 d after virus challenge. Subsequent amplification of MMTV infection between days 2 and 6 requires SAg expression and occurs in the absence of any further requirement for reverse transcription. We therefore conclude that clonal expansion of infected B cells via cognate interaction with SAg-reactive T cells is the predominant mechanism for increasing the level of MMTV infection. Since infected B cells display a memory (surface IgG+) phenotype, both clonal expansion and possibly longevity of the virus carrier cells may contribute to stable MMTV infection.

The impact of fibrosis and steatosis on early viral kinetics in HCV genotype 1-infected patients treated with Peg-IFN-alfa-2a and ribavirin.

Summary.  Hepatitis C viral (HCV) kinetics after initiation of interferon-based therapy provide valuable insights for understanding virus pathogenesis, evaluating treatment antiviral effectiveness and predicting treatment outcome. Adverse effects of liver fibrosis and steatosis on sustained virological response have been frequently reported, yet their impacts on the early viral kinetics remain unclear. In this study, associations between histology status and early viral kinetics were assessed in 149 HCV genotype 1-infected patients treated with pegylated interferon alfa-2a and ribavirin (DITTO trial). In multivariate analyses adjusted for critical factors such as IL28B genotype and baseline viral load, presence of significant fibrosis (Ishak stage > 2) was found to independently reduce the odds of achieving an initial reduction (calculated from day 0 to day 4) in HCV RNA of ≥2 logIU/mL (adjusted OR 0.03, P = 0.004) but was not associated with the second-phase slope of viral decline (calculated from day 8 to day 29). On the contrary, presence of liver steatosis was an independent risk factor for not having a rapid second-phase slope, that is, ≥0.3 logIU/mL/week (adjusted OR 0.22, P = 0.012) but was not associated with the first-phase decline. Viral kinetic modelling theory suggests that significant fibrosis primarily impairs the treatment antiviral effectiveness in blocking viral production by infected cells, whereas the presence of steatosis is associated with a lower net loss of infected cells. Further studies will be necessary to identify the biological mechanisms underlain by these findings.

Oxygen uptake kinetics and middle distance swimming performance.

OBJECTIVE: The aim of this study was to determine whether V˙O(2) kinetics and specifically, the time constant of transitions from rest to heavy (τ(p)H) and severe (τ(p)S) exercise intensities, are related to middle distance swimming performance. DESIGN: Fourteen highly trained male swimmers (mean ± SD: 20.5 ± 3.0 yr; 75.4 ± 12.4 kg; 1.80 ± 0.07 m) performed an discontinuous incremental test, as well as square wave transitions for heavy and severe swimming intensities, to determine V˙O(2) kinetics parameters using two exponential functions. METHODS: All the tests involved front-crawl swimming with breath-by-breath analysis using the Aquatrainer swimming snorkel. Endurance performance was recorded as the time taken to complete a 400 m freestyle swim within an official competition (T400), one month from the date of the other tests. RESULTS: T400 (Mean ± SD) (251.4 ± 12.4 s) was significantly correlated with τ(p)H (15.8 ± 4.8s; r=0.62; p=0.02) and τ(p)S (15.8 ± 4.7s; r=0.61; p=0.02). The best single predictor of 400 m freestyle time, out of the variables that were assessed, was the velocity at V˙O(2max)vV˙O(2max), which accounted for 80% of the variation in performance between swimmers. However, τ(p)H and V˙O(2max) were also found to influence the prediction of T400 when they were included in a regression model that involved respiratory parameters only. CONCLUSIONS: Faster kinetics during the primary phase of the V˙O(2) response is associated with better performance during middle-distance swimming. However, vV˙O(2max) appears to be a better predictor of T400.

3

A Knudsen flow reactor has been used to quantify functional groups on the surface of seven different types of combustion particle samples: 3 amorphous carbons (FS 101, Printex 60, FW 2), 2 flame soots (hexane soot generated from a rich and a lean diffusion flame), and 2 Diesel particles (SRM 2975, Diesel soot recovered from a Diesel particulate filter). The technique is based on a heterogeneous titration reaction between a probe gas and a specific functional group on the particle surface. Six probe gases have been selected for the quantification of important functional groups: N(CH3)3 for the titration of acidic sites, NH2OH for carbonyl functions of aldehydes and ketones, CF3COOH and HCl for basic sites of different strength, O3 and NO2 for oxidizable groups. The limit of detection was generally well below 1% of a formal monolayer of adsorbed probe gas. Results obtained with N(CH3)3 were higher for the FW 2 amorphous carbon (post-oxidized sample, according to the manufacturer) and the Diesel particles (between 5.2·10 13 and 5.8·10 13 molecule/cm2), indicating a higher state of oxidation than for the other samples (between 1.3·10 12 and 3.7·10 12 molecule/cm2). The ratio of uptakes of CF3COOH and HCl inferred the presence of basic oxides on the particle surface, owing to the larger stability of the acetate compared to the chloride counter ion in the resulting pyrylium salt. The reactivity of the FS 101 amorphous carbon (3.7·10 15 molecule/cm2) and the hexane flame soot (between 1.9·10 15 and 2.7·10 15 molecule/cm2) towards O3 was very high, indicating the presence of a huge amount of oxidizable or reduced groups on the surface of these samples. Besides the quantification of surface functional groups, the kinetics of reactions between particles and probe gases has also been studied. The uptake coefficient γ0 was roughly correlated with the amount of probe gas taken up by the samples. Indeed, the presence of a high density of functional groups led to fast uptake of the probe gas. These different findings indicate that the particle surface appeared multi-functional, with the simultaneous presence of antagonistic functional groups which do not undergo internal chemical reactions, such as acid-base neutralization. Results also point to important differences in the surface reactivity of the samples, depending on the combustion conditions. The relative distribution of the surface functional groups may be a useful indicator for the state of oxidation and the reactivity of the particle surface.