The MODEL Procedure
Example 26.10 Systems of Differential Equations
The following is a simplified reaction scheme for the competitive inhibitors with recombinant human renin (Morelock et al. 1995).
Output 26.10.1: Competitive Inhibition of Recombinant Human Renin
In Output 26.10.1, E=enzyme, D=probe, and I=inhibitor.
The differential equations that describe this reaction scheme are as follows:
![\[ \frac{d\mi{D} }{dt} = k1r{\ast }\mi{ED} - k1f {\ast }\mi{E} {\ast }\mi{D} \]](images/etsug_model0661.png)
![\[ \frac{d\mi{ED} }{dt} = k1f{\ast } \mi{E} {\ast }\mi{D} - k1r{\ast }\mi{ED} \]](images/etsug_model0662.png)
![\[ \frac{d\mi{E} }{dt} = k1r{\ast } \mi{ED} - k1f{\ast } \mi{E} {\ast } \mi{D} + k2r {\ast } \mi{EI} - k2f {\ast } \mi{E} {\ast }\mi{I} \]](images/etsug_model0663.png)
![\[ \frac{d\mi{EI} }{dt} = k2f{\ast } \mi{E} {\ast }\mi{I} - k2r {\ast } \mi{EI} \]](images/etsug_model0664.png)
![\[ \frac{d\mi{I} }{dt} = k2r{\ast } \mi{EI} - k2f {\ast } \mi{E} {\ast } \mi{I} \]](images/etsug_model0665.png)
For this system, the initial values for the concentrations are derived from equilibrium considerations (as a function of parameters) or are provided as known values.
The experiment used to collect the data was carried out in two ways; preincubation (type='disassoc') and no preincubation (type='assoc'). The data also contain repeated measurements. The data contain values for fluorescence F, which is a function of concentration. Since there are no direct data for the concentrations, all the differential equations are simulated dynamically.
The SAS statements used to fit this model are as follows:
title1 'Systems of Differential Equations Example'; proc sort data=fit; by type time; run; %let k1f = 6.85e6 ; %let k1r = 3.43e-4 ; %let k2f = 1.8e7 ; %let k2r = 2.1e-2 ; %let qf = 2.1e8 ; %let qb = 4.0e9 ; %let dt = 5.0e-7 ; %let et = 5.0e-8 ; %let it = 8.05e-6 ;
proc model data=fit;
parameters qf = 2.1e8
qb = 4.0e9
k2f = 1.8e5
k2r = 2.1e-3
l = 0;
k1f = 6.85e6;
k1r = 3.43e-4;
/* Initial values for concentrations */
control dt 5.0e-7
et 5.0e-8
it 8.05e-6;
/* Association initial values --------------*/
if type = 'assoc' and time=0 then do;
ed = 0;
/* solve quadratic equation ----------*/
a = 1;
b = -(&it+&et+(k2r/k2f));
c = &it*&et;
ei = (-b-(((b**2)-(4*a*c))**.5))/(2*a);
d = &dt-ed;
i = &it-ei;
e = &et-ed-ei;
end;
/* Disassociation initial values ----------*/
if type = 'disassoc' and time=0 then do;
ei = 0;
a = 1;
b = -(&dt+&et+(&k1r/&k1f));
c = &dt*&et;
ed = (-b-(((b**2)-(4*a*c))**.5))/(2*a);
d = &dt-ed;
i = &it-ei;
e = &et-ed-ei;
end;
if time ne 0 then do;
dert.d = k1r* ed - k1f *e *d;
dert.ed = k1f* e *d - k1r*ed;
dert.e = k1r* ed - k1f* e * d + k2r * ei - k2f * e *i;
dert.ei = k2f* e *i - k2r * ei;
dert.i = k2r * ei - k2f* e *i;
end;
/* L - offset between curves */
if type = 'disassoc' then
F = (qf*(d-ed)) + (qb*ed) -L;
else
F = (qf*(d-ed)) + (qb*ed);
fit F / method=marquardt;
run;
This estimation requires the repeated simulation of a system of 41 differential equations (5 base differential equations and 36 differential equations to compute the partials with respect to the parameters).
The results of the estimation are shown in Output 26.10.2.
Output 26.10.2: Kinetics Estimation