A phase 2 trial of personalized cytotoxic therapy based on tumor immunohistochemistry in previously treated metastatic pancreatic cancer patients

Background

The choice of a regimen in metastatic pancreatic cancer patients following progression on 1st line therapy is empiric and outcomes are unsatisfactory. This phase II study was performed to evaluate the efficacy of therapy selected by immunohistochemistry (IHC) in these patients following progression after one or more therapies.

Methods

Eligible patients underwent a percutaneous biopsy of a metastatic lesion and treatment selection was determined by IHC. The study required 35 evaluable patients (power of 86%) for detecting a true 1-year survival rate of >20%.

Results

A tumor biopsy was performed in 48 of 49 accrued patients. Study therapy was not given (n=13) either due to insufficient tumor on biopsy (n=8) or due to worsening cancer related symptoms after biopsy (n=5). The demographics of evaluable patients (n=35) are male/female (59%/41%), with age range 34-78 years (median 63 years). Patients had 1-6 prior regimens (median of 2). The most common IHC targets were topoisomerase 1 or 2, thymidylate synthase, excision repair cross-complementation group 1 protein (ERCC1), and osteonectin secreted protein acidic and rich in cysteine (SPARC). Commercially available treatment regimens prescribed included FOLFIRI, FOLFOX, irinotecan, and doxorubicin. The response (RECIST) was 9%, the median survival was 5.6 months (94% CI, 3.8-8.2), and the 1-year survival was 20% (95% CI, 7-33%).

Conclusions

In all patients, IHC assays resulted in identification of at least two targets for therapy and a non-cross resistant regimen could be prescribed for therapy with evidence of some benefit. An IHC based treatment strategy is feasible and needs validation in larger studies.

Differential mast cell outcomes are sensitive to FcεRI-Syk binding kinetics

Cross-linking of immunoglobulin E-bound FcεRI triggers multiple cellular responses, including degranulation and cytokine production. Signaling is dependent on recruitment of Syk via docking of its dual SH2 domains to phosphorylated tyrosines within the FcεRI immunoreceptor tyrosine-based activation motifs. Using single-molecule imaging in live cells, we directly visualized and quantified the binding of individual mNeonGreen-tagged Syk molecules as they associated with the plasma membrane after FcεRI activation. We found that Syk colocalizes transiently to FcεRI and that Syk-FcεRI binding dynamics are independent of receptor aggregate size. Substitution of glutamic acid for tyrosine between the Syk SH2 domains (Syk-Y130E) led to an increased Syk-FcεRI off-rate, loss of site-specific Syk autophosphorylation, and impaired downstream signaling. Genome edited cells expressing only Syk-Y130E were deficient in antigen-stimulated calcium release, degranulation, and production of some cytokines (TNF-a, IL-3) but not others (MCP-1, IL-4). We propose that kinetic discrimination along the FcεRI signaling pathway occurs at the level of Syk-FcεRI interactions, with key outcomes dependent upon sufficiently long-lived Syk binding events.

Allergen Valency, Dose, and FcεRI Occupancy Set Thresholds for Secretory Responses to Pen a 1 and Motivate Design of Hypoallergens

Ag-mediated crosslinking of IgE-FcεRI complexes activates mast cells and basophils, initiating the allergic response. Of 34 donors recruited having self-reported shrimp allergy, only 35% had significant levels of shrimp-specific IgE in serum and measurable basophil secretory responses to rPen a 1 (shrimp tropomyosin). We report that degranulation is linked to the number of FcεRI occupied with allergen-specific IgE, as well as the dose and valency of Pen a 1. Using clustered regularly interspaced palindromic repeat-based gene editing, human RBL^rαKO cells were created that exclusively express the human FcεRIα subunit. Pen a 1-specific IgE was affinity purified from shrimp-positive plasma. Cells primed with a range of Pen a 1-specific IgE and challenged with Pen a 1 showed a bell-shaped dose response for secretion, with optimal Pen a 1 doses of 0.1-10 ng/ml. Mathematical modeling provided estimates of receptor aggregation kinetics based on FcεRI occupancy with IgE and allergen dose. Maximal degranulation was elicited when ∼2700 IgE-FcεRI complexes were occupied with specific IgE and challenged with Pen a 1 (IgE epitope valency of ≥8), although measurable responses were achieved when only a few hundred FcεRI were occupied. Prolonged periods of pepsin-mediated Pen a 1 proteolysis, which simulates gastric digestion, were required to diminish secretory responses. Recombinant fragments (60-79 aa), which together span the entire length of tropomyosin, were weak secretagogues. These fragments have reduced dimerization capacity, compete with intact Pen a 1 for binding to IgE-FcεRI complexes, and represent a starting point for the design of promising hypoallergens for immunotherapy.

BioNetFit: a fitting tool compatible with BioNetGen, NFsim and distributed computing environments

Rule-based models are analyzed with specialized simulators, such as those provided by the BioNetGen and NFsim open-source software packages. Here, we present BioNetFit, a general-purpose fitting tool that is compatible with BioNetGen and NFsim. BioNetFit is designed to take advantage of distributed computing resources. This feature facilitates fitting (i.e. optimization of parameter values for consistency with data) when simulations are computationally expensive.

AVAILABILITY AND IMPLEMENTATION

BioNetFit can be used on stand-alone Mac, Windows/Cygwin, and Linux platforms and on Linux-based clusters running SLURM, Torque/PBS, or SGE. The BioNetFit source code (Perl) is freely available (http://bionetfit.nau.edu).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CONTACT

bionetgen.help@gmail.com.

Use of mechanistic models to integrate and analyze multiple proteomic datasets

Proteins in cell signaling networks tend to interact promiscuously through low-affinity interactions. Consequently, evaluating the physiological importance of mapped interactions can be difficult. Attempts to do so have tended to focus on single, measurable physicochemical factors, such as affinity or abundance. For example, interaction importance has been assessed on the basis of the relative affinities of binding partners for a protein of interest, such as a receptor. However, multiple factors can be expected to simultaneously influence the recruitment of proteins to a receptor (and the potential of these proteins to contribute to receptor signaling), including affinity, abundance, and competition, which is a network property. Here, we demonstrate that measurements of protein copy numbers and binding affinities can be integrated within the framework of a mechanistic, computational model that accounts for mass action and competition. We use cell line-specific models to rank the relative importance of protein-protein interactions in the epidermal growth factor receptor (EGFR) signaling network for 11 different cell lines. Each model accounts for experimentally characterized interactions of six autophosphorylation sites in EGFR with proteins containing a Src homology 2 and/or phosphotyrosine-binding domain. We measure importance as the predicted maximal extent of recruitment of a protein to EGFR following ligand-stimulated activation of EGFR signaling. We find that interactions ranked highly by this metric include experimentally detected interactions. Proteins with high importance rank in multiple cell lines include proteins with recognized, well-characterized roles in EGFR signaling, such as GRB2 and SHC1, as well as a protein with a less well-defined role, YES1. Our results reveal potential cell line-specific differences in recruitment.

Enhanced dimerization drives ligand-independent activity of mutant epidermal growth factor receptor in lung cancer

Epidermal growth factor receptor kinase mutations drive oncogenesis, but the molecular mechanism of pathological signal initiation is poorly understood. Using high-resolution microscopy methods, the authors reveal that these kinase mutations induce structural changes in the receptor ectodomain that lead to enhanced, ligand-independent dimerization.

Mutations within the epidermal growth factor receptor (EGFR/erbB1/Her1) are often associated with tumorigenesis. In particular, a number of EGFR mutants that demonstrate ligand-independent signaling are common in non–small cell lung cancer (NSCLC), including kinase domain mutations L858R (also called L834R) and exon 19 deletions (e.g., ΔL747-P753insS), which collectively make up nearly 90% of mutations in NSCLC. The molecular mechanisms by which these mutations confer constitutive activity remain unresolved. Using multiple subdiffraction-limit imaging modalities, we reveal the altered receptor structure and interaction kinetics of NSCLC-associated EGFR mutants. We applied two-color single quantum dot tracking to quantify receptor dimerization kinetics on living cells and show that, in contrast to wild-type EGFR, mutants are capable of forming stable, ligand-independent dimers. Two-color superresolution localization microscopy confirmed ligand-independent aggregation of EGFR mutants. Live-cell Förster resonance energy transfer measurements revealed that the L858R kinase mutation alters ectodomain structure such that unliganded mutant EGFR adopts an extended, dimerization-competent conformation. Finally, mutation of the putative dimerization arm confirmed a critical role for ectodomain engagement in ligand-independent signaling. These data support a model in which dysregulated activity of NSCLC-associated kinase mutants is driven by coordinated interactions involving both the kinase and extracellular domains that lead to enhanced dimerization.

Modeling for (physical) biologists: an introduction to the rule-based approach

Models that capture the chemical kinetics of cellular regulatory networks can be specified in terms of rules for biomolecular interactions. A rule defines a generalized reaction, meaning a reaction that permits multiple reactants, each capable of participating in a characteristic transformation and each possessing certain, specified properties, which may be local, such as the state of a particular site or domain of a protein. In other words, a rule defines a transformation and the properties that reactants must possess to participate in the transformation. A rule also provides a rate law. A rule-based approach to modeling enables consideration of mechanistic details at the level of functional sites of biomolecules and provides a facile and visual means for constructing computational models, which can be analyzed to study how system-level behaviors emerge from component interactions.

Modeling biomolecular site dynamics in immunoreceptor signaling systems

The immune system plays a central role in human health. The activities of immune cells, whether defending an organism from disease or triggering a pathological condition such as autoimmunity, are driven by the molecular machinery of cellular signaling systems. Decades of experimentation have elucidated many of the biomolecules and interactions involved in immune signaling and regulation, and recently developed technologies have led to new types of quantitative, systems-level data. To integrate such information and develop nontrivial insights into the immune system, computational modeling is needed, and it is essential for modeling methods to keep pace with experimental advances. In this chapter, we focus on the dynamic, site-specific, and context-dependent nature of interactions in immunoreceptor signaling (i.e., the biomolecular site dynamics of immunoreceptor signaling), the challenges associated with capturing these details in computational models, and how these challenges have been met through use of rule-based modeling approaches.

Computational analysis of an autophagy/translation switch based on mutual inhibition of MTORC1 and ULK1

We constructed a mechanistic, computational model for regulation of (macro)autophagy and protein synthesis (at the level of translation). The model was formulated to study the system-level consequences of interactions among the following proteins: two key components of MTOR complex 1 (MTORC1), namely the protein kinase MTOR (mechanistic target of rapamycin) and the scaffold protein RPTOR; the autophagy-initiating protein kinase ULK1; and the multimeric energy-sensing AMP-activated protein kinase (AMPK). Inputs of the model include intrinsic AMPK kinase activity, which is taken as an adjustable surrogate parameter for cellular energy level or AMP:ATP ratio, and rapamycin dose, which controls MTORC1 activity. Outputs of the model include the phosphorylation level of the translational repressor EIF4EBP1, a substrate of MTORC1, and the phosphorylation level of AMBRA1 (activating molecule in BECN1-regulated autophagy), a substrate of ULK1 critical for autophagosome formation. The model incorporates reciprocal regulation of mTORC1 and ULK1 by AMPK, mutual inhibition of MTORC1 and ULK1, and ULK1-mediated negative feedback regulation of AMPK. Through analysis of the model, we find that these processes may be responsible, depending on conditions, for graded responses to stress inputs, for bistable switching between autophagy and protein synthesis, or relaxation oscillations, comprising alternating periods of autophagy and protein synthesis. A sensitivity analysis indicates that the prediction of oscillatory behavior is robust to changes of the parameter values of the model. The model provides testable predictions about the behavior of the AMPK-MTORC1-ULK1 network, which plays a central role in maintaining cellular energy and nutrient homeostasis.

Phosphorylation site dynamics of early T-cell receptor signaling

In adaptive immune responses, T-cell receptor (TCR) signaling impacts multiple cellular processes and results in T-cell differentiation, proliferation, and cytokine production. Although individual protein-protein interactions and phosphorylation events have been studied extensively, we lack a systems-level understanding of how these components cooperate to control signaling dynamics, especially during the crucial first seconds of stimulation. Here, we used quantitative proteomics to characterize reshaping of the T-cell phosphoproteome in response to TCR/CD28 co-stimulation, and found that diverse dynamic patterns emerge within seconds. We detected phosphorylation dynamics as early as 5 s and observed widespread regulation of key TCR signaling proteins by 30 s. Development of a computational model pointed to the presence of novel regulatory mechanisms controlling phosphorylation of sites with central roles in TCR signaling. The model was used to generate predictions suggesting unexpected roles for the phosphatase PTPN6 (SHP-1) and shortcut recruitment of the actin regulator WAS. Predictions were validated experimentally. This integration of proteomics and modeling illustrates a novel, generalizable framework for solidifying quantitative understanding of a signaling network and for elucidating missing links.

Optimal aggregation of FcεRI with a structurally defined trivalent ligand overrides negative regulation driven by phosphatases

To investigate why responses of mast cells to antigen-induced IgE receptor (FcεRI) aggregation depend nonlinearly on antigen dose, we characterized a new artificial ligand, DF3, through complementary modeling and experimentation. This ligand is a stable trimer of peptides derived from bacteriophage T4 fibritin, each conjugated to a hapten (DNP). We found low and high doses of DF3 at which degranulation of mast cells sensitized with DNP-specific IgE is minimal, but ligand-induced receptor aggregation is comparable to aggregation at an intermediate dose, optimal for degranulation. This finding makes DF3 an ideal reagent for studying the balance of negative and positive signaling in the FcεRI pathway. We find that the lipid phosphatase SHIP and the protein tyrosine phosphatase SHP-1 negatively regulate mast cell degranulation over all doses considered. In contrast, SHP-2 promotes degranulation. With high DF3 doses, relatively rapid recruitment of SHIP to the plasma membrane may explain the reduced degranulation response. Our results demonstrate that optimal secretory responses of mast cells depend on the formation of receptor aggregates that promote sufficient positive signaling by Syk to override phosphatase-mediated negative regulatory signals.

Recruitment of the adaptor protein Grb2 to EGFR tetramers

Adaptor protein Grb2 binds phosphotyrosines in the epidermal growth factor (EGF) receptor (EGFR) and thereby links receptor activation to intracellular signaling cascades. Here, we investigated how recruitment of Grb2 to EGFR is affected by the spatial organization and quaternary state of activated EGFR. We used the techniques of image correlation spectroscopy (ICS) and lifetime-detected Förster resonance energy transfer (also known as FLIM-based FRET or FLIM–FRET) to measure ligand-induced receptor clustering and Grb2 binding to activated EGFR in BaF/3 cells. BaF/3 cells were stably transfected with fluorescently labeled forms of Grb2 (Grb2–mRFP) and EGFR (EGFR–eGFP). Following stimulation of the cells with EGF, we detected nanometer-scale association of Grb2–mRFP with EGFR–eGFP clusters, which contained, on average, 4 ± 1 copies of EGFR–eGFP per cluster. In contrast, the pool of EGFR–eGFP without Grb2–mRFP had an average cluster size of 1 ± 0.3 EGFR molecules per punctum. In the absence of EGF, there was no association between EGFR–eGFP and Grb2–mRFP. To interpret these data, we extended our recently developed model for EGFR activation, which considers EGFR oligomerization up to tetramers, to include recruitment of Grb2 to phosphorylated EGFR. The extended model, with adjustment of one new parameter (the ratio of the Grb2 and EGFR copy numbers), is consistent with a cluster size distribution where 2% of EGFR monomers, 5% of EGFR dimers, <1% of EGFR trimers, and 94% of EGFR tetramers are associated with Grb2. Together, our experimental and modeling results further implicate tetrameric EGFR as the key signaling unit and call into question the widely held view that dimeric EGFR is the predominant signaling unit.

An Interaction Library for the FcεRI Signaling Network

Antigen receptors play a central role in adaptive immune responses. Although the molecular networks associated with these receptors have been extensively studied, we currently lack a systems-level understanding of how combinations of non-covalent interactions and post-translational modifications are regulated during signaling to impact cellular decision-making. To fill this knowledge gap, it will be necessary to formalize and piece together information about individual molecular mechanisms to form large-scale computational models of signaling networks. To this end, we have developed an interaction library for signaling by the high-affinity IgE receptor, FcεRI. The library consists of executable rules for protein–protein and protein–lipid interactions. This library extends earlier models for FcεRI signaling and introduces new interactions that have not previously been considered in a model. Thus, this interaction library is a toolkit with which existing models can be expanded and from which new models can be built. As an example, we present models of branching pathways from the adaptor protein Lat, which influence production of the phospholipid PIP₃ at the plasma membrane and the soluble second messenger IP₃. We find that inclusion of a positive feedback loop gives rise to a bistable switch, which may ensure robust responses to stimulation above a threshold level. In addition, the library is visualized to facilitate understanding of network circuitry and identification of network motifs.

Mitochondrial morphological features are associated with fission and fusion events

Mitochondria are dynamic organelles that undergo constant remodeling through the regulation of two opposing processes, mitochondrial fission and fusion. Although several key regulators and physiological stimuli have been identified to control mitochondrial fission and fusion, the role of mitochondrial morphology in the two processes remains to be determined. To address this knowledge gap, we investigated whether morphological features extracted from time-lapse live-cell images of mitochondria could be used to predict mitochondrial fate. That is, we asked if we could predict whether a mitochondrion is likely to participate in a fission or fusion event based on its current shape and local environment. Using live-cell microscopy, image analysis software, and supervised machine learning, we characterized mitochondrial dynamics with single-organelle resolution to identify features of mitochondria that are predictive of fission and fusion events. A random forest (RF) model was trained to correctly classify mitochondria poised for either fission or fusion based on a series of morphological and positional features for each organelle. Of the features we evaluated, mitochondrial perimeter positively correlated with mitochondria about to undergo a fission event. Similarly mitochondrial solidity (compact shape) positively correlated with mitochondria about to undergo a fusion event. Our results indicate that fission and fusion are positively correlated with mitochondrial morphological features; and therefore, mitochondrial fission and fusion may be influenced by the mechanical properties of mitochondrial membranes.

Rule-based modeling: a computational approach for studying biomolecular site dynamics in cell signaling systems

Rule-based modeling was developed to address the limitations of traditional approaches for modeling chemical kinetics in cell signaling systems. These systems consist of multiple interacting biomolecules (e.g., proteins), which themselves consist of multiple parts (e.g., domains, linear motifs, and sites of phosphorylation). Consequently, biomolecules that mediate information processing generally have the potential to interact in multiple ways, with the number of possible complexes and posttranslational modification states tending to grow exponentially with the number of binary interactions considered. As a result, only large reaction networks capture all possible consequences of the molecular interactions that occur in a cell signaling system, which is problematic because traditional modeling approaches for chemical kinetics (e.g., ordinary differential equations) require explicit network specification. This problem is circumvented through representation of interactions in terms of local rules. With this approach, network specification is implicit and model specification is concise. Concise representation results in a coarse graining of chemical kinetics, which is introduced because all reactions implied by a rule inherit the rate law associated with that rule. Coarse graining can be appropriate if interactions are modular, and the coarseness of a model can be adjusted as needed. Rules can be specified using specialized model-specification languages, and recently developed tools designed for specification of rule-based models allow one to leverage powerful software engineering capabilities. A rule-based model comprises a set of rules, which can be processed by general-purpose simulation and analysis tools to achieve different objectives (e.g., to perform either a deterministic or stochastic simulation).

Exploring higher-order EGFR oligomerisation and phosphorylation--a combined experimental and theoretical approach

The epidermal growth factor receptor (EGFR) kinase is generally considered to be activated by either ligand-induced dimerisation or a ligand-induced conformational change within pre-formed dimers. Ligand-induced higher-order EGFR oligomerisation or clustering has been reported but it is not clear how EGFR oligomers, as distinct from EGFR dimers, influence signaling outputs. To address this question, we combined measures of receptor clustering (microscopy; image correlation spectroscopy) and phosphorylation (Western blots) with modelling of mass-action chemical kinetics. A stable BaF/3 cell-line that contains a high proportion (>90%) of inactive dimers of EGFR-eGFP but no secreted ligand and no other detectable ErbB receptors was used as the model cell system. EGF at concentrations of greater than 1 nM was found to cluster EGFR-eGFP dimers into higher-order complexes and cause parallel increases in EGFR phosphorylation. The kinetics of EGFR clustering and phosphorylation were both rapid, plateauing within 2 minutes after stimulation with 30 nM EGF. A rule-based model was formulated to interpret the data. This model took into account ligand binding, ligand-induced conformational changes in the cytosolic tail, monomer-dimer-trimer-tetramer transitions via ectodomain- and kinase-mediated interactions, and phosphorylation. The model predicts that cyclic EGFR tetramers are the predominant phosphorylated species, in which activated receptor dimers adopt a cyclic side-by-side orientation, and that receptor kinase activation is stabilised by the intramolecular interactions responsible for cyclic tetramerization.

Single-cell measurements of IgE-mediated FcεRI signaling using an integrated microfluidic platform

Heterogeneity in responses of cells to a stimulus, such as a pathogen or allergen, can potentially play an important role in deciding the fate of the responding cell population and the overall systemic response. Measuring heterogeneous responses requires tools capable of interrogating individual cells. Cell signaling studies commonly do not have single-cell resolution because of the limitations of techniques used such as Westerns, ELISAs, mass spectrometry, and DNA microarrays. Microfluidics devices are increasingly being used to overcome these limitations. Here, we report on a microfluidic platform for cell signaling analysis that combines two orthogonal single-cell measurement technologies: on-chip flow cytometry and optical imaging. The device seamlessly integrates cell culture, stimulation, and preparation with downstream measurements permitting hands-free, automated analysis to minimize experimental variability. The platform was used to interrogate IgE receptor (FcεRI) signaling, which is responsible for triggering allergic reactions, in RBL-2H3 cells. Following on-chip crosslinking of IgE-FcεRI complexes by multivalent antigen, we monitored signaling events including protein phosphorylation, calcium mobilization and the release of inflammatory mediators. The results demonstrate the ability of our platform to produce quantitative measurements on a cell-by-cell basis from just a few hundred cells. Model-based analysis of the Syk phosphorylation data suggests that heterogeneity in Syk phosphorylation can be attributed to protein copy number variations, with the level of Syk phosphorylation being particularly sensitive to the copy number of Lyn.

Binding of nucleoid-associated protein fis to DNA is regulated by DNA breathing dynamics

Physicochemical properties of DNA, such as shape, affect protein-DNA recognition. However, the properties of DNA that are most relevant for predicting the binding sites of particular transcription factors (TFs) or classes of TFs have yet to be fully understood. Here, using a model that accurately captures the melting behavior and breathing dynamics (spontaneous local openings of the double helix) of double-stranded DNA, we simulated the dynamics of known binding sites of the TF and nucleoid-associated protein Fis in Escherichia coli. Our study involves simulations of breathing dynamics, analysis of large published in vitro and genomic datasets, and targeted experimental tests of our predictions. Our simulation results and available in vitro binding data indicate a strong correlation between DNA breathing dynamics and Fis binding. Indeed, we can define an average DNA breathing profile that is characteristic of Fis binding sites. This profile is significantly enriched among the identified in vivo E. coli Fis binding sites. To test our understanding of how Fis binding is influenced by DNA breathing dynamics, we designed base-pair substitutions, mismatch, and methylation modifications of DNA regions that are known to interact (or not interact) with Fis. The goal in each case was to make the local DNA breathing dynamics either closer to or farther from the breathing profile characteristic of a strong Fis binding site. For the modified DNA segments, we found that Fis-DNA binding, as assessed by gel-shift assay, changed in accordance with our expectations. We conclude that Fis binding is associated with DNA breathing dynamics, which in turn may be regulated by various nucleotide modifications.

Cellular transcription factors (TFs) are proteins that regulate gene expression, and thereby cellular activity and fate, by binding to specific DNA segments. The physicochemical determinants of protein-DNA binding specificity are not completely understood. Here, we report that the propensity of transient opening and re-closing of the double helix, resulting from thermal fluctuations, aka “DNA breathing” or “DNA bubbles,” can be associated with binding affinity in the case of Fis, a well-studied nucleoid-associated protein in Escherichia coli. We found that a particular breathing profile is characteristic of high-affinity Fis binding sites and that DNA fragments known to bind Fis in vivo are statistically enriched for this profile. Furthermore, we used simulations of DNA breathing dynamics to guide design of gel-shift experiments aimed at testing the idea that local breathing influences Fis binding. As a result, we show that via nucleotide modifications but without modifying nucleotides that directly contact Fis, we were able to transform a low-affinity Fis binding site into a high-affinity site and vice versa. The nucleotide modifications were designed only based on DNA breathing simulations. Our study suggests that strong Fis-DNA binding depends on DNA breathing - a novel physicochemical characteristic that could be used for prediction and rational design of TF binding sites.

Computational model for autophagic vesicle dynamics in single cells

Macroautophagy (autophagy) is a cellular recycling program essential for homeostasis and survival during cytotoxic stress. This process, which has an emerging role in disease etiology and treatment, is executed in four stages through the coordinated action of more than 30 proteins. An effective strategy for studying complicated cellular processes, such as autophagy, involves the construction and analysis of mathematical or computational models. When developed and refined from experimental knowledge, these models can be used to interrogate signaling pathways, formulate novel hypotheses about systems, and make predictions about cell signaling changes induced by specific interventions. Here, we present the development of a computational model describing autophagic vesicle dynamics in a mammalian system. We used time-resolved, live-cell microscopy to measure the synthesis and turnover of autophagic vesicles in single cells. The stochastically simulated model was consistent with data acquired during conditions of both basal and chemically-induced autophagy. The model was tested by genetic modulation of autophagic machinery and found to accurately predict vesicle dynamics observed experimentally. Furthermore, the model generated an unforeseen prediction about vesicle size that is consistent with both published findings and our experimental observations. Taken together, this model is accurate and useful and can serve as the foundation for future efforts aimed at quantitative characterization of autophagy.

Improved predictions of transcription factor binding sites using physicochemical features of DNA

Typical approaches for predicting transcription factor binding sites (TFBSs) involve use of a position-specific weight matrix (PWM) to statistically characterize the sequences of the known sites. Recently, an alternative physicochemical approach, called SiteSleuth, was proposed. In this approach, a linear support vector machine (SVM) classifier is trained to distinguish TFBSs from background sequences based on local chemical and structural features of DNA. SiteSleuth appears to generally perform better than PWM-based methods. Here, we improve the SiteSleuth approach by considering both new physicochemical features and algorithmic modifications. New features are derived from Gibbs energies of amino acid-DNA interactions and hydroxyl radical cleavage profiles of DNA. Algorithmic modifications consist of inclusion of a feature selection step, use of a nonlinear kernel in the SVM classifier, and use of a consensus-based post-processing step for predictions. We also considered SVM classification based on letter features alone to distinguish performance gains from use of SVM-based models versus use of physicochemical features. The accuracy of each of the variant methods considered was assessed by cross validation using data available in the RegulonDB database for 54 Escherichia coli TFs, as well as by experimental validation using published ChIP-chip data available for Fis and Lrp.