Collaboration & Service Projects

MMBioS will engage in multiple, diverse Collaboration and Service projects (C&SP).  These C&SPs  will benefit from the tools and technologies developed by MMBioS.  Current CS&Ps are listed below; this list will grow over time.

 More details are available for Collaboration & Service Projects listed below.


1) Relationship between structural Dynamics and Sequence Evolution; Application to molecular machined such as bacterial chaperonons and HSP70 chaperones

PI: I. Bahar, University of Pittsburgh  Collaborators: Amnon Horovitz, Weizmann Institute; R. Altman, Stanford University; L. Gierasch, University of Massachusetts.

aanm We have expanded the scope of this project to include Dr. Amnon Horovitz from Weizmann Institute, Israel, an expert in the allostery of chaperonins as well as sequence co-evolution. We are exploring the sequence->structure->dynamics->function mapping of allosteric proteins. Our overarching goal is to establish the computational methodology for simulating the machinery of biomolecular systems on the order of Megadaltons. The molecular machines that we selected for developing and implementing the methodology are bacterial chaperonins (GroEL), molecular chaperones (E. coli Hsp70 DnaK) and their co-chaperones. We study the ATP-regulated mechanisms of intersubunit (GroEL) or interdomain (DnaK) communication, as well as the role of co-chaperonins or co-chaperones in modulating the conformational changes occurring along the allosteric cycle of these machines.

In Year 3 we made progress in quantitative analysis of the evolution of proteins, and in particular in assessing the relationship, if any, between their structural dynamics and sequence evolution. In a recent study performed toward this goal (Mao et al., 2015), we carried out a systematic comparative analysis of current methods for analyzing sequence co-evolution and we proposed a hybrid model that yields high performance. The two metrics for performance were the elimination of false positives between non-interacting proteins, and the verification of tertiary contacts between co-evolving pairs. The methods used in this newly published (that acknowledged MMBioS support) are currently being implemented into our ProDy API.

Publications Resulting from This Work

    • General IJ, Liu Y, Blackburn M, Mao W, Gierasch LM, Bahar I (2014) ATPase subdomain IA is a mediator of interdomain allostery in Hsp70 molecular chaperones PLoS Comp Bio 10: e1003624 PMID: 24831085, PMC4022485
    • Mao W, Kaya C, Dutta A, Horovitz A, Bahar I (2015) Comparative Study of the Effectiveness and Limitations of Current Methods for Detecting Sequence Coevolution Bioinformatics 31: 1929-37 PMID: 25697822; PMC4481699



2) Membrane Proteins Structure and Dynamics Consortium (MPSDC) Computational Core, and AMPA receptors structural dynamics

PI: I. Bahar, University of Pittsburgh   Collaborators: B. Roux and E. Perozo, University of Chicago; K. Schulten and E. Tajkhorshid, University of Illinois at Urbana-Champaign;  I. Greger, University of Cambridge, UK; H. Weinstein, Cornell University

To assess whether ANM-predicted AMPAR NTD pivoting motions can occur in cellulo, the Greger lab mutated K262 to cysteine (K262C), hypothesizing that this mutation would trap the NTD tetramer in a closed conformer via formation of a disulfide bridge. GluA2 wild type (WT) and K262C mutant were expressed in HEK293 cells and protein extracts were analyzed by Western blot. The results showed that our hypothesis was valid, thus lending support to computational predictions. We further examined the dynamics of AMPAR and NMDAR to observe that the differ along a few slow modes only (Figure 1). Overall, our data provided a first glimpse into the dynamic spectrum of AMPAR and NMDARs and delineated conserved mechanisms underlying allosteric communication in iGluRs.

csp2Figure 1. AMPAR to NMDAR conformations differ by rearrangements along a few soft modes. The dark blue bars show the overlap (correlation cosine) between ANM modes of AMPAR and the structural difference vector between AMPAR and NMDAR; the red curve displays the cumulative overlap; the dotted green curve shows the expected cumulative overlap if the modes were equally contributing each. The inset shows the initial structures of AMPAR (left) and NMDAR (right).

Publications Resulting from This Work

    • Krieger J, Bahar I, Greger IH (2015) Structure, Dynamics, and Allosteric Potential of Ionotropic Glutamate Receptor N-Terminal Domains Biophys J109: 1136-48 PMID: 26255587, PMC4576161
    • Dutta A, Krieger J, Lee JY, Garcia-Nafria J, Greger IH, Bahar I (2015) Cooperative Dynamics of Intact AMPA and NMDA Glutamate Receptors: Similarities and Subfamily-Specific Differences Structure23: 1692-1704 PMID: 26256538, PMC4558295
    • Dutta A, Shrivastava IH, Sukumaran M, Greger IH, Bahar I (2012) Comparative dynamics of NMDA- and AMPA-Glutamate receptor N-terminal domains Structure20: 1838-49 PMID: 22959625, PMC3496038




3) Integration, prediction, and generation of mixed mode information using graphical models, with applications to protein-protein interactions

PIs: C. J. Langmead, Carnegie Mellon University Collaborators: C. Bailey-Kellog, Dartmouth University; N. Ramakrishnan, Virginia Tech;  A. Friedman, Purdue University

This C&SP was an NSF-funded project titled: “Integration, prediction, and generation of mixed mode information using graphical models, with applications to protein-protein interactions” (NSF IIS-0905193). It was collaboration between Carnegie Mellon University, Dartmouth College, Virginia Tech, and Purdue. The project concluded in August, 2013. The focus of this C&SP was to develop novel probabilistic graphical models for modeling protein-protein interactions. Three broad classes of techniques were developed to integrate attribute-value and relational information, integrate statistical and physical information, and utilize probabilistic models generatively. The research has resulted in new graphical models, new algorithms, and the application of those algorithms to the modeling of protein structures and dynamics, including protein-protein interactions. Specific accomplishments include: the first optimal algorithm for learning regularized undirected graphical models of protein sequences; the first method for performing binding free energy calculations for protein-protein interactions via graphical models; the first probabilistic graphical models of molecular dynamics; the first algorithm for learning the parameters of molecular mechanics force fields by minimizing functionals over Boltzmann distributions; the first undirected model of distributions on the hypersphere (for modeling distributions over angles); a game-theoretic method for modeling the emergence of drug resistance-causing mutations in proteins; new algorithms for learning semi- and non-parametric distributions; and a new class of regression models for predicting binding free energies. The project resulted in seventeen publications, three Ph.D. dissertations, and one patent. Primary applications of this research include computer-aided drug design and computer-aided protein engineering. Secondary applications of this work include techniques for identifying functionally important residues mediating binding and allostery.

Publications Resulting from This Work

    • Kamisetty H, Ghosh B, Langmead CJ, Bailey-Kellogg C (2014) Learning Sequence Determinants of Protein: protein Interaction Specificity with Sparse Graphical Models Res Comput Mol Biol8394:129-143 PMID: 25414914, PMC4235964.
    • Kamisetty H, Ghosh B, Langmead CJ, Bailey-Kellogg C (2015) Learning sequence determinants of protein: protein interaction specificity with sparse graphical models J Comput Biol22:474- 86 PMID: 25973864, PMC4449715



4) Calcium entry and transmitter release at the neuromuscular junction

PIs: T. Bartol, T. Sejnowski, Salk Institute; R. Laghaei, G. Hood, M. Dittrich, Pittsburgh Supercomputing Center  Collaborator: S. Meriney, University of Pittsburgh

Communication between cells in the nervous system (synaptic transmission) underlies all complex behaviors, is often disrupted in neurological disease, and is a focus for therapeutic intervention. Synapses work by releasing chemical transmitter from a region called the active zone. In this project we synergistically combine computer simulation using MCell (Dittrich lab) and synaptic anatomy, physiology, and Ca2+ imaging (Meriney lab) to investigate the structure and function of synapses. We have previously completed and published a baseline model of frog neuromuscular junction (NMJ) active zone structure and function (Dittrich et al., 2013). Based on this model we have recently investigated several possible mechanisms underlying short-term synaptic facilitation at the frog NMJ (Ma et al., 2015). Our study showed that a vesicle release mechanism in which a second set of Ca2+ sensor sites was responsible for facilitation provided good overall agreement with our experimental constraints. In addition, we used our MCell modeling approach combined with Ca2+ imaging, pharmacological Ca2+ channel block, and postsynaptic recording to show that release of individual synaptic vesicles at the frog NMJ is predominately triggered by Ca2+ ions entering the nerve terminal through the nearest open calcium channel. This work has recently been published (Luo et al., 2015). Building on these insights, we are currently expanding our investigation to the role of detailed active zone structure in synapse function with a focus on neuromuscular diseases such as Lambert-Eaton Myasthenic Syndrome.

Publications Resulting from This Work

    • Tarr TB, Dittrich M, Meriney SD (2013) Are unreliable release mechanisms conserved from NMJ to CNS? Trends Neurosci36(1):14-22. PMID: 23102681, PMC4076818
    • Meriney SD, Dittrich M (2013) Organization and function of transmitter release sites at the neuromuscular junction J Physiol591(13):3159-65. PMID: 23613535, PMC3717219
    • Dittrich M, Pattillo JM, King JD, Cho S, Stiles JR, Meriney SD (2013) An excess-calcium-binding-site model predicts neurotransmitter release at the neuromuscular junction Biophys J104(12):2751-63 PMID: 23790384, PMC3686347
    • Ma J, Kelly L, Ingram J, Price TJ, Meriney SD, Dittrich M (2015) New insights into short-term synaptic facilitation at the frog neuromuscular junction J Neurophysiol113:71-87 PMID: 25210157, PMC4294575
    • Luo F, Dittrich M, Cho S, Stiles JR, Meriney SD (2015) Transmitter release is evoked with low probability predominately by calcium flux through single channel openings at the frog neuromuscular junction J Neurophysiol113:2480-9 PMID: 25652927, PMC4416571



5) Modeling immunoreceptor signaling, autophagy, and endocytosis, Hardening Software for Rule-Based Modeling

PI:  J. Faeder, University of Pittsburgh  Collaborator: W. S. Hlavacek, Los Alamos National Laboratory

The fundamental events in FcεRI signaling initiation were considered during the pioneering work of Goldstein, Faeder, Hlavacek and colleagues to establish the utility of rule-based methods for mechanistic modeling of the dynamics of molecular interactions in signaling networks. This led to the creation of BioNetGen, already introduced in the BTRC proposal as a computational framework for modeling biochemical networks. Here, we propose to take advantage of successful initiatives at the BTRC for integration of the BioNetGen rule-based modeling approach into the mesh-based architecture that underlies MCell spatial simulations. In particular, we will focus on the proposed new capabilities for adaptive meshing, which can superimpose well-characterized changes in surface geometry over time during simulations. This unique feature will allow us to develop novel computational models to describe simultaneously slow endocytic vesicle dyanamics and fast signaling dynamics at vesicle surfaces.

Aim 1. To utilize the CellBlender interface to create a framework for cell shape changes that occur during mast cell signaling, such as microvillar-to-ruffled transformation, appearance of signaling patches, and changes associated with enhanced endocytosis.

Aim 2. To evaluate the effects of a changing 3D geometry on signaling output, by incorporating existing Rule-Based Models for FcεRI chemical reaction networks into MCell-BioNetGen and refining the model with new parameters for spatial distributions, rate constants and diffusion characteristics collected through funded initiatives.

During the project we developed MCell-R, which now enables network-free simulations to be performed directly in MCell based on a BioNetGen model specification. We have validated that the simulator gives correct results for a range of FceRI signaling models and also for the bivalent-ligand bivalent receptor, a simple example of unrestricted aggregation. In principle the simulator can now be used to simulate spatial signaling dynamics of any model that can be described in BNGL with any geometry that can be specified in MCell, including a geometry that changes in time. We plan to implement several such models in the future.

Publications Resulting from This Work

    • Chylek LA, Harris LA, Tung CS, Faeder JR, Lopez CF, Hlavacek WS (2013) Rule-based modeling: a computational approach for studying biomolecular site dynamics in cell signaling systems Wiley Interdiscip Rev Syst Biol Med6: 13-36 PMID: 24123887, PMC3947470
    • Chylek LA, Harris LA, Faeder JR, Hlavacek WS (2015) Modeling for (physical) biologists: an introduction to the rule-based approach Phys Biol12: 045007 PMID: 26178138, PMC4526164

In progress


6) Distance-dependent structure and function of neuronal dendrites

PIs: T. Sejnowski and T. Bartol, The Salk Institute Collaborator: K. Harris, University of Texas at Austin

The purpose of this CS&P is to collaborate with Kristen Harris' lab to create realistic MCell models of dendritic structure and function.

MCell can simulate the diffusion and interaction of molecules involved in biochemical signaling pathways within the 3D subcellular structure of cells. To do so, surface meshes used to represent cell membranes and subcellular structures must meet very strict geometric standards (e.g. water-tight, non-intersecting, manifold). We help create realistic MCell models of CA1 dendritic spines, in particular to learn how more the accurate serial section tomographic reconstructions of core organelles impact simulations of molecular signaling within dendritic spines.

Recently two papers were published relating to this C&SP. The first paper, published in Frontiers in Synaptic Neuroscience, reported on reconstitution of calcium dynamics in dendritic spines. Nine parameters of the model were optimized within realistic experimental limits by a process that compared results of simulations to published data. Simulations in the optimized model reproduce the timing and amplitude of Ca(2+) transients measured experimentally in intact neurons. Thus, the characteristics of individual isolated proteins determined in vitro could accurately reproduce the dynamics of experimentally measured Ca(2+) transients in spines. The second paper, published in eLife, found that dendritic spines that receive input from the same axon are the same size. This finding allowed estimation of the variability of synaptic plasticity and we found that the amount of information stored at synapses is approximately 4.7 bits, which is an order of magnitude larger than previous estimates.

Publications Resulting from This Work

    • Edwards J, Daniel E, Kinney J, Bartol T, Sejnowski T, Johnston D, Harris K, Bajaj C (2014) VolRoverN: enhancing surface and volumetric reconstruction for realistic dynamical simulation of cellular and subcellular function Neuroinformatics12:277-89 PMID: 24100964, PMC4033674
    • Kinney JP, Spacek J, Bartol TM, Bajaj CL, Harris KM, Sejnowski TJ (2013) Extracellular sheets and tunnels modulate glutamate diffusion in hippocampal neuopil J Comp Neurol521: 448-64 PMID: 22740128, PMC3540825
    • Bartol TM, Keller DX, Kinney JP, Bajaj CL, Harris KM, Sejnowski TJ, Kennedy MB (2015) Computational reconstitution of spine calcium transients from individual proteins Front Synaptic Neurosci7: 17 PMID: 26500546, PMC4595661
    • Bartol TM, Bromer C, Kinney J, Chirillo MA, Bourne JN, Harris KM, Sejnowski TJ (2015) Nanoconnectomic upper bound on the variability of synaptic plasticity Elife4.pii: e10778 PMID: 26618907, PMC4737657



7) Interfacing image-derived generative models with cell simulation engines

PIs: R. Murphy and G. Rohde, Carnegie Mellon University; M. Dittrich, Pittsburgh Supercomputing Center; J. Faeder, University of PittsburghCollaborators: T. Bartol, The Salk Institute; S. Andrews and R. Brent, Fred Hutchison Cancer Research Center; I. Moraru, J. Schaff and L. Loew, University of Connecticut

The primary goal of this collaborative project is to develop methods to exchange generative models of the spatiotemporal organization of cells between various cell simulation engines. This will enable those engines to have a greatly expanded source of cell geometries, to take full advantage of dynamic spatial models, and to combine models constructed by various means. The vision is for CellOrganizer to be able to work synergistically with Virtual Cell, Smoldyn, MCell, and any other cell simulator. A secondary, but important, goal will be to establish a de facto standard by which future systems for constructing generative cell models can interface with simulation engines. The interfacing would go well beyond the SBML spatial package, consisting of modules to generate various types of models from different image sources, to save them in exchangeable files, and most importantly, plugins for simulators to generate instances from the models to use in simulations and specify how those instances change over time.

In progress

8) Using generative models of cell organization to investigate tumor cell heterogeneity

PIs: R. Murphy and G. Rohde, Carnegie Mellon University Collaborators: S. Altschuler and L. Wu, University of Texas Southwestern Medical Center

Tumor cell populations have long been known to be complex mixtures of multiple phenotypes. Given the precision with which one is able to measure properties of cellular structure and function, it is often the case that the measure of two cells will differ. Paramount to deriving useful information from a potentially large set of measurements is to be able to determine which differences are statistically significant as well as their biological meaning. The established framework for analyzing the heterogeneity of tumor cell populations involves analyzing these measurements in a high dimensional feature space. Features can be simple and easily interpretable (such as cell size), or complex and difficult to interpret (such as texture features). The goal of this subproject is to try to represent cell images using generative models that are more easily interpretable and could potentially offer more insights into the underlying phenomena that affect the distribution of phenotypes in tumor cells.

In progress


9) Large-scale electron microscopy of calcium-imaged neuron populations

PIs: A. Wetzel and G. Hood, Pittsburgh Supercomputing Center Collaborator: D. Bock, Janelia Farm Research Campus, Howard Hughes Medical Institute

During the last year of this project, the biological focus of Bock's work has shifted from localized mouse brain circuits to the mapping of complete drosophila brains with an emphasis on visual circuits. This change did not affect the relation of our collaboration which is focused on the technical aspects of maximum scale serial-TEM data assembly.

Bock's group has improved the throughput of their automated 4 camera TEM system. Although the cameras and their raw data rates are unchanged new mechanisms have been tested in the last year for automatic sample cassette loading. This has substantially improved the ratio of imaging to sample exchange time. The long term goal for this new sample loading method is to enable continuous runs up to 2 weeks for datasets of 15,000 sections and >50 TBytes. Individual sections are imaged in arrays of ~20,000 5 MPixel tiles, typically 112 columns by 196 rows, with a unique stepping pattern to accommodate constraints on the spacing of the optical cameras and lenses. One difficult aspect of these data is the use of the small 5 MPixel tiles which, particularly in the case of hierarchical registration methods, become limited by poor performance of large filesystems when processing small files.

Regular production of datasets at this scale would be difficult to register using our previous method of raw data transmission to the PSC followed by the return of aligned results to Janelia Farm. Therefore we are preparing for Bock's team to run our alignment codes on the large 6,144 core Janelia Farm computing cluster. Wetzel and Hood were part of an April 1-11, 2015 Janelia Farm image registration Hackathon organized by Stephan Saafeld's group. During that trip we were able to test both the AlignTK and SWiFT codes on the Janelia systems using data from Bock's work as well as many others from Janelia and elsewhere (i.e. Winfried Denk and our C&SP10 collaborator Jeff Lichtman). We were also able to adapt a key component of the SWiFT method for use within one of the Janelia registration codes and found that it greatly improved the reliability of registration on difficult regions of Bock's data.

Publications Resulting from This Work

    • Bock DD, Lee WC, Kerlin AM, Andermann ML, Hood G, Wetzel AW, Yurgenson S, Soucy ER, Kim HS, Reid RC (2011) Network anatomy in vivo physiology of visual cortical neurons Nature471: 177-82 PMID: 21390124, PMC3095821



10) Advancing high-throughput thin-section scanning EM to study relationships between neuronal circuit structure and function

PIs: A. Wetzel and G. Hood, Pittsburgh Supercomputing Center Collaborator: J. Lichtman, Harvard University

We have continued work on automated alignment of serial SEM image sets produced using tape collecting ultramicrotome sectioning and wafer mounted imaging procedures developed at Harvard. The primary datasets have been the 16,000 section zebrafish stack, which is the focus of a new CS&P project with the Engert lab and the 114 TB lateral geniculate nucleus (LGN) volume that was previously acquired and is the focus of this C&SP.

The different requirements and resolutions of these datasets have highlighted the need for a variety of approaches to handle different types of biological content and different types of image characteristics. The large number of tape cracks in the LGN dataset has been particularly difficult since the cracks from any section to its neighbors interfere with one another during the alignment process. Even though the detailed cracking pattern is different on each section the fact that most of the cracks have similar sizes and orientations means that they tend to pull the alignment such that we have seen errors on the order of 1 per 1000 sections. We are currently working on a strategy to produce a low-resolution map of the cracks, 1/4th scale, that can then be used to constrain the SWiFT alignment process.

We are entering a new and much higher-throughput phase of EM data capture. Previous datasets, including the zebrafish and LGN, were captured over several months using a Zeiss Merlin microscope and pixel rates up to 20MHz. Lichtman's team recently installed the first Zeiss 61-beam 1.2 Gpixel/sec scanning microscope. This microscope is now being tested with large single sections in the 400GB/section range, which are captured in a hexagonal mosaic pattern. These data will require a substantial number of registration code adaptations. We will test speed and accuracy of these modifications once the new equipment is producing serial spans of at least 100 sections, which is where the power of the SWiFT approach is most useful.

Publications Resulting from This Work

    • Morgan JL, Berger DR, Wetzel AW, Lichtman JW (2016) The fuzzy logic of network connectivity in mouse visual thalamus Cell165: 192-206 PMID: 27015312, PMC4808248
    • Hildebrand DGC, Torres RM, Choi W, Tran Minh Quan, Arthur Willis Wetzel, George ScottPlummer, Ruben Portugues, Isaac Henry Bianco, Owen Randlett, Stephan Saalfeld, Alex Baden, Kunal Lillaney, Randal Burns,Joshua Tzvi Vogelstein, Won-Ki Jeong, Jeff William Lichtman, Florian Engert (2016) Whole-brain serial-section electron microscopy in larval zebrafish Nature, revised version in preparation.



11) Morphological and regulatory models of neuronal differentiation

PIs: H Busch and M Böerries, University of Freiburg Collaborators: R Murphy and G Rohde, Carnegie Mellon University

The goal of this project is to develop a spatiotemporal generative model of the changes in cell size, shape and subcellular organization that occur during the differentiation of PC12 cells induced by nerve growth factor. This model will be related to a regulatory model constructed from parallel measurements of RNA expression over the time course of differentiation. The results are expected to lead to the identification of gene expression changes that lead to specific morphological changes and to test their importance for the differentiation process.

Publications Resulting from This Work

    • Johnson GR, Bierschenk I, Nitschke R, Boerries M, Busch H, and Murphy RF. (2017) Image-derived Models of Cell Organization Changes During Differentiation of PC12 Cells, submitted.


12) Generative models of plant organelle distribution and differentiation

PI: K Palme, University of Freiburg Collaborator: R Murphy, Carnegie Mellon University

Plants exhibit the ability to undergo dramatic dedifferentiation and redifferentiation such that both protoplasts (cells extracted from mature plants) and microspores (cells that arise during gametogenesis) can give rise to mature plants under the right circumstances. The goal of this project is to identify the specific changes in protein expression and localization that are associated with such processes, in part by building spatiotemporal generative models directly from images and movies of cells obtained by light microscopy.

Publications Resulting from This Work

    • Johnson GR, Kangas JD, Dovzhekno A, Trojok R, Voigt K, Majarian TD, Palme K and Murphy RF. (2017) A Method for Characterizing Phenotypic Changes in Highly Variable Cell Populations and its Application to High Content Screening of Arabidopsis thaliana Protoplasts. Cytometry Part A, in press.


13) Microtubule pattern analysis and drug sensitivity

PI: P Giannakakou, Cornell University Collaborator: R Murphy, Carnegie Mellon University

The goal of this new project is to determine whether and how the distributions of microtubules differ between tumor cells sensitive and resistant to drugs that act on microtubules. Initial work has focused on training systems to recognize these patterns, with a complicating factor being the large differences in cell size, shape and unperturbed microtubule pattern between cells from different tumors. Preliminary results indicate that three broad but distinct patterns can be distinguished, corresponding to untreated cells, drug-treated resistant cells, and drug-treated sensitive cells. Future work will focus on using generative models to remove variation due to cell size and shape. Ultimately, this would be used to determine whether a particular tumor is likely to be resistant to a particular drug and thereby to choose an appropriate drug for that tumor.



14) Actin filament patterns and cell motility

PI: J Theriot, Stanford University Collaborators: G Rohde and D Slepčev, Carnegie Mellon University

This is a new project, the goal of which is to quantitatively describe the patterns of actin filaments and related proteins as they relate to cell motility in neutrophils. The project aims to use recently developed methods for measuring the similarity between such patterns, as well as transport-based method for finding correlation between intensity patterns, to decode relationships between different proteins as cells move. Results are expected is to help characterize cell motion as a function of the subcellular protein patterns.


15) Reconstructing zebrafish neural circuits controlling visually induced behaviors

PI: Art Wetzel and Greg Hood, Pittsburgh Supercomputing Center Collaborators: Florian Engert and David Hildebrand, Harvard University

The general goal of this new C&SP is the development of the larval zebrafish as a model system for the comprehensive identification and examination of neural circuits controlling visually induced behaviors. To that end our collaborators plan to establish and quantify a series of visually induced behaviours and analyze the individual resulting motor components. Using these assays they will monitor neuronal activity throughout the fish brain in an awake and intact preparation.

The underlying data are highly accurate zebrafish reconstructions from serial section SEM. The primary data is an 18,000 section series prepared by our collaborator David Hildebrand of a 6 day post fertilization specimen with full body cross section coverage extending from the front of the animal through the entire brain and into the initial spinal cord region. This SEM data consists of 3 resolutions 1) a very low resolution overview, ~1 micron/pixel, of the full18,000 section series that was used for ROI selection and targeting 2) a moderate resolution 60nm per pixel scan of full body cross sections over a 16,000 section ROI used for tracing myelinated pathways and 3) a 20 nm per pixel scan of a 12,000 section sub-region covering the entire brain that is used for detailed circuit tracing. As part of this collaboration we have applied our SWIFT approach to help produce highly accurate alignments of these image stacks. This required a number of significant improvements to the underlying core algorithms in our SWIFT software. Using SWIFT, we were already able to produce a very well aligned stack of the 60 nm image data shown in Figure 1 as a greatly reduced cut plane view. We aren currently using the SWIFT approach align the higher resolution 20nm stack, which is needed for accurate segmentation of the zebrafish brain area, into a geometry that is consistent with the lower resolution volumes.

The longer-term goal of correlating in-vivo neural activity with brain structure requires multi-modal alignment of EM data with 2-photon imaging as the animal reacts to visual stimuli. A preliminary example of the 2-photon data, acquired from a different specimen, is shown in Figure 2 with the primary forebrain ROI for in-vivo study indicated by the green square. After detailed reconstruction of pathways in the current dataset we will prepare a series of dual SEM plus 2-photo datasets for this forebrain region that will require accurate and efficient registration of the EM data for structural tracing along with alignment of the 2-photon data onto the EM data for correlation to in-vivo data recordings.

Publications Resulting from This Work

      • Hildebrand DGC, Torres RM, Choi W, Tran Minh Quan, Arthur Willis Wetzel, George ScottPlummer, Ruben Portugues, Isaac Henry Bianco, Owen Randlett, Stephan Saalfeld, Alex Baden, Kunal Lillaney, Randal Burns,Joshua Tzvi Vogelstein, Won-Ki Jeong, Jeff William Lichtman, Florian Engert (2016) Whole-brain serial-section electron microscopy in larval zebrafish Nature, revised version in preparation.

In progress


16) Latent Factor Models for Identification of Novel Neuroprotectives for Huntington's Disease

PI: Ivet Bahar and Bing Liu, University of Pittsburgh;  Collaborators: Robert Friedlander, Mark Shurdak, Andrew Stern, D. Lancing Taylor, University of Pittsburgh

This C&SP is supported by NIH award 5R01NS077748-04 (PI: Friedlander) entitled “Functional role of microRNAs in Huntington's disease pathogenesis”. Our focus is the identification of novel neuroprotectives using latent factor modeling (LFM). Firstly, we built a set of known neuroprotectives based on the data in the literature by collecting compounds that were identified in a large scale screen for neuroprotectives in Huntington’s disease (HD) as well as compounds that have been tested in clinical trials for their neuroprotective effect. Then, the LFM-based algorithm learns the most successful latent factor model of the drug-target interactions in the STITCH dataset after assessing the performance of multiple LFM learning algorithms. The LFM of STITCH allowed us to quantitatively assess the interaction profiles of the known neuroprotectives, which we then used to find new compounds toward (i) the discovery of new neuroprotectives by target-based diversification of known neuroprotectives, (ii) the discovery of mechanism of action hypotheses by identifying targets of interest. Specifically, for each target of known neuroprotectives, we identified other chemicals known to be interacting with high confidence and sorted them for maximal dissimilarity to the known neuroprotectives interacting with that target. This method identifies compounds that share the target of interest, but are otherwise as diverse in their interaction profile as possible. Compounds selected will be active if the target of interest is important for the HD etiology, therefore they represent new potential treatment candidates. Furthermore if one or more targets yield many neuroprotectives, that target is implicated as being central to the disease pathophysiology; therefore this method allows us to infer about the proteins regulating the disease phenotype by testing compounds alone.

Publications Resulting from This Work

      • Stern, A. M.; Schurdak, M. E.; Bahar, I.; Berg, J. M.; Taylor, D. L.(2016) A Perspective on Implementing a Quantitative Systems Pharmacology Platform for Drug Discovery and the Advancement of Personalized Medicine. J Biomol. Screen., in press. PMID: 26962875

In progress


17) Structure and Function of Voltage-Gated Calcium Channels

PI: Ivet Bahar and Mary Cheng, University of Pittsburgh; Tom Bartol, Salk Institute; Collaborators: Peter Wipf and Steve Meriney, University of Pittsburgh

This C&SP is a joint Charles E. Kaufman Foundation Grant proposal titled “Structure and Function of Voltage-Gated Calcium Channels (VGCCs)”, submitted on June 30 2015. The goal is to establish a chemistry - computational biology - neuroscience collaboration and launch a novel combination of strategic chemical synthesis (Wipf), VGCC mutagenesis and biophysical measurements of VGCC function (Meriney), and computational modeling (Bahar) toward the development a model of Cav2 VGCC gating. This work will fill a critical gap in our understanding of VGCC structure and function and significantly impact future studies in the field. VGCC are critical regulators of cellular function. We are interested in an important family of VGCCs (Cav2) that regulate communication between neurons in the brain. Despite its critical importance in cellular function, surprisingly little is known about VGCC structure and function, primarily because a crystal structure has proven very difficult to generate for this large membrane protein. In preliminary studies, the Wipf and Meriney labs have developed new tools targeted selectively to voltage-gated calcium channels. With the expertise of a new member of our team, Dr. Bahar, we propose to deploy novel chemical synthesis, calcium channel mutagenesis, biophysical measurements of calcium channel function, and computational modeling of channel movement to significantly increase our understanding of the structure-function relationships of voltage-gated calcium channels. This work will be transformative since it will impact the study of communication between cells, especially in the brain, and likely lead to the development of novel drugs that may be useful for a variety of diseases. Bahar lab will focus on structural modeling and simulations of VGCC gating and drug binding: a) modeling the structures and dynamics of the Cav2 VGCC in different conformational states, with a focus on key residues and interactions that control gating; and b) determination of the mechanism of action of the agonist gating modifiers of Cav2 VGCCs .

In progress


18) Effect of RNA-editing of ADAR1 on the activity of ionotropic glutamate receptors

PI: Ivet Bahar, Mary H Cheng, University of Pittsburgh  Collaborators: QingDe Wang, University of Pittsburgh

This is a C&SP, supported by NIH 5R21CA158650-02 (PI: Wang). The goal of the project is to define the function of RNA editing enzyme adenosine deaminase acting on RNA 1 (ADAR1) in hematopoietic and leukemia stem cells. ADAR1 is an essential protein for embryonic and adult hematopoiesis, while the leukemia cells are more susceptible to the gene deletion that codes ADAR1. The impact of ADAR1 on normal and leukemia stem cells (LSC) remain unknown, and there is a need to develop a therapeutic strategy to eliminate LSCs by targeting ADAR1. Although ADAR1 is identified as an RNA editing enzyme and RNA editing has been shown to play critical roles in stem cells and other biological processes, multiple attempts to date have been unable to identify an editing target that accounts for the death of normal hematopoietic and leukemia cells of ADAR1 knockouts. The R21 seeks a definitive answer whether the editing activity of ADAR1 is necessary for the proliferation and differentiation of LSCs.

ADAR has been found to modify the AMPA receptor GluR2 and G-protein coupling receptors, i.e, SR2C and GABA receptor. The knockout animal models (Wang lab) have shown that this molecule is essential for embryonic survival and tissue homeostasis. However the molecular mechanism is largely unknown. The Bahar lab has initiated a new collaboration with the Wang Lab, to help investigate (a) the effects of the mutations induced upon ADAR1 editing on calcium channeling properties of iGluRs, (ii the effect of specific mutations on the functions of serotonin receptors and potassium channels, using molecular modeling, and (c) mathematical modeling and quantitative assessment of the effect of suppressing cytosolic RNA pathway by ADAR1 on relevant cellular pathways networks.

In progress


19) Modeling Dendritic Sparks at Dendritic Branch Points

PI: Terry Sejnowski and Tom Bartol, Salk Institute  Collaborators: Mark Ellisman, University of California at San Diego

The purpose of this C&SP is to collaborate with Mark Ellisman's lab to create realistic MCell models of calcium sparks that occur at dendritic branch points.

Neuronal calcium events similar to calcium sparks in cardiac myocytes have been observed to occur at dendritic branch points (Manita and Ross, J Neurosci, 2009). The origin of these events has been attributed to ryanodine receptors (RyRs) and IP3 receptors (IP3Rs). The Ellisman lab is working to reconstruct the subcellular ultrastructure of dendritic branch points and to identify the localization and relative concentrations of channels and molecules involved in calcium signaling and processing, including: ryanodine receptors (RyR1, RyR2, RyR3), IP3R), metabotropic glutamate receptors (mGluR), various voltage-gated calcium channels (T-Type, L-Type, and N-Type), various calcium-activated potassium channels (BK, SK, Kv4.2), and sodium channels. From these data we will help the Ellisman lab create a realistic model of dendritic branch points that incorporates both voltage dynamics at the neuronal level and calcium and other molecular dynamics at single branch points with via a combined NEURON and MCell system where both are using realistic morphologies from the same sample.

In progress


20) Structural plasticity of chaperonins determined by cryo-electron microscopy: Modeling the machinery of GroEL and TRiC

Collaborating Investigators: Ivet Bahar, University of Pittsburgh, Wah Chiu, Baylor College of Medicine

Chaperonins are essential mediators of many functions in the cell, including protein folding. These are oligomeric machines, of ~ 1 megadalton, with two back-to-back rings enclosing a central cavity that accommodates polypeptide substrates. The Chiu lab has solved near-atomic resolution cryoEM structures of GroEL from bacteria, Mm-Cpn from archaea, and TRiC from eukaryotes. In addition, they have studied them in different nucleotide-binding states and with a variety of substrates. These studies show that chaperonins assume various conformations under different conditions, i.e. structural plasticity is an inherent property of chaperonins.

Our ongoing activity is to determine GroEL cryoEM structure with GroES and carboxylase substrate. In our newly funded program project, we will also pursue the cryoEM characterization of TRiC complexed with mutant huntingtin exon extracted from Huntington's disease (HD) cell model and mouse neuron model. These experiments will provide a mechanistic understanding of how TRiC prevents the progression of aggregation in polyQ, and mutant huntingtin exon as part of our effort to develop a TRiC-like reagent for HD therapeutics. The Bahar lab has experience on the machinery and conformational variability of GroEL-GroES, which will be further analyzed in the light of the new data generated by the Chiu lab, and with the help of novel extensions of elastic network models (ENMs). The computationally predicted dynamics will be compared to the structural variability observed in cryoEM, and the interactions that modulate GroEL allosteric machinery will be identified. Similar methods will be used for TRiC, toward elucidating its spectrum of motions and how its dynamics relates to its chaperoning activities. Target sites whose perturbation may alter/control TRiC machinery will be identified using ProDy tools, toward assisting in the design of new therapeutic strategies against HD. We anticipate a productive interaction between this project and C&SP16 on the identification of neuroprotectives for HD with Dr. Friedlander (UPMC), an expert in the etiology of HD as well as discovery of neuroprotectives against HD.

In progress


21) Functional significance of the dynamics of AMPAR extracellular region

Collaborating Investigators: Ingo H. Greger, MRC Lab of Molecular Biology, Ivet Bahar, QingDe Wang, University of Pittsburgh, Tom M. Bartol, Terry J. Sejnowski, Salk Institute

Ionotropic glutamate receptor (iGluRs) are ligand-gated ion channels that allow for the flow of cations into the postsynaptic cell in response to glutamate binding, thus regulating neurotransmission upon depolarization of the cell membrane. Among iGluR subfamilies, AMPAR and NMDAR play a key role in learning and memory, and in particular the AMPAR is essential to rapid neurotransmission and synaptic plasticity. The Greger and Bahar labs have been productively collaborating in recent years on AMPAR dynamics, first using the NTD dimer structures and more recently the intact tetrameric structures. These studies demonstrated that the NTD domains exhibit structural flexibilities comparable to those of AMPAR NTDs. Furthermore, the global modes of motions predicted by ANM (or ProDy) revealed the propensity of homotetrameric AMPAR to assume more compact forms similar to NMDARs. The validity of these modes of motions were confirmed by cross-linking experiments between NTD sites predicted by ANM to come into close proximity. In the new term, we will first adopt ANM-based analysis to characterize the mode spectrum of the heterotetrameric AMPAR. ProDy analysis already revealed that the O ↔ N transition is enabled by a global ANM mode. We will characterize thoroughly the whole spectrum of motions and generate the energy landscape of Glu2/3 heterotetramer, using the recently introduced extension of coMD. Then we will focus on the ECR motions that induce a pore opening (or cooperative twisting) at the TMD and analyze the conformational events that enable the allosteric coupling between the ECR and the TMD with the help of accelerated MD simulations. In the next phase, we plan to examine the significance of GluA2/3 ECR flexibility in adapting to its interactions with auxiliary proteins such as cornichon homologs, TARPs or in forming clusters, which will be further tested/validated with structural and single-particle tracking methods in the Greger lab.

Publications Resulting from This Work

    • Krieger J, Bahar I, Greger IH (2015) Structure, Dynamics, and Allosteric Potential of Ionotropic Glutamate Receptor N-Terminal Domains Biophys J 109: 1136-48 PMID: 26255587, PMC45761612
    • Dutta A, Krieger J, Lee JY, Garcia-Nafria J, Greger IH, Bahar I (2015) Cooperative Dynamics of Intact AMPA and NMDA Glutamate Receptors: Similarities and Subfamily-Specific Differences Structure 23: 1692-1704 PMID: 26256538, PMC4558295
    • Dutta A, Shrivastava IH, Sukumaran M, Greger IH, Bahar I (2012) Comparative dynamics of NMDA- and AMPA-Glutamate receptor N-terminal domains Structure 20: 1838-49 PMID: 22959625, PMC3496038

In progress


22) Modeling T cell fate decisions

Collaborating Investigators: Penelope Morel, Robin E. C. Lee, James R. Faeder, University of Pittsburgh

The T cell receptor (TCR) is an exquisitely sensitive molecular machine that can translate small differences in peptide/MHC ligands into profoundly different outcomes. Our long-term goal is to understand how T cells perceive signals of differing strengths in order to elucidate how this determines T helper (Th) cell fate. Antigen (Ag) dose affects Th differentiation; and several in vivo models have shown that low Ag dose favors T regulatory (Treg) and Th2 differentiation, whereas high Ag dose induces inflammatory Th1 cells. On the one hand, Treg cells are critical for the maintenance of self-tolerance and the prevention of autoimmune diseases such as diabetes; on the other, tumors induce Treg cells to suppress immune surveillance. Despite the powerful effect on T cell outcome, how these signals contribute to determining Th cell fate remains poorly understood.

By coupling mathematical modeling with detailed in vitro activation profiling we have identified several important feedback loops that control the degree of Akt/mTOR activation. Based on these and other biochemical mechanisms we have developed a rule-based model in BioNetGen and calibrated it to biochemical data for low and high dose stimulation using Bayesian parameter estimation with the ptempest software. This model predicts a sharp threshold for Akt activation with respect to antigen dose that depends on several key variables, most notably PTEN expression level, which exhibit considerable cell-to-cell variation. We are currently conducting experiments in mice that express a fluorescent reporter of TCR activation to determine whether there is a sharp activation threshold as predicted by the model, and in mice heterozygous for PTEN expression to examine how PTEN levels affect this threshold. We will use these data to refine the model and identify additional factors that can affect antigen dose thresholds. In addition, we will use the Amnis ImageStream technology available to the Morel lab and the single-cell live imaging pipeline being developed in the Lee lab to spatially resolve key signaling components, including Akt, PIP3, and NFκB. We will use these data to confirm temporal oscillations observed in preliminary data and to calibrate a coarse-grained spatial model.

Publications Resulting from This Work

    • Hawse WF, Sheehan RP, Miskov-Zinanov N, Menk AV, Kane LP, Faeder JR, Morel PA (2015) Cutting edge: Differential regulation of PTEN by TCR, Akt, and Fox01 Controls CD4+ T Cell fate decisions J Immunol 194: 4615-9 PMID: 25855357, PMC4418530
    • Morel PA, Faeder JR, Hawse WF, Miskov-Zivanov N (2014) Modeling the T cell immune response: a fascinating challenge J Pharmacokinet Pharmacodyn 41: 401-13 PMID: 25155903, PMC4210366
    • Miskov-Zivanov N, Turner MS, Kane LP, Morel PA, Faeder JR (2013) The duration of T cell stimulation is a critical determinant of cell fate and plasticity Sci Signal 6: ra97 PMID: 24194584, PMC4074924

In progress


23) CaMKII Structural Dynamics

Collaborating Investigators: John Kuriyan, University of California-Berkeley, Mary Kennedy, Caltech, Ivet Bahar, James R. Faeder, University of Pittsburgh, Tom M. Bartol, Terry J. Sejnowski, Salk Institute

This C&SP will focus on studying mechanisms of function of the complex biomolecular assembly of CaMKII and its dynamics. We recently showed that the human alpha isoform, CaMKIIalpha, literally acts as a device for propagating the activation by calcium signaling. This is achieved through subunit exchanges between activated CaMKII and holoenzyme assemblies that have not been activated. The spread of activation through these exchanges now suggests that activation can be maintained long after the disappearance of the initial signal. This is remarkable because it reveals a possible mechanism by which CaMKII assists in the induction of long-term memory in the brain. It is clear that the phosphorylation of the calmodulin-recognition element triggers the exchange. But in addition to localized interactions, there is something unique about the overall architecture of CaMKII so that it accommodates different oligomerization states. Starting from early crystal structures and EM images, we have seen that the holoenzyme could assume both dodecameric and tetradecameric forms. Our new observations now show that they coexist. An intriguing observation was: why does the vertical dimer dissociate, rather than a single subunit, or a horizontal dimer? ANM results unambiguously show that the dodecamer structure indeed favors the vertical dimers to be released. The type of analysis performed by ProDy can open the way to a better understanding of CaMKII dynamics. Additionally, coevolution analysis may give us additional information on interfacial interactions that modulate subunit exchange. These studies will help further understanding the dynamics and interactions of CaMKII.

In progress


24) Spatio-temporal cell biology

Collaborating Investigators: Alan Horwitz, Gregory Johnson, Allen Institute, Robert Murphy, Carnegie Mellon University

The goal of the Allen Institute for Cell Science (AICS) is to develop predictive models of cell behavior. The initial goal is to assemble a map of the shapes and positions of all major molecular machines and signaling complexes in human induced Pluripotent Stem (hiPS) cells and capture how they change during their differentiation into cardiomyocytes. hiPS cells will be genome-edited to express multiple fluorescent-tagged markers, and 3D images and short movies will be collected for 1000s of cells at various time points after initiation of differentiation, using a large fluorescence microscopy pipeline. We will combine images with various combinations of markers into a comprehensive 3D model of the undifferentiated stem cell and its changes during differentiation. To achieve this goal, several computational tools must be developed: 1) Models of causal spatio-temporal relationships of subcellular structures. 2) Integrated multi-protein subcellular organization models from images where only a small subset of markers are visible, and 3) Cellular population models.

This project takes advantage of the new CellOrganizer capabilities being developed for point process models and causal inference. It will also draw on the existing capabilities for learning point process models and for constructing models for a full differentiation process from short movies of parts of that process that were previously developed in conjunction with completed C&SP11 and similar to those in our recent study. The collaboration will involve both large-scale testing of CellOrganizer in the context of the AICS images and joint development and refinement of methods.

In progress


25) Efficient parallel sampling at multiple scales using the weighted ensemble strategy

Collaborating Investigators: Lillian T. Chong, University of Pittsburgh, Daniel Zuckerman, Oregon Health and Science University

There is a "silicon ceiling" that ultimately limits many, if not most, types of dynamical biological simulations. That is, even the world's most powerful computers cannot generate sufficiently long simulations, whether for atomistic models of proteins or for realistic models of cell behavior. In many cases, the key events may occur beyond simulation timescales - such as protein folding, conformational transitions of proteins, assembly of protein complexes, or transitions of cell behavior from healthy to pathological states. The work in this C&SP will continue the development of WESTPA, a tool for controlling other software tools: it orchestrates up to thousands of trajectories run natively by other software at any scale (e.g., Gromacs, Amber, BioNetGen, MCell) using a "weighted ensemble" strategy. Not only does WESTPA parallelize the use of dynamics engines - but because of the statistical process by which trajectories are added and removed, WESTPA can obtain estimates of key kinetic as well as equilibrium observables in significantly less computing time than would be required in ordinary parallelization. Our aims are to improve the ease of use and interoperability of WESTPA; to improve its performance and reliability; to demonstrate the effectiveness of WESTPA through a series of "showcase" examples from molecular to cellular scale using a variety of dynamics engines; and to improve instructional materials based on the showcase examples.

Publications Resulting from This Work

    • Zwier MC, Adelman JL, Kaus JW, Pratt AJ, Wong KF, Rego NB, Suárez E, Lettieri S, Wang DW, Grabe M, Zuckerman DM, Chong LT (2015) WESTPA: an interoperable, highly scalable software package for weighted ensemble simulation and analysis J Chem Theory Comput 11: 800-809. PMID: 26392815; PMC4573570
    • Suárez E, Pratt AJ, Chong LT, Zuckerman DM (2016) Estimating First Passage Time Distributions from Weighted Ensemble Simulations and non-Markovian Analyses Protein Science 25: 67-78. PMID: 26131764; PMC4815309
    • Suárez E, Pratt AJ, Chong LT, Zuckerman DM (2016) Estimating First Passage Time Distributions from Weighted Ensemble Simulations and non-Markovian Analyses Protein Science 25: 67-78. PMID: 26131764; PMC4815309

In progress


26) GPCR signaling and interactions. Learning from SILAC/proteomics

Collaborating Investigators: Kunhong (Kevin) Xiao, Ivet Bahar, James R. Faeder, University of Pittsburgh

This focus of this C&SP will focus on the development and application of quantitative proteomics and systems biology approaches to study GPCR (G protein couple receptor) structures, their signaling networks, and regulation. The approach will be to use mass spectrometry-based quantitative proteomics, in combination with systems, chemical, and structural biology to study the receptor structure-function relationship, macromolecular interactions, and the expression and posttranslational modifications of signaling molecules downstream of GPCR. The ultimate goal of this research is to provide a better understanding of structure-function relationship of GPCRs and their related proteins and provide tools and resources for structure-based drug design.

In progress


27) Circuit reconstruction of association cortex

Collaborating Investigators: Wei-Chung Allen Lee, Harvard, Art Wetzel, Greg Hood, Pittsburgh Supercomputing Center

This C&SP with Wei Chung Allen Lee at Harvard Medical School builds upon the successful DBP5 with Clay Reid of Harvard and the Allen Brain Institute during the 2012-2017 funding period. Lee is a former member of Reid's lab, and has continued this line of connectomics research as an independent researcher at HMS using both the original TEMCA microscope, and an improved successor, TEMCA-GT. Our previous collaborative work with Reid and Lee has focused on studying the pattern of anatomical connections among functionally characterized neurons within mouse visual cortex. Lee plans to study next a new cortical area in the mouse. This region of the cortex is hypothesized to be involved in key cognitive functions including decision-making, movement planning, attention, and reward. Recent work has demonstrated that groups of neurons in this brain region are activated sequentially. Neurons fire selectively during discrete epochs of a navigational task, and the resulting sequences of neuronal activity accurately predict the behavioral choice of the animal. Two of the unique aspects of this project are that we will, for the first time, be able to compare an association area with a sensory area and identify conserved and distinguishing features of connectivity motifs underlying their different network dynamics.

Publications Resulting from This Work

    • Lee, W-C. A, Vincent B, Reed M, Graham BJ, Hood G, Glattfelder K, Reid RC (2016) Anatomy and function of an excitatory network in the visual cortex Nature PMID:27018655, PMC4844839

In progress


28) Dynamic modulation of interferon binding affinity as a mechanism to regulate interferon receptor signaling

Collaborating Investigators: Gideon Schreiber, Weizmann Institute, Ivet Bahar, James R. Faeder, University of Pittsburgh

Type I interferons (IFNs) are multifunctional cytokines that mediate/induce diverse cellular responses, including both innate and adaptive immune responses, stimulation of antiviral responses, and cancer surveillance, upon forming a ternary complex with two surface receptors, IFNAR1 and IFNAR2. The activities of IFN-a subtypes correlate with their affinities to bind to IFNAR1 and IFNAR2. While the Schreiber lab made seminal contributions to understanding the molecular basis of IFNARs, the mechanism of regulation of differential IFN activities through interactions with IFNAR1 and 2, remains unclear. Our integrated computational (TR&D1) and experimental preliminary studies point to the significance of the intrinsic dynamics in modulating binding affinity. We adopted a closely integrated computational/experimental strategy that yielded promising results, which we are currently further pursuing and that illustrate the adaptability of proteins to different bound states or to sequence variations/mutations via their softest modes of motion. Further cross-linking, fluorescence quenching and gene induction experiments will be conducted in the Schreiber lab, in close coordination with TR&D1 computational studies at the Bahar lab.

In progress


29)  Spatiotemporal models of autophagy

PI: Joern Dengjel, University of Frieberg Collaborator: R. Murphy, Carnegie Mellon University

Macroautophagy is thought to be mainly cytoprotective, although it also has been linked to autophagic cell death (Feng et al., 2015). Autophagosomes are de-novo formed double membrane vesicles, which enwrap cytoplasm and target it for lysosomal degradation. Their biogenesis is complex and still holds many open questions (Lamb et al., 2013). We will generate wide-field as well as confocal fluorescent images of forming autophagosomes using different stress conditions such as amino acid and growth factor starvation. We will use cells stably expressing the autophagosomal marker MAP1LC3B (LC3) coupled to mRFP-GFP which allows the analysis of autophagosome turnover (Kimura et al., 2007). Also, we will generate images of all six fluorescently tagged isoforms of LC3 to determine if different isoforms have an influence on site of formation of autophagosomes and size/dynamics of vesicles. We will use these to construct image-derived models to enable comparison of different conditions.

Literature Cited

    • Feng, Y., Z. Yao, and D.J. Klionsky. (2015) How to control self-digestion: transcriptional, post-transcriptional, and post-translational regulation of autophagy. Trends Cell Biol.25:354-363.
    • Kimura, S., T. Noda, and T. Yoshimori. (2007) Dissection of the autophagosome maturation process by a novel reporter protein, tandem fluorescent-tagged LC3. Autophagy. 3:452-460.
    • Lamb, C.A., T. Yoshimori, and S.A. Tooze. 2013. The autophagosome: origins unknown, biogenesis complex. Nat Rev Mol Cell Biol.14:759-774.

Publications Resulting from This Work

    • C. Gretzmeier, S. Eiselein, G. R. Johnson, R. Engelke, H. Nowag, M. Zarei, V. Küttner, A. C. Becker, K. T. G. Rigbolt, M. Høyer-Hansen, J. S. Andersen, C. Münz, R. F. Murphy, and J. Dengjel (2017) Degradation of protein translation machinery by amino acid starvation-induced macroautophagy. Autophagy, in press.

In progress


30) Modeling of Tetrahymena basal body dynamics

Collaborating Investigators: Chad G. Pearson, University of Colorado, Robert Murphy, Carnegie Mellon University

Mapping basal body assembly in Tetrahymena cells.
The Pearson laboratory works on understanding how basal bodies assemble through the cell cycle in Tetrahymena cells. While most of this is directed at specific molecules in the assembly process, the Pearson lab is also interested in mapping where, spatially within Tetrahymena cells, basal bodies assemble and to understand where they go / separate upon assembly. This will allow us to explore the pattern that basal bodies organize into for normal ciliary beating and hydrodynamic flow. Therefore, the goals of this collaboration are the following: 1) Analyze the data previously acquired by the Pearson lab to map the stereotypical organization of basal bodies relative to a) the cell cycle and b) the cell size, 2) Modeling localized basal body separation to establish the temporal dynamics of new, daughter basal body separation from old, mother basal body separation, and 3) Live imaging of basal body separation to visualize “parts of” basal body separation during assembly and maturation and use CellOrganizer (TR&D4) to reassemble our “parts” to the full view.

Mapping the Tetrahymena ciliary beat pattern.
We want to measure/reconstruct the ciliary beat pattern in Tetrahymena cells. Using DIC imaging and high temporal resolution imaging (~500fps), we have low resolution imaging of the Tetrahymena ciliary beat pattern. We will marry these data with fixed timepoint fluorescence data to reconstruct a generative model of the Tetrahymena cilia, basal body and associated structures during the ciliary beat stroke.

There are two basic methods in CellOrganizer that will be used for this project. The first is to automatically assemble collections of static images of fixed cells that are at various points in the cell cycle into a model of the changes in basal body arrangement through the cycle. This will be done by measuring pairwise distances between images and assembling them into a shape space as we have done previously for cell and nuclear shape analysis. The second is to take short movies of living cells and assemble them into a dynamic shape space as has been done previously for the dynamics of H1299 cell shape changes. The main difference is that the distance measure will be based on basal body positions rather than overall cell shape.

In progress


31) Receptor Clustering and the Kinetics of T Cell-Mediated Killing

Collaborating Investigators: Yuri Sykulev, Thomas Jefferson University, James R. Faeder, University of Pittsburgh

The goal of the project is to develop the model Membrane systems that can recapitulate clustering of immune receptors in a controlled manner. We propose to utilize biodegradable nanolipoprotein particles (NLPs) as a universal platform to mimic molecular Membrane clustering. NLPs are self-assembled in solution to form discoidal Nanostructures containing Lipid Bilayers stabilized at the perimeter by apolipoprotein molecules. The size of the NLP ranges from 8 to 30 nm that allow capturing up to 50 molecules of soluble Ligands and enable us to achieve model cluster size and density close to physiological. We will use the NLPs to assemble pMHC and other MembraneLigands into model Membrane patches. We will study how changes in the Ligands density and composition of these patches Affect Binding of the model Membrane patches to live T cells and the kinetics and magnitude of TCR-mediated Signaling. This will provide a basis for the engineering of aAPC bearing the model Membrane patches incorporated into Lipid Bilayers covering the surface of glass beads. Such aAPCs will allow us to calibrate the strength of T cell stimulation. We will utilize these novel aAPCs to vary the strength of stimulation of naïve CD8+ T cells derived from OT-1 TCR Transgenic Mice in order to induce different Subsets of activated T cells with the same Specificity. Building of aAPCs is expected to enable us to expand T cells with instructional programs that allow T cells to persist, function, and migrate in a desired fashion after Adoptive Transfer. The experimental Data will also provide the basis for building a mathematical model to characterize how clustering of Ligands on an APC surface determines speed, sensitivity and Discrimination of the pMHC I ligand by activated and naïve CD8 T cells.

In progress


32) Integrating compartmental rule-based modeling into VCell

Collaborating Investigators: James C. Shcaff, Mikhail L. Blinov, Ion I. Moraru, Lew Loew, University of Connecticut, James R. Faeder, University of Pittsburgh

This C&SP is a collaboration with The National Resource for Cell Analysis and Modeling (NRCAM), which develops new technologies for modeling cell biological processes. The technologies are integrated through Virtual Cell (VCell), a problem solving environment built on a central database and disseminated as a client-server application. It will develop methods to simulate multi-molecular interactions at a sub-cellular scale for problems where molecular shape and cellular geometry influence the system behavior. MMBioS investigators will collaborate to integrate compartmental rule-based modeling into the VCell platform.

In progress


33) Integration of rule-based modeling capabilities with pySB modeling platform

Collaborating Investigators: Carlos Lopez, Vanderbilt, James R. Faeder, University of Pittsburgh

This C&SP is a collaboration that will integrate rule-base modeling into the PySB modeling platform. PySB is a framework for building mathematical models of biochemical systems as Python programs. PySB abstracts the complex process of creating equations describing interactions among multiple proteins or other biomolecules into a simple and intuitive domain specific programming language, which is internally translated into BioNetGen or Kappa rules and from there into systems of equations. PySB makes it straightforward to divide models into modules and to call libraries of reusable elements (macros) that encode standard biochemical actions. These features promote model transparency, reuse and accuracy. PySB also interoperates with standard scientific Python libraries such as NumPy, SciPy and SymPy, enabling model simulation and analysis.

In progress


34) Development of a high-level, rule-based whole-cell modeling language

Collaborating Investigators: Jonathan Karr, Mt. Sinai, James R. Faeder, University of Pittsburgh

Despite decades of research and the growing wealth of data, we do not have a unified understanding of cell biology. Whole-cell (WC) models that represent every gene function and that predict the cell cycle dynamics of every molecular species have the potential to unify our understanding of cell biology. Recently, Karr and others developed the first WC model. This was achieved by combining multiple mathematically distinct submodels of individual cellular pathways, and by developing custom software to simulate the model. However, this ad hoc approach is not easily reusable, expandable, or reproducible. New WC modeling tools are needed to build better models that can guide bioengineering and medicine.

This project will develop a high-level rule-based WC modeling language and a rule-based multi-algorithm simulator. The language will enable researchers to compactly describe biological processes, such as the translation of each amino acid, using high-level rule patterns and sequence, genomic, and biochemical data. WC-Lang will provide high-level data-based rules to succinctly describe large reactions networks in terms of reaction patterns and sequence, genomic, and other data, that will be encoded in BioNetGen. This high-level language will make it easier for researchers to develop and communicate WC models, and the BioNetGen implementation will make it possible to rapidly develop a reusable and efficient WC model simulator that will take advantage of BioNetGen’s network generation as well as the network-free simulation capabilities of NFsim. WC-Lang will be implemented as a Python library which encodes the high-level rules into lower-level BioNetGen rules. Use of rules to encode molecular interactions will open many possibilities for the visualization of WC models using contact maps and regulatory graphs. The new language along with an accompanying expanded simulator will enable researchers to build much larger and more accurate WC models that can guide bioengineering and precision medicine.

In progress


Copyright © 2018 National Center for Multiscale Modeling of Biological Systems. All Rights Reserved.