Search for: All records

ABSTRACT Many microscopic images and simulations of cells give results in different kinds of formats, making it difficult for people lacking computational skills to visualize and interact with them. Minecraft—known for its three-dimensional, open-world, voxel-based environment—offers a unique solution by allowing the direct insertion of voxel-based cellular structures from light microscopy and simulations into its worlds without modification. This integration enables Minecraft players to explore the ultrastructure of cells in a highly immersive and interactive environment. Here, we demonstrate several workflows that can convert images and simulation results into Minecraft worlds. Using the workflows, students can easily import and interact with a variety of cellular content, including bacteria, yeast, and cancer cells. This approach not only opens new avenues for science education but also demonstrates the potential of combining scientific visualization with interactive gaming platforms for facilitating research and improving appreciation of cellular structure for a broad audience. 

Stochastic nature of gene expression leads to the complex formation of the bacterial transcriptome and proteome. In contrast to typical transcriptome studies, we employ a near wild-type, Syn1.0, of the naturally genome-reducedMycoplasmas, and the dramatically further genome-reduced JCVI-syn3A thus avoiding additional contributions from many non-essential cellular functions. To aid in profiling the transcriptional landscape within these bacteria, we present a bioinformatic analysis of the genetic sequence motifs implicated in modulating the stochastic gene expression events, coupled with genome-wide short-read (Illumina) and long-read (Oxford Nanopore Technologies and Pacfic Biosciences) RNA sequencing. The bioinformatic analysis coupled with information from structural studies assigns strengths of the Shine-Dalgarno signatures and identifies both transcription initiation and termination sites, leading to predictions of RNA isoforms in Syn1.0 (and related organisms). The long-read and short-read RNA sequencing characterized the predicted transcriptional activity, and the long-read methods provide direct insight into the RNA isoform complexity within Syn1.0. Comparison of the RNA sequencing results with that of the bioinformatic analysis highlights the inability of bioinformatics alone to capture the results of bacterial transcription without including effects of RNA degradation. This study emphasizes the need for comparative analysis and potential dangers of genome reduction, exemplified through the discovery of altered gene expression patterns of JCVI-syn1.0 and JCVI-syn3A, achieved via the union of our transcriptome study with their proteomics data. Analysis of the transcriptomics data sets through a Jupyter notebook allows any genomic region to be easily examined. Table of Content Image 

Computational models of cells cannot be considered complete unless they include the most fundamental process of life, the replication and inheritance of genetic material. By creating a computational framework to model systems of replicating bacterial chromosomes as polymers at 10 bp resolution with Brownian dynamics, we investigate changes in chromosome organization during replication and extend the applicability of an existing whole-cell model (WCM) for a genetically minimal bacterium, JCVI-syn3A, to the entire cell-cycle. To achieve cell-scale chromosome structures that are realistic, we model the chromosome as a self-avoiding homopolymer with bending and torsional stiffnesses that capture the essential mechanical properties of dsDNA in Syn3A. In addition, the conformations of the circular DNA must avoid overlapping with ribosomes identitied in cryo-electron tomograms. While Syn3A lacks the complex regulatory systems known to orchestrate chromosome segregation in other bacteria, its minimized genome retains essential loop-extruding structural maintenance of chromosomes (SMC) protein complexes (SMC-scpAB) and topoisomerases. Through implementing the effects of these proteins in our simulations of replicating chromosomes, we find that they alone are sufficient for simultaneous chromosome segregation across all generations within nested theta structures. This supports previous studies suggesting loop-extrusion serves as a near-universal mechanism for chromosome organization within bacterial and eukaryotic cells. Furthermore, we analyze ribosome diffusion under the influence of the chromosome and calculatein silicochromosome contact maps that capture inter-daughter interactions. Finally, we present a methodology to map the polymer model of the chromosome to a Martini coarse-grained representation to prepare molecular dynamics models of entire Syn3A cells, which serves as an ultimate means of validation for cell states predicted by the WCM. 

The ultimate microscope, directed at a cell, would reveal the dynamics of all the cell’s components with atomic resolution. In contrast to their real-world counterparts, computational microscopes are currently on the brink of meeting this challenge. In this perspective, we show how an integrative approach can be employed to model an entire cell, the minimal cell, JCVI-syn3A, at full complexity. This step opens the way to interrogate the cell’s spatio-temporal evolution with molecular dynamics simulations, an approach that can be extended to other cell types in the near future. 

Abstract Eukaryotic cells are spatially organized into functionally-distinct compartments. This three-dimensional(3D) organization generates intracellular heterogeneities that can modulate regulatory dynamics. Despite this knowledge of subcellular organization, most quantitative gene-regulation models still assume a well-mixed environment in which molecules can react regardless of their spatial positions. Here, we use the well-established galactose switch in budding yeast (Saccharomyces cerevisiae) to develop spatially-resolved models that integrate experimentally-derived intracellular architectures, including chromosome organization, the endoplasmic reticulum (ER) and spatially distinct ribosome populations. We implement a hybrid stochastic–deterministic framework in which gene expression is modeled using a reaction–diffusion master equation that enforces locality (i.e., reactions occur only when molecules are in physical proximity), while metabolic and transport processes are captured by ordinary differential equations. Guided by electron microscopy and biochemical constraints, we quantify how accounting for intracellular spatial organization alters regulatory predictions in the galactose switch. We show that chromosome geometry has little effect on Gal2p output, whereas ER-associated translation reduces Gal2p delivery to the plasma membrane; the largest decrease of Gal2p abundance occurs when translation ofGAL2mRNA is restricted to a population of ribosomes physically bound to the ER. Together, these results demonstrate that more realistic 3D cellular architectures and local reaction rules can qualitatively change regulatory predictions, motivating integration of intracellular organization in future whole-cell models. Author SummaryEukaryotic cells are highly organized spaces with distinct subcellular compartments and heterogeneous distributions of molecules which influence how cells function. Yet, most computational models of gene regulation assume a spatially homogeneous intracellular environment. To quantify how intracellular architecture influences regulatory predictions, we developed spatially resolved models of the galactose switch in budding yeast, a well-characterized gene regulatory system controlling the response to extracellular galactose. We compared conventional well-stirred simulations with models that explicitly incorporate experimentally informed intracellular organization, including chromosome positioning, endoplasmic reticulum geometry, and functionally distinct ribosome populations. Incorporating these spatial features substantially altered the dynamics of predicted gene activity, protein production, and intracellular sugar levels. Our results demonstrate that the three-dimensional organization of eukaryotic cells can significantly change regulatory outcomes of computational models, underscoring the need for integrating realistic spatial architectures in future models of eukaryotic gene regulation. 

JCVI-syn3A is a genetically minimal bacterial cell, consisting of 493 genes and only a single 543 kbp circular chromosome. Syn3A’s genome and physical size are approximately one-tenth those of the model bacterial organism Escherichia coli ’s, and the corresponding reduction in complexity and scale provides a unique opportunity for whole-cell modeling. Previous work established genome-scale gene essentiality and proteomics data along with its essential metabolic network and a kinetic model of genetic information processing. In addition to that information, whole-cell, spatially-resolved kinetic models require cellular architecture, including spatial distributions of ribosomes and the circular chromosome’s configuration. We reconstruct cellular architectures of Syn3A cells at the single-cell level directly from cryo-electron tomograms, including the ribosome distributions. We present a method of generating self-avoiding circular chromosome configurations in a lattice model with a resolution of 11.8 bp per monomer on a 4 nm cubic lattice. Realizations of the chromosome configurations are constrained by the ribosomes and geometry reconstructed from the tomograms and include DNA loops suggested by experimental chromosome conformation capture (3C) maps. Using ensembles of simulated chromosome configurations we predict chromosome contact maps for Syn3A cells at resolutions of 250 bp and greater and compare them to the experimental maps. Additionally, the spatial distributions of ribosomes and the DNA-crowding resulting from the individual chromosome configurations can be used to identify macromolecular structures formed from ribosomes and DNA, such as polysomes and expressomes.