Synthesis of biomolecular circuits for controlling molecular-scale processes is an important goal of synthetic biology with a wide range of in vitro and in vivo applications, including biomass maximization, nanoscale drug delivery, and many others. In this paper, we present new results on how abstract chemical reactions can be used to implement commonly used system theoretic operators such as the polynomial functions, rational functions and Hill-type nonlinearity. We first describe how idealised versions of multi-molecular reactions, catalysis, annihilation, and degradation can be combined to implement these operators. We then show how such chemical reactions can be implemented using enzyme-free, entropy-driven DNA reactions. Our results are illustrated through three applications: (1) implementation of a Stan-Sepulchre oscillator, (2) the computation of the ratio of two signals, and (3) a PI+antiwindup controller for regulating the output of a static nonlinear plant.