We describe the theoretical and computational framework for the Dynamic Signatures for Genetic Regulatory Network ( ) database. The motivation stems from urgent need to understand the global dynamics of biologically relevant signal transduction/gene regulatory networks that have at least 5 to 10 nodes, involve multiple interactions, and decades of parameters. The input to the database computations is a regulatory network, i.e. a directed graph with edges indicating up or down regulation. A computational model based on switching networks is generated from the regulatory network. The phase space dimension of this model equals the number of nodes and the associated parameter space consists of one parameter for each node (a decay rate), and three parameters for each edge (low level of expression, high level of expression, and threshold at which expression levels change). Since the nonlinearities of switching systems are piece-wise constant, there is a natural decomposition of phase space into cells from which the dynamics can be described combinatorially in terms of a state transition graph. This in turn leads to a compact representation of the global dynamics called an annotated Morse graph that identifies recurrent and nonrecurrent dynamics. The focus of this paper is on the construction of a natural computable finite decomposition of parameter space into domains where the annotated Morse graph description of dynamics is constant. We use this decomposition to construct an SQL database that can be effectively searched for dynamical signatures such as bistability, stable or unstable oscillations, and stable equilibria. We include two simple 3-node networks to provide small explicit examples of the type of information stored in the database. To demonstrate the computational capabilities of this system we consider a simple network associated with p53 that involves 5 nodes and a 29-dimensional parameter space.

The blood stage of the infection of the malaria parasite Plasmodium falciparum exhibits a 48-hour developmental cycle that culminates in the synchronous release of parasites from red blood cells, which triggers 48-hour fever cycles in the host. This cycle could be driven extrinsically by host circadian processes or by a parasite-intrinsic oscillator. To distinguish between these hypotheses, we examine the P. falciparum cycle in an in vitro culture system and show that the parasite has molecular signatures associated with circadian and cell cycle oscillators. Each of the four strains examined has a different period, which indicates strain-intrinsic period control. Finally, we demonstrate that parasites have low cell-to-cell variance in cycle period, on par with a circadian oscillator. We conclude that an intrinsic oscillator maintains Plasmodium’s rhythmic life cycle.

Many arthropods use filiform hairs as mechanoreceptors to detect air motion. In common house crickets (Acheta domestica) the hairs cover two antenna-like appendages called cerci at the rear of the abdomen. The biomechanical stimulus-response properties of individual filiform hairs have been investigated and modeled extensively in several earlier studies. However, only a few previous studies have considered viscosity-mediated coupling between pairs of hairs, and only in particular configurations. Here we present a model capable of calculating hair-to-hair coupling in arbitrary configurations. We simulate the coupled motion of a small group of mechanosensory hairs on a cylindrical section of cercus. We have found that the coupling effects are non-negligible, and likely constrain the operational characteristics of the cercal sensory array.

Mosquito host-seeking behavior and heterogeneity in host distribution are important factors in predicting the transmission dynamics of mosquito-borne infections such as dengue fever, malaria, chikungunya, and West Nile virus. We develop and analyze a new mathematical model to describe the effect of spatial heterogeneity on the contact rate between mosquito vectors and hosts. The model includes odor plumes generated by spatially distributed hosts, wind velocity, and mosquito behavior based on both the prevailing wind and the odor plume. On a spatial scale of meters and a time scale of minutes, we compare the effectiveness of different plume-finding and plume-tracking strategies that mosquitoes could use to locate a host. The results show that two different models of chemotaxis are capable of producing comparable results given appropriate parameter choices and that host finding is optimized by a strategy of flying across the wind until the odor plume is intercepted. We also assess the impact of changing the level of host aggregation on mosquito host-finding success near the end of the host-seeking flight. When clusters of hosts are more tightly associated on smaller patches, the odor plume is narrower and the biting rate per host is decreased. For two host groups of unequal number but equal spatial density, the biting rate per host is lower in the group with more individuals, indicative of an attack abatement effect of host aggregation. We discuss how this approach could assist parameter choices in compartmental models that do not explicitly model the spatial arrangement of individuals and how the model could address larger spatial scales and other probability models for mosquito behavior, such as Lévy distributions.

We present a new modeling and computational tool that computes rigorous summaries of network dynamics over large sets of parameter values. These summaries, organized in a database, can be searched for observed dynamics, e.g., bistability and hysteresis, to discover parameter regimes over which they are supported. We illustrate our approach on several networks underlying the restriction point of the cell cycle in humans and yeast. We rank networks by how robustly they support hysteresis, which is the observed phenotype. We find that the best 6-node human network and the yeast network share similar topology and robustness of hysteresis, in spite of having no homology between the corresponding nodes of the network. Our approach provides a new tool linking network structure and dynamics.