Biochemical oscillations, such as glycolytic oscillations, are often believed to be caused by a single so-called ‘oscillophore’. The main characteristics of yeast glycolytic oscillations, such as frequency and amplitude, are however controlled by several enzymes. In this paper, we develop a method to quantify to which extent any enzyme determines the occurrence of oscillations. Principles extrapolated from metabolic control analysis are applied to calculate the control exerted by individual enzymes on the real and imaginary parts of the eigenvalues of the Jacobian matrix. We propose that the control exerted by an enzyme on the real part of the smallest eigenvalue, in terms of absolute value, quantifies to which extent that enzyme contributes to the emergence of instability. Likewise the control exerted by an enzyme on the imaginary part of complex eigenvalues may serve to quantify the extent to which that enzyme contributes to the tendency of the system to oscillate. The method was applied both to a core model and to a realistic model of yeast glycolytic oscillations. Both the control over stability and the control over oscillatory tendency were distributed among several enzymes, of which glucose transport, pyruvate decarboxylase and ATP utilization were the most important. The distributions of control were different for stability and oscillatory tendency, showing that control of instability does not imply control of oscillatory tendency nor vice versa. The control coefficients summed up to 1, suggesting the existence of a new summation theorem. These results constitute proof that glycolytic oscillations in yeast are not caused by a single oscillophore and provide a new, subtle, definition for the oscillophore strength of an enzyme.
Read full abstract