PEER REVIEWED Evaluation of Validation Dissolution Test Results Using Lower Probability Bound Distribution Charts Pramote Cholayudth INTRODUCTION The concept of using probability charts in the article entitled “Evaluation of Validation Content Uniformity Test Results Using Probability Distribution Charts,” published in the Journal of Validation Technology (JVT) (1), is applicable by its title to the content uniformity test—one of the critical quality attributes (CQAs) in process validation of oral solid dosage forms (e.g., tablets and capsules). Another important CQA in such a validation exercise is the dissolution test. The approach to compute the probability of passing the test was introduced, in addition to computing the probability for content uniformity test, by J.S. Bergum in 1990 (2) using the SAS program and developed into a Microsoft Excel (MS Excel) program by the author in 2006 (3). To provide an alternative to computation of the probability of passing the dissolution test by the MS Excel program, this article introduces an approach of using the lower probability bound distribution charts, constructed using Bergum’s method, for evaluation of validation dissolution test results. The approach is similar to that described in the previously-mentioned JVT article. USP TEST ACCEPTANCE CRITERIA The USP dissolution test is summarized in Table I. [ For more Author information, go to ivthome.com 72 JOURNAL OF PROBABILITY OF PASSING THE DISSOLUTION TEST Based on Bergum’s method described in “Using the Bergum Method and MS Excel to Determine the Probability of Passing the USP Dissolution Test” (3), the probability formulae are illustrated in Tables II and III. Readers can either copy the formulae in the tables to construct a two-sheet MS Excel file for computing the probability or download the MS Excel file from www.ivthome.com/JVTSummer. The computed probability values (in Table II) are the lower bounds, because the upper bound (UB) for the standard deviation (SD) and the lower bound (LB) for the mean are used for the computation. Using the MS Excel program, the lower probability bound distribution charts for various sample means and standard deviations (SDs) can be constructed, and their example is illustrated in Figure 1. Figure 2 is identical to Figure 1, except that the difference values between the sample means and the Q value are plotted instead of the mean values (the formula is derived in Appendix 1). Now the new chart can be applied to any Q value, which becomes the advantage for construction of the probability charts with “Mean-Q” values on the X-axis in this article. From the two figures, all the curves are overlapping at the 100% probability level. To reduce the difficulty of identifying the probability results at about 100%, another chart is constructed for initial evaluation, as illustrated in Figure 3 (reproduced from page 4 of Appendix 2). ABOUT THE AUTHOR Pramote Cholayudth is a validation consultant to Biolab in Thailand and executive director of Valitech, a wellestablished GMP and validation service firm to the pharmaceutical industry. He is a guest speaker on process validation to the industry organized by local FDA. He can be contacted by e-mail at [email protected]. VALIDATION TECHNOLOGY [SUMMER 2008] ivthome.com P R A M O T E C H O L AY U D T H Table I: Dissolution test acceptance criteria. Stage # Test Sample Size Acceptance Criteria Stage 1 (S1) Test 6 units. Pass if each unit is not less than Q+5%. Stage 2 (S2) Test 6 additional units. Pass if mean of 12 units (S1+S2) is not less than Q% and no unit is less than Q-15%. Stage 3 (S3) Test 12 additional units. Pass if mean of 24 units (S1+S2+S3) is not less than Q% and not more than two units are less than Q-15% and no unit is less than Q-25%. Where Q is the amount of dissolved active ingredient specified in the individual monograph, expressed as a percentage of the labeled content (LC) of the dosage unit, e.g., Q = 80% LC. Table II: Microsoft Excel formula sheet1. Excel Sheet1 A B C 80 1 Q Value (% LC) 80 2 Joint Confidence Level (%) 95 3 Sample Size (n) 20 4 Sample Mean (% LC) 84.00 5 Sample SD (% LC) Lower (0.95^0.5)x100 % Conf. Limit – Mean 3.00 6 7 Upper (0.95^0.5)x100 % Conf. Limit - SD 8 Probability of Passing USP Test 9 Sampling Plan # D 82.09 86.27 4.38 4.33 99.03% 100.00% Plan 1 Plan 2 E For sampling plan 1, follow column B. B1: Q Value (entry data) B2: Joint Confidence Level (entry data) = (0.95^0.5)x(0.95^0.5) = 0.95 = 95% B3: Sample Size (entry data) B4: Sample Mean (entry data) B5: Sample SD (entry data) B6 =B4-B7*NORMSINV(1-(1-(B2/100)^0.5))/B3^0.5 B7 =B5*((B3-1)/(CHIINV((B2/100)^0.5,B3-1)))^0.5 B8 =MAX((1-NORMSDIST(((B1+5)-B6)/B7))^6,(1-NORMSDIST(((B1-15)-B6)/B7))^12+(1-NORMSDIST((12^0.5)*(B1-B6)/B7))-1,(1NORMS DIST(((B1-15)-B6)/B7))^24+COMBIN(24,1)*(1-NORMSDIST(((B1-15)-B6)/B7))^23*(NORMSDIST(((B1-15)-B6)/B7)NORMSDIST(((B1-25)-B6)/B7))+COMBIN(24,2)*(1-NORMSDIST(((B1-15)-B6)/B7))^22*(NORMSDIST(((B1-15)-B6)/B7)-NORMSDIST(((B1-25)-B6)/B7))^2+(1-NORMSDIST((24^0.5)*(B1-B6)/B7))-1) (Reproduced from references [3, 4]) For sampling plan 2, follow column C and Table III. Plan 2 is, for example, 6x3 (6 units/locationx3 locations). C6 =Sheet2!D9-NORMSINV(1-(1-0.95^(1/3))/2)*Sheet2!D12 C7 =Sheet2!D13 The formula for C8 is the same as B8, replace B1, B6 & B7 with C1, C6 & C7, respectively. JOURNAL OF VALIDATION TECHNOLOGY [SUMMER 2008] 73 PEER REVIEWED Table III: Microsoft Excel formula sheet2. Excel Sheet2 A B 1 C D 6 units/location x 3 locations E F n 6 G H Unit # 2 Loc # 1 Loc # 2 Loc # 3 L 3 nxL 18 3 1 100.30 99.00 103.80 4 2 100.90 101.40 100.50 5 3 101.30 100.60 97.50 6 4 101.90 103.20 100.60 7 5 101.60 100.80 99.00 8 6 100.90 103.20 102.90 9 Sample Mean 101.08 10 Within-location SD (SE) 1.68 11 Between-location SD (SM) 0.81 12 95% Upper Bound – SM 6.20 13 95% Upper Bound – SD 4.33 B3-D8 is the entry data area. D9 =AVERAGE(B3:D8) D10 =((SUM(DEVSQ(B3:B8),DEVSQ(C3:C8),DEVSQ(D3:D8)))/((F1-1)*F2))^0.5 D11 =((DEVSQ(B3:D8)-(SUM(DEVSQ(B3:B8),DEVSQ(C3:C8),DEVSQ(D3:D8))))/(F2-1))^0.5 D12 =(D11^2*(F2-1)/CHIINV(0.95^(1/3),(F2-1)))^0.5 (linked to Sheet1) D13 =(IF((D10*((F1-1)*F2)/CHIINV(0.95^(1/3),((F1-1)*F2)))>=(D11*(F2-1)/CHIINV(0.95^(1/3),(F2-1))),(D10*((F1-1)*F2)/CHIINV(0.95^(1/3),((F11)*F2))),(2/3)*(D10*((F1-1)*F2)/CHIINV(0.95^(1/3),((F1-1)*F2)))+(1/3)*(D11*(F2-1)/CHIINV(0.95^(1/3),(F2-1)))))^0.5 (linked to Sheet1) 74 JOURNAL OF VALIDATION TECHNOLOGY [SUMMER 2008] ivthome.com P R A M O T E C H O L AY U D T H Figure 1: Probability distributions for dissolution test (with Q value). Figure 2: Probability distributions for dissolution test (without Q value). JOURNAL OF VALIDATION TECHNOLOGY [SUMMER 2008] 75 PEER REVIEWED Figure 3: Probability distributions at 100% level (Page 4 of Appendix 2). About 2.7 Figure 4: Probability distributions (Page 6 of Appendix 2). Approx. 100% SD= 1.4 76 JOURNAL OF VALIDATION TECHNOLOGY [SUMMER 2008] ivthome.com P R A M O T E C H O L AY U D T H HOW TO USE THE PROBABILITY CHARTS Several probability charts are provided in Appendix 2. The table of contents for Appendix 2 assists practitioners in determining the probability result. Once the validation dissolution test results (sampling plan 1: meanQ, SD, & n and sampling plan 2: mean-Q, SE, SM, & n) are obtained, one may determine if the probability result is 100% using the initial evaluation charts in the early pages of the Appendix 2. If the actual mean-Q result is greater than the determined result, the answer is 100% probability. If the actual mean-Q result is less than the determined result, continue with determining on the probability charts in the following pages of Appendix 2. For example, if one has a dissolution result – meanQ = 2.4 and SD = 1.4 for sampling plan 1 data of n = 1x20 = 20 (5), the initial evaluation using Figure 3 (page 4 of Appendix 2) shows that at the point SD = 1.4, the mean-Q value is about 2.7. This implies that the dissolution results with SD = 1.4 and mean-Q at 2.7 or greater will generate the probability result exactly at 100% (i.e., determining on the other charts is not necessary). In the given example the mean-Q is 2.4; therefore, one has to continue to determine on Figure 4 (reproduced from page 6 of Appendix 2) and the probability result is found to be approximately 100%. Page 6 may be printed out, marked as illustrated in Figure 4, and attached to the validation report as supporting documentation. Sampling plan 2 data are required to compute the sample mean, SE, and SM using the formulae in Table III and use the three computed parameters (including mean-Q) for determination of the probability result. is accomplished. Please note that it is a requirment of the US Food and Drug Administration that any software that is used to determine quality attributes of a medical device or drug, even off-the-shelf software such as MS Excel and SAS, must be validated for its intended use prior to implementation. In other words, the company that chooses to use the calculations in this paper within an MS Excel program must independently verify the inputs and outputs, and document this activity, prior to implementation. REFERENCES 1. Cholayudth, P., “Evaluation of Validation Content Uniformity Test Results Using Probability Distribution Charts,” Journal of Validation Technology, Volume 14, No. 3, Spring 2008. 2. J. S. Bergum, “Constructing Acceptance Limits for Multiple Stage Tests,” Drug Development and Industrial Pharmacy, 16 (14), Marcel Dekker, Inc, 1990, pp. 2153-2166. 3. Cholayudth, P., “Using the Bergum Method and MS Excel to Determine the Probability of Passing the USP Dissolution Test,” Pharmaceutical Technology, Volume 30, Number 1, January 2006, www.pharmtech.com. 4. Cholayudth, P., “Application of Probability of Passing Multiple Stage Tests in Benchmarking and Validation of Processes” Journal of Validation Technology, Volume 13, No. 4, August 2007. 5. Pluta, P., “Keep in Step with Confidence; Keep in Step through Benchmarking” Journal of Validation Technology, Volume 12, No. 1, November 2005. 6. “Chapter <711> Dissolution,” in United States Pharmacopeia 30–National Formulary 25 (US Pharmacopeial Convention, Rockville, MD, 2007). JVT CONCLUSION The concept of using probability charts will help validation practitioners and QA/QC personnel involved visualize the integrated views of the lower probability bound distributions, and also understand the influence of sample mean location and variability (SD, SE, and SM) on the probability of passing the USP dissolution test. To provide a complete validation documentation, a completed (drawing & notes are taken) copy of the page (in Appendix 2) with the probability chart matching the validation dissolution result should be prepared and attached. To use the probability charts accurately one should realize which sampling plan (1 or 2) is being used and then use the correct chart. Sampling plan 1 (e.g., 1x20 or 1x30) is recommended so that a more simple evaluation of the validation data ARTICLE ACRONYM LISTING CQAs LB LC MS Q SD SE Sigma SM UB USP Critical Quality Attributes Lower Bound Label Claim Microsoft Q value (The amount of dissolved active ingredient specified in the individual monograph, expressed as a percentage of the labeled content of the dosage unit) [6] Standard Deviation (for a Sample) Square Root of Mean Square Error Population Standard Deviation Standard Deviation of Location Means Upper Bound United States Pharmacopeia JOURNAL OF VALIDATION TECHNOLOGY [SUMMER 2008] 77 PEER REVIEWED APPENDIX 1: EXCEL FORMULAE Table IV: Probability of passing dissolution test (with Q value). Probability as a function of Mu and Sigma =MAX((1-NORMSDIST(((Q+5)-Mu)/Sigma))^6,(1-NORMSDIST(((Q-15)-Mu)/Sigma))^12+(1-NORMSDIST((12^0.5)*(QMu)/Sigma))-1,(1-NORMSDIST(((Q-15)-Mu)/Sigma))^24+COMBIN(24,1)*(1-NORMSDIST(((Q-15)Mu)/Sigma))^23*(NORMSDIST(((Q-15)-Mu)/Sigma)-NORMSDIST(((Q-25)-Mu)/Sigma))+COMBIN(24,2)*(1-NORMSDIST(((Q-15)-Mu)/Sigma))^22*(NORMSDIST(((Q-15)-Mu)/Sigma)-NORMSDIST(((Q-25)-Mu)/Sigma))^2+(1-NORMSDIST ((24^0.5)*(Q-Mu)/Sigma))-1) * * Reproduced from [3, 4], Mu (m) = Lower Bound (LB) for Mean, Sigma (s) = Upper Bound (UB) for standard deviation (SD) Table V: Probability of passing dissolution test (with Mean-Q value). Probability as a function of Mean-Q, Sigma, and n (modified from Table IV) =MAX((1-NORMSDIST((-(Mu-Q)+5)/Sigma))^6,(1-NORMSDIST((-(Mu-Q)-15)/Sigma))^12+(1-NORMSDIST((12^0.5)*((Mu-Q))/Sigma))-1,(1-NORMSDIST((-(Mu-Q)-15)/Sigma))^24+COMBIN(24,1)*(1-NORMSDIST((-(Mu-Q)15)/Sigma))^23*(NORMSDIST((-(Mu-Q)-15)/Sigma)-NORMSDIST((-(Mu-Q)-25)/Sigma))+COMBIN(24,2)*(1-NORMSDIST((-(Mu-Q)-15)/Sigma))^22*(NORMSDIST((-(Mu-Q)-15)/Sigma)-NORMSDIST((-(Mu-Q)-25)/Sigma))^2+(1-NORMSDI ST((24^0.5)*(-(Mu-Q))/Sigma))-1) =MAX((1-NORMSDIST((-(Mean-Sigma*Z/n^0.5-Q)+5)/Sigma))^6,(1-NORMSDIST((-(Mean-Sigma*Z/n^0.5-Q)15)/Sigma))^12+(1-NORMSDIST((12^0.5)*(-(Mean-Sigma*Z/n^0.5-Q))/Sigma))-1,(1-NORMSDIST((-(MeanSigma*Z/n^0.5-Q)-15)/Sigma))^24+COMBIN(24,1)*(1-NORMSDIST((-(Mean-Sigma*Z/n^0.5-Q)15)/Sigma))^23*(NORMSDIST((-(Mean-Sigma*Z/n^0.5-Q)-15)/Sigma)-NORMSDIST((-(Mean-Sigma*Z/n^0.5-Q)-25)/ Sigma))+COMBIN(24,2)*(1-NORMSDIST((-(Mean-Sigma*Z/n^0.5-Q)-15)/Sigma))^22*(NORMSDIST((-(MeanSigma*Z/n^0.5-Q)-15)/Sigma)-NORMSDIST((-(Mean-Sigma*Z/n^0.5-Q)-25)/Sigma))^2+(1-NORMSDIST((24^0.5)*((Mean-Sigma*Z/n^0.5-Q))/Sigma))-1) =MAX((1-NORMSDIST((-(Mean-Q)-Sigma*Z/n^0.5+5)/Sigma))^6,(1-NORMSDIST((-(Mean-Q)-Sigma*Z/n^0.515)/Sigma))^12+(1-NORMSDIST((12^0.5)*(-(Mean-Q)-Sigma*Z/n^0.5)/Sigma))-1,(1-NORMSDIST((-(Mean-Q)Sigma*Z/n^0.5-15)/Sigma))^24+COMBIN(24,1)*(1-NORMSDIST((-(Mean-Q)-Sigma*Z/n^0.5-15)/Sigma))^23*(NORMSDIST((-(Mean-Q)-Sigma*Z/n^0.5-15)/Sigma)-NORMSDIST((-(Mean-Q)-Sigma*Z/n^0.5-25)/Sigma))+COMBIN(24,2)*(1NORMSDIST((-(Mean-Q)-Sigma*Z/n^0.5-15)/Sigma))^22*(NORMSDIST((-(Mean-Q)-Sigma*Z/n^0.5-15)/Sigma)NORMSDIST((-(Mean-Q)-Sigma*Z/n^0.5-25)/Sigma))^2+(1-NORMSDIST((24^0.5)*(-(Mean-Q)Sigma*Z/n^0.5)/Sigma))-1) Mu = Lower Bound for Mean = (Sample) Mean-Sigma*NORMSINV(1-(1-(0.95)^0.5))/n^0.5 = Mean-Sigma*Z/n^0.5, Sigma = Upper Bound for SD, n = Sample Size APPENDIX 2: PROBABILITY CHARTS Appendix 2 is available for download at www.ivthome.com/JVTSummer. Probability charts provided in Appendix 2 are presented in order of sampling plan types and overall sample sizes. All the charts are constructed on 95% joint confidence level. The charts on pages 4 and 5 are recommended for initial evaluation and followed by the following pages (as appropriate) as described in the article. For test results outside the provided charts, please contact [email protected]. 78 JOURNAL OF VALIDATION TECHNOLOGY [SUMMER 2008] ivthome.com
© Copyright 2024 Paperzz