A model‐driven approach towards rational microbial bioprocess optimization

Due to sustainability concerns, bio‐based production capitalizing on microbes as cell factories is in demand to synthesize valuable products. Nevertheless, the nonhomogenous variations of the extracellular environment in bioprocesses often challenge the biomass growth and the bioproduction yield. To enable a more rational bioprocess optimization, we have established a model‐driven approach that systematically integrates experiments with modeling, executed from flask to bioreactor scale, and using ferulic acid to vanillin bioconversion as a case study. The impacts of mass transfer and aeration on the biomass growth and bioproduction performances were examined using minimal small‐scale experiments. An integrated model coupling the cell factory kinetics with the three‐dimensional computational hydrodynamics of bioreactor was developed to better capture the spatiotemporal distributions of bioproduction. Full‐factorial predictions were then performed to identify the desired operating conditions. A bioconversion yield of 94% was achieved, which is one of the highest for recombinant Escherichia coli using ferulic acid as the precursor.     

Origins of regulated cell-to-cell variability

Single-cell measurements and lineage-tracing experiments are revealing that phenotypic cell-to-cell variability is often the result of deterministic processes, despite the existence of intrinsic noise in molecular networks. In most cases, this determinism represents largely uncharacterized molecular regulatory mechanisms, which places the study of cell-to-cell variability in the realm of molecular cell biology. Further research in the field will be important to advance quantitative cell biology because it will provide new insights into the mechanisms by which cells coordinate their intracellular activities in the spatiotemporal context of the multicellular environment.     

Modeling spatiotemporal abundance of mobile wildlife in highly variable environments using boosted GAMLSS hurdle models

Modeling organism distributions from survey data involves numerous statistical challenges, including accounting for zero‐inflation, overdispersion, and selection and incorporation of environmental covariates. In environments with high spatial and temporal variability, addressing these challenges often requires numerous assumptions regarding organism distributions and their relationships to biophysical features. These assumptions may limit the resolution or accuracy of predictions resulting from survey‐based distribution models. We propose an iterative modeling approach that incorporates a negative binomial hurdle, followed by modeling of the relationship of organism distribution and abundance to environmental covariates using generalized additive models (GAM) and generalized additive models for location, scale, and shape (GAMLSS). Our approach accounts for key features of survey data by separating binary (presence‐absence) from count (abundance) data, separately modeling the mean and dispersion of count data, and incorporating selection of appropriate covariates and response functions from a suite of potential covariates while avoiding overfitting. We apply our modeling approach to surveys of sea duck abundance and distribution in Nantucket Sound (Massachusetts, USA), which has been proposed as a location for offshore wind energy development. Our model results highlight the importance of spatiotemporal variation in this system, as well as identifying key habitat features including distance to shore, sediment grain size, and seafloor topographic variation. Our work provides a powerful, flexible, and highly repeatable modeling framework with minimal assumptions that can be broadly applied to the modeling of survey data with high spatiotemporal variability. Applying GAMLSS models to the count portion of survey data allows us to incorporate potential overdispersion, which can dramatically affect model results in highly dynamic systems. Our approach is particularly relevant to systems in which little a priori knowledge is available regarding relationships between organism distributions and biophysical features, since it incorporates simultaneous selection of covariates and their functional relationships with organism responses.     

Coregionalized single- and multiresolution spatially varying growth curve modeling with application to weed growth

Modeling of longitudinal data from agricultural experiments using growth curves helps understand conditions conducive or unconducive to crop growth. Recent advances in Geographical Information Systems (GIS) now allow geocoding of agricultural data that help understand spatial patterns. A particularly common problem is capturing spatial variation in growth patterns over the entire experimental domain. Statistical modeling in these settings can be challenging because agricultural designs are often spatially replicated, with arrays of subplots, and interest lies in capturing spatial variation at possibly different resolutions. In this article, we develop a framework for modeling spatially varying growth curves as Gaussian processes that capture associations at single and multiple resolutions. We provide Bayesian hierarchical models for this setting, where flexible parameterization enables spatial estimation and prediction of growth curves. We illustrate using data from weed growth experiments conducted in Waseca, Minnesota, that recorded growth of the weed Setaria spp. in a spatially replicated design.     

In channelrhodopsin-2 Glu-90 is crucial for ion selectivity and is deprotonated during the photocycle

The light-activated microbial ion channel channelrhodopsin-2 (ChR2) is a powerful tool to study cellular processes with high spatiotemporal resolution in the emerging field of optogenetics. To customize the channel properties for optogenetic experiments, a detailed understanding of its molecular reaction mechanism is essential. Here, Glu-90, a key residue involved in the gating and selectivity mechanism of the ion channel is characterized in detail. The deprotonation of Glu-90 during the photocycle is elucidated by time-resolved FTIR spectroscopy, which seems to be part of the opening mechanism of the conductive pore. Furthermore, Glu-90 is crucial to ion selectivity as also revealed by mutation of this residue combined with voltage clamp experiments. By dynamic homology modeling, we further hypothesized that the conductive pore is flanked by Glu-90 and located between helices A, B, C, and G.     

Integration of software tools for integrative modeling of biomolecular systems

Integrative modeling computes a model based on varied types of input information, be it from experiments or prior models. Often, a type of input information will be best handled by a specific modeling software package. In such a case, we desire to integrate our integrative modeling software package, Integrative Modeling Platform (IMP), with software specialized to the computational demands of the modeling problem at hand. After several attempts, however, we have concluded that even in collaboration with the software’s developers, integration is either impractical or impossible. The reasons for the intractability of integration include software incompatibilities, differing modeling logic, the costs of collaboration, and academic incentives. In the integrative modeling software ecosystem, several large modeling packages exist with often redundant tools. We reason, therefore, that the other development groups have similarly concluded that the benefit of integration does not justify the cost. As a result, modelers are often restricted to the set of tools within a single software package. The inability to integrate tools from distinct software negatively impacts the quality of the models and the efficiency of the modeling. As the complexity of modeling problems grows, we seek to galvanize developers and modelers to consider the long-term benefit that software interoperability yields. In this article, we formulate a demonstrative set of software standards for implementing a model search using tools from independent software packages and discuss our efforts to integrate IMP and the crystallography suite Phenix within the Bayesian modeling framework.     

Hemodynamic traveling waves in human visual cortex

Functional MRI (fMRI) experiments rely on precise characterization of the blood oxygen level dependent (BOLD) signal. As the spatial resolution of fMRI reaches the sub-millimeter range, the need for quantitative modelling of spatiotemporal properties of this hemodynamic signal has become pressing. Here, we find that a detailed physiologically-based model of spatiotemporal BOLD responses predicts traveling waves with velocities and spatial ranges in empirically observable ranges. Two measurable parameters, related to physiology, characterize these waves: wave velocity and damping rate. To test these predictions, high-resolution fMRI data are acquired from subjects viewing discrete visual stimuli. Predictions and experiment show strong agreement, in particular confirming BOLD waves propagating for at least 5-10 mm across the cortical surface at speeds of 2-12 mm s-1. These observations enable fundamentally new approaches to fMRI analysis, crucial for fMRI data acquired at high spatial resolution.     

Odontocete spatial patterns and temporal drivers of detection at sites in the Hawaiian islands

Successful conservation and management of marine top predators rely on detailed documentation of spatiotemporal behavior. For cetacean species, this information is key to defining stocks, habitat use, and mitigating harmful interactions. Research focused on this goal is employing methodologies such as visual observations, tag data, and passive acoustic monitoring (PAM) data. However, many studies are temporally limited or focus on only one or few species. In this study, we make use of an existing long-term (2009–2019), labeled PAM data set to examine spatiotemporal patterning of at least 10 odontocete (toothed whale) species in the Hawaiian Islands using compositional analyses and modeling techniques. Species composition differs among considered sites, and this difference is robust to seasonal movement patterns. Temporally, hour of day was the most significant predictor of detection across species and sites, followed by season, though patterns differed among species. We describe long-term trends in species detection at one site and note that they are markedly similar for many species. These trends may be related to long-term, underlying oceanographic cycles that will be the focus of future study. We demonstrate the variability of temporal patterns even at relatively close sites, which may imply that wide-ranging models of species presence are missing key fine-scale movement patterns. Documented seasonal differences in detection also highlights the importance of considering season in survey design both regionally and elsewhere. We emphasize the utility of long-term, continuous monitoring in highlighting temporal patterns that may relate to underlying climatic states and help us predict responses to climate change. We conclude that long-term PAM records are a valuable resource for documenting spatiotemporal patterns and can contribute many insights into the lives of top predators, even in highly studied regions such as the Hawaiian Islands.     

Wavelet bicoherence analysis as a method for investigating coherent structures in an electron beam with an overcritical current

Results are presented from numerical modeling of the effect of the inhomogeneity of the ion background on the complicated spatiotemporal dynamics of an electron beam with a virtual cathode in plane geometry. The possibility is demonstrated of increasing the generation frequency without changing the beam current. The spatiotemporal structures that form in the beam and govern the complicated stochastic dynamics of the nonuniform electron-plasma system under consideration are investigated by the methods of wavelet bicoherence and by analyzing the calculated electron trajectories on the space-time diagrams.     

Oscillatory NAD(P)H Waves and Calcium Oscillations in Neutrophils? A Modeling Study of Feasibility

The group of Howard Petty has claimed exotic metabolic wave phenomena together with mutually phase-coupled NAD(P)H- and calcium-oscillations in human neutrophils. At least parts of these phenomena are highly doubtful due to extensive failure of reproducibility by several other groups and hints that unreliable data from the Petty lab are involved in publications concerning circular calcium waves. The aim of our theoretical spatiotemporal modeling approach is to propose a possible and plausible biochemical mechanism which would, in principle, be able to explain metabolic oscillations and wave phenomena in neutrophils. Our modeling suggests the possibility of a calcium-controlled glucose influx as a driving force of metabolic oscillations and a potential role of polarized cell geometry and differential enzyme distribution for various NAD(P)H wave phenomena. The modeling results are supposed to stimulate further controversial discussions of such phenomena and potential mechanisms and experimental efforts to finally clarify the existence and biochemical basis of any kind of temporal and spatiotemporal patterns of calcium signals and metabolic dynamics in human neutrophils. Independent of Petty's observations, they present a general feasibility study of such phenomena in cells.     

Basic science through engineering? Synthetic modeling and the idea of biology-inspired engineering

Synthetic biology is often understood in terms of the pursuit for well-characterized biological parts to create synthetic wholes. Accordingly, it has typically been conceived of as an engineering dominated and application oriented field. We argue that the relationship of synthetic biology to engineering is far more nuanced than that and involves a sophisticated epistemic dimension, as shown by the recent practice of synthetic modeling. Synthetic models are engineered genetic networks that are implanted in a natural cell environment. Their construction is typically combined with experiments on model organisms as well as mathematical modeling and simulation. What is especially interesting about this combinational modeling practice is that, apart from greater integration between these different epistemic activities, it has also led to the questioning of some central assumptions and notions on which synthetic biology is based. As a result synthetic biology is in the process of becoming more “biology inspired.”     

Bridging experiments, models and simulations: an integrative approach to validation in computational cardiac electrophysiology

Computational models in physiology often integrate functional and structural information from a large range of spatiotemporal scales from the ionic to the whole organ level. Their sophistication raises both expectations and skepticism concerning how computational methods can improve our understanding of living organisms and also how they can reduce, replace, and refine animal experiments. A fundamental requirement to fulfill these expectations and achieve the full potential of computational physiology is a clear understanding of what models represent and how they can be validated. The present study aims at informing strategies for validation by elucidating the complex interrelations among experiments, models, and simulations in cardiac electrophysiology. We describe the processes, data, and knowledge involved in the construction of whole ventricular multiscale models of cardiac electrophysiology. Our analysis reveals that models, simulations, and experiments are intertwined, in an assemblage that is a system itself, namely the model-simulation-experiment (MSE) system. We argue that validation is part of the whole MSE system and is contingent upon 1) understanding and coping with sources of biovariability; 2) testing and developing robust techniques and tools as a prerequisite to conducting physiological investigations; 3) defining and adopting standards to facilitate the interoperability of experiments, models, and simulations; 4) and understanding physiological validation as an iterative process that contributes to defining the specific aspects of cardiac electrophysiology the MSE system targets, rather than being only an external test, and that this is driven by advances in experimental and computational methods and the combination of both.     

Value of models for membrane budding

The budding of membranes and curvature generation is common to many forms of trafficking in cells. Clathrin-mediated endocytosis, as a prototypical example of trafficking, has been studied in great detail using a variety of experimental systems and methods. Recently, advances in experimental methods have led to great strides in insights on the molecular mechanisms and the spatiotemporal dynamics of the protein machinery associated with membrane curvature generation. These advances have been ably supported by computational models, which have given us insights into the underlying mechanical principles of clathrin-mediated endocytosis. On the other hand, targeted experimental perturbation of membranes has lagged behind that of proteins in cells. In this area, modeling is especially critical to interpret experimental measurements in a mechanistic context. Here, we discuss the contributions made by these models to our understanding of endocytosis and identify opportunities to strengthen the connections between models and experiments.     

Spatially resolved in silico modeling of NKG2D signaling kinetics suggests a key role of NKG2D and Vav1 Co-clustering in generating natural killer cell activation

Natural Killer (NK) cells provide key resistance against viral infections and tumors. A diverse set of activating and inhibitory NK cell receptors (NKRs) interact with cognate ligands presented by target host cells, where integration of dueling signals initiated by the ligand-NKR interactions determines NK cell activation or tolerance. Imaging experiments over decades have shown micron and sub-micron scale spatial clustering of activating and inhibitory NKRs. The mechanistic roles of these clusters in affecting downstream signaling and activation are often unclear. To this end, we developed a predictive in silico framework by combining spatially resolved mechanistic agent based modeling, published TIRF imaging data, and parameter estimation to determine mechanisms by which formation and spatial movements of activating NKG2D microclusters affect early time NKG2D signaling kinetics in a human cell line NKL. We show co-clustering of NKG2D and the guanosine nucleotide exchange factor Vav1 in NKG2D microclusters plays a dominant role over ligand (ULBP3) rebinding in increasing production of phospho-Vav1(pVav1), an activation marker of early NKG2D signaling. The in silico model successfully predicts several scenarios of inhibition of NKG2D signaling and time course of NKG2D spatial clustering over a short (~3 min) interval. Modeling shows the presence of a spatial positive feedback relating formation and centripetal movements of NKG2D microclusters, and pVav1 production offers flexibility towards suppression of activating signals by inhibitory KIR ligands organized in inhomogeneous spatial patterns (e.g., a ring). Our in silico framework marks a major improvement in developing spatiotemporal signaling models with quantitatively estimated model parameters using imaging data.     

Spatial clustering of average risks and risk trends in Bayesian disease mapping

Spatiotemporal disease mapping focuses on estimating the spatial pattern in disease risk across a set of nonoverlapping areal units over a fixed period of time. The key aim of such research is to identify areas that have a high average level of disease risk or where disease risk is increasing over time, thus allowing public health interventions to be focused on these areas. Such aims are well suited to the statistical approach of clustering, and while much research has been done in this area in a purely spatial setting, only a handful of approaches have focused on spatiotemporal clustering of disease risk. Therefore, this paper outlines a new modeling approach for clustering spatiotemporal disease risk data, by clustering areas based on both their mean risk levels and the behavior of their temporal trends. The efficacy of the methodology is established by a simulation study, and is illustrated by a study of respiratory disease risk in Glasgow, Scotland.     

Visual pattern discrimination by population retinal ganglion cells’ activities during natural movie stimulation

In the visual system, neurons often fire in synchrony, and it is believed that synchronous activities of group neurons are more efficient than single cell response in transmitting neural signals to down-stream neurons. However, whether dynamic natural stimuli are encoded by dynamic spatiotemporal firing patterns of synchronous group neurons still needs to be investigated. In this paper we recorded the activities of population ganglion cells in bullfrog retina in response to time-varying natural images (natural scene movie) using multi-electrode arrays. In response to some different brief section pairs of the movie, synchronous groups of retinal ganglion cells (RGCs) fired with similar but different spike events. We attempted to discriminate the movie sections based on temporal firing patterns of single cells and spatiotemporal firing patterns of the synchronous groups of RGCs characterized by a measurement of subsequence distribution discrepancy. The discrimination performance was assessed by a classification method based on Support Vector Machines. Our results show that different movie sections of the natural movie elicited reliable dynamic spatiotemporal activity patterns of the synchronous RGCs, which are more efficient in discriminating different movie sections than the temporal patterns of the single cells’ spike events. These results suggest that, during natural vision, the down-stream neurons may decode the visual information from the dynamic spatiotemporal patterns of the synchronous group of RGCs’ activities.