Simulating Quantitative Cellular Responses Using Asynchronous Threshold Boolean Network Ensembles

<p>Abstract</p> <p>Background</p> <p>With increasing knowledge about the potential mechanisms underlying cellular functions, it is becoming feasible to predict the response of biological systems to genetic and environmental perturbations. Due to the lack of homogeneity in living tissues it is difficult to estimate the physiological effect of chemicals, including potential toxicity. Here we investigate a biologically motivated model for estimating tissue level responses by aggregating the behavior of a cell population. We assume that the molecular state of individual cells is independently governed by discrete non-deterministic signaling mechanisms. This results in noisy but highly reproducible aggregate level responses that are consistent with experimental data.</p> <p>Results</p> <p>We developed an asynchronous threshold Boolean network simulation algorithm to model signal transduction in a single cell, and then used an ensemble of these models to estimate the aggregate response across a cell population. Using published data, we derived a putative crosstalk network involving growth factors and cytokines - i.e., Epidermal Growth Factor, Insulin, Insulin like Growth Factor Type 1, and Tumor Necrosis Factor α - to describe early signaling events in cell proliferation signal transduction. Reproducibility of the modeling technique across ensembles of Boolean networks representing cell populations is investigated. Furthermore, we compare our simulation results to experimental observations of hepatocytes reported in the literature.</p> <p>Conclusion</p> <p>A systematic analysis of the results following differential stimulation of this model by growth factors and cytokines suggests that: (a) using Boolean network ensembles with asynchronous updating provides biologically plausible noisy individual cellular responses with reproducible mean behavior for large cell populations, and (b) with sufficient data our model can estimate the response to different concentrations of extracellular ligands. Our results suggest that this approach is both quantitative, allowing statistical verification and calibration, and extensible, allowing modification and revision as guided by experimental evidence. The simulation methodology is part of the US EPA Virtual Liver, which is investigating the effects of everyday contaminants on living tissues. Future models will incorporate additional crosstalk surrounding proliferation as well as the putative effects of xenobiotics on these signaling cascades within hepatocytes.</p

httk: R Package for High-Throughput Toxicokinetics

Thousands of chemicals have been profiled by high-throughput screening programs such as ToxCast and Tox21; these chemicals are tested in part because most of them have limited or no data on hazard, exposure, or toxicokinetics. Toxicokinetic models aid in predicting tissue concentrations resulting from chemical exposure, and a "reverse dosimetry" approach can be used to predict exposure doses sufficient to cause tissue concentrations that have been identified as bioactive by high-throughput screening. We have created four toxicokinetic models within the R software package httk. These models are designed to be parameterized using high-throughput in vitro data (plasma protein binding and hepatic clearance), as well as structure-derived physicochemical properties and species-specific physiological data. The package contains tools for Monte Carlo sampling and reverse dosimetry along with functions for the analysis of concentration vs. time simulations. The package can currently use human in vitro data to make predictions for 553 chemicals in humans, rats, mice, dogs, and rabbits, including 94 pharmaceuticals and 415 ToxCast chemicals. For 67 of these chemicals, the package includes rat-specific in vitro data. This package is structured to be augmented with additional chemical data as they become available. Package httk enables the inclusion of toxicokinetics in the statistical analysis of chemicals undergoing high-throughput screening

Hubungan Penggunaan Dan Penanganan Pestisida Pada Petani Bawang Merah Terhadap Residu Pestisida Dalam Tanah Di Lahan Pertanian Desa Wanasari Kecamatan Wanasari Kabupaten Brebes

Excessive use of pesticides causing pollution and environmental damage agriculture. Examination in Brebes on 31 samples of fruits and vegetables, found 22% of samples contain detectable residues of organophosphate and found two soil samples (10%) contained residues organochlorin. The purpose of this study was to determine the relationship of the use and handling of pesticides on their onion farmers against pesticide residues in the soil on agricultural land Wanasari Village, District Wanasari, Brebes. This study is observational method with cross sectional approach. The population in this study were all farmers in the Wanasari conducting spraying. Collecting data using the tool Banu questionnaire and examination of pesticide residues in soil using GC-MS Gas Chromatography - Mass Spectrometry. The results of this study are of 55 69.1 onion farmers use pesticides are not good. The use of pesticides covering 80% is not good in mixing pesticides, 87.3% use a smaller dose, 49.1% use pesticides that are not registered with the Ministry of Agriculture, 87.3% is not good in the way of spraying and 87.3 does well in frequency spraying. Handling pesticides in agricultural land is not good 59.1%, ie 74.5% is not good in handling pesticide containers, 90.9% is not good in storage of pesticides, 89.1% is not good in handling a spill and 87.3% did not either in place to clean pesticide containers. The research result is negative soil samples pesticide residues. The conclusion was that no pesticide residue class organochlorin

Exploring consumer exposure pathways and patterns of use for chemicals in the environment

AbstractHumans are exposed to thousands of chemicals in the workplace, home, and via air, water, food, and soil. A major challenge in estimating chemical exposures is to understand which chemicals are present in these media and microenvironments. Here we describe the Chemical/Product Categories Database (CPCat), a new, publically available (http://actor.epa.gov/cpcat) database of information on chemicals mapped to “use categories” describing the usage or function of the chemical. CPCat was created by combining multiple and diverse sources of data on consumer- and industrial-process based chemical uses from regulatory agencies, manufacturers, and retailers in various countries. The database uses a controlled vocabulary of 833 terms and a novel nomenclature to capture and streamline descriptors of chemical use for 43,596 chemicals from the various sources. Examples of potential applications of CPCat are provided, including identifying chemicals to which children may be exposed and to support prioritization of chemicals for toxicity screening. CPCat is expected to be a valuable resource for regulators, risk assessors, and exposure scientists to identify potential sources of human exposures and exposure pathways, particularly for use in high-throughput chemical exposure assessment

Toxicokinetic Triage for Environmental Chemicals

Toxicokinetic (TK) models link administered doses to plasma, blood, and tissue concentrations. High-throughput TK (HTTK) performs in vitro to in vivo extrapolation to predict TK from rapid in vitro measurements and chemical structure-based properties. A significant toxicological application of HTTK has been “reverse dosimetry,” in which bioactive concentrations from in vitro screening studies are converted into in vivo doses (mg/kg BW/day). These doses are predicted to produce steady-state plasma concentrations that are equivalent to in vitro bioactive concentrations. In this study, we evaluate the impact of the approximations and assumptions necessary for reverse dosimetry and develop methods to determine whether HTTK tools are appropriate or may lead to false conclusions for a particular chemical. Based on literature in vivo data for 87 chemicals, we identified specific properties (eg, in vitro HTTK data, physico-chemical descriptors, and predicted transporter affinities) that correlate with poor HTTK predictive ability. For 271 chemicals we developed a generic HT physiologically based TK (HTPBTK) model that predicts non-steady-state chemical concentration time-courses for a variety of exposure scenarios. We used this HTPBTK model to find that assumptions previously used for reverse dosimetry are usually appropriate, except most notably for highly bioaccumulative compounds. For the thousands of man-made chemicals in the environment that currently have no TK data, we propose a 4-element framework for chemical TK triage that can group chemicals into 7 different categories associated with varying levels of confidence in HTTK predictions. For 349 chemicals with literature HTTK data, we differentiated those chemicals for which HTTK approaches are likely to be sufficient, from those that may require additional data

Identifying xenobiotic metabolites with in silico prediction tools and LCMS suspect screening analysis

Understanding the metabolic fate of a xenobiotic substance can help inform its potential health risks and allow for the identification of signature metabolites associated with exposure. The need to characterize metabolites of poorly studied or novel substances has shifted exposure studies towards non-targeted analysis (NTA), which often aims to profile many compounds within a sample using high-resolution liquid-chromatography mass-spectrometry (LCMS). Here we evaluate the suitability of suspect screening analysis (SSA) liquid-chromatography mass-spectrometry to inform xenobiotic chemical metabolism. Given a lack of knowledge of true metabolites for most chemicals, predictive tools were used to generate potential metabolites as suspect screening lists to guide the identification of selected xenobiotic substances and their associated metabolites. Thirty-three substances were selected to represent a diverse array of pharmaceutical, agrochemical, and industrial chemicals from Environmental Protection Agency’s ToxCast chemical library. The compounds were incubated in a metabolically-active in vitro assay using primary hepatocytes and the resulting supernatant and lysate fractions were analyzed with high-resolution LCMS. Metabolites were simulated for each compound structure using software and then combined to serve as the suspect screening list. The exact masses of the predicted metabolites were then used to select LCMS features for fragmentation via tandem mass spectrometry (MS/MS). Of the starting chemicals, 12 were measured in at least one sample in either positive or negative ion mode and a subset of these were used to develop the analysis workflow. We implemented a screening level workflow for background subtraction and the incorporation of time-varying kinetics into the identification of likely metabolites. We used haloperidol as a case study to perform an in-depth analysis, which resulted in identifying five known metabolites and five molecular features that represent potential novel metabolites, two of which were assigned discrete structures based on in silico predictions. This workflow was applied to five additional test chemicals, and 15 molecular features were selected as either reported metabolites, predicted metabolites, or potential metabolites without a structural assignment. This study demonstrates that in some–but not all–cases, suspect screening analysis methods provide a means to rapidly identify and characterize metabolites of xenobiotic chemicals

Simulating Microdosimetry in a Virtual Hepatic Lobule

The liver plays a key role in removing harmful chemicals from the body and is therefore often the first tissue to suffer potentially adverse consequences. To protect public health it is necessary to quantitatively estimate the risk of long-term low dose exposure to environmental pollutants. Animal testing is the primary tool for extrapolating human risk but it is fraught with uncertainty, necessitating novel alternative approaches. Our goal is to integrate in vitro liver experiments with agent-based cellular models to simulate a spatially extended hepatic lobule. Here we describe a graphical model of the sinusoidal network that efficiently simulates portal to centrilobular mass transfer in the hepatic lobule. We analyzed the effects of vascular topology and metabolism on the cell-level distribution following oral exposure to chemicals. The spatial distribution of metabolically inactive chemicals was similar across different vascular networks and a baseline well-mixed compartment. When chemicals were rapidly metabolized, concentration heterogeneity of the parent compound increased across the vascular network. As a result, our spatially extended lobule generated greater variability in dose-dependent cellular responses, in this case apoptosis, than were observed in the classical well-mixed liver or in a parallel tubes model. The mass-balanced graphical approach to modeling the hepatic lobule is computationally efficient for simulating long-term exposure, modular for incorporating complex cellular interactions, and flexible for dealing with evolving tissues