Stomata, microscopic pores in leaf surfaces through which water loss and carbon dioxide uptake occur, are closed in response to drought by the phytohormone abscisic acid (ABA). This process is vital for drought tolerance and has been the topic of extensive experimental investigation in the last decades. Although a core signaling chain has been elucidated consisting of ABA binding to receptors, which alleviates negative regulation by protein phosphatases 2C (PP2Cs) of the protein kinase OPEN STOMATA 1 (OST1) and ultimately results in activation of anion channels, osmotic water loss, and stomatal closure, over 70 additional components have been identified, yet their relationships with each other and the core components are poorly elucidated. We integrated and processed hundreds of disparate observations regarding ABA signal transduction responses underlying stomatal closure into a network of 84 nodes and 156 edges and, as a result, established those relationships, including identification of a 36-node, strongly connected (feedback-rich) component as well as its in- and out-components. The network’s domination by a feedback-rich component may reflect a general feature of rapid signaling events. We developed a discrete dynamic model of this network and elucidated the effects of ABA plus knockout or constitutive activity of 79 nodes on both the outcome of the system (closure) and the status of all internal nodes. The model, with more than 10 24 system states, is far from fully determined by the available data, yet model results agree with existing experiments in 82 cases and disagree in only 17 cases, a validation rate of 75%. Our results reveal nodes that could be engineered to impact stomatal closure in a controlled fashion and also provide over 140 novel predictions for which experimental data are currently lacking. Noting the paucity of wet-bench data regarding combinatorial effects of ABA and internal node activation, we experimentally confirmed several predictions of the model with regard to reactive oxygen species, cytosolic Ca 2+ (Ca 2+ c), and heterotrimeric G-protein signaling. We analyzed dynamics-determining positive and negative feedback loops, thereby elucidating the attractor (dynamic behavior) repertoire of the system and the groups of nodes that determine each attractor. Based on this analysis, we predict the likely presence of a previously unrecognized feedback mechanism dependent on Ca 2+ c. This mechanism would provide model agreement with 10 additional experimental observations, for a validation rate of 85%. Our research underscores the importance of feedback regulation in generating robust and adaptable biological responses. The high validation rate of our model illustrates the advantages of discrete dynamic modeling for complex, nonlinear systems common in biology.
Guard cells, located in pairs on the surface of plant leaves, circumscribe microscopic pores called stomata, through which plants take up gaseous carbon dioxide that will be fixed by photosynthesis into sugars. However, plants also inevitably lose water vapor to the atmosphere through open stomata. Under drought conditions, the plant hormone abscisic acid (ABA) causes volume changes in guard cells that result in stomatal closure, thereby restricting water loss. Given the paramount importance of drought tolerance for plant survival, it is important to understand the cellular mechanisms underlying guard cell response to ABA, and over 100 studies in the literature have addressed this topic. We synthesized this information into a network that contains 84 cellular components and 156 interactions between them and then applied logic-based analyses to predict how these components coordinately transduce the ABA signal. We identified several positive feedback loops and mutual inhibition loops that can lead to sustained activity of their constituent components in the presence, or absence, of ABA. Control of these loops, for example, by other stimuli present in the natural environment, may sensitize the system to ABA. We validated some of the novel predictions from our model with new experiments.