The impact of fluctuations on the dynamical behaviour of complex biological systems is a longstanding issue, whose understanding would elucidate how evolutionary pressure tends to modulate intrinsic noise. Using the Itō stochastic differential equation formalism, we performed analytic and numerical analyses of model systems containing different molecular species in contact with the environment and interacting with each other through mass-action kinetics. For networks of zero deficiency, which admit a detailed-or complex-balanced steady state, all molecular species are uncorrelated and their Fano factors are Poissonian. Systems of higher deficiency have nonequilibrium steady states and non-zero reaction fluxes flowing between the complexes. When they model homo-oligomerization, the noise on each species is reduced when the flux flows from the oligomers of lowest to highest degree, and amplified otherwise. In the case of hetero-oligomerization systems, only the noise on the highest-degree species shows this behaviour.& 2018 The Author(s) Published by the Royal Society. All rights reserved.transcription [11][12][13][14]. This mechanism has indeed been shown to lower noise while reducing the metabolic cost of protein production, and speeds up the rise-times of transcription units [15]. However, not all systems with negative feedback loops decrease in the intrinsic noise levels [16][17][18]. Similarly, although it is generally accepted that positive feedback loops tend to increase noise levels [19], some appear to decrease them [20]. Hence, the problem is far from being totally elucidated.Despite the many valuable advances in the field, the mechanisms used to amplify or to suppress the fluctuation levels need to be further understood and clarified. Indeed, the huge complexity of biological systems, their dependence on a large number of variables and the system-to-system variability make the unravelling of these issues, whether using experimental or computational approaches, a highly non-trivial task.More specifically, while noise control is relatively well understood for small and simple networks, it is still far from clear how the fluctuations propagate through more general and complicated networks and what is the link between network topology and complexity with noise buffering or amplification. Different investigations addressed these issues from various perspectives, for example by characterizing the stochastic properties of the chemical reaction networks (CRNs) and studying the propagation of the fluctuations [21][22][23][24]. From a physics-oriented perspective, other studies have analysed the connection between the nonequilibrium thermodynamic properties of the network and the noise level [25 -27]. It has furthermore been shown that the increase in network complexity tends to decrease intrinsic noise as well as to reduce the effect of extrinsic noise for some multistable model systems [28], whereas the dependence of noise reduction or amplification on the system parameters has been found in [29].This paper focu...