We introduce the basic concepts and develop a theory for nonequilibrium steady-state biochemical systems applicable to analyzing large-scale complex isothermal reaction networks. In terms of the stoichiometric matrix, we dem-onstrate both Kirchhoff’s flux lawR‘J‘ ¼0 over a bio-chemical species, and potential law R‘l‘ ¼0overa reaction loop. They reflect mass and energy conservation, respectively.