Open-Systems-Pharmacology / Forum

Discussion forum for the Open Systems Pharmacology Project
65 stars 19 forks source link

Unestimate the concentration of vemurafenib #424

Closed wangwei1619 closed 4 years ago

wangwei1619 commented 5 years ago

Hi, all

I am doing a vemurafenib (Log P 5.1, pka 7.1, MW 500, fu <0.01)simualtion, but I can't reach the Cmax like the observed data. image

The hepatic clearance rate is according to publication value, and the fu is fitted to be 2.70E-4 to get consistent with blood plasma ratio of 0.72. The distribution model is Poulin & Theil. It's an oral formulation, however, choosing dissolved formulation can only make littele difference I have condiser transporters like ABCG2, ABCB1, OATP1B3 and also the EHC. After parameters identification, the parameters for transporters have been set to make the distribution in brain be of small fraction, which is agreement with publication, and biliary clearance get very small so that EHC might get insignificant.

Is there any else factor to better the fitting? Moreover, published data says vemurafenib is mainly excreted unchanged by biliary pathway, especially within first 48 h, and I think it suggest the EHC for vemurafenib is important. Can EHC influence the Cmax?

Any comment is appreciated.

prvmalik commented 5 years ago

Try optimizing intestinal permeability (transcellular) in the parameter identification tool.

With pKa 7.1 a notable portion of your molecule will be charged at physiologic pH. Try switching the cell permeability calculation algorithm to charge-dependent schmitt. Of course this will change your blood:plasma ratio.

On opinion alone, I find that using a dummy enzyme inside hepatic intracellular space works better than using the total hepatic plasma clearance function, especially in the case when permeability may be limiting to elimination.

And of course, this drug is a major substrate for transporters. You'll have to read into the literature to find which ones and to which degrees the drug maybe transported. Try using literature-based Km for PgP and then optimizing Vmax according to the fraction excreted unchanged in bile/feces.

wangwei1619 commented 5 years ago

Dear @prvmalik Thanks for your help. I tried to find Km, Vmax for transporters on vemurafenib but failed. What I know is the fraction of unchanged vemurafenib excreted to feces should be 0.38 at 48h and 0.55 at 96h, but these data is for 960 mg vemurafenib mutiple dose after 14 days of treatment (on day 15 , 960 mg vemurafenib labled with isotope was given and fraction was then calculated )not for 160mg which I am simulating. The observed data is far above my result (with the difference of 7 magnitude). It seems I can never reach that degree.

wangwei1619 commented 5 years ago

Besides the problem above, I also find something very weird. Basing on the parameters for transporters, the plasma concentration curve would be like this even if the specific clearance is set to be 0: image I want to see where the drug accumulates in and I find concentration curve for every tissue is showing decreasing shape. But for small intestine and the intracellular space of duodenum, jejunum and ileum, the Cmax are especially high. image image How can drug concentration in every tissue decrease when the clearance is 0? And, I know there is a mechanism of scaling the permeability between intracellular and interstial space in mucosa to avoid drug accumulation in enterocytes, does the figure mean an accumulation in my simulation (I use the specific intestinal permeability and organ permeability calculated automatically )?

tobiasK2001 commented 5 years ago

Hi, It is indeed very hard to explain your results without having a chance to look at the model. Do you think you can upload it? Some general things. Your compound is very lipophil and a week base. Why do you use Poulin Theil distribution? Have you tried fitting Lipophilicity in a range of 4-6 with variation of distribution model calculation as a fist step?

Best, Tobias

wangwei1619 commented 5 years ago

Hi, @tobiasK2001 Thanks for your help. Here is my PK-Sim file, the simulation named "Shan" is what I am focusing on now, and only the individual "Shan" is loaded with transporters expression. I forgeted to mention that I changed the distribution model to Rodgers & Rowland. However, with fu is fitted to blood plasma ratio, results form two model seem to be similar. I have tried less Log P, and the result was truly better, but I don't know if there is another way to make it. Vemurafenib randoomly.zip

StephanSchaller commented 5 years ago

Hi @wangwei1619 ,

Besides the problem above, I also find something very weird. Basing on the parameters for transporters, the plasma concentration curve would be like this even if the specific clearance is set to be 0

If I switch off all your clearance processess in the "Shan" simulation, I get this: image

Maybe you changes the clearances in the Building Block, but of course you need to set it to zero in the simulation if you want it to affect the results:

image

StephanSchaller commented 5 years ago

To continue:

for every tissue is showing decreasing shape. But for small intestine and the intracellular space of duodenum, jejunum and ileum, the Cmax are especially high.

Yes, this is expected, as the drug gets firs absorbed in small intestine mucosal tissue. The concentrations (technically) always follow a gradient that decreases from source (i.e. the gut) to sink (i.e. the plasma in case of plasma clearance). Partitioning coefficients may offset this gradient, but not to go into details here.

How can drug concentration in every tissue decrease when the clearance is 0? As pointed aout above, you may have set the clearance --> 0 in the Building Block and not in the Simulation.

And, I know there is a mechanism of scaling the permeability between intracellular and interstial space in mucosa to avoid drug accumulation in enterocytes, does the figure mean an accumulation in my simulation (I use the specific intestinal permeability and organ permeability calculated automatically )? No, this is not indicative of an accumulation... as pointed out above, this is just the first entry-point of the drug into tissue where it still has a high concentration before distributing into the body.

Looking at your file and trying a fit, I didn't see a problem. I used the Zhang model and fitted the following: image

Varying both, calculation methods for partitioning and permeabilities, the curves could be fitted very well using the Zhang model (the one with only few processes)... I didn't check on the feasability of the end results for the parameters: image

Of course you had a lot of processes in the Shan model and you will need to inform them to some extent from literature. Just fitting is not the ideal way to go in PBPK.

Maybe this helps. Stephan

wangwei1619 commented 5 years ago

Dear @StephanSchaller ,

I think I come across some bug. image image Thanks for your help, I will folllow your guide .

wangwei1619 commented 5 years ago

Hi, @StephanSchaller

I feel I am a troublemaker. I thought I came across a bug, so I tried to unload and reload the clearance and transporter processes to solve it, but failed. Then I am trying to re-build this mode from the beginning. However, I can not load the transporters processes any more, like this: image I only define one individual and one compound, and associate the vemurafenib to the transporters, but I transporters are not showed in simulation. Is there anyway to solve this?

msevestre commented 5 years ago

Do you have the transporters in your individual?

wangwei1619 commented 5 years ago

Hi, @msevestre Yes, I firstly load transporters in individual, then I assocoate them to drug

msevestre commented 5 years ago

It does not look like it from the screenshot. Can you please go to the expression tab of your individual and make a screenshot?

wangwei1619 commented 5 years ago

Hi, @msevestre

Yes, I find I have made a mistake by adding transporters to enzymes when I check the individual building block. And from this new work I can get the result what @StephanSchaller showed.

Thanks for your help

prvmalik commented 5 years ago

Consider checking out this recent publication, which features a discussion of PBPK modeling of vemurafenib doi: 10.1002/cpt.1672

leeking-55 commented 4 years ago

Hi, all

I am doing a vemurafenib (Log P 5.1, pka 7.1, MW 500, fu <0.01)simualtion, but I can't reach the Cmax like the observed data. image

The hepatic clearance rate is according to publication value, and the fu is fitted to be 2.70E-4 to get consistent with blood plasma ratio of 0.72. The distribution model is Poulin & Theil. It's an oral formulation, however, choosing dissolved formulation can only make littele difference I have condiser transporters like ABCG2, ABCB1, OATP1B3 and also the EHC. After parameters identification, the parameters for transporters have been set to make the distribution in brain be of small fraction, which is agreement with publication, and biliary clearance get very small so that EHC might get insignificant.

Is there any else factor to better the fitting? Moreover, published data says vemurafenib is mainly excreted unchanged by biliary pathway, especially within first 48 h, and I think it suggest the EHC for vemurafenib is important. Can EHC influence the Cmax?

Any comment is appreciated.

hi,wangwei I have the same problem,have you solved this problem? Can I add your wechat? best wish, leeking

wangwei1619 commented 4 years ago

Hi, Ieeking,

Yes, I have finished this model. You can contact me via my email wangwei1619@yeah.net, and I give you ma wechat there.

Best regards, Wang Wei