skip to main content

Title: Mapping parameter spaces of biological switches
Since the seminal 1961 paper of Monod and Jacob, mathematical models of biomolecular circuits have guided our understanding of cell regulation. Model-based exploration of the functional capabilities of any given circuit requires systematic mapping of multidimensional spaces of model parameters. Despite significant advances in computational dynamical systems approaches, this analysis remains a nontrivial task. Here, we use a nonlinear system of ordinary differential equations to model oocyte selection in Drosophila , a robust symmetry-breaking event that relies on autoregulatory localization of oocyte-specification factors. By applying an algorithmic approach that implements symbolic computation and topological methods, we enumerate all phase portraits of stable steady states in the limit when nonlinear regulatory interactions become discrete switches. Leveraging this initial exact partitioning and further using numerical exploration, we locate parameter regions that are dense in purely asymmetric steady states when the nonlinearities are not infinitely sharp, enabling systematic identification of parameter regions that correspond to robust oocyte selection. This framework can be generalized to map the full parameter spaces in a broad class of models involving biological switches.  more » « less
Award ID(s):
1934924 1839294
Author(s) / Creator(s):
; ; ; ; ; ; ;
Finley, Stacey
Date Published:
Journal Name:
PLOS Computational Biology
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. Concerns about food security under climate change motivate efforts to better understand future changes in crop yields.Process-based crop models, which represent plant physiological and soil processes, are necessary tools for this purpose since they allow representing future climate and management conditions not sampled in the historical record and new locations to which cultivation may shift.However, process-based crop models differ in many critical details, and their responses to different interacting factors remain only poorly understood.The Global Gridded Crop Model Intercomparison (GGCMI) Phase 2 experiment, an activity of the Agricultural Model Intercomparison and Improvement Project (AgMIP), is designed to provide a systematic parameter sweep focused on climate change factors and their interaction with overall soil fertility, to allow both evaluating model behavior and emulating model responses in impact assessment tools.In this paper we describe the GGCMI Phase 2 experimental protocol and its simulation data archive.A total of 12 crop models simulate five crops with systematic uniform perturbations of historical climate, varying CO2, temperature, water supply, and applied nitrogen (“CTWN”) for rainfed and irrigated agriculture, and a second set of simulations represents a type of adaptation by allowing the adjustment of growing season length.We present some crop yield results to illustrate general characteristics of the simulations and potential uses of the GGCMI Phase 2 archive.For example, in cases without adaptation, modeled yields show robust decreases to warmer temperatures in almost all regions, with a nonlinear dependence that means yields in warmer baseline locations have greater temperature sensitivity.Inter-model uncertainty is qualitatively similar across all the four input dimensions but is largest in high-latitude regions where crops may be grown in the future. 
    more » « less
  2. Abstract

    In neuroscience, computational modeling is an effective way to gain insight into cortical mechanisms, yet the construction and analysis of large-scale network models—not to mention the extraction of underlying principles—are themselves challenging tasks, due to the absence of suitable analytical tools and the prohibitive costs of systematic numerical exploration of high-dimensional parameter spaces. In this paper, we propose a data-driven approach assisted by deep neural networks (DNN). The idea is to first discover certain input-output relations, and then to leverage this information and the superior computation speeds of the well-trained DNN to guide parameter searches and to deduce theoretical understanding. To illustrate this novel approach, we used as a test case a medium-size network of integrate-and-fire neurons intended to model local cortical circuits. With the help of an accurate yet extremely efficient DNN surrogate, we revealed the statistics of model responses, providing a detailed picture of model behavior. The information obtained is both general and of a fundamental nature, with direct application to neuroscience. Our results suggest that the methodology proposed can be scaled up to larger and more complex biological networks when used in conjunction with other techniques of biological modeling.

    more » « less
  3. Abstract

    Planets in synchronous rotation around low-mass stars are the most salient targets for current ground- and space-based missions to observe and characterize. Such model calculations can help to prioritize targets for observation with current and future missions; however, intrinsic differences in the complexity and physical parameterizations of various models can lead to different predictions of a planet’s climate state. Understanding model differences is necessary if such models are to guide target selection and aid in the analysis of observations. This paper presents a protocol to intercompare models of a hypothetical planet with a 15-day synchronous rotation period around a 3000 K blackbody star across a parameter space of surface pressure and incident instellation. We conduct a sparse sample of 16 cases from a previously published exploration of this parameter space with the ExoPlaSim model. By selecting particular cases across this broad parameter space, the SAMOSA intercomparison will identify areas where simpler models are sufficient, as well as areas where more complex GCMs are required. Our preliminary comparison using ExoCAM shows general consistency between the climate state predicted by ExoCAM and ExoPlaSim except in regions of the parameter space most likely to be in a steam atmosphere or incipient runaway greenhouse state. We use this preliminary analysis to define several options for participation in the intercomparison by models of all levels of complexity. The participation of other GCMs is crucial to understand how the atmospheric states across this parameter space differ with model capabilities.

    more » « less
  4. Abstract

    We consider steady state solutions of the massive, asymptotically flat, spherically symmetric Einstein–Vlasov system, i.e., relativistic models of galaxies or globular clusters, and steady state solutions of the Einstein–Euler system, i.e., relativistic models of stars. Such steady states are embedded into one-parameter families parameterized by their central redshift$$\kappa >0$$κ>0. We prove their linear instability when$$\kappa $$κis sufficiently large, i.e., when they are strongly relativistic, and prove that the instability is driven by a growing mode. Our work confirms the scenario of dynamic instability proposed in the 1960s by Zel’dovich & Podurets (for the Einstein–Vlasov system) and by Harrison, Thorne, Wakano, & Wheeler (for the Einstein–Euler system). Our results are in sharp contrast to the corresponding non-relativistic, Newtonian setting. We carry out a careful analysis of the linearized dynamics around the above steady states and prove an exponential trichotomy result and the corresponding index theorems for the stable/unstable invariant spaces. Finally, in the case of the Einstein–Euler system we prove a rigorous version of the turning point principle which relates the stability of steady states along the one-parameter family to the winding points of the so-called mass-radius curve.

    more » « less
  5. One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data. 
    more » « less