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.     

The CR3 motif of Rrp44p is important for interaction with the core exosome and exosome function

The 10-subunit RNA exosome is involved in a large number of diverse RNA processing and degradation events in eukaryotes. These reactions are carried out by the single catalytic subunit, Rrp44p/Dis3p, which is composed of three parts that are conserved throughout eukaryotes. The exosome is named for the 3′ to 5′ exoribonuclease activity provided by a large C-terminal region of the Rrp44p subunit that resembles other exoribonucleases. Rrp44p also contains an endoribonuclease domain. Finally, the very N-terminus of Rrp44p contains three Cys residues (CR3 motif) that are conserved in many eukaryotes but have no known function. These three conserved Cys residues cluster with a previously unrecognized conserved His residue in what resembles a metal-ion-binding site. Genetic and biochemical data show that this CR3 motif affects both endo- and exonuclease activity in vivo and both the nuclear and cytoplasmic exosome, as well as the ability of Rrp44p to associate with the other exosome subunits. These data provide the first direct evidence that the exosome-Rrp44p interaction is functionally important and also provides a molecular explanation for the functional defects when the conserved Cys residues are mutated.     

The cell cycle modulated glycoprotein GP115 is one of the major yeast proteins containing glycosylphosphatidylinositol

The cell cycle modulated protein gp115 (115 kDa, isoelectric point about 4.8-5) of Saccharomyces cerevisiae undergoes various post-translational modifications. It is N-glycosylated during its maturation along the secretory pathway where an intermediary precursor of 100 kDa (p100), dynamically related to the mature gp115 protein, is detected at the level of endoplasmic reticulum. Moreover, we have shown by the use of metabolic labeling with [35S]methionine, [3H]palmitic acid and myo-[3H]inositol combined with high resolution two-dimensional gel electrophoresis and immunoprecipitation with a specific antiserum, that gp115 is one of the major palmitate- and inositol-containing proteins in yeast. These results, and the susceptibility of gp115 to phosphatidylinositol-specific phospholipase C treatment strongly indicate that gp115 contains the glycosylphosphatidylinositol (GPI) structure as membrane anchor domain. The two-dimensional analysis of the palmitate- and inositol-labeled proteins has also allowed the characterization of other polypeptides which possibly contain a GPI structure.