1. Biochemical oscillators—a search procedure
- Author
-
Paul E. Rapp
- Subjects
Statistics and Probability ,General Immunology and Microbiology ,Differential equation ,Applied Mathematics ,Mathematical analysis ,Describing function ,Block (permutation group theory) ,General Medicine ,General Biochemistry, Genetics and Molecular Biology ,Reaction rate ,Nonlinear system ,Simultaneous equations ,Modeling and Simulation ,General Agricultural and Biological Sciences ,Constant (mathematics) ,Fourier series ,Mathematics ,Mathematical physics - Abstract
The paper shows how the describing function technique can be employed in a search for periodic solutions to large systems of nonlinear differential equations that describe chemical reaction networks of a type encountered in biological contexts. The equations studied describe a reaction scheme in which the first reaction is inhibited by the product of the last reaction. The differential equations take the form dS 1 ⧸dt = K⧸( 1 + α(S n ) ρ ) and dS j ⧸dt = g j− 1 S j− 1 −b j S j , j = 2, 3,⋯n . The parameter ϱ is a positive integer and represents the number of molecules of inhibitor required to block the S 1 to S 2 reaction. The equations can be interpreted as a simple model of inhibition of transcription, or they can be interpreted as an inhibited enzyme reaction system in which initial substrate concentration and enzyme concentrations are held constant. The system of simultaneous equations is reduced to a single equation in S n : (p + b 1 )(p + b 2 )⋯(p + b n )S n = d 1 /( 1 + d 2 S ρ n ) . Here p is the differentiation operator and d 1 and d 2 are positive constants expressible in terms of α, g j 's, and K . The transfer function block diagram equivalent to this equation is established, and is seen to have closed feedback loop form and to contain a linear block and a nonlinear block. The nonlinear block is replaced by the describing function which is related to the Fourier series of the nonlinear block's output. A study of the system containing the describing function is then developed. It is shown that if ϱ =1, then the describing function analysis indicates, for arbitrary n , that stable limit cycles will not appear if K , α, g j , and b j are positive. Since these constants represents the products of reaction rate constants and concentrations, this requirement is automatically met. The nonexistence of stable limit cycles for ρ =1 holds for all values of n ; it is also shown that for n ⩽4 unstable limit cycles are also necessarily absent. The extension of the method for larger values of ρ is outlined. Since the 1 (1 + z ρ ) nonlinearity frequently appears in biological problems, the describing function will be useful in several areas of research.
- Published
- 1975
- Full Text
- View/download PDF