Parallel implementation of stochastic simulation for large scale cellular processes

Experimental and theoretical studies have shown the importance of stochastic processes in genetic regulatory networks and cellular processes. Cellular networks and genetic circuits often involve small numbers of key proteins such as transcriptional factors and signaling proteins. In recent years stochastic models have been used successfully for studying noise in biological pathways, and stochastic modelling of biological systems has become a very important research field in computational biology. One of the challenge problems in this field is the reduction of the huge computing time in stochastic simulations. Based on the system of the mitogen-activated protein kinase cascade that is activated by epidermal growth factor, this work give a parallel implementation by using OpenMP and parallelism across the simulation. Special attention is paid to the independence of the generated random numbers in parallel computing, that is a key criterion for the success of stochastic simulations. Numerical results indicate that parallel computers can be used as an efficient tool for simulating the dynamics of large-scale genetic regulatory networks and cellular processes

Searching oligo sets of human chromosome 12 using evolutionary strategies

DNA microarray is a powerful tool to measure the level of a mixed population of nucleic acids at one time, which has great impact in many aspects of life sciences research. In order to distinguish nucleic acids with very similar composition by hybridization, it is necessary to design probes with high specificities, i.e. uniqueness, and also sensitivities, i.e., suitable melting temperature and no secondary structure. We make use of available biology tools to gain necessary sequence information of human chromosome 12, and combined with evolutionary strategy (ES) to find unique subsequences representing all predicted exons. The results are presented and discussed.

The generation of H-1-NMR-detectable mobile lipid in stimulated lymphocytes: Relationship ts cellular activation, the cell cycle, and phosphatidylcholine-specific phospholipase C

Mobile Lipids detected using H-1-NMR in stimulated lymphocytes were correlated with cell cycle phase, expression of the interleukin-2 receptor alpha and proliferation to assess the activation status of the lymphocytes. Mobile lipid levels, IL-2R alpha expression and proliferation increased after treatment with PMA and ionomycin. PMA or ionomycin stimulation alone induced increased IL-2R alpha expressiom but not proliferation, PMA- but not ionomycin-stimulation generated mobile lipid, Treatment with anti-CD3 antibody did not increase IL-2R alpha expression or proliferation but did generate increased amounts of mobile lipid, The cell cycle status of thymocytes treated with anti-CD3, PMA or ionomycin alone indicated an. accumulation of the cells in the G(1) phase of the cell cycle, The generation of mobile lipid was abrogated in anti-CD3 antibody-stimulated thymic lymphocytes but not in splenic lymphocytes, using a phosphatidylcholine-specific phospholipase C (PC-PLC) inhibitor which blocked cells in the G(1)/S phase of the cell cycle, This suggests that the H-1-NMR-detectable mobile Lipid may be generated in anti-CD3 antibody-stimulated thymic lymphocytes by the action of PC-PLC activity via the catabolism of PC, in the absence of classical signs of activation. (C) 1997 Academic Press.

Mechanics of cellular adhesion to artificial artery templates

We are using polymer templates to grow artificial artery grafts in vivo for the replacement of diseased blood vessels. We have previously shown that adhesion of macrophages to the template starts the graft formation. We present a study of the mechanics of macrophage adhesion to these templates on a single cell and single bond level with optical tweezers. For whole cells, in vitro cell adhesion densities decreased significantly from polymer templates polyethylene to silicone to Tygon (167, 135, and 65 cells/mm(2)). These cell densities were correlated with the graft formation success rate (50%, 25%, and 0%). Single-bond rupture forces at a loading rate of 450 pN/s were quantified by adhesion of trapped 2-mm spheres to macrophages. Rupture force distributions were dominated by nonspecific adhesion (forces, < 40 pN). On polystyrene, preadsorption of fibronectin or presence of serum proteins in the cell medium significantly enhanced adhesion strength from a mean rupture force of 20 pN to 28 pN or 33 pN, respectively. The enhancement of adhesion by fibronectin and serum is additive (mean rupture force of 43 pN). The fraction of specific binding forces in the presence of serum was similar for polystyrene and polymethyl-methacrylate, but specific binding forces were not observed for silica. Again, we found correlation to in vivo experiments, where the density of adherent cells is higher on polystyrene than on silica templates, and can be further enhanced by fibronectin adsorption. These findings show that in vitro adhesion testing can be used for template optimization and to substitute for in-vivo experiments.

Simulation of Brugada syndrome using cellular and three-dimensional whole-heart modeling approaches

Brugada syndrome (BS) is a genetic disease identified by an abnormal electrocardiogram ( ECG) ( mainly abnormal ECGs associated with right bundle branch block and ST-elevation in right precordial leads). BS can lead to increased risk of sudden cardiac death. Experimental studies on human ventricular myocardium with BS have been limited due to difficulties in obtaining data. Thus, the use of computer simulation is an important alternative. Most previous BS simulations were based on animal heart cell models. However, due to species differences, the use of human heart cell models, especially a model with three-dimensional whole-heart anatomical structure, is needed. In this study, we developed a model of the human ventricular action potential (AP) based on refining the ten Tusscher et al (2004 Am. J. Physiol. Heart Circ. Physiol. 286 H1573 - 89) model to incorporate newly available experimental data of some major ionic currents of human ventricular myocytes. These modified channels include the L-type calcium current (ICaL), fast sodium current (I-Na), transient outward potassium current (I-to), rapidly and slowly delayed rectifier potassium currents (I-Kr and I-Ks) and inward rectifier potassium current (I-Ki). Transmural heterogeneity of APs for epicardial, endocardial and mid-myocardial (M) cells was simulated by varying the maximum conductance of IKs and Ito. The modified AP models were then used to simulate the effects of BS on cellular AP and body surface potentials using a three-dimensional dynamic heart - torso model. Our main findings are as follows. (1) BS has little effect on the AP of endocardial or mid-myocardial cells, but has a large impact on the AP of epicardial cells. (2) A likely region of BS with abnormal cell AP is near the right ventricular outflow track, and the resulting ST-segment elevation is located in the median precordium area. These simulation results are consistent with experimental findings reported in the literature. The model can reproduce a variety of electrophysiological behaviors and provides a good basis for understanding the genesis of abnormal ECG under the condition of BS disease.

Sequence and tissue distribution of chicken cellular nucleic acid binding protein cDNA

We sequenced cDNAs coding for chicken cellular nucleic acid binding protein (CNBP). Two slightly different variations of the open reading frame were found, each of which translates into a protein with seven zinc finger domains. The longest transcript contains an in-frame insert of 3 bp. The sequence conservation between chick CNBP cDNAs with human, rat and mouse CNBP cDNAs is extreme, especially in the coding region, where the deduced amino acid sequence identity with human, rat and mouse CNBP is 99%. CNBP-like transcripts were also found in various tissues from insect, shrimp, fish and lizard. Regions with remarkable nucleotide conservation were also found in the 3' untranslated region, indicating important functions for these regions. Quantitative reverse transcription polymerase chain reaction (RT-PCR) indicated that in the chick, CNBP is present in all tissues examined in approximately equal ratios to total RNA. RT-PCR of total RNA isolated from different phyla indicate CNBP-like proteins art widespread throughout the animal kingdom. The extraordinary level of conservation suggests an important physiological role for CNBP. (C) 1997 Elsevier Science Inc.

A protein with protective properties against the cellular defense reactions in insects

The molecular mechanism of how insects recognize intruding microorganisms and parasites and distinguish them from own body structures is not well known. We explored evolutionary adaptations in an insect parasitoid host interaction to identify components that interfere with the recognition of foreign objects and cellular encapsulation. Because some parasitoids provide protection for the developing wasp in the absence of an overt suppression of the insect host defense, we analyzed the surface of eggs and symbiotic viruses for protective properties. Here we report on the molecular cloning of a 32-kDa protein (Crp32) that is one of the major protective components. It is produced in the calyx cells of the female wasp ovaries and attached to the surface of the egg and other particles including polydnaviruses. The recombinant protein confers protection to coated objects in a cellular encapsulation assay suggesting that a layer of Crp32 may prevent cellular encapsulation reactions by a local inactivation of the host defense system.

Cellular origin of chlorinated diketopiperazines in the dictyoceratid sponge Dysidea herbacea (Keller)

The tropical marine sponge Dysidea herbacea (Keller) contains the filamentous unicellular cyanobacterium Oscillatoria spongeliae (Schulze) Hauck as an endosymbiont, plus numerous bacteria, both intracellular and extracellular. Archaeocytes and choanocytes are the major sponge cell types present. Density gradient centrifugation of glutaraldehyde-fixed cells with Percoll as the support medium has been used to separate the cyanobacterial symbiont from the sponge cells on the basis of their differing densities. The protocol also has the advantage of separating broken from intact cells of O. spongeliae. The lighter cell preparations contain archaeocytes and choanocytes together with damaged cyanobacterial cells, whereas heavier cell preparations contain intact cyanobacterial cells, with less than 1% contamination by sponge cells. Gas chromatography/mass spectrometry analysis has revealed that the terpene spirodysin is concentrated in preparations containing archaeocytes and choanocytes, whereas nuclear magnetic resonance analysis of the symbiont cell preparations has shown that they usually contain the chlorinated diketopiperazines, dihydrodysamide C and didechlorodihydrodysamide C, which are the characteristic metabolites of the sponge/symbiont association. However, one symbiont preparation, partitioned by a second Percoll gradient, has been found to be devoid of chlorinated diketopiperazines. The capability to synthesize secondary metabolites may depend on the physiological state of the symbiont; alternatively, there may be two closely related cyanobacterial strains within the sponge tissue.

A sponge/dinoflagellate association in the haplosclerid sponge Haliclona sp.: cellular origin of cytotoxic alkaloids by Percoll density gradient fractionation

Light-microscopic and electron-microscopic studies of the tropical marine sponge Haliclona sp. (Or der: Haplosclerida Family: Haliclonidae) from Heron Island, Great Barrier Reef, have revealed that this sponge is characterized by the presence of dinoflagellates and by nematocysts. The dinoflagellates are 7-10 mu m in size, intracellular, and contain a pyrenoid with a single stalk, whereas the single chloroplast is branched, curved, and lacks grana. Mitochondria are present, and the nucleus is oval and has distinct chromosomal structure. The dinoflagellates are morphologically similar to Symbiodinium microadriaticum, the common intracellular symbiont of corals, although more detailed biochemical and molecular studies are required to provide a precise taxonomic assignment. The major sponge cell types found in Haliclona sp, are spongocytes, choanocytes, and archaeocytes; groups of dinoflagellates are enclosed within large vacuoles in the archaeocytes. The occurrence of dinoflagellates in marine sponges has previously been thought to be restricted to a small group of sponges including the excavating hadromerid sponges; the dinoflagellates in these sponges are usually referred to as symbionts. The role of the dinoflagellates present in Haliclona sp. as a genuine symbiotic partner requires experimental investigation. The sponge grows on coral substrates, from which it may acquire the nematocysts, and shows features, such as mucus production, which are typical of some excavating sponges. The cytotoxic alkaloids, haliclonacyclamines A and B, associated with Haliclona sp. are shown by Percoll density gradient fractionation to be localized within the sponge cells rather than the dinoflagellates. The ability to synthesize bioactive compounds such as the haliclonacyclamines may help Haliclona sp. to preserve its remarkable ecological niche.

Effect of estrogen on vascular smooth muscle cells is dependent upon cellular phenotype

To investigate the growth-regulating action of estrogen on vascular smooth muscle cells (SMC), effects of beta-17-estradiol (beta-E-2) on phenotypic modulation and proliferation of rabbit aortic SMC were observed in vitro. At 10(-8) M, beta-E-2 significantly slowed the decrease in volume fraction of myofilaments (V(v)myo) of freshly dispersed SMCs in primary culture, indicating an inhibitory effect of beta-E-2 On spontaneous phenotypic modulation of SMC from a contractile to a synthetic phenotype. Freshly dispersed SMCs treated with beta-E-2 also had a relatively longer quiescent phase than control cells before intense proliferation occurred. This was in contrast to SMCs in passage 2-3 (synthetic state), where beta-E-2-treated cells replicated significantly faster than untreated cells. beta-E-2 also markedly enhanced the serum-induced DNA synthesis of synthetic SMCs in a concentration-dependent manner within physiological range (10(-10) to 10-8 M). These findings indicate that the growth-regulating effect of estrogen on vascular SMC is dependent on the cell's phenotypic stare. It delays the cell cycle re-entry of the contractile SMCs by retarding their phenotypic modulation however, once cells have modulated to the synthetic phenotype, it promotes their replication. (C) 1998 Elsevier Science Ireland Ltd. All rights reserved.

Murine cytomegalovirus homologues of cellular immunomodulatory genes

The study of 'molecular mimicry' or 'genetic piracy', with respect to the utilisation of cellular genes captured and modified during the course of virus evolution, has been an area of increasing research with the expansion in virus genome sequencing. Examples of cellular immunomodulatory genes which have been captured from hosts have been identified in a number of viruses. This review concentrates upon studies of murine cytomegalovirus (MCMV), investigating the functions of viral genes homologous to G protein-coupled receptors, MHC class I and chemokines, The study of recombinant MCMV engineered with specific disruptions of these genes has revealed their significance during virus replication and dissemination within the host, In the case of the latter two classes of genes, evidence suggests they interfere with cellular immune responses, although the detailed mechanisms underlying this interference have yet to be delineated. Copyright (C) 2000 S. Karger AG, Basel.

Evolution of stress deficit and changing rates of seismicity in cellular automaton models of earthquake faults

We investigate the internal dynamics of two cellular automaton models with heterogeneous strength fields and differing nearest neighbour laws. One model is a crack-like automaton, transferring ail stress from a rupture zone to the surroundings. The other automaton is a partial stress drop automaton, transferring only a fraction of the stress within a rupture zone to the surroundings. To study evolution of stress, the mean spectral density. f(k(r)) of a stress deficit held is: examined prior to, and immediately following ruptures in both models. Both models display a power-law relationship between f(k(r)) and spatial wavenumber (k(r)) of the form f(k(r)) similar tok(r)(-beta). In the crack model, the evolution of stress deficit is consistent with cyclic approach to, and retreat from a critical state in which large events occur. The approach to criticality is driven by tectonic loading. Short-range stress transfer in the model does not affect the approach to criticality of broad regions in the model. The evolution of stress deficit in the partial stress drop model is consistent with small fluctuations about a mean state of high stress, behaviour indicative of a self-organised critical system. Despite statistics similar to natural earthquakes these simplified models lack a physical basis. physically motivated models of earthquakes also display dynamical complexity similar to that of a critical point system. Studies of dynamical complexity in physical models of earthquakes may lead to advancement towards a physical theory for earthquakes.

Accelerating seismic energy release and evolution of event time and size statistics: Results from two heterogeneous cellular automaton models

The evolution of event time and size statistics in two heterogeneous cellular automaton models of earthquake behavior are studied and compared to the evolution of these quantities during observed periods of accelerating seismic energy release Drier to large earthquakes. The two automata have different nearest neighbor laws, one of which produces self-organized critical (SOC) behavior (PSD model) and the other which produces quasi-periodic large events (crack model). In the PSD model periods of accelerating energy release before large events are rare. In the crack model, many large events are preceded by periods of accelerating energy release. When compared to randomized event catalogs, accelerating energy release before large events occurs more often than random in the crack model but less often than random in the PSD model; it is easier to tell the crack and PSD model results apart from each other than to tell either model apart from a random catalog. The evolution of event sizes during the accelerating energy release sequences in all models is compared to that of observed sequences. The accelerating energy release sequences in the crack model consist of an increase in the rate of events of all sizes, consistent with observations from a small number of natural cases, however inconsistent with a larger number of cases in which there is an increase in the rate of only moderate-sized events. On average, no increase in the rate of events of any size is seen before large events in the PSD model.

Identification of cellular changes associated with increased production of human growth hormone in a recombinant Chinese hamster ovary cell line

A proteomics approach was used to identify the proteins potentially implicated in the cellular response concomitant with elevated production levels of human growth hormone in a recombinant Chinese hamster ovary (CHO) cell line following exposure to 0.5 mM butyrate and 80 muM zinc sulphate in the production media. This involved incorporation of two-dimensional (2-D) gel electrophoresis and protein identification by a combination of N-terminal sequencing, matrix-assisted laser desorption/ionisation-time of flight mass spectrometry, amino acid analysis and cross species database matching. From these identifications a CHO 2-D reference,map and annotated database have been established. Metabolic labelling and subsequent autoradiography showed the induction of a number of cellular proteins in response to the media additives butyrate and zinc sulphate. These were identified as GRP75, enolase and thioredoxin. The chaperone proteins GRP78, HSP90, GRP94 and HSP70 were not up-regulated under these conditions.