• 1_2fa67f482133f1c934235b73c2a03954_m.jpg
  • 2_e0a70f72bdae9885bfc32d7cd19a26a1_m.jpg
  • 3_94d43e327d9303539cb1e2aac7032668_m.jpg

C&SP30: Modeling of Tetrahymena basal body dynamics

A. Collaborating Investigators: Chad G. Pearson, 1 Robert F. Murphy2

B. Institutions: 1University of Colorado and 2Carnegie Mellon University

C. Funding Status of Project: 5R01GM099820-04 Mechanisms of Centriole Assembly and Stability (Pearson) 9/24/12-7/31/17

 D. Biomedical Research Problem:

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 109 to map the stereotypical organization of basal bodies relative to a) the cell cycle and b) the cell size. This will then be the basis for future work to model in ciliary beating and the efficiency of cellular motility. Importantly, because larger cells during cell division swim more slowly, we predict that the slower motility reflects changes in basal body positioning. We will also use these analyses to establish a baseline for future mutant analyses.

2. Modeling localized basal body separation. This is a micro-version of aim 1 and will establish the temporal dynamics of new, daughter basal body separation from old, mother basal body separation. Such studies are currently ongoing in the lab to assess protein incorporation into basal bodies as new basal bodes separate from their mother basal body and mature. We currently make the assumption (likely incorrect) that as basal bodies separate from their mother that this is a linear proxy for their age. We will use CellOrganizer110 to define the temporal separation of basal bodies. These studies will also be correlated with ongoing studies following the dynamics of protein incorporation at newly forming basal bodies.

3. Live imaging of basal body separation. In a direct test of aim 2, we will visualize “parts of” basal body separation during assembly and maturation and use CellOrganizer110 (TR&D4) to reassemble our “parts” to the full view. This will reveal the rate of basal body separation relative to basal body protein incorporation and structural maturation.

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.

 

E. Methods and Procedures. 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.85 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.

 

C&SP28: Dynamic modulation of interferon binding affinity as a mechanism to regulate interferon receptor signaling

A. Collaborating Investigators: Gideon Schreiber,1 Ivet Bahar,2 James R Faeder2

B. Institutions: 1Weizmann Institute and 2University of Pittsburgh

C. Funding Status of Project: ISF-I-Core- Center for Integrated Structural Cell Biology 2013-2018; ISF grant: Protein binding and catalysis in crowded environments: from in vitro to in vivo (2014-2018)

 

Fig VIII.8. Intrinsic flexibility of interferon receptor IFNAR1, predicted by ProDy ANM, is (a) essential to facilitating the formation of ternary complex with IFN and IFNAR2, (b) explains the distance changes observed in FRET. Pairs of residues making interdomain contacts, identified by GNM to control IFNAR1 global motions.

D. Biomedical Research Problem:

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 (Fig VIII.8a).95-97 The activities of IFN-a subtypes correlate with their affinities to bind to IFNAR1 and IFNAR2.98 While the Schreiber lab made seminal contributions to understanding the molecular basis of IFNARs,95,96,98-105 the mechanism of regulation of differential IFN activities through interactions with IFNAR1 and 2,102 remains unclear. Our integrated computational (TR&D1) and experimental preliminary studies described below point to the significance of the intrinsic dynamics in modulating binding affinity.

E. Methods and Procedures: We adopted a closely integrated computational/experimental strategy that yielded promising results, which we are currently further pursuing: (i) we retrieved and analyzed existing structures, including a model constructed for human IFNAR1 ectodomain (EC) in the unbound state using data from two PDB structures,105,106 (ii) we determined, using the GNM core function, the hinge sites that control the relative movements of IFNAR1 EC four subdomains SD1-SD4, (iii) we identified 5 residue pairs near those regions, which exhibit large fluctuations in their inter-residue distances, <(DRij)2>, during global motions; we hypothesized that those at the interface between subdomains SD3 and SD4 would be critical to enabling subdomain rearrangements that modulate, if not optimize, IFN binding (Fig VIII.8b), and obstructing their adaptability by locking the distance between those pairs, and in particular the pair with the largest moment arm (e.g. G162-F267) with respective to the hinge site, could impair the function, (iv) in vitro binding experiments and gene induction activity assays by the Schreiber lab confirmed that hindering the conformational flexibility of IFNAR1 at those regions (by cysteine-trapping) decreased binding affinity, and suppressed activity, and the relative strengths of these effects matched our predicted rank-ordered list, (v) ANM-predicted global (energetically softest) mode of motion also showed a highly cooperative flexing movement of the entire IFNAR1 EC (Fig VIII.8c), with end-to-end distance changes in accord with FRET experiments107, (vi) ANM examination of the murine IFNAR1 structure resolved106 in IFNb-bound state, demonstrated that this structure has the intrinsic ability to readily reconfigure toward the human counterpart. This result is important: it shows that the seemingly different human and murine structures are simply alternative easily exchangeable conformers, stabilized by the different sequences and/or bound substrates, and that IFNAR1 EC favors such structural shifts as required for function, consistent with ProDy predictions. This is another illustration of the adaptability of proteins to different bound states or to sequence variations/mutations via their softest modes of motion.108 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.

 

C&SP27: Circuit reconstruction of association cortex

A. Collaborating Investigators: Wei-Chung Allen Lee,1 Art Wetzel,2 Greg Hood2

B. Institutions: 1Harvard Medical School (HMS) and 2Pittsburgh Supercomputing Center

C. Funding Status: NIDCD 5R01-DC013622-03" Network Anatomy of Olfactory Processing (Lee) 8/9/2013 - 8/31/2016; DP2 OD022472-01 'Network anatomy of behavioral choice' (PI: Lee) Pending

D. Biomedical Research Problem:

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,87 and an improved successor, TEMCA-GT. Our previous collaborative work87,88with 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.89-93 Recent work has demonstrated that groups of neurons in this brain region are activated sequentially.94 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.94 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.

 

 

Fig VIII.7. Preliminary acquisition of a10 gigapixel section from a cortical area of mouse brain involved in decision-making investigated by the Lee lab.

 

E. Methods and Procedures: We will trace a portion of the neural connectome within a region of tissue about 1mm3 in size. To do this will require ~1Petabyte of raw image data comprising ~25,000 sections. This volume will encompass 4 cortical columns, which are hypothesized to be modular, canonical, local circuits. Prior datasets have been up to Terabyte in size, so this will necessitate an order-of-magnitude improvement in our ability to acquire EM images and analyze them. Several technical hurdles must be overcome, and we focus here on the issues affecting alignment of these raw images into a form which can then be traced, either manually or in a semi-automated way. Specifically, we must reduce the amount of manual intervention so that the effort is much less per gigapixel than on prior datasets. Most of the manual intervention that has been necessary in aligning prior datasets has been associated with artifacts introduced during sectioning and handling (e.g., tissue folds and tears), so an important step in reducing labor has been the development of an improved automated pipeline for sample handling. This new workflow is centered around a novel substrate that allows automated collection of 1000s of serial thin sections from the microtome, and imaging them in the vacuum chamber of the TEMCA-GT transmission electron microscope. TEMCA-GT's improved stage allows for more accurate positioning of the samples, and this information can be fed into our AlignTK pipeline to constrain the model’s initial conditions and reduce the probability of registration errors (Fig VIII.7). This new approach will allow 10 times more data (5-10 TB) to be imaged per day, resulting in an increased workload for the alignment software. AlignTK is being converted to use GPUs for its most computationally-intensive tasks (TR&D2 Subaim 2.1). With current Harvard Medical School and MMBioS facilities, we estimate we will be able to register 64 Mpixels/sec of image data, with the ability to process a ~1PB dataset in roughly 6 months. The main concern will be to increase the robustness of the software to deal well with atypical image content (e.g. voids within capillaries), and to make it easy to incrementally process the data as it is generated. To increase robustness and throughput, we may also perform registration using MMBioS' SWiFT-IR software, which has heretofore been used primarily on SEM images, but should also perform well with TEM data.

 

C&SP24: Spatio-temporal cell biology

A. Collaborating Investigators: Alan Horwitz,1 Gregory Johnson,1 and Robert Murphy2

B. Institutions: 1Allen Institute for Cell Science and 2Carnegie Mellon University

C. Funding Status of Project: Data generation fully funded by $100 million gift from the Allen Foundation

 

D. Biomedical Research Problem:

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. Understanding how cell organization depends on the location and structure of other cellular components is an unknown and critical to understanding cell behavior and organization. These representations will allow others to model and understand the mechanisms underlying the ordered changes in cellular organization during differentiation. These models will be constructed from single images and short time series, and identify cellular spatio-temporal relationships at the mesoscale.

 

2) Integrated multi-protein subcellular organization models from images where only a small subset of markers are visible. Generating images of many labeled proteins in single cells is not feasible due to the limited number of non perturbing XFP colors. Therefore, it is necessary to build models of subcellular organization by combining information from images containing different, but overlapping, subsets of labeled structures.

 

3) Population modeling. Modeling the origin of the variation among cells in a population and how it may be based on position in cell cycle, partial early differentiation, and unknown effects of neighboring cells is necessary to understand how large populations of cells develop and mature. To parse out these relationships, we will need samples of large numbers of "similar" cells in the context of their differentiation, as well as how neighboring cells are related to cell organization.

 

E. Methods and Procedures:

This project takes advantage of the new CellOrganizer capabilities that will be developed in TR&D4, especially Aim 1 on point process models and causal inference. It will also draw on the existing capabilities for learning point process models84,85 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.86 The collaboration will involve both large-scale testing of CellOrganizer in the context of the AICS images and joint development and refinement of methods.

 

Research Highlights

 

Direct coupling of oligomerization and oligomerization-driven endocytosis of the dopamine transporter to its conformational mechanics and activity

The Bahar (TR&D1) and Sorkin (DBP3) labs published an article in the Journal of Biological Chemistry, selected as one of JBC's "Editors' Picks. Our results demonstrate a direct coupling between conformational dynamics of DAT, functional activity of the transporter and its oligomerization leading to endocytosis. The high specificity of such coupling for DAT makes the TM4-9 hub a new target for pharmacological modulation of DAT activity and subcellular localization. (Read more)

 

 

Differences in the intrinsic spatial dynamics of the chromatin contribute to cell differentiation

Comparison with RNA-seq expression data reveals a strong overlap between highly expressed genes and those distinguished by high mobilities in the present study, in support of the role of the intrinsic spatial dynamics of chromatin as a determinant of cell differentiation. (Read more)

 

 

Nanoscale co-organization and coactivation of AMPAR, NMDAR, and mGluR at excitatory synapses

Work by TR&D2 Investigators and collaborators provide insights into the nanometer scale organization of postsynaptic glutamate receptors using a combination of dual-color superresolution imaging, electrophysiology, and computational modeling. (Read more)

 

 

Parallel Tempering with Lasso for model reduction in systems biology

TR&D3 Investigators and collaborators develop PTLasso, a Bayesian model reduction approach that combines Parallel Tempering with Lasso regularization, to automatically extract minimal subsets of detailed models that are sufficient to explain experimental data. On both synthetic and real biological data, PTLasso is an effective method to isolate distinct parts of a larger signaling model that are sufficient for specific data. (Read more)

 

Image-derived models of cell organization changes during differentiation and drug treatments

Our work on modeling PC12 cells undergoing differentiation into neuron-like morphologies (under C&SP11, completed) has been published in Molecular Biology of the Cell. We have also made the large dataset of 3D images collected in that study available through Dryad. (Read more)

 

Monoamine transporters: structure, intrinsic dynamics and allosteric regulation

T&RD1 investigators Mary Cheng and Ivet Bahar published an invited review article in Nature Structural & Molecular Biology, addressing recent progress in the elucidation of the structural dynamics of MATs and their conformational landscape and transitions, as well as allosteric regulation mechanisms. (Read more)

Trimerization of dopamine transporter triggered by AIM-100 binding

The Bahar (TR&D1) and Sorkin (DBP3) labs explored the trimerization of dopamine transporter (DAT) triggered by a furopyrimidine, AIM-100, using a combination of computational and biochemical methods, and single-molecule live-cell imaging assays. (Read more)

Pre-post synaptic alignment through neuroligin-1 tunes synaptic transmission efficiency

TR&D2 investigators and collaborators describe organizing role of neuroligin-1 to align post-synaptic AMPA Receptors with pre-synaptic release sites into trans-synaptic “nano-columns” to enhance signaling.(Read more)

Inferring the Assembly Network of Influenza Virus

In an article in PLoS Computational Biology, MMBioS TR&D4 members Xiongto Ruan and Bob Murphy collaborated with Seema Lakdawala to address this question of the assembly network of the Influenza virus.(Read more)

PINK1 Interacts with VCP/p97 and Activates PKA to Promote NSFL1C/p47 Phosphorylation and Dendritic Arborization in Neurons

Our findings highlight an important mechanism by which proteins genetically implicated in Parkinson’s disease (PD; PINK1) and frontotemporal dementia (FTD; VCP) interact to support the health and maintenance of neuronal arbors.(Read more)

Improved methods for modeling cell shape

In a recent paper in Bioinformatics, Xiongtao Ruan and Bob Murphy of TR&D4 addressed the question of how best to model cell and nuclear shape.(Read more)

New tool to predict pathogenicity of missense variants based on structural dynamics: RHAPSODY

We demonstrated that the analysis of a protein’s intrinsic dynamics can be successfully used to improve the prediction of the effect of point mutations on a protein functionality. This method employs ANM/GNM tools (Read more)

New method for investigating chromatin structural dynamics.

By adapting the Gaussian Network Model (GNM) protein-modeling framework, we were able to model chromatin dynamics using Hi-C data, which led to the identification of novel cross-correlated distal domains (CCDDs) that were found to also be associated with increased gene co-expression.  (Read more)

 

Structural elements coupling anion conductance and substrate transport identified

We identified an intermediate anion channeling state (iChS) during the global transition from the outward facing (OF) to inward facing state (IFS). Our prediction was tested and validated by experimental study conducted in the Amara lab (NIMH). Critical residues and interactions were analyzed by SCAM, electrophysiology and substrate uptake experiments (Read more)

Integrating MMBioS technologies for multiscale discovery

TR&D teams driven by individual DBPs are naturally joining forces, integrating their tools to respond to the needs of the DBP, and creating integrative frameworks for combining structural and kinetic data and computing technologies at multiple scales.  (Read more)

 

Large scale visualization of rule-based models.

Signaling in living cells is mediated through a complex network of chemical interactions. Current predictive models of signal pathways have hundreds of reaction rules that specify chemical interactions, and a comprehensive model of a stem cell or cancer cell would be expected to have many more. Visualizations of rules and their interactions are needed to navigate, organize, communicate and analyze large signaling models.  (Read more)

Integration of MCellR into MCell/CellBlender

Using spatial biochemical models of SynGAP/PSD95, MMBioS investigators were able to merge the MCellR code-base with the MCell code-base and validate its utility and correctness of this sophisticated technology now easily accessible through the MCell/CellBlender GUI.  (Read more)

 

csp29

Causal relationships of spatial distributions of T cell signaling proteins

The idea is to identify a relationship in which a change in the concentration of one protein in one cell region consistently is associated with a change in the concentration of another protein in the same or a different region. We used the data from our Science Signaling paper reported last year to construct a model for T cells undergoing stimulation by both the T cell receptor and the costimulatory receptor. (Read more...)

T-Cell Receptor Signaling

BioNetGen modeling helps reveal immune system response decision

To attack or to let be is an important decision that our immune systems must make to protect our bodies from foreign invaders or protect bodily tissues from an immune attack. Using modeling and experiments, we have painted a sharper picture of how T cells make these critical decisions.  (Read more)

 

 

distancecell

Tools for determining the spatial relationships between different cell components

An important task for understanding how cells are organized is determining which components have spatial patterns that are related to each other.Read more

 

4d rtd

Pipeline for creation of spatiotemporal maps

Using a combination of diffeomorphic methods and improved cell segmentation, we developed a CellOrganizer pipeline for use in DPB4 to construct models of the 4D distributions of actin and 8 of its regulators during the response of T cells to antigen presentation. Read more

 

Multi-scale Hybrid Methodology

The hybrid methodology, coMD, that we have recently developed [1] has been recently extended to construct the energy landscape near the functional states of LeuT (Fig 1) [2]. This is the first energy landscape constructed for this NSS family member. Read more

 


Insights into the cooperative dynamics of AMPAR

Comparative analysis of AMPAR and NMDAR dynamics reveals striking similarities, opening the way to designing new modulators of allosteric interactions. Read more


Improved Sampling of Cell-Scale Models using the WE Strategy

The WE strategy for orchestrating a large set of parallel simulations has now been extended to spatially resolved cell-scale systems. The WESTPA implementation of WE has been used to control MCell simulations, including models built using a BioNetGen-CellOrganizer pipeline for situating complex biochemistry within spatially realistic cell models. Read more

Mouse visual cortex
Anatomy and Function of an Excitatory Network in the Visual Cortex

MMBioS researcher Greg Hood’s collaboration with Wei-Chung Allen Lee of Harvard University and R. Clay Reid of the Allen Institute for Brain Science concerning the reconstruction of an excitatory nerve-cell network in the mouse brain cortex at a subcellular level using the AlignTK software has been published in Nature. Read more

 

Molecular Mechanism of Dopamine Transport by hDAT

Dopamine transporters (DATs) control neurotransmitter dopamine (DA) homeostasis by reuptake of excess DA, assisted by sodium and chloride ions. The recent resolution of DAT structure (dDAT) from Drosophila permits us for the first time to directly view the sequence of events involved in DA reuptake in human DAT (hDAT). Read more

 

 

 

figure good 170Synaptic Facilitation Revealed

An investigation of several mechanisms of short-term facilitation at the frog neuromuscular junction concludes that the presence of a second class of calcium sensor proteins distinct from synaptotagmin can explain known properties of facilitation. Read more

 

langmead2 200Sparse Graphical Models of Protein:Protein Interactions

DgSpi is a new method for learning and using graphical models that explicitly represent the amino acid basis for interaction specificity and extend earlier classification-oriented approaches to predict ΔG of binding. Read more

 

Picture1 180Advancing Parallel Bio-simulations

A new non-Markovian analysis can eliminate bias in estimates of long-timescale behavior, such as the mean first-passage time for the dissociation of methane molecules in explicit solvent. Read more

 

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