We present the first-ever, fully discrete, stochastic model of triggered cardiac Ca2+ dynamics. Using anatomically accurate subcellular cardiac myocyte geometries, we simulate the molecular players involved in Ca2+ handling using high-resolution stochastic and explicit-particle methods at the level of an individual cardiac dyadic junction. Integrating data from multiple experimental sources, the model not only replicates the findings of traditional in silico studies and complements in vitro experimental data but also reveals new insights into the molecular mechanisms driving cardiac dysfunction under stress and disease conditions. We improve upon older, nondiscrete models using the same realistic geometry by incorporating molecular mechanisms for spontaneous, as well as triggered calcium-induced calcium release (CICR). Action potentials are used to activate L-type calcium channels (LTCC), triggering CICR through ryanodine receptors (RyRs) on the surface of the sarcoplasmic reticulum. These improvements allow for the specific focus on the couplon: the structure-function relationship between LTCC and RyR. We investigate the electrophysical effects of normal and diseased action potentials on CICR and interrogate the effects of dyadic junction deformation through detubulation and orphaning of RyR. Our work demonstrates the importance of the electrophysical integrity of the calcium release unit on CICR fidelity, giving insights into the molecular basis of heart disease. Finally, we provide a unique, detailed, molecular view of the CICR process using advanced rendering techniques. This easy-to-use model comes complete with tutorials and the necessary software for use and analysis to maximize usability and reproducibility. Our work focuses on quantifying, qualifying, and visualizing the behavior of the molecular species that underlie the function and dysfunction of subcellular cardiomyocyte systems.

A lattice-based method is presented for creating 3D models of entire bacterial nucleoids integrating ultrastructural information cryoelectron tomography, genomic and proteomic data, and experimental atomic structures of biomolecules and assemblies. The method is used to generate models of the minimal genome bacterium JCVI-Syn3A, producing a series of models that test hypotheses about transcription, condensation, and overall distribution of the genome. Lattice models are also used to generate atomic models of an entire JCVI-Syn3A cell.

The advent of cryo-electron microscopy (cryo-EM) and cryo-electron tomography (cryo-ET), coupled with computational modeling, has enabled the creation of integrative 3D models of viruses, bacteria, and cellular organelles. These models, composed of thousands of macromolecules and billions of atoms, have historically posed significant challenges for manipulation and visualization without specialized molecular graphics tools and hardware. With the recent advancements in GPU rendering power and web browser capabilities, it is now feasible to render interactively large molecular scenes directly on the web. In this work, we introduce Mesoscale Explorer, a web application built using the Mol* framework, dedicated to the visualization of large-scale molecular models ranging from viruses to cell organelles. Mesoscale Explorer provides unprecedented access and insight into the molecular fabric of life, enhancing perception, streamlining exploration, and simplifying visualization of diverse data types, showcasing the intricate details of these models with unparalleled clarity.

In January 2024, a targeted conference, ‘CellVis2’, was held at Scripps Research in La Jolla, USA, the second in a series designed to explore the promise, practices, roadblocks, and prospects of creating, visualizing, sharing, and communicating physical representations of entire biological cells at scales down to the atom.

Understanding structures and interactions of biomolecules benefits from the use of interactive 3D graphics, which are greatly facilitated by stereoscopy and intuitive user interfaces. The recent coronavirus disease 2019 (COVID-19) pandemic has brought remote learning to the fore and has exposed its limitations, especially regarding 3D structure-based topics. The ever-expanding computer gaming industry has brought virtual reality (VR) and remote multiplayer technology into the consumer realm, creating the currently most affordable stereoscopic viewing platform with immersive presence and natural 3D interactivity.

Models of insulin secretory vesicles from pancreatic beta cells have been created using the cellPACK suite of tools to research, curate, construct and visualise the current state of knowledge. The model integrates experimental information from proteomics, structural biology, cryoelectron microscopy and X-ray tomography, and is used to generate models of mature and immature vesicles. A new method was developed to generate a confidence score that reconciles inconsistencies between three available proteomes using expert annotations of cellular localisation. The models are used to simulate soft X-ray tomograms, allowing quantification of features that are observed in experimental tomograms, and in turn, allowing interpretation of X-ray tomograms at the molecular level.

Building structural models of entire cells has been a long-standing cross-discipline challenge for the research community, as it requires an unprecedented level of integration between multiple sources of biological data and enhanced methods for computational modeling and visualization. Here, we present the first 3D structural models of an entire Mycoplasma genitalium (MG) cell, built using the CellPACK suite of computational modeling tools. Our model recapitulates the data described in recent whole-cell system biology simulations and provides a structural representation for all MG proteins, DNA and RNA molecules, obtained by combining experimental and homology-modeled structures and lattice-based models of the genome. We establish a framework for gathering, curating and evaluating these structures, exposing current weaknesses of modeling methods and the boundaries of MG structural knowledge, and visualization methods to explore functional characteristics of the genome and proteome. We compare two approaches for data gathering, a manually-curated workflow and an automated workflow that uses homologous structures, both of which are appropriate for the analysis of mesoscale properties such as crowding and volume occupancy. Analysis of model quality provides estimates of the regularization that will be required when these models are used as starting points for atomic molecular dynamics simulations.

We present a new technique for the rapid modeling and construction of scientifically accurate mesoscale biological models. The resulting 3D models are based on a few 2D microscopy scans and the latest knowledge available about the biological entity, represented as a set of geometric relationships. Our new visual-programming technique is based on statistical and rule-based modeling approaches that are rapid to author, fast to construct, and easy to revise. From a few 2D microscopy scans, we determine the statistical properties of various structural aspects, such as the outer membrane shape, the spatial properties, and the distribution characteristics of the macromolecular elements on the membrane. This information is utilized in the construction of the 3D model. Once all the imaging evidence is incorporated into the model, additional information can be incorporated by interactively defining the rules that spatially characterize the rest of the biological entity, such as mutual interactions among macromolecules, and their distances and orientations relative to other structures. These rules are defined through an intuitive 3D interactive visualization as a visual-programming feedback loop. We demonstrate the applicability of our approach on a use case of the modeling procedure of the SARS-CoV-2 virion ultrastructure. This atomistic model, which we present here, can steer biological research to new promising directions in our efforts to fight the spread of the virus.

Coarse-grained models have long been considered indispensable tools in the investigation of biomolecular dynamics and assembly. However, the process of simulating such models is arduous because unconventional force fields and particle attributes are often needed, and some systems are not in thermal equilibrium. Although modern molecular dynamics programs are highly adaptable, software designed for preparing all-atom simulations typically makes restrictive assumptions about the nature of the particles and the forces acting on them. Consequently, the use of coarse-grained models has remained challenging. Moltemplate is a file format for storing coarse-grained molecular models and the forces that act on them, as well as a program that converts moltemplate files into input files for LAMMPS, a popular molecular dynamics engine. Moltemplate has broad scope and an emphasis on generality. It accommodates new kinds of forces as they are developed for LAMMPS, making moltemplate a popular tool with thousands of users in computational chemistry, materials science, and structural biology. To demonstrate its wide functionality, we provide examples of using moltemplate to prepare simulations of fluids using many-body forces, coarse-grained organic semiconductors, and the motor-driven supercoiling and condensation of an entire bacterial chromosome.

CellPAINT is an interactive digital tool that allows non-expert users to create illustrations of the molecular structure of cells and viruses. We present a new release with several key enhancements, including the ability to generate custom ingredients from structure information in the Protein Data Bank, and interaction, grouping, and locking functions that streamline the creation of assemblies and illustration of large, complex scenes. An example of CellPAINT as a tool for hypothesis generation in the interpretation of cryoelectron tomograms is presented. CellPAINT is freely available at http://ccsb.scripps.edu/cellpaint.

Interest is growing for 3D models of the biological mesoscale, the intermediate scale between the nanometer scale of molecular structure and micrometer scale of cellular biology. However, it is currently difficult to gather, curate and integrate all the data required to define such models. To address this challenge we developed Mesoscope (mesoscope.scripps.edu/beta), a web-based data integration and curation tool. Mesoscope allows users to begin with a listing of molecules (such as data from proteomics), and to use resources at UniProt and the PDB to identify, prepare and validate appropriate structures and representations for each molecule, ultimately producing a portable output file used by CellPACK and other modeling tools for generation of 3D models of the biological mesoscale. The availability of this tool has proven essential in several exploratory applications, given the high complexity of mesoscale models and the heterogeneity of the available data sources.

In this paper, we describe cellPACKexplorer, a system designed to help developers of cellPACK find errors in and improve their algorithm. cellPACKexplorer focuses on visualizing the effects of cellPACK recipe parameters on the final packing output. We found that the developers have two different methods for understanding the output, numerical and visual, depending on their background. We designed cellPACKexplorer with a flexible interface to support both types of users. We evaluated our tool through case studies and questionnaires. Novice users were able to create cell models with cellPACK and explore the behavior of different parameters. Further, expert users discovered an error in the code and were able to locate the problem quickly with our new analysis tool. We conclude with a discussion of the implications of our findings in the wider visualization community.

Visualization of biological mesoscale models provides a glimpse at the inner workings of living cells. One of the most complex components of these models is DNA, which is of fundamental importance for all forms of life. Modeling the 3D structure of genomes has previously only been attempted by sequential approaches. We present the first parallel approach for the instant construction of DNA structures. Traditionally, such structures are generated with algorithms like random walk, which have inherent sequential constraints. These algorithms result in the desired structure, are easy to control, and simple to formulate. Their execution, however, is very time‐consuming, as they are not designed to exploit parallelism. We propose an approach to parallelize the process, facilitating an implementation on the GPU.

Visualizations of hierarchical data can often be explored interactively. For example, in geographic visualization, there are continents, which can be subdivided into countries, states, counties and cities. Similarly, in models of viruses or bacteria at the highest level are the compartments, and below that are macromolecules, secondary structures (such as α-helices), amino-acids, and on the finest level atoms. Distinguishing between items can be assisted through the use of color at all levels. However, currently, there are no hierarchical and adaptive color mapping techniques for very large multi-scale visualizations that can be explored interactively. We present a novel, multi-scale, color-mapping technique for adaptively adjusting the color scheme to the current view and scale. Color is treated as a resource and is smoothly redistributed. The distribution adjusts to the scale of the currently observed detail and maximizes the color range utilization given current viewing requirements. Thus, we ensure that the user is able to distinguish items on any level, even if the color is not constant for a particular feature. The coloring technique is demonstrated for a political map and a mesoscale structural model of HIV. The technique has been tested by users with expertise in structural biology and was overall well received.

The small program Illustrate generates non-photorealistic images of biological molecules for use in dissemination, outreach, and education. The method has been used as part of the "Molecule of the Month,'' an ongoing educational column at the RCSB Protein Data Bank (http://rcsb.org). Insights from 20 years of application of the program are presented, and the program has been released both as open-source Fortran at GitHub and through an interactive web-based interface.

Information from proteomics, microscopy, and structural biology are integrated to create structural models of exosomes, small vesicles released from cells. Three visualization methods are employed and compared: 2D painting of a cross section using traditional media, manual creation of a cross section using the mesoscale 2.5D digital painting software cellPAINT, and generation of a 3D atomic model using the mesoscale modeling program cellPACK.

Mesoscale molecular modeling is providing a new window into the inner workings of living cells. Modeling of genomes, however, remains a technical challenge, due to their large size and complexity. We describe a lattice method for rapid generation of bacterial nucleoid models that integrates experimental data from a variety of biophysical techniques and provides a starting point for simulation and hypothesis generation. The current method builds models of a circular bacterial genome with supercoiled plectonemes, packed within the small space of the bacterial cell. Lattice models are generated for Mycoplasma genitalium and Escherichia coli nucleoids, and used to simulate interaction data. The method is rapid enough to allow generation of multiple models when analyzing structure/function relationships, and we demonstrate use of the lattice models in creation of an all-atom representation of an entire cell.

Computational studies on viruses can reveal a whole host of new information, including protein cooperativity, conformational changes and numerical binding data. Unfortunately, computational studies on viruses have been lacking in the literature due to their size and realistic time constraints. Thus, there are several holes in the current state-of-the-art to be filled. Currently, it is thought that influenza A has multiple conformations, including one where hemagglutinin (HA) potentially interacts with neuraminidase (NA) to facilitate binding of sialic acid (SA). However, this cooperativity is not clear as cytoskeletal proteins may fix the structure and influence the dynamics of HA and NA. One additional inconsistency is that SA preferentially binds to NA, even though binding to HA is needed for endocytosis. Decoding the reasons for this preference will require long molecular dynamics (MD) simulations; one way to realize this is through software engineering. Investigating viral characteristics as the defining science goal drives the development of cellPACK, a software program able to incorporate experimental and theoretical data to probabilistically pack proteins into a defined space. By creating cellPACK functionalities and linking them to simulation programs, accurate protein-protein interactions can be modeled. The first step in this process is building a model; here we built the world's first full viral model, of the influenza A, H1N1 Cali09 strain, which is needed to better understand intrinsic characteristics of viruses and how cross-species mutations can influence virulence. This model can then be used in MD and Brownian dynamics (BD) simulations to tease out how protein cooperativity affects influenza virulence.

CellPAINT allows nonexpert users to create interactive mesoscale illustrations that integrate a variety of biological data. Like popular digital painting software, scenes are created using a palette of molecular “brushes.” The current release allows creation of animated scenes with an HIV virion, blood plasma, and a simplified T-cell.

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers. L’archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d’enseignement et de recherche français ou étrangers, des laboratoires publics ou privés. In Silico Studies on DARC. Alexandre De Brevern, Ludovic Autin, Yves Colin, Olivier Bertrand, Catherine Etchebest