Trypanosoma cruzi: blood parasitism kinetics and their correlation with heart parasitism intensity during long-term infection of Beagle dogs

The goals of the present study were to evaluate the kinetics of blood parasitism by examination of fresh blood, blood culture (BC) and PCR assays and their correlation with heart parasitism during two years of infection in Beagle dogs inoculated with the Be-78, Y and ABC Trypanosoma cruzi strains. Our results showed that the parasite or its kDNA is easily detected during the acute phase in all infected animals. On the other hand, a reduced number of positive tests were verified during the chronic phase of the infection. The frequency of positive tests was correlated with T. cruzi strain. The percentage of positive BC and blood PCR performed in samples from animals inoculated with Be-78 and ABC strains were similar and significantly larger in relation to animals infected with the Y strain.Comparison of the positivity of PCR tests performed using blood and heart tissue samples obtained two years after infection showed two different patterns associated with the inoculated T. cruzi strain: (1) high PCR positivity for both blood and tissue was observed in animals infected with Be-78 or ABC strains; (2) lower and higher PCR positivity for the blood and tissue, respectively, was detected in animals infected with Y strains. These data suggest that the sensitivity of BC and blood PCR was T. cruzi strain dependent and, in contrast, the heart tissue PCR revealed higher sensitivity regardless of the parasite stock.

Urokinase receptor (uPAR, CD87) is a platelet receptor important for kinetics and TNF-induced endothelial adhesion in mice.

BACKGROUND: Urokinase plasminogen activator receptor (uPAR, CD87) is a widely distributed 55-kD, glycoprotein I-anchored surface receptor. On binding of its ligand uPA, it is known to increase leukocyte adhesion and traffic. Using genetically deficient mice, we explored the role of uPAR in platelet kinetics and TNF-induced platelet consumption. METHODS AND RESULTS: Anti-uPAR antibody stained platelets from normal (+/+) but not from uPAR-/- mice, as seen by fluorescence-activated cell sorter analysis. 51Cr-labeled platelets from uPAR-/- donors survived longer than those from +/+ donors when injected into a +/+ recipient. Intratracheal TNF injection induced thrombocytopenia and a platelet pulmonary localization, pronounced in +/+ but absent in uPAR-/- mice. Aprotinin, a plasmin inhibitor, decreased TNF-induced thrombocytopenia. TNF injection markedly reduced the survival and increased the pulmonary localization of 51Cr-labeled platelets from +/+ but not from uPAR-/- donors, indicating that it is the platelet uPAR that is critical for their response to TNF. As seen by electron microscopy, TNF injection increased the number of platelets and polymorphonuclear neutrophils (PMNs) in the alveolar capillaries of +/+ mice, whereas in uPAR-/- mice, platelet trapping was insignificant and PMN trapping was slightly reduced. Platelets within alveolar capillaries of TNF-injected mice were activated, as judged from their shape, and this was evident in +/+ but not in uPAR-/- mice. CONCLUSIONS: These results demonstrate for the first time the critical role of platelet uPAR for kinetics as well as for activation and endothelium adhesion associated with inflammation.

Infection kinetics of human adenovirus serotype 41 in HEK 293 cells

The purpose of this work was to acquire an overview of the infectious cycle of HAdV-41 in permissive HEK 293 cells and compare it to that observed with the prototype of the genus, Human adenovirus C HAdV-2. HEK 293 cells were infected with each virus separately and were harvested every 12 h for seven days. Infection kinetics were analysed using confocal and electronic microscopy. The results show that, when properly cultivated, HAdV-41 was not fastidious. It had a longer multiplication cycle, which resulted in the release of complete viral particles and viral stocks reached high titres. After 60 h of infection, the export of viral proteins from the infected cell to the extracellular milieu was observed, with a pattern similar to that previously described for HAdV-2 penton-base trafficking after 30 h of infection. HAdV-41 had a non-lytic cycle and the infection spread from the first infected cell to its neighbours. The release process of the viral particles is unknown. The results observed for HAdV-41 infection in HEK 293 cells show how different this virus is from the prototype HAdV-2 and provides information for the development of this vector for use in gene therapy.

Growth, cysts and kinetics of Borrelia garinii (Spirochaetales: Spirochaetacea) in different culture media

The aim of the present paper was to evaluate cyst formation and growth parameters of Borrelia garinii in a range of media differing in formulation and cost. A qualitative assessment of morphology and motility of B. garinii was conducted. All media were prepared aseptically and used in test tubes or Petri dishes. For each medium, the initial spirochete concentration was standardized to 10³ spirochets/mL. The following culture media were suitable to grow B. garinii: Barbour-Stoenner-Kelly, brain heart infusion and PMR. Growth was minimal at six weeks post-inoculation and maximum spirochete density was observed between 9-12 weeks. Often, the cultures developed cysts of different sizes, isolated or in groups, with a spiraled portion of variable sizes, mainly in unfavorable culture media. Brazilian Lyme disease-like illness, also known as Baggio-Yoshinari syndrome (BYS), is a new and interesting emerging tick-borne disease, caused by Borrelia burgdorferi sensu lato spirochetes, only during its cystic forms. It has been assumed that the peculiar clinical and laboratory features of BYS are consequential to the absence of a human sucker Ixodes ricinus complex tick at risk areas in Brazil, supporting the concept that the borrelia phenotypic expression pattern is modified as it is transmitted through the host.

Steady-state brain glucose transport kinetics re-evaluated with a four-state conformational model.

Glucose supply from blood to brain occurs through facilitative transporter proteins. A near linear relation between brain and plasma glucose has been experimentally determined and described by a reversible model of enzyme kinetics. A conformational four-state exchange model accounting for trans-acceleration and asymmetry of the carrier was included in a recently developed multi-compartmental model of glucose transport. Based on this model, we demonstrate that brain glucose (G(brain)) as function of plasma glucose (G(plasma)) can be described by a single analytical equation namely comprising three kinetic compartments: blood, endothelial cells and brain. Transport was described by four parameters: apparent half saturation constant K(t), apparent maximum rate constant T(max), glucose consumption rate CMR(glc), and the iso-inhibition constant K(ii) that suggests G(brain) as inhibitor of the isomerisation of the unloaded carrier. Previous published data, where G(brain) was quantified as a function of plasma glucose by either biochemical methods or NMR spectroscopy, were used to determine the aforementioned kinetic parameters. Glucose transport was characterized by K(t) ranging from 1.5 to 3.5 mM, T(max)/CMR(glc) from 4.6 to 5.6, and K(ii) from 51 to 149 mM. It was noteworthy that K(t) was on the order of a few mM, as previously determined from the reversible model. The conformational four-state exchange model of glucose transport into the brain includes both efflux and transport inhibition by G(brain), predicting that G(brain) eventually approaches a maximum concentration. However, since K(ii) largely exceeds G(plasma), iso-inhibition is unlikely to be of substantial importance for plasma glucose below 25 mM. As a consequence, the reversible model can account for most experimental observations under euglycaemia and moderate cases of hypo- and hyperglycaemia.

Modulation of expression and activity of cytochrome P450s and alteration of praziquantel kinetics during murine schistosomiasis

In this study, we investigated the expression and activity of liver cytochrome P450s (CYPs) and praziquantel (PZQ) kinetics in mice infected with Schistosoma mansoni. Swiss Webster (SW) mice of both genders were infected (100 cercariae) on postnatal day 10 and killed on post-infection days (PIDs) 30 or 55. Non-infected mice of the same age and sex served as controls. Regardless of mouse sex, infection depressed the activities of CYP1A [ethoxy/methoxy-resorufin-O-dealkylases (EROD/MROD)], 2B9/10 [pentoxy/benzyloxy-resorufin-O-dealkylases (PROD, BROD)], 2E1 [p-nitrophenol-hydroxylase (PNPH)] and 3A11 [erythromycin N-demethylase (END)] on PID 55 but not on PID 30. On PID 55, infection decreased liver CYP mRNA levels (real-time reverse transcription-polymerase chain reaction). On PID 30, whereas mRNA levels remained unaltered in males, they were depressed in females. Plasma PZQ (200 and 400 mg/kg body weight intraperitoneally) levels were measured (high-performance liquid chromatography) at different post-treatment intervals. In males and females, infection delayed the PZQ clearance on PID 55, but not on PID 30. Therefore, it can be concluded that schistosomiasis down-modulated CYP expression and activity and delayed PZQ clearance on PID 55, when a great number of parasite eggs were lodged in the liver. On PID 30, when egg-laying was initiated by the worms, no change of CYP expression and activity was found, except for a depression of CYP1A2 and 3A11 mRNAs in female mice.

A single nucleotide polymorphism, rs129679860, in the IL28B locus is associated with the viral kinetics and a sustained virological response in a chronic, monoinfected hepatitis C virus genotype-1 Brazilian population treated with pegylated interferon-ribavirin

Single nucleotide polymorphisms (SNPs) in the interleukin (IL)28B locus have been associated with a sustained virological response (SVR) in interferon-ribavirin (IFN-RBV)-treated chronic hepatitis C virus (HCV)-infected patients in European and African populations. In this study, the genotype frequency of two IL28B SNPs (rs129679860 and rs8099917) in a cohort of chronic HCV-monoinfected patients in Brazil was evaluated and the SNP sufficient to predict the treatment response outcome was determined. A total of 66 naïve genotype-1 chronic HCV-infected patients were genotyped and the associated viral kinetics and SVR were assessed. The overall SVR was 38%. Both the viral kinetics and SVR were associated with rs129679860 genotypes (CC = 62% vs. CT = 33% vs. TT = 18%, p = 0.016). However, rs8099917 genotypes were only associated with SVR (TT = 53% vs. TG = 33% vs. GG = 18%; p = 0.032). In this population, the analysis of a single SNP, rs12979860, successfully predicts SVR in the IFN-RBV treatment of HCV.

Mammalian genes are transcribed with widely different bursting kinetics.

In prokaryotes and eukaryotes, most genes appear to be transcribed during short periods called transcriptional bursts, interspersed by silent intervals. We describe how such bursts generate gene-specific temporal patterns of messenger RNA (mRNA) synthesis in mammalian cells. To monitor transcription at high temporal resolution, we established various gene trap cell lines and transgenic cell lines expressing a short-lived luciferase protein from an unstable mRNA, and recorded bioluminescence in real time in single cells. Mathematical modeling identified gene-specific on- and off-switching rates in transcriptional activity and mean numbers of mRNAs produced during the bursts. Transcriptional kinetics were markedly altered by cis-regulatory DNA elements. Our analysis demonstrated that bursting kinetics are highly gene-specific, reflecting refractory periods during which genes stay inactive for a certain time before switching on again.

Dissecting TCR-MHC/peptide complex interactions with A2/peptide multimers incorporating tumor antigen peptide variants: crucial role of interaction kinetics on functional outcomes.

Soluble peptide/MHC-class-I (pMHC) multimers have recently emerged as unique reagents for the study of specific interactions between the pMHC complex and the TCR. Here, we assessed the relative binding efficiency of a panel of multimers incorporating single-alanine-substituted variants of the tumor-antigen-derived peptide MAGE-A10(254-262) to specific CTL clones displaying different functional avidity. For each individual clone, the efficiency of binding of multimers incorporating MAGE-A10 peptide variants was, in most cases, in good although not linear correlation with the avidity of recognition of the corresponding variant. In addition, we observed two types of discrepancies between efficiency of recognition and multimer binding. First, for some peptide variants, efficient multimer binding was detected in the absence of measurable effector functions. Some of these peptide variants displayed antagonist activity. Second, when comparing different clones we found clear discrepancies between the dose of peptide required to obtain half-maximal lysis in CTL assays and the binding efficiency of the corresponding multimers. These discrepancies, however, were resolved when the differential stability of the TCR/pMHC complexes was determined. For individual clones, decreased recognition correlated with increased TCR/pMHC off-rate. TCR/pMHC complexes formed by antagonist ligands displayed off-rates faster than those of TCR/pMHC complexes formed with weak agonists. In addition, when comparing different clones, the efficiency of multimer staining correlated better with relative multimer off-rates than with half-maximal lysis values. Altogether, the data presented here reconcile and extend our previous results on the impact of the kinetics of interaction of TCR with pMHC complexes on multimer binding and underline the crucial role of TCR/pMHC off-rates for the functional outcome of such interactions.

A non-invasive assessment of hepatic glycogen kinetics and post-absorptive gluconeogenesis in man.

A novel approach to the study of hepatic glycogen kinetics and fractional gluconeogenesis in vivo is described. Ten healthy female subjects were fed an iso-caloric diet containing 55% carbohydrate energy with a 13C abundance of 1.083 atom percent for a 3-day baseline period; then, a diet of similar composition, but providing carbohydrate with a 13C abundance of 1.093 atom percent was started and continued for 5 days. Resting respiratory gas exchanges, urinary nitrogen excretion, breath 13CO2 and plasma 13C glucose were measured every morning in the fasting state. The enrichment in 13C of hepatic glycogen was calculated from these measured data. 13C glycogen enrichment increased after switching to a 13C enriched carbohydrate diet, and was identical to the 13C enrichment of dietary carbohydrates after 3 days. The time required to renew 50% of hepatic glycogen, as determined from the kinetics of 13C glycogen enrichment, was 18.9 +/- 3.6 h. Fractional gluconeogenesis, as determined from the difference between the enrichments of glucose oxidized originating from hepatic glycogen and plasma glucose 13C was 50.8 +/- 5.3%. This non-invasive method will allow the study of hepatic glycogen metabolism in insulin-resistant patients.

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.

Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps

Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Résumé grand public Comme toute matière, la roche est sensible à son environnement et cherche à s'adapter pour acquérir un état stable (état d'équilibre). Les changements des conditions physiques (température et pression) vont ainsi impliquer des modifications dans la roche. Le métamorphisme est l'étude de ces changements. Les minéraux qui constituent la roche peuvent modifier, leur structure, leur chimie ou être remplacer par d'autres minéraux plus stables. Il est ainsi crucial de déterminer les processus responsables et limitant de la croissance minérale. Trois processus permettent la croissance ; (1) la dissolution des éléments du réactant, (2) le transport de ces éléments vers le site de croissance, (3) l'incorporation de ces éléments dans la nouvelle structure. Cette thèse se focalise sur les structures des minéraux de haute pression (forme, zonation chimique, structure interne) pour essayer de déterminer les facteurs importants à l'origine de leur état final. Les zones d'étude se situent dans la zone de Sésia. La première partie traite de la problématique liée à l'incorporation d'un élément dans une structure minérale. A l'image de la croissance humaine, les irrégularités minéralogiques permettent de mettre en lumière un dysfonctionnement de la croissance due à un excès ou à une carence d'un élément. Bien dosé, cet élément est cependant essentiel à la croissance. Les zoisites (épidotes) des métabasites de la région de Cima di Bonze montrent une zonation chimique en sablier. Dans cette zonation la teneur en fer excède la capacité maximum que peut contenir la structure orthorhombique de la zoisite. Des défauts de structure permettent l'accommodation de cet excès. La zoisite peut ainsi adapter sa structure pour permettre l'incorporation d'une relativement grande quantité de fer. Les études précédentes montraient, pour des conditions similaires, la formation de deux épidotes distinctes. La deuxième partie se penche sur la compétition entre le minéral qui fait sa croissance et les minéraux (réactants) qui l'entourent. Les métapélites de la région du Monte Mucrone contiennent des grenats atollaires. Des études détaillées de la texture et de la zonation chimique du grenat ainsi qu'une modélisation thermodynamique ont permis de mieux cerner les facteurs importants responsables de la forme atollaire. Cette structure est obtenue par un changement du comportement de la croissance du grenat le long d'un chemin P-T hercynien. Dans un premier stade, le grenat croît rapidement et consume peu le quartz de la matrice. La croissance se fait ainsi le long des jointures des grains de quartz. Dans un second temps, les changements de conditions PT donnent une croissance lente du grenat et une forte consommation du quartz. Le grenat peut ainsi développer sa forme dodécaédrale classique. La troisième partie s'intéresse aux distances de transport par diffusion d'un élément (ici l'argon) durant la haute pression. Pour ce faire, un profile d'âges 40Ar/39Ar sur biotite a été mesuré depuis un veine de haute pression riche en argon jusque dans son encaissant (granitoïd du Monte Mucrone). Le profile montre une répartition des âges suivant une courbe de diffusion. Le transport se fait sur une longueur de deux centimètre avec l'aide d'un fluide. Il est réduit à une échelle millimétrique quand la phase fluide disparaît. Cette étude montre ainsi les difficultés de transport des éléments durant la haute pression ne permettant pas un rééquilibrage de la roche à grande échelle. Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Résumé de thèse Les processus de croissance (diffusion des éléments et les réactions d'interface) et les conditions dans lesquelles les minéraux grandissent (température, pression, fluide, composition chimique de la roche), déterminent la texture ainsi que la zonation des minéraux. Cette thèse se focalise, par le biais de textures peu communes, sur trois différents processus impliqués dans la croissance minérale à haute pression (Zone de Sésia, Alpes de l'Ouest, Italie). L'incorporation d'un élément dans une structure minérale ne peut se faire que dans des sites en accord avec la taille et la charge ionique de l'élément. De plus, la balance de charge doit être maintenue dans le minéral. La régularité de la structure cristalline fixe ainsi une limite maximum de concentration d'un élément donné. Les zoisites provenant des métabasites de la région de Cima di Bonze montrent des zonations en sablier caractérisées par une concentration anormale en fer. La zonation se marque par une différente teinte de biréfringence et par un plus grand angle d'extinction que le reste de la zoisite. Une inter-croissance de clinozoisite à l'intérieur de la structure orthorhombique de la zoisite peut ainsi être suspectée. Les analyses XRD (diffraction des rayons x) ainsi que les analyses Raman ne confirment pas cette suspicion. Seules les analyses TEM (microscope à électrons transmis) montrent des défauts de structure pouvant être interprétés comme des modules de clinozoisite. Ils ne peuvent cependant pas être considérés comme une phase thermodynamique. Un nouveau trou d'immiscibilité entre deux zoisite (X ep= 0.1 and Xep = 0.15) a ainsi pu être établi. Dans les métapélites la région du Monte Mucrone, des grenats fortement zonés montrent une évolution texturale singulière. Ils présentent une forme initiale de `champignon' qui se développe pour former une structure atollaire finale. L'étude conjuguée de la structure 3D et des zonations, ainsi que l'établissement d'un model thermodynamique, indiquent que ces structures proviennent de deux épisodes de croissances : (1) La croissance du grenat durant un chemin prograde hercynien (de 525 °C et 6.2 kbar à 640 °C et 9 kbar) permet la formation des textures atollaires. Elles sont le résultat d'une croissance poecilitique initiale suivie d'une croissance idiomorphique du grenat. (2) La structure est rendue plus complexe par la cristallisation d'un grenat homogène tout autour ainsi qu'à l'intérieur du grenat hercynien durant la haute pression alpine (550 °C and 20 kbar). L'arrivée de l'eau durant la haute pression facilite le transport d'éléments et permet une cristallisation rapide du grenat. La diffusion peut être un facteur limitant de la croissance minéralogique. Elle a aussi une grande importance pour la géochronologie. Une veine de haute pression à l'intérieur du granitoïde du Monte Mucrone a été étudiée dans le but de déterminer la distance de diffusion de l'argon. Le profile d'âges 40Ar/39Ar sur biotites, établi de la veine vers le métagranitoïde, suit une courbe de diffusion. Les âges sont élevés proche de la veine (800 Ma) puis décroissent jusqu'à des âges homogènes (170-150 Ma) à deux centimètres de la veine. La présence de fluide, marqué par de hautes concentrations en chlore, permet une diffusion centimétrique. Cependant, la distance est réduite à une échelle millimétrique quand le fluide est absent. Les très faibles distances de diffusion préservent les âges pré-alpins et impliquent un événement géologique pour les âges de 170-150 Ma. Kinetics of crystal growth and equilibrium domains in eclogite of the Sesia Zone, Western Alps Darbellay Bastien Institut de Minéralogie et Géochimie Thesis abstract Rock textures and zonings are the consequence of growth processes (element diffusion and interface reaction) steered by the environment in which they grew (pressure, differential stress, temperature, fluid and rock composition). The thesis presented here focuses on three different topics, each of it dealing with aspects of mineral growth processes during subduction, in a high-pressure environment. All studies were conducted in the Sesia Zone of the Western European Alps, Italy. The first study addresses the crystallography and geochemistry of element incorporation in zoisite, one of the major hydrous minerals found in subduction zone rocks. Elements can be incorporated into a mineral structure only on crystallographic sites that offer enough space for the ion and the overall charge balance has to be maintained. Element concentrations are hence limited. Incorporation of some elements produces complex zoning, including hourglass like patterns, which are the focus of the first contribution. Zoisites from Cima di Bonze (Sesia Zone) show spectacular hourglass zoning defined by Fe-content variations. The hourglass zones have a distinct birefringence and a different extinction angle than the regular part of the zoisite. We show by detailed XRD (X-ray diffraction) and confocal Raman analyses that the high Fe-zones are nevertheless zoisite, and not clinozoisite as one might expect. High resolution TEM (transmission electron microscopy) analyses show planar defects on (100) that can be interpreted as small-scale clinozoisite modules. However, these clinozoisites cannot be interpreted as a distinctive thermodynamic phase and the entire mineral has to be considered as zoisite. The miscibility gap between two zoisites (Xep = 0.1 and Xep = 0.15) can be then definite at 550 ± 50°C and 14 to 20 Kbar. Strongly zoned garnets in quartz rich metapelite from the Monte Mucrone area (Sesia Zone) show evolution form 3D mushroom to atoll structure. The second contribution presents textural investigations, garnet zoning and thermodynamic modeling that demonstrate that atoll garnets are the result of two distinctive growth events. (1) Garnet atoll structure is already formed during a prograde Hercynian path from 525 °C and 6.2 kbar to 640 °C and 9 kbar. It results in an initial poikilitic growth followed by a final idiomorphic growth event. (2) Alpine HP garnet are homogenous (550 °C and 20 kbar) and grew around and also inside the Hercynian garnet. Lack of prograde Alpine garnet and fast growth of the HP garnet is explained by the absence of water during much of the prograde path. Water saturation was only observed towards the end, close towards the peak metamorphic conditions. Diffusion could be a limiting factor for crystal growth. It has also a great importance in geochronology. HP vein inside the metagranitoide of the Monte Mucrone (~300 Ma) was investigated to determine argon diffusion scales during high-pressure metamorphism. 40Ar/39Ar biotite ages profile from the vein toward the metagranodiorite show a diffusion curve: old ages (800 Ma) located close to the vein decrease until homogenous 170-150 Ma ages are obtained, two centimeter away from the vein. Centimeter-scale diffusion occurs with help of a fluid phase marked by high chlorine concentrations. Argon diffusion is reduced to a millimeter scale when free fluid is absent. Very short diffusion distance permits to preserve pre-Alpine ages. The 170-150 Ma ages are considered to be geologic meaningful, probably resulting from the extensional tectonics linked to opening of the Tethian ocean.

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.

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.