Stochastic chemical kinetics and the quasi-steady-state assumption : application to the stochastic simulation algorithm and chemical master equation

Recently the application of the quasi-steady-state approximation (QSSA) to the stochastic simulation algorithm (SSA) was suggested for the purpose of speeding up stochastic simulations of chemical systems that involve both relatively fast and slow chemical reactions [Rao and Arkin, J. Chem. Phys. 118, 4999 (2003)] and further work has led to the nested and slow-scale SSA. Improved numerical efficiency is obtained by respecting the vastly different time scales characterizing the system and then by advancing only the slow reactions exactly, based on a suitable approximation to the fast reactions. We considerably extend these works by applying the QSSA to numerical methods for the direct solution of the chemical master equation (CME) and, in particular, to the finite state projection algorithm [Munsky and Khammash, J. Chem. Phys. 124, 044104 (2006)], in conjunction with Krylov methods. In addition, we point out some important connections to the literature on the (deterministic) total QSSA (tQSSA) and place the stochastic analogue of the QSSA within the more general framework of aggregation of Markov processes. We demonstrate the new methods on four examples: Michaelis–Menten enzyme kinetics, double phosphorylation, the Goldbeter–Koshland switch, and the mitogen activated protein kinase cascade. Overall, we report dramatic improvements by applying the tQSSA to the CME solver.

Development of a colorimetric assay for heparanase activity suitable for kinetic analysis and inhibitor screening

The role that heparanase plays during metastasis and angiogenesis in tumors makes it an attractive target for cancer therapeutics. Despite this enzyme’s significance, most of the assays developed to measure its activity are complex. Moreover, they usually rely on labeling variable preparations of the natural substrate heparan sulfate, making comparisons across studies precarious. To overcome these problems, we have developed a convenient assay based on the cleavage of the synthetic heparin oligosaccharide fondaparinux. The assay measures the appearance of the disaccharide product of heparanase-catalyzed fondaparinux cleavage colorimetrically using the tetrazolium salt WST-1. Because this assay has a homogeneous substrate with a single point of cleavage, the kinetics of the enzyme can be reliably characterized, giving a Km of 46 μM and a kcat of 3.5 s−1 with fondaparinux as substrate. The inhibition of heparanase by the published inhibitor, PI-88, was also studied, and a Ki of 7.9 nM was determined. The simplicity and robustness of this method, should, not only greatly assist routine assay of heparanase activity but also could be adapted for high-throughput screening of compound libraries, with the data generated being directly comparable across studies.

Evolving genetic regulatory networks performing as stochastic switches

Recent studies have shown that small genetic regulatory networks (GRNs) can be evolved in silico displaying certain dynamics in the underlying mathematical model. It is expected that evolutionary approaches can help to gain a better understanding of biological design principles and assist in the engineering of genetic networks. To take the stochastic nature of GRNs into account, our evolutionary approach models GRNs as biochemical reaction networks based on simple enzyme kinetics and simulates them by using Gillespie’s stochastic simulation algorithm (SSA). We have already demonstrated the relevance of considering intrinsic stochasticity by evolving GRNs that show oscillatory dynamics in the SSA but not in the ODE regime. Here, we present and discuss first results in the evolution of GRNs performing as stochastic switches.

A quasi stationary approach to drying kinetics of cylindrical food particulates in a heat pump asisted fluid bed dryer

Experiments were undertaken to study drying kinetics of moist cylindrical shaped food particulates during fluidised bed drying. Cylindrical particles were prepared from Green beans with three different length:diameter ratios, 3:1, 2:1 and 1:1. A batch fluidised bed dryer connected to a heat pump system was used for the experimentation. A Heat pump and fluid bed combination was used to increase overall energy efficiency and achieve higher drying rates. Drying kinetics, were evaluated with non-dimensional moisture at three different drying temperatures of 30, 40 and 50o C. Numerous mathematical models can be used to calculate drying kinetics ranging from analytical models with simplified assumptions to empirical models built by regression using experimental data. Empirical models are commonly used for various food materials due to their simpler approach. However problems in accuracy, limits the applications of empirical models. Some limitations of empirical models could be reduced by using semi-empirical models based on heat and mass transfer of the drying operation. One such method is the quasi-stationary approach. In this study, a modified quasi-stationary approach was used to model drying kinetics of the cylindrical food particles at three drying temperatures.

Prediction of fluidization behaviour and a quasi-stationary approach to drying kinetics of irregular particulate food materials

Changes in fluidization behaviour behaviour was characterised for parallelepiped particles with three aspect ratios, 1:1, 2:1 and 3:1 and spherical particles. All drying experiments were conducted at 500C and 15 % RH using a heat pump dehumidifier system. Fluidization experiments were undertaken for the bed heights of 100, 80, 60 and 40 mm and at 10 moisture content levels. Due to irregularities in shape minimum fluidisation velocity of parallelepiped particulates (potato) could not fitted to any empirical model. Also a generalized equation was used to predict minimum fluidization velocity. The modified quasi-stationary method (MQSM) has been proposed to describe drying kinetics of parallelepiped particulates at 30o C, 40o C and 50o C that dry mostly in the falling rate period in a batch type fluid bed dryer.

Effects of inoculating dose on the kinetics of Chlamydia muridarum genital infection in female mice

Chlamydia trachomatis infections have been implicated in problems such as pelvic inflammatory disease and infertility in females. Although there are some studies examining the kinetics of ascending infection, there is limited information on the kinetics of pathology development and cellular infiltrate into the reproductive tissues in relation to the effects of inoculating dose, and a better understanding of these is needed. The murine model of female genital tract Chlamydia muridarum infection is frequently used as a model of human C. trachomatis reproductive tract infection. To investigate the kinetics of ascending genital infection and associated pathology development, female BALB/c mice were intravaginally infected with C. muridarum at doses ranging from 5102 to 2.6106 inclusion forming units. We found that the inoculating dose affects the course of infection and the ascension of bacteria, with the highest dose ascending rapidly to the oviducts. By comparison, the lowest dose resulted in the greatest bacterial load in the lower reproductive tract. Interestingly, we found that the dose did not significantly affect inflammatory cell infiltrate in the various regions. Overall, this data show the effects of infectious dose on the kinetics of ascending chlamydial infection and associated inflammatory infiltration in BALB/c mice.

Effect of inaccuracies in anthropometric data and linked-segment inverse dynamic modeling on kinetics of gait in persons with partial foot amputation

The accuracy of data derived from linked-segment models depends on how well the system has been represented. Previous investigations describing the gait of persons with partial foot amputation did not account for the unique anthropometry of the residuum or the inclusion of a prosthesis and footwear in the model and, as such, are likely to have underestimated the magnitude of the peak joint moments and powers. This investigation determined the effect of inaccuracies in the anthropometric input data on the kinetics of gait. Toward this end, a geometric model was developed and validated to estimate body segment parameters of various intact and partial feet. These data were then incorporated into customized linked-segment models, and the kinetic data were compared with that obtained from conventional models. Results indicate that accurate modeling increased the magnitude of the peak hip and knee joint moments and powers during terminal swing. Conventional inverse dynamic models are sufficiently accurate for research questions relating to stance phase. More accurate models that account for the anthropometry of the residuum, prosthesis, and footwear better reflect the work of the hip extensors and knee flexors to decelerate the limb during terminal swing phase.

Characterisation of mesenchymal cells from osteophytes in osteoarthritis

Osteophytes form through the process of chondroid metamorphosis of fibrous tissue followed by endochondral ossification. Osteophytes have been found to consist of three different mesenchymal tissue regions including endochondral bone formation within cartilage residues, intra-membranous bone formation within fibrous tissue and bone formation within bone marrow spaces. All these features provide evidence of mesenchymal stem cells (MSC) involvement in osteophyte formation; nevertheless, it remains to be characterised. MSC from numerous mesenchymal tissues have been isolated but bone marrow remains the “ideal” due to the ease of ex vivo expansion and multilineage potential. However, the bone marrow stroma has a relatively low number of MSC, something that necessitates the need for long-term culture and extensive population doublings in order to obtain a sufficient number of cells for therapeutic applications. MSC in vitro have limited proliferative capacity and extensive passaging compromises differentiation potential. To overcome this barrier, tissue derived MSC are of strong interest for extensive study and characterisation, with a focus on their potential application in therapeutic tissue regeneration. To date, no MSC type cell has been isolated from osteophyte tissue, despite this tissue exhibiting all the hallmark features of a regenerative tissue. Therefore, this study aimed to isolate and characterise cells from osteophyte tissues in relation to their phenotype, differentiation potential, immuno-modulatory properties, proliferation, cellular ageing, longevity and chondrogenesis in in vitro defect model in comparison to patient matched bone marrow stromal cells (bMSC). Osteophyte derived cells were isolated from osteophyte tissue samples collected during knee replacement surgery. These cells were characterised by the expression of cell surface antigens, differentiation potential into mesenchymal lineages, growth kinetics and modulation of allo-immune responses. Multipotential stem cells were identified from all osteophyte samples namely osteophyte derived mesenchymal stem cells (oMSC). Extensively expanded cell cultures (passage 4 and 9 respectively) were used to confirm cytogenetic stability and study signs of cellular aging, telomere length and telomerase activity. Cultured cells at passage 4 were used to determine 84 pathway focused stem cell related gene expression profile. Micro mass pellets were cultured in chondrogenic differentiation media for 21 days for phenotypic and chondrogenic related gene expression. Secondly, cell pellets differentiated overnight were placed into articular cartilage defects and cultured for further 21 days in control medium and chondrogenic medium to study chondrogenesis and cell behaviour. The surface antigen expression of oMSC was consistent with that of mesenchymal stem cells, such as lacking the haematopoietic and common leukocyte markers (CD34, CD45) while expressing those related to adhesion (CD29, CD166, CD44) and stem cells (CD90, CD105, CD73). The proliferation capacity of oMSC in culture was superior to that of bMSC, and they readily differentiated into tissues of the mesenchymal lineages. oMSC also demonstrated the ability to suppress allogeneic T-cell proliferation, which was associated with the expression of tryptophan degrading enzyme indoleamine 2,3 dioxygenase (IDO). Cellular aging was more prominent in late passage bMSC than in oMSC. oMSC had longer telomere length in late passages compared with bMSC, although there was no significant difference in telomere lengths in the early passages in either cell type. Telomerase activity was detectable only in early passage oMSC and not in bMSC. In osteophyte tissues telomerase positive cells were found to be located peri vascularly and were Stro-1 positive. Eighty-four pathway-focused genes were investigated and only five genes (APC, CCND2, GJB2, NCAM and BMP2) were differentially expressed between bMSC and oMSC. Chondrogenically induced micro mass pellets of oMSC showed higher staining intensity for proteoglycans, aggrecan and collagen II. Differential expression of chondrogenic related genes showed up regulation of Aggrecan and Sox 9 in oMSC and collagen II in bMSC. The in vitro defect models of oMSC in control medium showed rounded and aggregated cells staining positively for proteoglycan and presence of some extracellular matrix. In contrast, defects with bMSC showed fragmentation and loss of cells, fibroblast-like cell morphology staining positively for proteoglycans. For defects maintained in chondrogenic medium, rounded, aggregated and proteoglycan positive cells were found in both oMSC and bMSC cultures. Extracellular matrix and cellular integration into newly formed matrix was evident only in oMSC defects. For analysis of chondrocyte hypertrophy, strong expression of type X collagen could be noticed in the pellet cultures and transplanted bMSC. In summary, this study demonstrated that osteophyte derived cells had similar properties to mesenchymal stem cells in the expression of antigen phenotype, differential potential and suppression of allo-immune response. Furthermore, when compared to bMSC, oMSC maintained a higher proliferative capacity due to a retained level of telomerase activity in vitro, which may account for the relatively longer telomeres delaying growth arrest by replicative senescence compared with bMSC. oMSC behaviour in defects supported chondrogenesis which implies that cells derived from regenerative tissue can be an alternative source of stem cells and have a potential clinical application for therapeutic stem cell based tissue regeneration.

Effects of pediatric obesity on joint kinematics and kinetics during 2 walking cadences

Objective: To determine whether differences existed in lower-extremity joint biomechanics during self-selected walking cadence (SW) and fast walking cadence (FW) in overweight- and normal-weight children.---------- Design: Survey.---------- Setting: Institutional gait study center.---------- Participants: Participants (N=20; mean age ± SD, 10.4±1.6y) from referred and volunteer samples were classified based on body mass index percentiles and stratified by age and sex. Exclusion criteria were a history of diabetes, neuromuscular disorder, or recent lower-extremity injury.---------- Main Outcome Measures: Sagittal, frontal, and transverse plane angular displacements (degrees) and peak moments (newton meters) at the hip, knee, and ankle joints.---------- Results: The level of significance was set at P less than .008. Compared with normal-weight children, overweight children had greater absolute peak joint moments at the hip (flexor, extensor, abductor, external rotator), the knee (flexor, extensor, abductor, adductor, internal rotator), and the ankle (plantarflexor, inverter, external/internal rotators). After including body weight as a covariate, overweight children had greater peak ankle dorsiflexor moments than normal-weight children. No kinematic differences existed between groups. Greater peak hip extensor moments and less peak ankle inverter moments occurred during FW than SW. There was greater angular displacement during hip flexion as well as less angular displacement at the hip (extension, abduction), knee (flexion, extension), and ankle (plantarflexion, inversion) during FW than SW.---------- Conclusions: Overweight children experienced increased joint moments, which can have long-term orthopedic implications and suggest a need for more nonweight-bearing activities within exercise prescription. The percent of increase in joint moments from SW to FW was not different for overweight and normal-weight children. These findings can be used in developing an exercise prescription that must involve weight-bearing activity.