Complex regulation of Hsf1-Skn7 activities by the catalytic subunits of PKA in Saccharomyces cerevisiae: experimental and computational evidences

BACKGROUND: The cAMP-dependent protein kinase regulatory network (PKA-RN) regulates metabolism, memory, learning, development, and response to stress. Previous models of this network considered the catalytic subunits (CS) as a single entity, overlooking their functional individualities. Furthermore, PKA-RN dynamics are often measured through cAMP levels in nutrient-depleted cells shortly after being fed with glucose, dismissing downstream physiological processes. RESULTS: Here we show that temperature stress, along with deletion of PKA-RN genes, significantly affected HSE-dependent gene expression and the dynamics of the PKA-RN in cells growing in exponential phase. Our genetic analysis revealed complex regulatory interactions between the CS that influenced the inhibition of Hsf1/Skn7 transcription factors. Accordingly, we found new roles in growth control and stress response for Hsf1/Skn7 when PKA activity was low (cdc25Δ cells). Experimental results were used to propose an interaction scheme for the PKA-RN and to build an extension of a classic synchronous discrete modeling framework. Our computational model reproduced the experimental data and predicted complex interactions between the CS and the existence of a repressor of Hsf1/Skn7 that is activated by the CS. Additional genetic analysis identified Ssa1 and Ssa2 chaperones as such repressors. Further modeling of the new data foresaw a third repressor of Hsf1/Skn7, active only in theabsence of Tpk2. By averaging the network state over all its attractors, a good quantitative agreement between computational and experimental results was obtained, as the averages reflected more accurately the population measurements. CONCLUSIONS: The assumption of PKA being one molecular entity has hindered the study of a wide range of behaviors. Additionally, the dynamics of HSE-dependent gene expression cannot be simulated accurately by considering the activity of single PKA-RN components (i.e., cAMP, individual CS, Bcy1, etc.). We show that the differential roles of the CS are essential to understand the dynamics of the PKA-RN and its targets. Our systems level approach, which combined experimental results with theoretical modeling, unveils the relevance of the interaction scheme for the CS and offers quantitative predictions for several scenarios (WT vs. mutants in PKA-RN genes and growth at optimal temperature vs. heat shock). ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12918-015-0185-8) contains supplementary material, which is available to authorized users.