Pharmacokinetic Simulations Involving Convolution Approaches
Total Page:16
File Type:pdf, Size:1020Kb
PPhhaarrmmaaccookkiinneettiicc SSiimmuullaattiioonnss IInnvvoollvviinngg CCoonnvvoolluuttiioonn AApppprrooaacchheess By Abdul Hakim Abdullah Ahmed Khaled A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Mathematics May 2014 Department of Mathematics The Islamia University of Bahawalpur Bahawalpur 63100, Pakistan Pharmacokinetic Simulations Involving Convolution Approaches By Abdul Hakim Abdullah Ahmed khaled A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Mathematics May 2014 Supervised by Dr. Khalid Pervaiz Akhter Dr. Ghulam Murtaza Department of Mathematics The Islamia University of Bahawalpur, Pakistan Bahawalpur 63100, Pakistan In The Name Of Allah, The Most Beneficent, The Most Merciful Dedicated To My respectful mother, my lovely wife, sons, and my daughter, who always pray, love, support and encourage me. I II Approval It is hereby certified that the work presented by Abdul Hakim Abdullah Ahmed Khaled S/O Abdullah Ahmed in the dissertation entitled “Pharmacokinetic Simulations Involving Convolution Approaches” has been successfully carried out under our supervision in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Mathematics, at the Department of Mathematics, Faculty of Science, The Islamia University of Bahawalpur, Pakistan. [email protected] Mobile: +92-3468783832 [email protected] Mobile: + 92-3142082826; fax: + 92-62925556 III IV Acknowledgement All praises and gratitude to my Lord Allah, the most gracious, the most merciful for guiding me out of darkness, helping in difficulties and bestowing upon me with the courage to accomplish this dissertation. All respect to the holy prophet Muhammad peace be upon him for showing the right path to humanity. I wish to express my profound gratitude to my supervisors, Dr. Khalid Pervaiz Akhter, Associate Professor, and Dr. Ghulam Murtaza, Assistant Professor, for their professional guidance, help, and scientific advice. I would like to express my sincere thanks to Prof. Dr. Tahir Mahmood, Chairman Department of Mathematics, The Islamia University of Bahawalpur, whose guidance, advice, and great help during the different phases in my study. Hearty appreciations to my teachers, and class fellows, I would like to thank my dearest friends, Dr. Abdul Gaffar, Dr. Aqeel Khan, Saddam Al-Adwar, and Amir Ahmed who never failed in helping me whenever I needed their assistance. My last few words devoted to my family: Mother, brothers and sisters, to my wife I would like to say that I owe my success to her moral support. My deepest gratitude and thanks to the two expert evaluators: Prof. Dr. Muhammad Saleem, Department of Mathematics, San Jose State University, U.S.A. and Prof. Dr. Doron Levy, Department of Mathematics, University of Maryland, U.S.A. who have accepted to evaluate my thesis. I am very thankful to external examiners Prof. Dr. Muhammad Ozair Ahmed, Chairman Department of Mathematic, UET Lahore, and Dr. Muhammad Nawaz Naeem, Associate Professor, Department of Mathematics, the University of Agriculture, Faisalabad , all of your notes and highlights will be really helpful and worth to be added. Finally, thanks to the Amran University which sent me to Pakistan in order to study for a Ph. D. the peace and the safe to my beloved homeland Yemen. V List of symbols and abbreviations PK Pharmacokinetic ADME Absorption, Distribution, Metabolism, Excretion BCS Bio-pharmaceutics classification system FDA Food and drug administration USP United States Pharmacopeia AUC Area under the concentration time curve 퐴푈퐶0−∞ Total area under the curve 퐴푈퐶0−푡 The area under the plasma drug concentration-time curve from zero to t 퐶푚푎푥 Maximum concentration IVIVC In vitro in vivo correlation IVIVP In vitro to in vivo profiling ∗ Symbol of convolution ∕∕ Symbol of deconvolution 푘푎 First – order absorption rate constant 푘푒푙 First – order elimination rate constant 푘12 Distribution rate constant for transfer of drug from compartment one to two 푘21 Distribution Rate Constant for transfer of drug from Compartment two to one UIR Unit impulse response (퐶훿) 푇푚푎푥 Time needed to reach maximum blood drug concentration PH Potential of Hydrogen SD Standard deviation VI SPSS Statistical Product and Service Solution ANOVA Analysis of variance IRF Input response function 퐹푇 Fraction of absorption at any time T 푓2 Similarity factor (퐶푃)0 Concentration at time zero 퐾퐻 Higuchi model Γ(. ) Gamma function VII Abstract This thesis is a part of a research project on pharmacokinetic modeling of enteric coated microparticulate formulations of Metoprolol tartrate. First part of this study dealt with pharmaceutical aspects i.e. formulation development, while second part of study dealt with mathematical modeling of pharmacokinetics. Firstly, this study was aimed to develop in vitro in vivo correlation (IVIVC) level A, B and C for encapsulated Metoprolol tartrate (T1, T2 and T3 having Metoprolol tartrate/polymer ratio of 1:1, 1:1.5 and 1:2 by weight/weight). The in vitro data was correlated with in vivo data. For IVIVC level A, drug absorption data was calculated using Wagner-Nelson method. In addition, convolution approach was used to approximate plasma drug levels from in vitro dissolution data. The coefficient of determination (R2) for level A was 0.720, 0.905, 0.928 and 0.878 for Mepressor®, T, T2 and T3 formulations, respectively with acceptable percent error (<15%). The value of (R2) for level B and C was 0.231 and 0.714, respectively. It is also concluded that IVIVC level A is a proficient mathematical model for biowaiver studies involving study parameters as those implemented for T1S (T1 formulation tested for dissolution in the presence of sodium lauryl sulphate) revealing that IVIVC level A is dosage form specific, rather than to be drug specific. Secondly, the aimed of this study was to assess and apply the in vitro to in vivo profiling (IVIVP), a new biowaiver approach, in designing a product with specific release pattern. The IVIVC was established by plotting the observed and predicted plasma drug concentrations. For IVIVC, convolution approach was employed to estimate plasma drug concentrations from in vitro dissolution profiles. The IVIVC for T1S exhibited a good correlation coefficient (R2 = 0.963) followed by the T2 (R2 = 0.682), T3 (R2 = 0.665), T1 (R2 = 0.616), and Mepressor® (R2 = 0.345). Establishing an IVIVP, based on VIII the convolution approach, can be more useful and practicable in the biowaiver studies, rather than present complicated practice of IVIVC estimated via deconvolution approach. This study also elaborates that there is good correlation between the IVIV profiles of the developed Metoprolol tartrate formulations, particularly for T1S. Keywords: Metoprolol tartrate, Eudragit® FS, Microparticles, Convolution, Deconvolution, IVIVC, IVIVP. IX List of content Serial No. Title Page No. Title page Bismillah Dedication I Declaration II Approval III Certificate V Acknowledgement IV List of abbreviations VI Abstracts VIII List of contents X List of tables XII List of figures XII Chapter 1 Introduction 1 1.1 Introduction 2 1.2 Basic definitions 5 1.3 Convolution /Deconvolution 10 1.4 Pharmacokinetic 18 1.4.1 Absorption 21 1.4.2 Distribution 22 1.4.3 Metabolism 22 1.4.4 Excretion 23 1.5 Simulation 23 1.6 Administration drugs and pharmacokinetic process 24 1.7 Classifications of drugs 25 1.7.1 BCS class I drugs 25 X 1.8 Metoprolol 26 Chapter 2 Literature review 27 2.1 Introduction 28 2.2 Pharmacokinetics 29 2.2.1 Pharmacokinetic modeling 32 2.2.2 Compartment models 37 2.2.2.1 One - compartment model 40 2.2.2.1.1 One - compartment model intravascular administration 40 2.2.2.1.2 One – compartment model extravascular administration 42 2.2.2.1.3 Application of one compartment model (Wagner – Nelson method) 49 2.2.3 Two – compartment model 51 2.2.3.1 Analytical solution for two compartments 54 2.2.3.2 Application of two - compartment model (Loo-Riegelman method) 57 2.2.4 푛- compartment Analysis 68 Chapter 3 Convolution/Deconvution and their applications 63 3.1 Introduction 64 3.2 Linear system analysis 65 3.2.1 Definition of convolution 66 3.2.1.1 Properties of convolution 67 3.2.1.2 Definition of Dirac Delta function 68 3.2.1.3 Properties of the Dirac Delta function 68 3.2.1.4 Properties of convolution of derivative 69 3.2.1.5 Convolution of Properties of integration 70 3.3.2 Definition Delta function (Kronecker function) 70 3.3.2.1 Properties of Delta function 훿 70 3.4 Convolution modeling and application 71 3.4.3 Plasma input by using convolution 76 3.6 Deconvolution modeling 77 3.6.1 Methods of deconvolution 80 XI Chapter 4 Materials and methods 90 4.1 Preparation of investigational tablets 91 4.1.1 Compatibility analysis 92 4.2 Drug release kinetics 94 4.3 Convolution of in vitro dissolution data to approximate plasma 96 4.4 Experimental protocols for in vivo study 103 4.5 Deconvolution of in vivo data IVIVC development 103 4.5.1 Computation of absorption data and IVIVC development 106 4.6 Mathematical and statistical analysis 106 Chapter 5 Results and discussion 107 Results and discussion 108 5.1 Compatibility analysis 108 5.2 Drug Release kinetics 110 5.3 In vivo study 112 5.4 Development of in vitro in vivo correlation 113 5.5 Development of in vitro in vivo profiling 120 Conclusion 122 References 125 Original publications from thesis 138 XII List of tables Serial No. Title Page No. 3.1 Convolution integral for different functions 88 4.1 Convolution of dissolution data Mepressor® of formulation 99 4.2 Convolution of dissolution data of formulation T1 100 4.3 Convolution of dissolution data of formulation TIS 101 4.4 Convolution of dissolution data of formulation T2 102 u 4.5 Convolution of dissolution data of formulation T3 103 5.1 Pharmacokinetic parameters for all tablets obtained from in vivo 115 experiments and convolution method 5.2 In vitro in vivo correlation data 120 5.3 In vitro to in vivo profiling data 120 List of figures Serial No.