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.

Fine-Tuning Translation Kinetics Selection as the Driving Force of Codon Usage Bias in the Hepatitis A Virus Capsid

Hepatitis A virus (HAV), the prototype of genus Hepatovirus, has several unique biological characteristics that distinguish it from other members of the Picornaviridae family. Among these, the need for an intact eIF4G factor for the initiation of translation results in an inability to shut down host protein synthesis by a mechanism similar to that of other picornaviruses. Consequently, HAV must inefficiently compete for the cellular translational machinery and this may explain its poor growth in cell culture. In this context of virus/cell competition, HAV has strategically adopted a naturally highly deoptimized codon usage with respect to that of its cellular host. With the aim to optimize its codon usage the virus was adapted to propagate in cells with impaired protein synthesis, in order to make tRNA pools more available for the virus. A significant loss of fitness was the immediate response to the adaptation process that was, however, later on recovered and more associated to a re-deoptimization rather than to an optimization of the codon usage specifically in the capsid coding region. These results exclude translation selection and instead suggest fine-tuning translation kinetics selection as the underlying mechanism of the codon usage bias in this specific genome region. Additionally, the results provide clear evidence of the Red Queen dynamics of evolution since the virus has very much evolved to re-adapt its codon usage to the environmental cellular changing conditions in order to recover the original fitness.

Comment on "Kinetics of spinodal decomposition in the ising model with vacancy diffusion"

In a recent paper [Phys. Rev. B 50, 3477 (1994)], P. Fratzl and O. Penrose present the results of the Monte Carlo simulation of the spinodal decomposition problem (phase separation) using the vacancy dynamics mechanism. They observe that the t1/3 growth regime is reached faster than when using the standard Kawasaki dynamics. In this Comment we provide a simple explanation for the phenomenon based on the role of interface diffusion, which they claim is irrelevant for the observed behavior.

Kinetics of martensitic transitions in Cu-Al-Mn under thermal cycling: Analysis at multiple length scales

In this paper we study the evolution of the kinetic features of the martensitic transition in a Cu-Al-Mn single crystal under thermal cycling. The use of several experimental techniques including optical microscopy, calorimetry, and acoustic emission, has enabled us to perform an analysis at multiple scales. In particular, we have focused on the analysis of avalanche events (associated with the nucleation and growth of martensitic domains), which occur during the transition. There are significant differences between the kinetics at large and small length scales. On the one hand, at small length scales, small avalanche events tend to sum to give new larger events in subsequent loops. On the other hand, at large length scales the large domains tend to split into smaller ones on thermal cycling. We suggest that such different behavior is the necessary ingredient that leads the system to the final critical state corresponding to a power-law distribution of avalanches.

Accuracy in the experimental calorimetric study of the crystallization kinetics and predicitve transformation diagrams: Application to a Ga-Te amorphous alloy

The uncertainties inherent to experimental differential scanning calorimetric data are evaluated. A new procedure is developed to perform the kinetic analysis of continuous heating calorimetric data when the heat capacity of the sample changes during the crystallization. The accuracy of isothermal calorimetric data is analyzed in terms of the peak-to-peak noise of the calorimetric signal and base line drift typical of differential scanning calorimetry equipment. Their influence in the evaluation of the kinetic parameters is discussed. An empirical construction of the time-temperature and temperature heating rate transformation diagrams, grounded on the kinetic parameters, is presented. The method is applied to the kinetic study of the primary crystallization of Te in an amorphous alloy of nominal composition Ga20Te80, obtained by rapid solidification.

Adsorption kinetics in the presence of external fields

In this article we present a detailed analysis of the kinetics of a class of sequential adsorption models that take into account the effect of externally applied fields (as an electric field, or a shear rate) on the adsorption. The excluded volume interactions related to the finite size of the adsorbing particles are modified by the external fields. As a result, new adsorption mechanisms appear with respect to the ones used to describe the kinetics in a quiescent fluid. In particular, if the adsorbing particles are allowed to roll over preadsorbed ones, adsorption becomes non local even in the simplest geometry. An exact analytic theory cannot be developed, but we introduce a self-consistent theory that turns out to agree with the simulation results over all the range of the parameters.

Adsorption kinetics in the presence of external fields

In this article we present a detailed analysis of the kinetics of a class of sequential adsorption models that take into account the effect of externally applied fields (as an electric field, or a shear rate) on the adsorption. The excluded volume interactions related to the finite size of the adsorbing particles are modified by the external fields. As a result, new adsorption mechanisms appear with respect to the ones used to describe the kinetics in a quiescent fluid. In particular, if the adsorbing particles are allowed to roll over preadsorbed ones, adsorption becomes non local even in the simplest geometry. An exact analytic theory cannot be developed, but we introduce a self-consistent theory that turns out to agree with the simulation results over all the range of the parameters.

Inhibition of glutathione efflux in the perfused rat liver and isolated hepatocytes by organic anions and bilirubin. Kinetics, sidedness, and molecular forms.

Using isolated, in situ, single-pass perfused rat livers, incubations of freshly isolated hepatocytes, and sinusoidal membrane-enriched vesicles, we and others have shown the saturability of transport (efflux) of hepatic glutathione (GSH). These observations have implicated a carrier mechanism. Our present studies were designed to provide further evidence in support of a carrier mechanism for hepatic GSH efflux by demonstrating competition by liver-specific ligands which are taken up by hepatocytes. Perfusing livers with different substances, we found that: (a) sulfobromophthalein-GSH (BSP-GSH) had a dose-dependent and fully reversible inhibitory effect on GSH efflux, while GSH alone did not have any effect; (b) taurocholate had no inhibitory effect; (c) all of the organic anions studied, i.e., BSP, rose bengal, indocyanine green, and unconjugated bilirubin (UCB), manifested potent, dose-dependent inhibitory effects, with absence of toxic effects and complete reversibility of inhibition in the case of UCB. The inhibitory effects of UCB could be overcome partially by raising (CoCl2-induced) hepatic GSH concentration. Because of the physiological importance of UCB, we conducted a detailed study of its inhibitory kinetics in the isolated hepatocyte model in the range of circulating concentrations of UCB. Studies with Cl- -free media, to inhibit the uptake of UCB by hepatocytes, showed that the inhibition of GSH efflux by UCB is apparently from inside the cell. This point was confirmed by showing that the inhibition is overcome only when bilirubin-loaded cells are cleared of bilirubin (incubation with 5% bovine serum albumin). Using Gunn rat hepatocytes and purified bilirubin mono- and diglucuronides, we found that both UCB and glucuronide forms of bilirubin inhibit GSH efflux in a dose-dependent manner. We conclude that the organic anions, although taken up by a mechanism independent of GSH, may competitively inhibit the carrier for GSH efflux from inside the hepatocyte.

Fine-Tuning Translation Kinetics Selection as the Driving Force of Codon Usage Bias in the Hepatitis A Virus Capsid

Hepatitis A virus (HAV), the prototype of genus Hepatovirus, has several unique biological characteristics that distinguish it from other members of the Picornaviridae family. Among these, the need for an intact eIF4G factor for the initiation of translation results in an inability to shut down host protein synthesis by a mechanism similar to that of other picornaviruses. Consequently, HAV must inefficiently compete for the cellular translational machinery and this may explain its poor growth in cell culture. In this context of virus/cell competition, HAV has strategically adopted a naturally highly deoptimized codon usage with respect to that of its cellular host. With the aim to optimize its codon usage the virus was adapted to propagate in cells with impaired protein synthesis, in order to make tRNA pools more available for the virus. A significant loss of fitness was the immediate response to the adaptation process that was, however, later on recovered and more associated to a re-deoptimization rather than to an optimization of the codon usage specifically in the capsid coding region. These results exclude translation selection and instead suggest fine-tuning translation kinetics selection as the underlying mechanism of the codon usage bias in this specific genome region. Additionally, the results provide clear evidence of the Red Queen dynamics of evolution since the virus has very much evolved to re-adapt its codon usage to the environmental cellular changing conditions in order to recover the original fitness.