Hi,
If I understand your question correctly, your indirect response model is in derivative form as: dProtein/dt = kin - kout*(1+Imax*Conc/(Conc+IC50))*Protein, is that right? If not, can you elaborate on your question?
If so and if Protein is defined in units of mass, amount, or concentration, in other words if Protein can be described using a Species, you can use reactions to define input/production (kin) and loss (kout) processes. See the screenshot below.
If you'd like to define Protein as a measured dimensionless quantitiy or with units other than mass, amount, or concentration (that cannot be descibed by a Species) you can use a Rate Rule instead to define it as Protein = kin - kout*(1+Imax*Conc/(Conc+IC50))*Protein
where Protein is a Parameter.
To address your question about different Protein baseline/initial values: I assume you're interested in applying different initial values to Protein for each individual during parameter estimation. Is that right? If that's the case, you can store these initial values in different Variants and apply these variants in a Fit Data program using Variants and Dose Setup table. Please see the screenshot below. Fulden