Opportunities and Challenges in Building a Spatiotemporal Multi-scale Model of the Human Pancreatic ß Cell.

The construction of a predictive model of an entire eukaryotic cell that describes its dynamic structure from atomic to cellular scales is a grand challenge at the intersection of biology, chemistry, physics, and computer science. Having such a model will open new dimensions in biological research and accelerate healthcare advancements. Developing the necessary experimental and modeling methods presents abundant opportunities for a community effort to realize this goal. Here, we present a vision for creation of a spatiotemporal multi-scale model of the pancreatic ß-cell, a relevant target for understanding and modulating the pathogenesis of diabetes.

Characterization of haloacid dehalogenase superfamily acid phosphatase from Staphylococcus lugdunensis.

The haloacid dehalogenase superfamily implicated in bacterial pathogenesis comprises different enzymes having roles in many metabolic pathways. Staphylococcus lugdunensis, a Gram-positive bacterium, is an opportunistic human pathogen causing infections in the central nervous system, urinary tract, bones, peritoneum, systemic conditions and cutaneous infection. The haloacid dehalogenase superfamily proteins play a significant role in the pathogenicity of certain bacteria, facilitating invasion, survival, and proliferation within host cells. The genome of S. lugdunensis encodes more than ten proteins belonging to this superfamily. However, none of them have been characterized. The present work reports the characterization of one of the haloacid dehalogenase superfamily proteins (SLHAD1) from Staphylococcus lugdunensis. The functional analysis revealed that SLHAD1 is a metal-dependent acid phosphatase, which catalyzes the dephosphorylation of phosphorylated metabolites of cellular pathways, including glycolysis, gluconeogenesis, nucleotides, and thiamine metabolism. Based on the substrate specificity and genomic analysis, the physiological function of SLHAD1 in thiamine metabolism has been tentatively assigned. The crystal structure of SLHAD1, lacking 49 residues at the C-terminal, was determined at 1.7 Å resolution with a homodimer in the asymmetric unit. It was observed that SLHAD1 exhibited time-dependent cleavage at a specific point, occurring through a self-initiated process. A combination of bioinformatics, biochemical, biophysical, and structural studies explored unique features of SLHAD1. Overall, the study revealed a detailed characterization of a critical enzyme of the human pathogen Staphylococcus lugdunensis, associated with several life-threatening infections.

Bayesian metamodeling of complex biological systems across varying representations.

Comprehensive modeling of a whole cell requires an integration of vast amounts of information on various aspects of the cell and its parts. To divide and conquer this task, we introduce Bayesian metamodeling, a general approach to modeling complex systems by integrating a collection of heterogeneous input models. Each input model can in principle be based on any type of data and can describe a different aspect of the modeled system using any mathematical representation, scale, and level of granularity. These input models are 1) converted to a standardized statistical representation relying on probabilistic graphical models, 2) coupled by modeling their mutual relations with the physical world, and 3) finally harmonized with respect to each other. To illustrate Bayesian metamodeling, we provide a proof-of-principle metamodel of glucose-stimulated insulin secretion by human pancreatic ß-cells. The input models include a coarse-grained spatiotemporal simulation of insulin vesicle trafficking, docking, and exocytosis; a molecular network model of glucose-stimulated insulin secretion signaling; a network model of insulin metabolism; a structural model of glucagon-like peptide-1 receptor activation; a linear model of a pancreatic cell population; and ordinary differential equations for systemic postprandial insulin response. Metamodeling benefits from decentralized computing, while often producing a more accurate, precise, and complete model that contextualizes input models as well as resolves conflicting information. We anticipate Bayesian metamodeling will facilitate collaborative science by providing a framework for sharing expertise, resources, data, and models, as exemplified by the Pancreatic ß-Cell Consortium.

Assessment of scoring functions to rank the quality of 3D subtomogram clusters from cryo-electron tomography.

Cryo-electron tomography provides the opportunity for unsupervised discovery of endogenous complexes in situ. This process usually requires particle picking, clustering and alignment of subtomograms to produce an average structure of the complex. When applied to heterogeneous samples, template-free clustering and alignment of subtomograms can potentially lead to the discovery of structures for unknown endogenous complexes. However, such methods require scoring functions to measure and accurately rank the quality of aligned subtomogram clusters, which can be compromised by contaminations from misclassified complexes and alignment errors. Here, we provide the first study to assess the effectiveness of more than 15 scoring functions for evaluating the quality of subtomogram clusters, which differ in the amount of structural misalignments and contaminations due to misclassified complexes. We assessed both experimental and simulated subtomograms as ground truth data sets. Our analysis showed that the robustness of scoring functions varies largely. Most scores were sensitive to the signal-to-noise ratio of subtomograms and often required Gaussian filtering as preprocessing for improved performance. Two scoring functions, Spectral SNR-based Fourier Shell Correlation and Pearson Correlation in the Fourier domain with missing wedge correction, showed a robust ranking of subtomogram clusters without any preprocessing and irrespective of SNR levels of subtomograms. Of these two scoring functions, Spectral SNR-based Fourier Shell Correlation was fastest to compute and is a better choice for handling large numbers of subtomograms. Our results provide a guidance for choosing an accurate scoring function for template-free approaches to detect complexes from heterogeneous samples.

Auto-segmentation and time-dependent systematic analysis of mesoscale cellular structure in ß-cells during insulin secretion.

The mesoscale description of the subcellular organization informs about cellular mechanisms in disease state. However, applications of soft X-ray tomography (SXT), an important approach for characterizing organelle organization, are limited by labor-intensive manual segmentation. Here we report a pipeline for automated segmentation and systematic analysis of SXT tomograms. Our approach combines semantic and first-applied instance segmentation to produce separate organelle masks with high Dice and Recall indexes, followed by analysis of organelle localization based on the radial distribution function. We demonstrated this technique by investigating the organization of INS-1E pancreatic ß-cell organization under different treatments at multiple time points. Consistent with a previous analysis of a similar dataset, our results revealed the impact of glucose stimulation on the localization and molecular density of insulin vesicles and mitochondria. This pipeline can be extended to SXT tomograms of any cell type to shed light on the subcellular rearrangements under different drug treatments.

An intensity-based post-processing tool for 3D instance segmentation of organelles in soft X-ray tomograms.

Investigating the 3D structures and rearrangements of organelles within a single cell is critical for better characterizing cellular function. Imaging approaches such as soft X-ray tomography have been widely applied to reveal a complex subcellular organization involving multiple inter-organelle interactions. However, 3D segmentation of organelle instances has been challenging despite its importance in organelle characterization. Here we propose an intensity-based post-processing tool to identify and separate organelle instances. Our tool separates sphere-like (insulin vesicle) and columnar-shaped organelle instances (mitochondrion) based on the intensity of raw tomograms, semantic segmentation masks, and organelle morphology. We validate our tool using synthetic tomograms of organelles and experimental tomograms of pancreatic ß-cells to separate insulin vesicle and mitochondria instances. As compared to the commonly used connected regions labeling, watershed, and watershed + Gaussian filter methods, our tool results in improved accuracy in identifying organelles in the synthetic tomograms and an improved description of organelle structures in ß-cell tomograms. In addition, under different experimental treatment conditions, significant changes in volumes and intensities of both insulin vesicle and mitochondrion are observed in our instance results, revealing their potential roles in maintaining normal ß-cell function. Our tool is expected to be applicable for improving the instance segmentation of other images obtained from different cell types using multiple imaging modalities.

Soft X-ray tomography to map and quantify organelle interactions at the mesoscale.

Inter-organelle interactions are a vital part of normal cellular function; however, these have proven difficult to quantify due to the range of scales encountered in cell biology and the throughput limitations of traditional imaging approaches. Here, we demonstrate that soft X-ray tomography (SXT) can be used to rapidly map ultrastructural reorganization and inter-organelle interactions in intact cells. SXT takes advantage of the naturally occurring, differential X-ray absorption of the carbon-rich compounds in each organelle. Specifically, we use SXT to map the spatiotemporal evolution of insulin vesicles and their co-localization and interaction with mitochondria in pancreatic ß cells during insulin secretion and in response to different stimuli. We quantify changes in the morphology, biochemical composition, and relative position of mitochondria and insulin vesicles. These findings highlight the importance of a comprehensive and unbiased mapping at the mesoscale to characterize cell reorganization that would be difficult to detect with other existing methodologies.

A new visual design language for biological structures in a cell.

As part of a project to build a spatiotemporal model of the pancreatic ß-cell, we are creating an immersive experience called "World in a Cell" that can be used to integrate and create new educational tools. To do this, we have developed a new visual design language that uses tetrahedral building blocks to express the structural features of biological molecules and organelles in crowded cellular environments. The tetrahedral language enables more efficient animation and user interaction in an immersive environment.

Visualizing insulin vesicle neighborhoods in ß cells by cryo-electron tomography.

Subcellular neighborhoods, comprising specific ratios of organelles and proteins, serve a multitude of biological functions and are of particular importance in secretory cells. However, the role of subcellular neighborhoods in insulin vesicle maturation is poorly understood. Here, we present single-cell multiple distinct tomogram acquisitions of ß cells for in situ visualization of distinct subcellular neighborhoods that are involved in the insulin vesicle secretory pathway. We propose that these neighborhoods play an essential role in the specific function of cellular material. In the regions where we observed insulin vesicles, a measurable increase in both the fraction of cellular volume occupied by vesicles and the average size (diameter) of the vesicles was apparent as sampling moved from the area near the nucleus toward the plasma membrane. These findings describe the important role of the nanometer-scale organization of subcellular neighborhoods on insulin vesicle maturation.

Visualizing subcellular rearrangements in intact ß cells using soft x-ray tomography.

Characterizing relationships between cell structures and functions requires mesoscale mapping of intact cells showing subcellular rearrangements following stimulation; however, current approaches are limited in this regard. Here, we report a unique application of soft x-ray tomography to generate three-dimensional reconstructions of whole pancreatic ß cells at different time points following glucose-stimulated insulin secretion. Reconstructions following stimulation showed distinct insulin vesicle distribution patterns reflective of altered vesicle pool sizes as they travel through the secretory pathway. Our results show that glucose stimulation caused rapid changes in biochemical composition and/or density of insulin packing, increased mitochondrial volume, and closer proximity of insulin vesicles to mitochondria. Costimulation with exendin-4 (a glucagon-like peptide-1 receptor agonist) prolonged these effects and increased insulin packaging efficiency and vesicle maturation. This study provides unique perspectives on the coordinated structural reorganization and interactions of organelles that dictate cell responses.

De Novo Structural Pattern Mining in Cellular Electron Cryotomograms.

Electron cryotomography enables 3D visualization of cells in a near-native state at molecular resolution. The produced cellular tomograms contain detailed information about a plethora of macromolecular complexes, their structures, abundances, and specific spatial locations in the cell. However, extracting this information in a systematic way is very challenging, and current methods usually rely on individual templates of known structures. Here, we propose a framework called "Multi-Pattern Pursuit" for de novo discovery of different complexes from highly heterogeneous sets of particles extracted from entire cellular tomograms without using information of known structures. These initially detected structures can then serve as input for more targeted refinement efforts. Our tests on simulated and experimental tomograms show that our automated method is a promising tool for supporting large-scale template-free visual proteomics analysis.

Understanding the functional impact of copy number alterations in breast cancer using a network modeling approach.

UNLABELLED: Copy number alterations (CNAs) are thought to account for 85% of the variation in gene expression observed among breast tumours. The expression of cis-associated genes is impacted by CNAs occurring at proximal loci of these genes, whereas the expression of trans-associated genes is impacted by CNAs occurring at distal loci. While a majority of these CNA-driven genes responsible for breast tumourigenesis are cis-associated, trans-associated genes are thought to further abet the development of cancer and influence disease outcomes in patients. Here we present a network-based approach that integrates copy-number and expression profiles to identify putative cis- and trans-associated genes in breast cancer pathogenesis. We validate these cis- and trans-associated genes by employing them to subtype a large cohort of breast tumours obtained from the METABRIC consortium, and demonstrate that these genes accurately reconstruct the ten subtypes of breast cancer. We observe that individual breast cancer subtypes are driven by distinct sets of cis- and trans-associated genes. Among the cis-associated genes, we recover several known drivers of breast cancer (e.g. CCND1, ERRB2, MDM2 and ZNF703) and some novel putative drivers (e.g. BRF2 and SF3B3). siRNA-mediated knockdown of BRF2 across a panel of breast cancer cell lines showed significant reduction in cell viability for ER-/HER2+ (MDA-MB-453) cells, but not in normal (MCF10A) cells thereby indicating that BRF2 could be a viable therapeutic target for estrogen receptor-negative/HER2-enriched (ER-/HER2+) cancers. Among the trans-associated genes, we identify modules of immune response (CD2, CD19, CD38 and CD79B), mitotic/cell-cycle kinases (e.g. AURKB, MELK, PLK1 and TTK), and DNA-damage response genes (e.g. RFC4 and FEN1). siRNA-mediated knockdown of RFC4 significantly reduced cell proliferation in ER-negative normal breast and cancer lines, thereby indicating that RFC4 is essential for both normal and cancer cell survival but could be a useful biomarker for aggressive (ER-negative) breast tumours. AVAILABILITY: under NetStrat.

Inferring synthetic lethal interactions from mutual exclusivity of genetic events in cancer.

BACKGROUND: Synthetic lethality (SL) refers to the genetic interaction between two or more genes where only their co-alteration (e.g. by mutations, amplifications or deletions) results in cell death. In recent years, SL has emerged as an attractive therapeutic strategy against cancer: by targeting the SL partners of altered genes in cancer cells, these cells can be selectively killed while sparing the normal cells. Consequently, a number of studies have attempted prediction of SL interactions in human, a majority by extrapolating SL interactions inferred through large-scale screens in model organisms. However, these predicted SL interactions either do not hold in human cells or do not include genes that are (frequently) altered in human cancers, and are therefore not attractive in the context of cancer therapy. RESULTS: Here, we develop a computational approach to infer SL interactions directly from frequently altered genes in human cancers. It is based on the observation that pairs of genes that are altered in a (significantly) mutually exclusive manner in cancers are likely to constitute lethal combinations. Using genomic copy-number and gene-expression data from four cancers, breast, prostate, ovarian and uterine (total 3980 samples) from The Cancer Genome Atlas, we identify 718 genes that are frequently amplified or upregulated, and are likely to be synthetic lethal with six key DNA-damage response (DDR) genes in these cancers. By comparing with published data on gene essentiality (~16000 genes) from ten DDR-deficient cancer cell lines, we show that our identified genes are enriched among the top quartile of essential genes in these cell lines, implying that our inferred genes are highly likely to be (synthetic) lethal upon knockdown in these cell lines. Among the inferred targets are tousled-like kinase 2 (TLK2) and the deubiquitinating enzyme ubiquitin-specific-processing protease 7 (USP7) whose overexpression correlates with poor survival in cancers. CONCLUSION: Mutual exclusivity between frequently occurring genetic events identifies synthetic lethal combinations in cancers. These identified genes are essential in cell lines, and are potential candidates for targeted cancer therapy. Availability: http://bioinformatics.org.au/tools-data/underMutExSL