MATLAB Answers


Issues fitting ka parameter

Asked by Brett Fleisher on 30 Apr 2019
Latest activity Answered by Brett Fleisher on 2 May 2019
I have simbiology 2019a. I am trying to fit ka using a one compartment PK model (I have Vd and Cl fixed).
I had previously scanned my model for ka and found that ka should be approximate 0.5. When I compared simulations to observed data, 0.5 allowed for the model to fit the data well.
I next tried fitting the data using the same concentrations in the model simulations. The ka's were much larger than expected. Upon further expection, the generated simulated concentration vs time profile doesn't produce a curve. It appears as if the model fittings are not generating a concentration profile. I checked if I was making a mistake in the dosing information or response information, but my fitting screen is similar to the example (
Any input or additional references to compare would be greatly appreciated.


Show 1 older comment
Sure, I can provide you with the data. I've attached the file below. My goal is to fit a 1 compartment oral comparment. I've run NCA on the data set and generated an apparent Volume of distribution (661.89 L) and apparent CL (23.39 L/hr). I need to account for bioavailability (F=0.45). Effect compartment and Gut compartments were set to 1 L. The drug concentration I'm estimating is in the central compartment.
My equations to explain this so far are the following:
Absorption = (ka*Drug_Oral_Dose*F)*Gut
Elimination = (CL/Vd*Drug_Central)*Effect
d(Drug_Central)/dt = 1/Effect*(Absorption - Elimination)
d(Drug_Oral_Dose)/dt = 1/Gut*(-Absorption)
Thanks in advance for the help
Just a quick comment that may help you. The way you are modeling bioavailability is essentially just changing the rate of absorption.
Dose*F should be incorporated as the initial condition of the gut compartment.
Great, thanks for help.
I had uploaded the wrong data set. I've attached for the corrected version for your reference.
Thanks again!

Sign in to comment.


More Answers in the  SimBiology Community

2 Answers

Answer by Sietse Braakman on 2 May 2019
 Accepted Answer

Hi Brett,
There are three ways you can approach this:
  1. by assuming that the dose is reduced by a factor F, such as Peter suggested above. However, keep in mind that your absorption rate is dependent on the concentration/amount of the drug_depot species. So by reducing the dose amount by a factor F, you also reduce the initial absorption rate.
  2. by assuming that a fraction F of the drug_depot species makes its way to the drug_central species, while a factor (1-F) makes is passed through the intestine without being absorbed. This ensures that the first order rate is not affected by the reduced dose amount. We have an example on bioavailability on our website as well:
  3. Get both IV and Oral data and approximate F by F = AUC_Oral/AUC_IV, where AUC is the area under the curve of the drug concentration in the central compartment.


Sign in to comment.

Answer by Brett Fleisher on 2 May 2019

Hi Sietse,
Thank you again for the help. I've uploaded the model you had helped me with for future reference.
Much appreciated!


Sign in to comment.