Challenges of Integrating Stochastic Dynamics and Cryo-Electron Tomograms in Whole-Cell Simulations

Harris S, Ramanathan A. Intelligent resolution: Integrating Cryo-EM with AI-driven multi-resolution simulations to observe the severe acute respiratory syndrome coronavirus-2 replication-transcription machinery in action

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) replication transcription complex (RTC) is a multi-domain protein responsible for replicating and transcribing the viral mRNA inside a human cell. Attacking RTC function with pharmaceutical compounds is a pathway to treating COVID-19. Conventional tools, e.g., cryo-electron microscopy and all-atom molecular dynamics (AAMD), do not provide sufficiently high resolution or timescale to capture important dynamics of this molecular machine. Consequently, we develop an innovative workflow that bridges the gap between these resolutions, using mesoscale fluctuating finite element analysis (FFEA) continuum simulations and a hierarchy of AI-methods that continually learn and infer features for maintaining consistency between AAMD and FFEA simulations. We leverage a multi-site distributed workflow manager to orchestrate AI, FFEA, and AAMD jobs, providing optimal resource utilization across HPC centers. Our study provides unprecedented access to study the SARS-CoV-2 RTC machinery, while providing general capability for AI-enabled multi-resolution simulations at scale.

Integrative modeling of the cell

A whole-cell model represents certain aspects of the cell structure and/or function. Due to the high complexity of the cell, an integrative modeling approach is often taken to utilize all available information including experimental data, prior knowledge and prior models. In this review, we summarize an emerging workflow of whole-cell modeling into five steps: (i) gather information; (ii) represent the modeled system into modules; (iii) translate input information into scoring function; (iv) sample the whole-cell model; (v) validate and interpret the model. In particular, we propose the integrative modeling of the cell by combining available (whole-cell) models to maximize the accuracy, precision, and completeness. In addition, we list quantitative predictions of various aspects of cell biology from existing whole-cell models. Moreover, we discuss the remaining challenges and future directions, and highlight the opportunity to establish an integrative spatiotemporal multi-scale whole-cell model based on a community approach.

The correlation between cell and nucleus size is explained by an eukaryotic cell growth model

In eukaryotes, the cell volume is observed to be strongly correlated with the nuclear volume. The slope of this correlation depends on the cell type, growth condition, and the physical environment of the cell. We develop a computational model of cell growth and proteome increase, incorporating the kinetics of amino acid import, protein/ribosome synthesis and degradation, and active transport of proteins between the cytoplasm and the nucleoplasm. We also include a simple model of ribosome biogenesis and assembly. Results show that the cell volume is tightly correlated with the nuclear volume, and the cytoplasm-nucleoplasm transport rates strongly influence the cell growth rate as well as the cell/nucleus volume ratio (C/N ratio). Ribosome assembly and the ratio of ribosomal proteins to mature ribosomes also influence the cell volume and the cell growth rate. We find that in order to regulate the cell growth rate and the cell/nucleus volume ratio, the cell must optimally control groups of kinetic and transport parameters together, which could explain the quantitative roles of canonical growth pathways. Finally, although not explicitly demonstrated in this work, we point out that it is possible to construct a detailed proteome distribution using our model and RNAseq data, provided that a quantitative cell division mechanism is known.

Generating Chromosome Geometries in a Minimal Cell From Cryo-Electron Tomograms and Chromosome Conformation Capture Maps

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.

[Cryo-electron microcopy for a new vision of the cell and its components]

Cryo-electron microscopy (cryo-EM) is a technique for imaging biological samples that plays a central role in structural biology, with high impact on research fields such as cell and developmental biology, bioinformatics, cell physics and applied mathematics. It allows the determination of structures of purified proteins within cells. This review describes the main recent advances in cryo-EM, illustrated by examples of proteins of biomedical interest, and the avenues for future development.

Visualizing Molecular Architectures of Cellular Condensates: Hints of Complex Coacervation Scenarios

In the last decade, liquid-liquid phase separation has emerged as a fundamental principle in the organization of crowded cellular environments into functionally distinct membraneless compartments. It is now established that biomolecules can condense into various physical phases, traditionally defined for simple polymer systems, and more recently elucidated by techniques employed in life sciences. We review pioneering cryo-electron tomography studies that have begun to unravel a wide spectrum of molecular architectures, ranging from amorphous to crystalline assemblies, that underlie cellular condensates. These observations bring into question current interpretations of microscopic phase behavior. Furthermore, by examining emerging concepts of non-classical phase separation pathways in small-molecule crystallization, we draw parallels with biomolecular condensation that highlight aspects not yet fully explored. In particular, transient and metastable intermediates that might be challenging to capture experimentally inside cells could be probed through computational simulations and enable a multi-scale understanding of the subcellular organization governed by distinct phases.

Art and Science of the Cellular Mesoscale

Experimental information from microscopy, structural biology, and bioinformatics may be integrated to build structural models of entire cells with molecular detail. This integrative modeling is challenging in several ways: the intrinsic complexity of biology results in models with many closely packed and heterogeneous components; the wealth of available experimental data is scattered among multiple resources and must be gathered, reconciled, and curated; and computational infrastructure is only now gaining the capability of modeling and visualizing systems of this complexity. We present recent efforts to address these challenges, both with artistic approaches to depicting the cellular mesoscale, and development and application of methods to build quantitative models.

Securing the future of research computing in the biosciences

Improvements in technology often drive scientific discovery. Therefore, research requires sustained investment in the latest equipment and training for the researchers who are going to use it. Prioritising and administering infrastructure investment is challenging because future needs are difficult to predict. In the past, highly computationally demanding research was associated primarily with particle physics and astronomy experiments. However, as biology becomes more quantitative and bioscientists generate more and more data, their computational requirements may ultimately exceed those of physical scientists. Computation has always been central to bioinformatics, but now imaging experiments have rapidly growing data processing and storage requirements. There is also an urgent need for new modelling and simulation tools to provide insight and understanding of these biophysical experiments. Bioscience communities must work together to provide the software and skills training needed in their areas. Research-active institutions need to recognise that computation is now vital in many more areas of discovery and create an environment where it can be embraced. The public must also become aware of both the power and limitations of computing, particularly with respect to their health and personal data.

Simulating biological processes: stochastic physics from whole cells to colonies

The last few decades have revealed the living cell to be a crowded spatially heterogeneous space teeming with biomolecules whose concentrations and activities are governed by intrinsically random forces. It is from this randomness, however, that a vast array of precisely timed and intricately coordinated biological functions emerge that give rise to the complex forms and behaviors we see in the biosphere around us. This seemingly paradoxical nature of life has drawn the interest of an increasing number of physicists, and recent years have seen stochastic modeling grow into a major subdiscipline within biological physics. Here we review some of the major advances that have shaped our understanding of stochasticity in biology. We begin with some historical context, outlining a string of important experimental results that motivated the development of stochastic modeling. We then embark upon a fairly rigorous treatment of the simulation methods that are currently available for the treatment of stochastic biological models, with an eye toward comparing and contrasting their realms of applicability, and the care that must be taken when parameterizing them. Following that, we describe how stochasticity impacts several key biological functions, including transcription, translation, ribosome biogenesis, chromosome replication, and metabolism, before considering how the functions may be coupled into a comprehensive model of a 'minimal cell'. Finally, we close with our expectation for the future of the field, focusing on how mesoscopic stochastic methods may be augmented with atomic-scale molecular modeling approaches in order to understand life across a range of length and time scales.

Multiscale Methods in Drug Design Bridge Chemical and Biological Complexity in the Search for Cures

Drug action is inherently multiscale: it connects molecular interactions to emergent properties at cellular and larger scales. Simulation techniques at each of these different scales are already central to drug design and development, but methods capable of connecting across these scales will extend understanding of complex mechanisms and the ability to predict biological effects. Improved algorithms, ever-more-powerful computing architectures and the accelerating growth of rich datasets are driving advances in multiscale modeling methods capable of bridging chemical and biological complexity from the atom to the cell. Particularly exciting is the development of highly detailed, structure-based, physical simulations of biochemical systems, which are now able to access experimentally relevant timescales for large systems and, at the same time, achieve unprecedented accuracy. In this Perspective, we discuss how emerging data-rich, physics-based multiscale approaches are of the cusp of realizing long-promised impact in the discovery, design and development of novel therapeutics. We highlight emerging methods and applications in this growing field, and outline how different scales can be combined in practical modelling and simulation strategies.

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.

Molecular simulations of cellular processes

It is, nowadays, possible to simulate biological processes in conditions that mimic the different cellular compartments. Several groups have performed these calculations using molecular models that vary in performance and accuracy. In many cases, the atomistic degrees of freedom have been eliminated, sacrificing both structural complexity and chemical specificity to be able to explore slow processes. In this review, we will discuss the insights gained from computer simulations on macromolecule diffusion, nuclear body formation, and processes involving the genetic material inside cell-mimicking spaces. We will also discuss the challenges to generate new models suitable for the simulations of biological processes on a cell scale and for cell-cycle-long times, including non-equilibrium events such as the co-translational folding, misfolding, and aggregation of proteins. A prominent role will be played by the wise choice of the structural simplifications and, simultaneously, of a relatively complex energetic description. These challenging tasks will rely on the integration of experimental and computational methods, achieved through the application of efficient algorithms. Graphical abstract.