RuleMonkey: software for stochastic simulation of rule-based models

<p>Abstract</p> <p>Background</p> <p>The system-level dynamics of many molecular interactions, particularly protein-protein interactions, can be conveniently represented using reaction rules, which can be specified using model-specification languages, such as the BioNetGen language (BNGL). A set of rules implicitly defines a (bio)chemical reaction network. The reaction network implied by a set of rules is often very large, and as a result, generation of the network implied by rules tends to be computationally expensive. Moreover, the cost of many commonly used methods for simulating network dynamics is a function of network size. Together these factors have limited application of the rule-based modeling approach. Recently, several methods for simulating rule-based models have been developed that avoid the expensive step of network generation. The cost of these "network-free" simulation methods is independent of the number of reactions implied by rules. Software implementing such methods is now needed for the simulation and analysis of rule-based models of biochemical systems.</p> <p>Results</p> <p>Here, we present a software tool called RuleMonkey, which implements a network-free method for simulation of rule-based models that is similar to Gillespie's method. The method is suitable for rule-based models that can be encoded in BNGL, including models with rules that have global application conditions, such as rules for intramolecular association reactions. In addition, the method is rejection free, unlike other network-free methods that introduce null events, i.e., steps in the simulation procedure that do not change the state of the reaction system being simulated. We verify that RuleMonkey produces correct simulation results, and we compare its performance against DYNSTOC, another BNGL-compliant tool for network-free simulation of rule-based models. We also compare RuleMonkey against problem-specific codes implementing network-free simulation methods.</p> <p>Conclusions</p> <p>RuleMonkey enables the simulation of rule-based models for which the underlying reaction networks are large. It is typically faster than DYNSTOC for benchmark problems that we have examined. RuleMonkey is freely available as a stand-alone application <url>http://public.tgen.org/rulemonkey</url>. It is also available as a simulation engine within GetBonNie, a web-based environment for building, analyzing and sharing rule-based models.</p

MRI radiomic features are independently associated with overall survival in soft tissue sarcoma

Purpose: Soft tissue sarcomas (STS) represent a heterogeneous group of diseases, and selection of individualized treatments remains a challenge. The goal of this study was to determine whether radiomic features extracted from magnetic resonance (MR) images are independently associated with overall survival (OS) in STS. Methods and Materials: This study analyzed 2 independent cohorts of adult patients with stage II-III STS treated at center 1 (N = 165) and center 2 (N = 61). Thirty radiomic features were extracted from pretreatment T1-weighted contrast-enhanced MR images. Prognostic models for OS were derived on the center 1 cohort and validated on the center 2 cohort. Clinical-only (C), radiomics-only (R), and clinical and radiomics (C+R) penalized Cox models were constructed. Model performance was assessed using Harrell\u27s concordance index. Results: In the R model, tumor volume (hazard ratio [HR], 1.5) and 4 texture features (HR, 1.1-1.5) were selected. In the C+R model, both age (HR, 1.4) and grade (HR, 1.7) were selected along with 5 radiomic features. The adjusted c-indices of the 3 models ranged from 0.68 (C) to 0.74 (C+R) in the derivation cohort and 0.68 (R) to 0.78 (C+R) in the validation cohort. The radiomic features were independently associated with OS in the validation cohort after accounting for age and grade (HR, 2.4; Conclusions: This study found that radiomic features extracted from MR images are independently associated with OS when accounting for age and tumor grade. The overall predictive performance of 3-year OS using a model based on clinical and radiomic features was replicated in an independent cohort. Optimal models using clinical and radiomic features could improve personalized selection of therapy in patients with STS

Management of anteriorly located C1-C2 neurofibromata

The authors discuss their recent experience with anteriorly located Cl-C2 neurofibromata in five patients with cervical myelopathy and magnetic resonance scans consistent with intradural extramedullary masses in this region. Surgery was performed using a posterolateral approach with microscopic intradural exploration. Gross total intradural tumor removal was achieved in all cases. Improvement in cervical myelopathy occurred in all patients.This report concludes that Cl-C2 neurofibromata located anterior to the spinal cord can be totally and safely removed using a posterolateral approach. Improvement in neurologic dysfunction accompanies posterior decompression and gross total intradural tumor removal.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/28766/1/0000598.pd

Outcome of ATP-based tumor chemosensitivity assay directed chemotherapy in heavily pre-treated recurrent ovarian carcinoma

BACKGROUND: We wished to evaluate the clinical response following ATP-Tumor Chemosensitivity Assay (ATP-TCA) directed salvage chemotherapy in a series of UK patients with advanced ovarian cancer. The results are compared with that of a similar assay used in a different country in terms of evaluability and clinical endpoints. METHODS: From November 1998 to November 2001, 46 patients with pre-treated, advanced ovarian cancer were given a total of 56 courses of chemotherapy based on in-vitro ATP-TCA responses obtained from fresh tumor samples or ascites. Forty-four patients were evaluable for results. Of these, 18 patients had clinically platinum resistant disease (relapse < 6 months after first course of chemotherapy). There was evidence of cisplatin resistance in 31 patients from their first ATP-TCA. Response to treatment was assessed by radiology, clinical assessment and tumor marker level (CA 125). RESULTS: The overall response rate was 59% (33/56) per course of chemotherapy, including 12 complete responses, 21 partial responses, 6 with stable disease, and 15 with progressive disease. Two patients were not evaluable for response having received just one cycle of chemotherapy: if these were excluded the response rate is 61%. Fifteen patients are still alive. Median progression free survival (PFS) was 6.6 months per course of chemotherapy; median overall survival (OAS) for each patient following the start of TCA-directed therapy was 10.4 months (95% confidence interval 7.9-12.8 months). CONCLUSION: The results show similar response rates to previous studies using ATP-TCA directed therapy in recurrent ovarian cancer. The assay shows high evaluability and this study adds weight to the reproducibility of results from different centre

Spatial Guilds in the Serengeti Food Web Revealed by a Bayesian Group Model

Food webs, networks of feeding relationships among organisms, provide fundamental insights into mechanisms that determine ecosystem stability and persistence. Despite long-standing interest in the compartmental structure of food webs, past network analyses of food webs have been constrained by a standard definition of compartments, or modules, that requires many links within compartments and few links between them. Empirical analyses have been further limited by low-resolution data for primary producers. In this paper, we present a Bayesian computational method for identifying group structure in food webs using a flexible definition of a group that can describe both functional roles and standard compartments. The Serengeti ecosystem provides an opportunity to examine structure in a newly compiled food web that includes species-level resolution among plants, allowing us to address whether groups in the food web correspond to tightly-connected compartments or functional groups, and whether network structure reflects spatial or trophic organization, or a combination of the two. We have compiled the major mammalian and plant components of the Serengeti food web from published literature, and we infer its group structure using our method. We find that network structure corresponds to spatially distinct plant groups coupled at higher trophic levels by groups of herbivores, which are in turn coupled by carnivore groups. Thus the group structure of the Serengeti web represents a mixture of trophic guild structure and spatial patterns, in contrast to the standard compartments typically identified in ecological networks. From data consisting only of nodes and links, the group structure that emerges supports recent ideas on spatial coupling and energy channels in ecosystems that have been proposed as important for persistence.Comment: 28 pages, 6 figures (+ 3 supporting), 2 tables (+ 4 supporting

Clinical Profiling of BCL-2 Family Members in the Setting of BRAF Inhibition Offers a Rationale for Targeting De Novo Resistance Using BH3 Mimetics

While response rates to BRAF inhibitiors (BRAFi) are high, disease progression emerges quickly. One strategy to delay the onset of resistance is to target anti-apoptotic proteins such as BCL-2, known to be associated with a poor prognosis. We analyzed BCL-2 family member expression levels of 34 samples from 17 patients collected before and 10 to 14 days after treatment initiation with either vemurafenib or dabrafenib/trametinib combination. The observed changes in mRNA and protein levels with BRAFi treatment led us to hypothesize that combining BRAFi with a BCL-2 inhibitor (the BH3-mimetic navitoclax) would improve outcome. We tested this hypothesis in cell lines and in mice. Pretreatment mRNA levels of BCL-2 negatively correlated with maximal tumor regression. Early increases in mRNA levels were seen in BIM, BCL-XL, BID and BCL2-W, as were decreases in MCL-1 and BCL2A. No significant changes were observed with BCL-2. Using reverse phase protein array (RPPA), significant increases in protein levels were found in BIM and BID. No changes in mRNA or protein correlated with response. Concurrent BRAF (PLX4720) and BCL2 (navitoclax) inhibition synergistically reduced viability in BRAF mutant cell lines and correlated with down-modulation of MCL-1 and BIM induction after PLX4720 treatment. In xenograft models, navitoclax enhanced the efficacy of PLX4720. The combination of a selective BRAF inhibitor with a BH3-mimetic promises to be an important therapeutic strategy capable of enhancing the clinical efficacy of BRAF inhibition in many patients that might otherwise succumb quickly to de novo resistance. Trial Registrations: ClinicalTrials.gov NCT01006980; ClinicalTrials.gov NCT01107418; ClinicalTrials.gov NCT01264380; ClinicalTrials.gov NCT01248936; ClinicalTrials.gov NCT00949702; ClinicalTrials.gov NCT0107217

Redesign of plates by large admissible perturbations

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76482/1/AIAA-12089-360.pd

Stratospheric AOD after the 2011 eruption of Nabro volcano measured by lidars over the Northern Hemisphere

International audienceNabro volcano (13.37°N, 41.70°E) in Eritrea erupted on 13 June 2011 generating a layer of sulfate aerosols that persisted in the stratosphere for months. For the first time we report on ground-based lidar observations of the same event from every continent in the Northern Hemisphere, taking advantage of the synergy between global lidar networks such as EARLINET, MPLNET and NDACC with independent lidar groups and satellite CALIPSO to track the evolution of the stratospheric aerosol layer in various parts of the globe. The globally averaged aerosol optical depth (AOD) due to the stratospheric volcanic aerosol layers was of the order of 0.018 ± 0.009 at 532 nm, ranging from 0.003 to 0.04. Compared to the total column AOD from the available collocated AERONET stations, the stratospheric contribution varied from 2% to 23% at 532 nm

DNA methylation-based classification of central nervous system tumours.

Accurate pathological diagnosis is crucial for optimal management of patients with cancer. For the approximately 100 known tumour types of the central nervous system, standardization of the diagnostic process has been shown to be particularly challenging-with substantial inter-observer variability in the histopathological diagnosis of many tumour types. Here we present a comprehensive approach for the DNA methylation-based classification of central nervous system tumours across all entities and age groups, and demonstrate its application in a routine diagnostic setting. We show that the availability of this method may have a substantial impact on diagnostic precision compared to standard methods, resulting in a change of diagnosis in up to 12% of prospective cases. For broader accessibility, we have designed a free online classifier tool, the use of which does not require any additional onsite data processing. Our results provide a blueprint for the generation of machine-learning-based tumour classifiers across other cancer entities, with the potential to fundamentally transform tumour pathology

Rib fracture after stereotactic radiotherapy on follow-up thin-section computed tomography in 177 primary lung cancer patients

<p>Abstract</p> <p>Background</p> <p>Chest wall injury after stereotactic radiotherapy (SRT) for primary lung cancer has recently been reported. However, its detailed imaging findings are not clarified. So this study aimed to fully characterize the findings on computed tomography (CT), appearance time and frequency of chest wall injury after stereotactic radiotherapy (SRT) for primary lung cancer</p> <p>Materials and methods</p> <p>A total of 177 patients who had undergone SRT were prospectively evaluated for periodical follow-up thin-section CT with special attention to chest wall injury. The time at which CT findings of chest wall injury appeared was assessed. Related clinical symptoms were also evaluated.</p> <p>Results</p> <p>Rib fracture was identified on follow-up CT in 41 patients (23.2%). Rib fractures appeared at a mean of 21.2 months after the completion of SRT (range, 4 -58 months). Chest wall edema, thinning of the cortex and osteosclerosis were findings frequently associated with, and tending to precede rib fractures. No patients with rib fracture showed tumors > 16 mm from the adjacent chest wall. Chest wall pain was seen in 18 of 177 patients (10.2%), of whom 14 patients developed rib fracture. No patients complained of Grade 3 or more symptoms.</p> <p>Conclusion</p> <p>Rib fracture is frequently seen after SRT for lung cancer on CT, and is often associated with chest wall edema, thinning of the cortex and osteosclerosis. However, related chest wall pain is less frequent and is generally mild if present.</p