Biological systems contain a large number of molecules that have diverse interactions. A fruitful path to understanding these systems is to represent them with interaction networks, and then describe flow processes in the network with a dynamic model. Boolean modeling, the simplest discrete dynamic modeling framework for biological networks, has proven its value in recapitulating experimental results and making predictions. A first step and major roadblock to the widespread use of Boolean networks in biology is the laborious network inference and construction process. Here we present a streamlined network inference method that combines the discovery of a parsimonious network structure and the identification of Boolean functions that determine the dynamics of the system. This inference method is based on a causal logic analysis method that associates a logic type (sufficient or necessary) to node-pair relationships (whether promoting or inhibitory). We use the causal logic framework to assimilate indirect information obtained from perturbation experiments and infer relationships that have not yet been documented experimentally. We apply this inference method to a well-studied process of hormone signaling in plants, the signaling underlying abscisic acid (ABA)-induced stomatal closure. Applying the causal logic inference method significantly reduces the manual work typically required for network and Boolean model construction. The inferred model agrees with the manually curated model. We also test this method by re-inferring a network representing epithelial to mesenchymal transition based on a subset of the information that was initially used to construct the model. We find that the inference method performs well for various likely scenarios of inference input information. We conclude that our method is an effective approach toward inference of biological networks and can become an efficient step in the iterative process between experiments and computations.

Stomatal pores play a central role in the control of carbon assimilation and plant water status. The guard cell pair that borders each pore integrates information from environmental and endogenous signals and accordingly swells or deflates, thereby increasing or decreasing the stomatal aperture. Prior research shows that there is a complex cellular network underlying this process. We have previously constructed a signal transduction network and a Boolean dynamic model describing stomatal closure in response to signals including the plant hormone abscisic acid (ABA), calcium or reactive oxygen species (ROS). Here, we improve the Boolean network model such that it captures the biologically expected response of the guard cell in the absence or following the removal of a closure-inducing signal such as ABA or external Ca2+. The expectation from the biological system is reversibility, i.e., the stomata should reopen after the closing signal is removed. We find that the model's reversibility is obstructed by the previously assumed persistent activity of four nodes. By introducing time-dependent Boolean functions for these nodes, the model recapitulates stomatal reopening following the removal of a signal. The previous version of the model predicts ∼20% closure in the absence of any signal due to uncertainty regarding the initial conditions of multiple network nodes. We systematically test and adjust these initial conditions to find the minimally restrictive combinations that appropriately result in open stomata in the absence of a closure signal. We support these results by an analysis of the successive stabilization of feedback motifs in the network, illuminating the system's dynamic progression toward the open or closed stomata state. This analysis particularly highlights the role of cytosolic calcium oscillations in causing and maintaining stomatal closure. Overall, we illustrate the strength of the Boolean network modeling framework to efficiently capture cellular phenotypes as emergent outcomes of intracellular biological processes.

The first mitigation response to the Covid-19 pandemic was to limit person-to-person interaction as much as possible. This was implemented by the temporary closing of many workplaces and people were required to follow social distancing. Networks are a great way to represent interactions among people and the temporary severing of these interactions. Here, we present a network model of human-human interactions that could be mediators of disease spread. The nodes of this network are individuals and different types of edges denote family cliques, workplace interactions, interactions arising from essential needs, and social interactions. Each individual can be in one of four states: susceptible, infected, immune, and dead. The network and the disease parameters are informed by the existing literature on Covid-19. Using this model, we simulate the spread of an infectious disease in the presence of various mitigation scenarios. For example, lockdown is implemented by deleting edges that denote non-essential interactions. We validate the simulation results with the real data by matching the basic and effective reproduction numbers during different phases of the spread. We also simulate different possibilities of the slow lifting of the lockdown by varying the transmission rate as facilities are slowly opened but people follow prevention measures like wearing masks etc. We make predictions on the probability and intensity of a second wave of infection in each of these scenarios.

The plant hormone abscisic acid (ABA) promotes stomatal closure via multifarious cellular signaling cascades. Our previous comprehensive reconstruction of the stomatal closure network resulted in an 81-node network with 153 edges. Discrete dynamic modeling utilizing this network reproduced over 75% of experimental observations but a few experimentally supported results were not recapitulated. Here we identify predictions that improve the agreement between model and experiment. We performed dynamics-preserving network reduction, resulting in a condensed 49 node and 113 edge stomatal closure network that preserved all dynamics-determining network motifs and reproduced the predictions of the original model. We then utilized the reduced network to explore cases in which experimental activation of internal nodes in the absence of ABA elicited stomatal closure in wet bench experiments, but not in our in silico model. Our simulations revealed that addition of a single edge, which allows indirect inhibition of any one of three PP2C protein phosphatases (ABI2, PP2CA, HAB1) by cytosolic Ca2+ elevation, resolves the majority of the discrepancies. Consistent with this hypothesis, we experimentally show that Ca2+ application to cellular lysates at physiological concentrations inhibits PP2C activity. The model augmented with this new edge provides new insights into the role of cytosolic Ca2+ oscillations in stomatal closure, revealing a mutual reinforcement between repeated increases in cytosolic Ca2+ concentration and a self-sustaining feedback circuit inside the signaling network. These results illustrate how iteration between model and experiment can improve predictions of highly complex cellular dynamics.

BACKGROUND:Cellular behaviors are governed by interaction networks among biomolecules, for example gene regulatory and signal transduction networks. An often used dynamic modeling framework for these networks, Boolean modeling, can obtain their attractors (which correspond to cell types and behaviors) and their trajectories from an initial state (e.g. a resting state) to the attractors, for example in response to an external signal. The existing methods however do not elucidate the causal relationships between distant nodes in the network.RESULTS:In this work, we propose a simple logic framework, based on categorizing causal relationships as sufficient or necessary, as a complement to Boolean networks. We identify and explore the properties of complex subnetworks that are distillable into a single logic relationship. We also identify cyclic subnetworks that ensure the stabilization of the state of participating nodes regardless of the rest of the network. We identify the logic backbone of biomolecular networks, consisting of external signals, self-sustaining cyclic subnetworks (stable motifs), and output nodes. Furthermore, we use the logic framework to identify crucial nodes whose override can drive the system from one steady state to another. We apply these techniques to two biological networks: the epithelial-to-mesenchymal transition network corresponding to a developmental process exploited in tumor invasion, and the network of abscisic acid induced stomatal closure in plants. We find interesting subnetworks with logical implications in these networks. Using these subgraphs and motifs, we efficiently reduce both networks to succinct backbone structures.CONCLUSIONS:The logic representation identifies the causal relationships between distant nodes and subnetworks. This knowledge can form the basis of network control or used in the reverse engineering of networks.

We present a simplistic approach to determine the electronic-bandgap (Eg) and the exciton binding energy (Ex) of hybrid perovskite (ABX3 where, A= CH3NH3+, B = Pb2+, X = Cl, Br or I ion) semiconductor film. In this study, we show the role of the halide anion (X-) on Eg and Ex using a simple room temperature UV-Vis absorption spectrum measurements and model our results using Elliot's theory of Wannier exciton model to determine the Eg and Ex. We also found that the disorder primarily arises from structure and not composition disorder. This is just a consequence of the fact the Urbach tail, which is a measure of disorder is reasonably independent of the halide anion among three perovskite films.

A description of Orthogonal Tensor Hermite Polynomials in 3-D is presented. These polynomials, as introduced by Grad in 1949 [1], can be used to obtain a series solution to the Boltzmann Transport Equation. The properties that are explored are scaling, translation and rotation. Order 6 Hermite Tensors are studied while obtaining the rotation relations. From the scaling of the independent variables of particle velocities, a criterion on temperature is obtained which implies that the equation can be applied to binary gas mixtures only if the temperature of the hotter constituent is less than four times that of the cooler one. This criterion and other properties of the tensor hermite polynomials obtained in this paper can be used to study gas dynamics in the thermosphere.