I want to estimate the parameters of my own model. The model contains several reactions, most of which obey Uni-Uni Michaelis-Menten kinetic law. There are also additional two reactions that are actually catalyzed by the same enzyme, and thus must have the same parameter values. % Add custom-made kinetic law ratelawexpression1 = 'kcat*e0*(S1*S2/Km1*Km2)/((1+S1/Km1)*(1+S2/Km2))'; mm_1 = sbioabstractkineticlaw('CustomMadeLaw',ratelawexpression1); set(mm_1,'SpeciesVariables',{'S1','S2'}); set(mm_1,'ParameterVariables',{'kcat','e0','Km1','Km2'}); sbioaddtolibrary(mm_1); % Add second reaction to the model r2 = addreaction(model,'B1 + C -> D'); k2 = addkineticlaw(r2,'CustomMadeLaw'); p2a = addparameter(k2,'kcat',110); p2b = addparameter(k2,'e0',0.6); p2c = addparameter(k2,'Km_B1',0.01); p2d = addparameter(k2,'Km_C',0.01); set(k2,'ParameterVariableNames',{'kcat','e0','Km_B1','Km_C'}); set(k2,'SpeciesVariableNames',{'B1','C'}); set(r2,'Name','r2'); % Add third reaction to the model r3 = addreaction(model,'B2 + C -> F'); k3 = addkineticlaw(r3,'CustomMadeLaw'); p3a = addparameter(k3,'kcat',110); p3b = addparameter(k3,'e0',0.6); p3c = addparameter(k3,'Km_B2',0.01); p3d = addparameter(k3,'Km_C',0.01); set(k3,'ParameterVariableNames',{'kcat','e0','Km_B2','Km_C'}); set(k3,'SpeciesVariableNames',{'B2','C'}); set(r3,'Name','r3'); From what I understand when using sbiofit, when performing parameter estimation, the kcat, e0 and Km_C parameters of both r2 and r3 are considered different. But I would like to configure the parameter estimation settings so that those parameters are considered the same (r2.kcat should be considered the same parameter as r3.kcat, r2.e0 is r3.e0, r2.Km_C is r3.Km_C, etc.). How should I write the model or how do I change the settings of sbiofit so that I can do the previously-mentioned task?