Complex network dynamics have been analyzed with models of systems of coupled switches or systems of coupled oscillators. However, many complex systems are composed of components with diverse dynamics whose interactions drive the system's evolution. We, therefore, introduce a new modeling framework that describes the dynamics of networks composed of both oscillators and switches. Both oscillator synchronization and switch stability are preserved in these heterogeneous, coupled networks. Furthermore, this model recapitulates the qualitative dynamics for the yeast cell cycle consistent with the hypothesized dynamics resulting from decomposition of the regulatory network into dynamic motifs. Introducing feedback into the cell-cycle network induces qualitative dynamics analogous to limitless replicative potential that is a hallmark of cancer. As a result, the proposed model of switch and oscillator coupling provides the ability to incorporate mechanisms that underlie the synchronized stimulus response ubiquitous in biochemical systems.Similarly, models of coupled oscillators along networks based upon the Kuramoto model [4] have been used to model synchronization of oscillators in diverse systems reviewed in [5]. In biochemical systems, in vivo oscillator synchronization has been observed in synthetic oscillatory fluorescent bacteria [6,7], yeast gene regulatory networks [8,9], and human cell fate decisions [10]. Such spontaneous synchronization has also been attributed to the development of the mammalian cardiac pacemaker cells (reviewed in [11]) and cortical systems (reviewed in [12]) including notably the circadian pacemaker (e.g., [13]). More recently, these network models have been found to be insufficient to model more complex dynamics in neuronal information transfer [12,[14][15][16][17] and cardiac arrhythmias [18][19][20][21]. These limitations extend to physical systems, such as the coupled lasers studied in [22]. Therefore, numerous studies have modified these network models to account for evolving networks [15,[23][24][25][26][27][28], dynamic frequencies [15, 29, 30], or phase delays [16,[31][32][33]. However, these mathematical modifications typically do not encode the mechanism underlying the limitations in the Kuramoto and Glass network models.We hypothesize that the observed limitations in the standard Kuramoto and Glass models arise from their exclusion of coupling components with qualitatively different dynamics. Several studies have inferred that biochemical systems contain "network motifs" with both oscillatory and switch-like dynamics [34,35]. The dynamics of these motifs are inferred from the topology of subgraphs in the networks of these systems. Their structures are statistically overrepresented in biochemical networks [36,37] such as intracellular regulatory networks [38], implicating evolutionary preservation (and thus utility) of these network motifs [39]. The dynamics of these motifs have been used to model yeast cell cycle regulation [40] and have been further confirmed in synthetic, desi...