Central functions in the cell are often linked to complex dynamic behaviours, such as sustained oscillations and multistability, in a biochemical reaction network. Determination of the specific mechanisms underlying such behaviours is important, e.g. to determine sensitivity, robustness, and modelling requirements of given cell functions. In this work we adopt a systems approach to the analysis of complex behaviours in intracellular reaction networks, described by ordinary differential equations with known kinetic parameters. We propose to decompose the overall system into a number of low complexity subsystems, and consider the importance of interactions between these in generating specific behaviours. Rather than analysing the network in a state corresponding to the complex non-linear behaviour, we move the system to the underlying unstable steady state, and focus on the mechanisms causing destabilisation of this steady state. This is motivated by the fact that all complex behaviours in unforced systems can be traced to destabilisation (bifurcation) of some steady state, and hence enables us to use tools from linear system theory to qualitatively analyse the sources of given network behaviours. One important objective of the present study is to see how far one can come with a relatively simple approach to the analysis of highly complex biochemical networks. The proposed method is demonstrated by application to a model of mitotic control in Xenopus frog eggs, and to a model of circadian oscillations in Drosophila. In both examples we are able to identify the subsystems, and the related interactions, which are instrumental in generating the observed complex non-linear behaviours.