g it is an open question whether CDP might up regulate Th2 speci

g. it is an open question whether CDP might up regulate Th2 specific genes or down regulate the genes in Th0 and Th1 lineages. The three TF hits having enriched predicted binding sites among the Th2 down regulated genes were the in terferon references regulatory factor family of TFs, IFN stimulated genes factor 3 and STAT6. IRF family consists of IRF1 to IRF10 and has been shown to be essential in ex pression of type I interferon genes, IFN stimulated genes and other pro inflammatory response related cyto kines. These genes are maintained down regulated during Th2 proliferation and therefore, the results are in line with the Th2 effector cells characteristics. More over, IFN�� induced expression of IRF1 and IRF2 has been shown to directly down regulate IL 4 production by repressing IL 4 promoter sites.

Opposing to other IRF family proteins, IRF4 has been shown to directly activate IL 4 promoter and IL 10 regulatory elements and be essential in Th2 cell differentiation by influencing the expression of GFI1, a transcriptional repressor in Th2 cells. However, the analysis relying on known TF bin ding specificities will not allow segregation of individual members of the IRF family. Further, an essential regulator of most ISGs is ISGF3 that is composed of STAT1, STAT2 and IRF9 complex and works in conjunction with IRFs. Identification of STAT6 as a regulator among the Th2 down regulated genes is well in line with our previ ously published results, although its effect was observed to be less profound within Th2 down regulated genes than among Th2 up regulated target genes.

Compari son analysis of the predicted STAT6 target genes and Th2 up regulated and down regulated genes gave 16 and 19 overlapping genes, respectively. The full lists of overlap ping genes are in Additional file 3, Table S2. We further analyzed the correlation between predicted STAT6 target promoters and experimentally observed promoter asso ciated binding sites, and observed signifi cant correlation between the target sites. The full list of predicted STAT6 target genes and promoter asso ciated STAT6 binding sites identified by ChIP seq as well as the overlapping genes are listed in the Additional file 3, Table S2. The overlapping binding sites included promo ters for C14orf177, CISH, HMMR, INO80, MGAT1, NUDCD2, SOCS1, SPINT2 and ZNF570 genes.

Discussion Identification of the key T helper cell regulators Batimastat provides possible targets for modulation of immune response. To reveal T cell subset specific genes and their often subtle differences in expression, we developed a novel compu tational method, LIGAP. Traditional ways of identifying differentially expressed genes, such as the t test, are pro blematic in studying time series data since there is a need to carry out hypothesis tests on individual time points. On the other hand, commonly used statistical tests for whole time course, including e. g.

As many of the small repeated sequences are highly helical in pre

As many of the small repeated sequences are highly helical in predicted structure, one could suggest they are involved in DNA binding thenthereby and regulation. Further work is needed to determine when they are expressed and at what stage of the life cycle. When analysis of the Pt and Pst genomes has been concluded, it can be determined if the repeated nature of these predicted genes is maintained within the wheat rust fungi. Methods Pt BAC library Total genomic DNA for the BAC library construction was isolated from P. triticina Race1, BBBD urediniospores collected from susceptible wheat cultivar Thatcher. Spores were increased on plants spray inoculated with a urediniospore suspension in light mineral oil.

The oil was allowed to evaporate for 30 min, then plants were moved to a dark dew chamber at 20 C and 100% relative humidity for 24 hrs for uredinios pore germination and appressorium formation. Plants were grown in a growth chamber under 16 hour day at 20 C. After 10 days, urediniospores were collected and germi nated by densely dusting them over sterile water in dishes for 8 hrs using a volatile nonanol solution, 1 ml acetone, 19 ml of ddH2O spotted on filter paper which was suspended in the lids to stimulate urediniospore germination under crowded conditions. The BAC library was constructed by BioS T. In brief, nu clei were isolated from collected germinated urediniospores and embedded in 1% low melting point agarose plugs. Total genomic DNA embedded in the plugs was partially digested with HindIII, separated by electrophoresis by pulse field gel electrophoresis, and the 100 200 kb region was isolated.

After electro elution and dialysis, the DNA fragments were cloned into the HindIII site of BAC vector pIndogoBAC5 and propagated in E. coli DH10B. BAC clone selection and sequencing The resulting BAC library of 15,360 individual clones was arrayed on nylon membranes. After colony lysis, DNA was bound to the membranes using standard procedures. BAC filters were probed to identify clones for sequencing. Several candidate fragments were selected as probes. The Sfi1 insert from a Pt cDNA clone, PT0313. J16. C21 was labeled with P32 dCTP using a random primer labeling kit. Selected BAC clones were sent as a stab culture to the Genome Center at Washington University, St. Louis, MO. BAC clones were cultured, subcloned, shot gun sequenced, and assembled.

Gene calls were made using FGENESH with gene models specific to Puccinia. BAC clone gene AV-951 predictions were compared to Pgt, Mlp and Um genomic resources using the BLASTN and BLASTX algorithms with settings of E value 1e 3, Matrix BLOSUM62, and gapped alignment. Repeats were identified using fungaldb of RepBase 17. 04, containing the repeats of Pt, Pgt and Pst. Long terminal repeats were determined by LTR Finder.

Moreover, to evaluate HOXB1 epigenetic regulation by the

Moreover, to evaluate HOXB1 epigenetic regulation by the selleck inhibitor histones acetylation deacetylation mechanisms, we treated the HL60 cells with 100 or 600 ng of the histone deacetylase inhibitor Trichostatin A for 48 and 72 hr. Following all the above mentioned treatments, we searched for HOXB1 mRNA re expression in HL60 cells by RT PCR. Statistical analysis All the experiments were repeated at least three times, unless otherwise stated. Reported values represent mean standard errors. The significance of differences between experimental variables was determined using parametric Students t test with P 0. 05 deemed statisti cally significant. P values relative to HOXB1 transduced cells were always referred to LXSN transduced cells.

Results HOXB1 is downregulated in leukemic cells We evaluated the endogenous expression of HOXB1 in a panel of representative primary acute myeloid leukemia cells, staged from M1 to M6, and some stabilized leukemic cell lines. As normal controls, we utilized termin ally differentiated cells, including granulocytes, monocytes, macrophages, erythroblasts and lymphocytes, as well as CD34 progenitors from peripheral blood. As determined by qReal Time and traditional RT PCR, HOXB1 was barely or not expressed in all the examined neoplastic cells, even after 40 cycles of amplification, whereas it was detectable, at RNA and protein levels, in normal cells purified from peripheral blood and in CD34 progenitors. Among the AMLs the exceptions, showing HOXB1 expression, were the M6 staged erythroleukemias and the K562 cell line, possibly in agreement with their predominant erythro blastic cells component.

In all the exper iments a 9 days ATRA induced teratocarcinoma NT2/D1 sample was included as a positive control. HOXB1 restored expression induces apoptosis and cell death in HL60 To investigate the functional role of HOXB1, we selected the AML193, U937, NB4 and HL60 cell lines as models for gene transduction. To this end was utilized the retro viral vector LB1SN and the correct transcription and translation of HOXB1 mRNA and protein were con firmed by qReal Time RT PCR and Western blot ana lysis. Unfortunately, as the enforced expression of HOXB1 resulted quickly lost in AML193, U937 and NB4, the sole HL60 cell line was exploitable to deter mine whether HOXB1 overexpression might actually affect the biological properties of HL60 cells.

We then performed some representative in vitro func tional assays in high and low serum condi tions. In order to evaluate the proliferative rate, cells were initially seeded Carfilzomib at 1��105/ml and monitored up to 7 days when a significant reduction of cell growth was visible in HOXB1 expressing cells, regard less of serum concentration. Looking for the cause of such reduction, we compared the total apoptotic rates detectable in HOXB1 and LXSN transduced cells.

Western

Western further blot analysis using an antibody specific for phosphorylated Ser473 residue in Akt demonstrated that the phosphorylation level of Akt protein started to decline after 16 hours of butyrate treatment, and again valproic acid had a weaker effect. Quantification of the Western blots verified that the levels of Akt phos phorylation were reduced by about 60% and 40% after 24 hours of butyrate and valproic acid treatment respectively, as compared to untreated control cells. Taken together, these data suggest that the down regulation of Akt activity by HDAC inhibitors, such as valproic acid and butyrate, is achieved through inhibition of gene expres sion of Akt1 and Akt2. Valproic acid and butyrate induce caspase 3 activation Programmed cell death can be mediated through either caspase dependent or independent mechanisms.

To determine the molecular pathway of cell apoptotic death induced by HDAC inhibitors, such as valproic acid and butyrate, we examined the active status of caspase 3 fol lowing the treatment of HeLa cells. Western blot analysis with an antibody specifically against caspase 3 showed that the HeLa cells contain an abundance of the precursor caspase 3. However, the cleaved or activated form of caspase 3 was only observed following valproic acid or butyrate treatment. Again, butyrate had a stronger effect than valproic acid on inducing caspase 3 activity, correlating directly with its higher efficacy than valproic acid on the deactivation of Akt and on the induc tion of cell apoptotic death.

Taken together, these data suggest that apoptotic cell death induced by HDAC inhibitors, such as valproic acid and butyrate, is mediated through the caspase dependent pathway. Valproic acid and butyrate activate both the caspase 8 and caspase 9 Caspase 3 is an executor caspase that can be cleaved or activated by either caspase 9, the initiator caspase of the mitochondrial pathway, or caspase 8, the initiator caspase of the death receptor pathway. Affymetrix microarray analysis showed that butyrate treatment increased the level of caspase 9 mRNA to a small extent. Quantitative real time RT PCR analysis con firmed that the capase 9 mRNA of the HeLa cells was indeed increased to about 1. 5 fold upon valproic acid or butyrate treatment. Next, we assessed the enzy matic activity of caspase 9 by using a fluorescence based assay in the presence or absence of LEHD CHO, a specific inhibitor for the caspase 9.

As shown in Fig. 3C, valproic acid or butyrate treatment induced the caspase 9 activity and again butyrate induced activity of the caspase 9 more robustly than valproic acid, about 9 fold versus 4 fold. Caspase 8 can be activated by HDAC inhibitors through the death receptors in several cancer cell lines. To determine the role of the death receptor pathway in HDAC Entinostat inhibitor induced HeLa cell death, we evaluated the activity of caspase 8 in the HeLa cells by using a fluo rescence based assay.

Briefly, the STE was analyzed by Hitachi L 6200 with an L 4500

Briefly, the STE was analyzed by Hitachi L 6200 with an L 4500 selleck chem KPT-330 Diode Array detector with a PE Sciex Qstar Pulsar ESI TOF mass spectrometer. Samples were injected onto a Merck LiChrospher 100 RP 18 column. The column was equilibrated in 0. 05% acetic acid water and elution of the compo nents was achieved by increasing the concentration of solution B from 0 to 100% in 30 min at a flow rate of 1 ml min. Absorbance was monitored at 254 nm. The molecular masses of the peaks were deter mined from electro spray ionization mass spectra using a multiply charged ion profile based on the modified method of Chang et al. For subsequent experi ments, the STE powder was dissolved in dimethyl sulfate to achieve designed concentrations.

Cell and cell culture A human nasopharyngeal carcinoma cell line from ATCC, HONE 1 cells, was cultured in RPMI 1640 medium, 10% fetal bovine serum, 2 mM glutamine, 100 U ml penicillin, and 100 ug ml streptomycin. All cell cultures were maintained at 37 C in a humidified atmosphere of 5% CO2. For STE treatment, appropriate amounts of stock solution of STE were added into the culture medium to achieve the indicated concentrations. The cells were then incubated for the indicated time periods. Dimethyl sulfox ide solution without STE was used as blank reagent. Analysis of cell viability To evaluate the cytotoxicity of STE, an MTT colorimet ric assay was performed to determine cell viability. Cells were seeded in 24 well plates at a density of 1��105 cells per well and treated with 0, 25, 50, 75, 100, 150 and 200 ug mL of STE at 37 C in 5% CO2 for 24 h and 48 h.

At the end of the exposure period, the cells were washed with PBS and incubated with 0. 8 mL of MTT per well at 37 C in 5% CO2 for 4 h. The viable cell number was directly proportional to the pro duction of formazan following solubilization with iso propanol, which was measured spectrophotometrically at 563 nm. Cell migration and invasion assays Cell migration and invasion were assayed according to the methods described by Chu et al. After treatment with STE for 24 h, the surviving HONE 1 cells were harvested and seeded to a Boyden chamber at 104 cells per well in serum free medium, and then incubated for 24 h at 37 C. To deter mine cell migration, the cells were seeded into the Boyden chamber on membrane filters that were not coated with Matrigel.

The filters were then air dried for 5 h in a lam inar flow hood. The migrating cells were fixed with methanol and stained with Giemsa. The GSK-3 cell numbers were counted by light microscopy. For the invasion assay, 10 uL Matrigel was applied to 8 um pore size polycarbonate membrane filters. The bottom cham ber contained standard medium. The invasion of cells treated or untreated with STE was measured as in the migration assay.

These functions include endocrine activities and intrauterine inv

These functions include endocrine activities and intrauterine invasion and modulation of the maternal vasculature and immune cells. Among the differentiation associated genes was a subgroup of genes encoding selleck chem inhibitor transcriptional regulators. Mouse mutagenesis experimentation has implicated a few of these genes as regulators of placen tal development. However, the specific roles of FOSL1, JUNB. CITED2, and the other transcriptional regulators in the regulation of trophoblast differentiation are yet to be determined. Some may participate in the regulation or maintenance of the differentiated tropho blast cell phenotype. There is a connection between the differentiation associated genes and the PI3K AKT signaling pathway. As trophoblast stem cells differentiate, the PI3K AKT signaling pathway becomes constitutively activated.

IGF2 and GRN are candidate autocrine activators of the PI3K AKT signaling pathway. Trb3 and Msn were also classified as differentiation associated genes. They encode proteins with potential roles downstream of PI3K AKT signaling pathway. PI3K signaling sensitive genes PI3K regulates the phenotype of differentiating tropho blast cells. Endoreduplication and or survival of tro phoblast giant cells are influenced by PI3K signaling. An active PI3K pathway favors trophoblast giant cells with lower ploidy levels. These cells may be more motile and phenotypically resemble midgestation trophoblast lining uterine spiral arteries. PI3K signaling also possesses dramatic effects on gene expression patterns. Overall, the functions of the PI3K sensitive genes are biologically less diverse.

Most interestingly, they include genes encoding proteins potentially impacting tropho blast invasion, directed to the maternal uterine environment influencing immune and vascular cells, and also regulating androgen bio synthesis. Cgm4 is one of the most abundant genes expressed by differentiating trophoblast cells. It encodes a member of the expanded pregnancy specific glycoprotein family called PSG16. PSGs act on immune cells, poten tially through CD9, to influence cytokine production, they also target the vasculature and modulate endothelial cell function. The presence of Cd9 in differentiating trophoblast cells implies that PSGs may also possess autocrine paracrine actions on trophoblast development, which may include regulating the tropho blast invasive phenotype.

FAS ligand, PRL like protein A, adrenomedullin, and interleukin 17f are cytokines produced by differentiating tro phoblast that are exquisitely sensitive to PI3K regulation. FASLG binds to the FAS receptor and can initiate cell death. Trophoblast derived FASLG has been implicated as a modulator of intraplacental immune Cilengitide cell trafficking and is hypothesized to be a key participant in uterine spiral arteriole remodeling.

Subsequently, protein concentrations were determined for the frac

Subsequently, protein concentrations were determined for the fractions using the BCA protein assay. Equivalent amounts of protein were loaded in each lane. Following trans fer to nitrocellulose, the blots http://www.selleckchem.com/products/INCB18424.html were first stained with Pon ceau Red to assess protein loading, and subsequently probed with the following antibodies phospho Hsp27S15 and Hsp27. Blots were cut and reprobed sequentially, and visualized with ECL reagents and exposure to X ray film. Developed films were subsequently digitized and densitometrically analyzed with a cyclone ChemiImager and AlphaEase software. Digital images of the blots were used to make composite figures with Adobe Photoshop graphics software. Background Gliomas, the most common brain tumor, are currently classified as astrocytic, ependymal, oligodendroglial and choroid plexus tumors.

Among astrocytic tumors, gliob lastoma is the most lethal primary malignant brain tumor. Although considerable progress has been made in its treatment, the clinical prognosis associated with this tumor remains poor. Histone deacetylases have recently become rec ognized as a promising target for cancer therapy, includ ing for the treatment of glioblastomas. Together with histone acetyltransferases, HDACs are responsible for chromatin packaging, which influences the transcrip tion process. In general, increased levels of acetylation are associated with increased transcrip tional activity, whereas decreased acetylation levels are associated with repression of transcrip tion. HDACs are classified into 4 major categories based on their homology to yeast HDACs, including structure and cellular localization.

Class I and class II HDAC proteins share a common enzy matic mechanism that is the Zn catalyzed hydrolysis of the acetyl lysine amide bond. Human class I HDACs includes HDAC1, 2, 3, and 8, which are enzymes simi lar to the yeast transcriptional regulator Rpd3, generally localized to the nucleus. These enzymes are ubiqui tously expressed and seems to function as a complex with other proteins. HDAC1 and 2 only show activity within a protein complex, which consists of proteins necessary for modulating their deacetylase activity and DNA binding, and the recruit ment of HDACs to gene promoters. Wilson AJ et al. have suggested that multiple class I HDAC members are also involved in repressing p21 and that the growth inhib itory and apoptotic effects induced by HDAC inhibitors are probably mediated through the inhibition of multiple HDACs.

Class II HDACs includes HDAC4, 5, 6, 7, 9a, Brefeldin_A 9b, and 10, which are homologous to yeast Hda1. These class II enzymes can be found in the nucleus and cytoplasm, sug gesting potential extranuclear functions by regulating the acetylation status of nonhistone substrates. HDAC members of class II are abundantly expressed in skeletal muscle, heart, brain, tissues with low levels of mitotic activity.

We found that the results were rela tively insensitive to the cut

We found that the results were rela tively insensitive to the cut off value and we set this to Palbociclib solubility be 10% of the average expression value. All sample expression profiles within a series were scaled to the, where sk is the expression level of the kth probe set database to be searchable with cross platform response profiles and gene lists it has to be rewritten as a data base of expression profiles over non redundant gene lists. The EF profiles across the probe sets were there fore mapped onto expression profiles for a non redun dant gene list. In general each gene is represented by multiple probe sets. For each platform we generated the EF statistics for each probe set across the totality of samples. The probe set with the most robust response across the samples was chosen to represent the gene.

Explicitly, the probe set with the highest root mean square deviation form zero was chosen to represent the given gene. The number of genes defined on each plat form were as follows GPL96 11,807, GPL570 15,983 genes, GPL1261 13,202 genes, GPL85 chip with 3,844 genes, GPL1355 chip with 6,341 genes. The database totals 106,101 samples and is searchable on a reasonably fast desktop PC in 10 minutes per query. Searching the database The query profile is a statistically thresholded non redun dant list of genes and associated fold values. Statistical significance is assigned to a fold change based on a sim ple Students t test between multiple control and treat ment sample expression values. This is compared to each profile in the database by means of a simple Pearson regression analysis, with a correlation coefficient r.

The experiments are ranked according to the significance. The significance is measured by scaling the correlation to the normal by a Fisher transformation and measuring the number of standard deviations from the mean. The tion coefficient and N is the number of genes making up the correlation. The final ranking score is CMAP combined profiles The CMAP contains ranked lists of probes for 6,100 separate perturbagen treatments of four different human cell lines, with the ranking based on response level rela tive to control. The treatments are various multiples of 1,306 different drug like compounds. To generate responder sets that can be used to search SPIED we combined rankings for each separate compound treat ment and converted these into pseudo fold values with associated statistics.

The pseudo fold value is defined by gene and min max are the minimal maximal ranks. Remembering that the highest rank corresponds to the most up regulated gene. The SPIED was searched with CMAP profiles corresponding Drug_discovery to folds with a p 0. 05 threshold and with at least three replicates. This left 1,218 separate perturbagen probes. We sought to cluster the perturbagens based on predicted target and response profile similarity.

NF B is a protein transcription factor that functions to enhance

NF B is a protein transcription factor that functions to enhance the transcription of www.selleckchem.com/products/PF-2341066.html a variety of genes, including cytokines and growth factors, adhesion molecules, immunoreceptors, and acute phase proteins. Upon activation by LPS, NF B is required for maximal tran scription of many cytokines, including tumor necrosis factor a, interleukin 1b, IL 6, and IL 8, which are thought to be important in the generation of ALI. These cytokines and chemokines contribute to the vigorous recruitment of neutrophils in lung. Therefore, ALI is substantially caused by excessive neutrophil and cytokine mediated inflammation. Despite advancement in understanding the pathophysiology of ALI ARDS and improved therapy methods, however, mortality rates of ALI ARDS are around 40%. Histone deacetylases regulate gene expres sion.

In general, inhibitors of HDACs result in a non specific increase in gene expression. Therefore, they are considered as a new class of therapeutic agents for the treatment of tumor. Agents such as trichostatin A or suberoylanilide hydroxamic acid induce differentiation and or apoptosis of transformed cells in vitro and inhibit tumor growth in vivo. An unexpected effect of HDAC inhibitors, however, was revealed by recent studies indicating that they are able to suppress transcription and reduce inflammatory cyto kines in models of autoimmune and inflammatory dis eases. Butyrate, a HDAC inhibitor, is a short chain fatty acid derived from bacterial metabolism of dietary fibers in the colon and produces cell cycle arrest, differ entiation and or apoptosis of colorectal cancer cells in vitro.

Previous study has shown that butyrate reduced inflammation in Crohns disease through NF B inhibition. To date, unfortunately, the protective role of HDAC inhibitors in ALI is not well characterized. The aims of this study were to investigate whether butyrate reduces inflammation in LPS induced ALI in mice and to determine whether the protective effect is produced by suppression of inflammatory cytokines production and NF B activation. Materials and methods Animals and Reagents Male BALB C mice weighing 20 25 g were purchased from the Animal Center of the Fourth Military Medical University. All animals were allowed to take food and tap water ad libitum. All procedures were in accordance with the Declaration of Helsinki of the World Medical Association.

The protocols were also approved by the Institutional Animal Care and Use Committee of the Fourth Military Medical University Tangdu Hospital. LPS and butyrate were obtained from Sigma Chemical Company, and were respectively dissolved in saline. Enzyme linked immunosorbent assay kits of TNF a, IL 1b, myeloperoxidase Anacetrapib and nitric oxide were purchased from R D Corporation. Antibodies speci fic for total NF B p65, Lamin B and b actin were obtained from the Wuhan Boster Biological Technology, Ltd.

RasGTP pull downs were performed as described previously Mechani

RasGTP pull downs were performed as described previously. Mechanistic model simulations MATLAB and the function ode15s was used to simulate a previously developed, ordinary differential equation based ERK cascade model, which is described in de tail in Tables 1 and 2. The function gamrnd was used to generate realizations of www.selleckchem.com/products/dorsomorphin-2hcl.html peak RasGTP, Raf, MEK, and ERK levels for individual cells in the stochastic simula tions according to the gamma distribution where N specifies a protein level, k is the shape param eter, and �� is the scale parameter. We specified the k parameter of each gamma distribution as 5. 4, as was measured for total ERK, assuming roughly similar expression regulation. Since the mean of a gamma distribution is equal to k��, the �� parameter of each gamma distribution was changed as needed to at tain the desired distribution mean.

To estimate the parameters for the RasGTP dynamics, which are described by a simple exponential rise and decay model, we used least squares optimization to ensure that desired initial magni tude, peak magnitude, time to peak, time to inflection, time to steady state, and steady state magnitude of the RasGTP dynamics matches well to that which the model prescribes. Additional file 1 Figure S6 describes these RasGTP dynamics metrics graphically. As there are four unknown parameters in the RasGTP dynam ics model, we need four equations, which we take as the following e5T where wi corresponds to a weight for optimization purposes. Eq. 2 specifies the proper steady state magnitude. Eq. 3 specifies that the 1st derivative at the time to peak is zero.

Eq. 4 specifies the proper magnitude at the time to steady state . and Eq. 5 specifies the proper peak magnitude. The following constraints are placed on this optimization problem d? Eq. 6 specifies that there is a maximum at the time to peak and Eq. 7 specifies that the 1st derivative is negative at the inflection point Ks and s are determined as described in Methods. Mean peak RasGTP levels were increased to simu late increasing input, and were linearly spaced between 10 nM and 200 nM using 6 points, which correspond to EGF doses of 0. 01, 0. 1, 0. 5, 1, 5, and 10. Following the trends of the experimental data in Additional file 1 Figure 2A and, peak times for RasGTP were sampled linearly between 7 min and 2 min, with 7 min corresponding to the lowest peak RasGTP level.

Also, we took ss as 10 min, Iss as 15% of Imax realizations, Io as 0, and infl as /2. All code is available upon request. Parameter Batimastat sensitivity analysis Five hundred different parameter sets were generated via latin hypercube sampling from a 23 dimensional uniform distribution http://www.selleckchem.com/products/Rapamycin.html that spans 1 order of magnitude around each nominal parameter value. For each of these parameter sets stochastic simu lations were performed as described above.