1
|
Molins E, Labes D, Schütz H, Cobo E, Ocaña J. An iterative method to protect the type I error rate in bioequivalence studies under two-stage adaptive 2×2 crossover designs. Biom J 2020; 63:122-133. [PMID: 33000873 DOI: 10.1002/bimj.201900388] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2019] [Revised: 04/20/2020] [Accepted: 06/22/2020] [Indexed: 11/10/2022]
Abstract
Bioequivalence studies are the pivotal clinical trials submitted to regulatory agencies to support the marketing applications of generic drug products. Average bioequivalence (ABE) is used to determine whether the mean values for the pharmacokinetic measures determined after administration of the test and reference products are comparable. Two-stage 2×2 crossover adaptive designs (TSDs) are becoming increasingly popular because they allow making assumptions on the clinically meaningful treatment effect and a reliable guess for the unknown within-subject variability. At an interim look, if ABE is not declared with an initial sample size, they allow to increase it depending on the estimated variability and to enroll additional subjects at a second stage, or to stop for futility in case of poor likelihood of bioequivalence. This is crucial because both parameters must clearly be prespecified in protocols, and the strategy agreed with regulatory agencies in advance with emphasis on controlling the overall type I error. We present an iterative method to adjust the significance levels at each stage which preserves the overall type I error for a wide set of scenarios which should include the true unknown variability value. Simulations showed adjusted significance levels higher than 0.0300 in most cases with type I error always below 5%, and with a power of at least 80%. TSDs work particularly well for coefficients of variation below 0.3 which are especially useful due to the balance between the power and the percentage of studies proceeding to stage 2. Our approach might support discussions with regulatory agencies.
Collapse
Affiliation(s)
- Eduard Molins
- Department of Statistics and Operations Research, Universitat Politecnica de Catalunya, Barcelona, Catalunya, Spain
| | | | | | - Erik Cobo
- Department of Statistics and Operations Research, Universitat Politecnica de Catalunya, Barcelona, Catalunya, Spain
| | - Jordi Ocaña
- Department of Genetics, Microbiology and Statistics - Statistics Section, Universitat de Barcelona, Barcelona, Catalunya, Spain
| |
Collapse
|
2
|
Fuglsang A. A Three-Treatment Two-Stage Design for Selection of a Candidate Formulation and Subsequent Demonstration of Bioequivalence. AAPS JOURNAL 2020; 22:109. [PMID: 32803519 DOI: 10.1208/s12248-020-00492-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2020] [Accepted: 07/24/2020] [Indexed: 11/30/2022]
Abstract
This paper introduces a two-stage bioequivalence design involving the selection of one out of two candidate formulations at an initial stage and quantifies the overall power (chance of ultimately showing bioequivalence) in a range of scenarios with CVs ranging from 0.1 to 1. The methods introduced are derivates of the methods invented in 2008 by Diane Potvin and co-workers (Pharm Stat. 7(4): 245-262, 2008). The idea is to test the two candidate formulations independently in an initial stage, making a selection of one of these formulations basis of the observed point estimates, and to run, when necessary, a second stage of the trial with pooling of data. Alpha levels are identified which are shown to control the maximum type I error at 5%. Results, expressed as powers and sample sizes, are also published for scenarios where the two formulations are far apart in terms of the match against the reference (one GMR being 0.80, the other GMR being 0.95) and in scenarios where the two test formulations have an actual better match (one GMR being 0.90, the other GMR being 0.95). The methods seem to be compliant with wording of present guidelines from EMA, FDA, WHO, and Health Canada. Therefore the work presented here may be useful for companies developing drugs whose approval hinges on in vivo proof of bioequivalence and where traditional in vitro screening methods (such as dissolution trials) may have poor ability to predict the best candidate.
Collapse
|
3
|
Kaza M, Sokolovskyi A, Rudzki PJ. 10th Anniversary of a Two-Stage Design in Bioequivalence. Why Has it Still Not Been Implemented? Pharm Res 2020; 37:140. [PMID: 32661944 PMCID: PMC7359142 DOI: 10.1007/s11095-020-02871-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2020] [Accepted: 07/02/2020] [Indexed: 11/05/2022]
Abstract
PURPOSE In 2010 the European Medicines Agency allowed a two-stage design in bioequivalence studies. However, in the public domain there are mainly articles describing the theoretical and statistical base for the application of the two-stage design. One of the reasons seems to be the lack of practical guidance for the Sponsors on when and how the two-stage design can be beneficial in bioequivalence studies. METHODS Different variants with positive and negative outcomes have been evaluated, including a pivotal study, pilot + pivotal study and two-stage study. The scientific perspective on the two-stage bioequivalence study has been confronted with the industrial one. RESULTS Key information needed to conduct a bioequivalence study - such as in vitro data and pharmacokinetics - have been listed and organized into a decision scheme. Advantages and disadvantages of the two-stage design have been summarized. CONCLUSION The use of the two-stage design in bioequivalence studies seems to be a beneficial alternative to the 2 × 2 crossover study. Basic information on the properties of the active substance and the characteristics of the drug form are needed to make an initial decision to carry out the two-stage study.
Collapse
Affiliation(s)
- Michał Kaza
- Pharmacokinetics Department, Łukasiewicz Research Network - Pharmaceutical Research Institute, 8 Rydygiera Str., 01-793, Warsaw, Poland.
| | | | - Piotr J Rudzki
- Pharmacokinetics Department, Łukasiewicz Research Network - Pharmaceutical Research Institute, 8 Rydygiera Str., 01-793, Warsaw, Poland
| |
Collapse
|
4
|
Rasmussen HE, Ma R, Wang JJ. Controlling type 1 error rate for sequential, bioequivalence studies with crossover designs. Pharm Stat 2018; 18:96-105. [PMID: 30370634 DOI: 10.1002/pst.1911] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2017] [Revised: 08/31/2018] [Accepted: 09/22/2018] [Indexed: 11/06/2022]
Abstract
Sample size reestimation in a crossover, bioequivalence study can be a useful adaptive design tool, particularly when the intrasubject variability of the drug formulation under investigation is not well understood. When sample size reestimation is done based on an interim estimate of the intrasubject variability and bioequivalence is tested using the pooled estimate of intrasubject variability, type 1 error inflation will occur. Type 1 error inflation is caused by the pooled estimate being a biased estimator of the intrasubject variability. The type 1 error inflation and bias of the pooled estimator of variability are well characterized in the setting of a two-arm, parallel study. The purpose of this work is to extend this characterization to the setting of a crossover, bioequivalence study with sample size reestimation and to propose an estimator of the intrasubject variability that will prevent type 1 error inflation.
Collapse
Affiliation(s)
| | - Rick Ma
- Amgen, Thousand Oaks, California
| | - Jessie J Wang
- University of North Carolina, Chapel Hill, North Carolina
| |
Collapse
|
5
|
Yan F, Zhu H, Liu J, Jiang L, Huang X. Design and inference for 3-stage bioequivalence testing with serial sampling data. Pharm Stat 2018; 17:458-476. [PMID: 29726096 PMCID: PMC6146059 DOI: 10.1002/pst.1865] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2017] [Revised: 03/21/2018] [Accepted: 03/29/2018] [Indexed: 11/07/2022]
Abstract
A bioequivalence test is to compare bioavailability parameters, such as the maximum observed concentration (Cmax ) or the area under the concentration-time curve, for a test drug and a reference drug. During the planning of a bioequivalence test, it requires an assumption about the variance of Cmax or area under the concentration-time curve for the estimation of sample size. Since the variance is unknown, current 2-stage designs use variance estimated from stage 1 data to determine the sample size for stage 2. However, the estimation of variance with the stage 1 data is unstable and may result in too large or too small sample size for stage 2. This problem is magnified in bioequivalence tests with a serial sampling schedule, by which only one sample is collected from each individual and thus the correct assumption of variance becomes even more difficult. To solve this problem, we propose 3-stage designs. Our designs increase sample sizes over stages gradually, so that extremely large sample sizes will not happen. With one more stage of data, the power is increased. Moreover, the variance estimated using data from both stages 1 and 2 is more stable than that using data from stage 1 only in a 2-stage design. These features of the proposed designs are demonstrated by simulations. Testing significance levels are adjusted to control the overall type I errors at the same level for all the multistage designs.
Collapse
Affiliation(s)
- Fangrong Yan
- Research Center of Biostatistics and Computational Pharmacy, China Pharmaceutical University, Nanjing 210009, P.R. China
| | - Huihong Zhu
- Research Center of Biostatistics and Computational Pharmacy, China Pharmaceutical University, Nanjing 210009, P.R. China
| | - Junlin Liu
- Research Center of Biostatistics and Computational Pharmacy, China Pharmaceutical University, Nanjing 210009, P.R. China
| | - Liyun Jiang
- Research Center of Biostatistics and Computational Pharmacy, China Pharmaceutical University, Nanjing 210009, P.R. China
| | - Xuelin Huang
- Department of Biostatistics, The University of Texas MD Anderson Cancer Center, Houston, Texas, 77030, USA
| |
Collapse
|
6
|
Knahl SIE, Lang B, Fleischer F, Kieser M. A comparison of group sequential and fixed sample size designs for bioequivalence trials with highly variable drugs. Eur J Clin Pharmacol 2018; 74:549-559. [PMID: 29362819 DOI: 10.1007/s00228-018-2415-7] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2017] [Accepted: 01/09/2018] [Indexed: 10/18/2022]
Abstract
PURPOSE A drug is defined as highly variable if its intra-individual coefficient of variation (CV) is greater than or equal to 30%. In such a case, bioequivalence may be assessed by means of methods that take the (high) variability into account. The Scaled Average Bioequivalence (SABE) approach is such a procedure and represents the recommendations of FDA. The aim of this investigation is to compare the performance characteristics of classical group sequential designs (GSD) and fixed design settings for three-period crossover bioequivalence studies with highly variable drugs, where the SABE procedure is utilized. METHODS Monte Carlo simulations were performed to assess type I error rate, power, and average sample size for GSDs with Pocock's and O'Brien-Fleming's stopping rules and various timings of the interim analysis and for fixed design settings. RESULTS Based on our investigated scenarios, the GSDs show comparable properties with regard to power and type I error rate as compared to the corresponding fixed designs. However, due to an advantage in average sample size, the most appealing design is Pocock's approach with interim analysis after 50% information fraction. CONCLUSIONS Due to their favorable performance characteristics, two-stage GSDs are an appealing alternative to fixed sample designs when assessing bioequivalence in highly variable drugs.
Collapse
Affiliation(s)
- Sophie I E Knahl
- Boehringer Ingelheim Pharma GmbH & Co. KG, Birkendorfer Straße 65, 88397, Biberach, Germany
| | - Benjamin Lang
- Boehringer Ingelheim Pharma GmbH & Co. KG, Birkendorfer Straße 65, 88397, Biberach, Germany
| | - Frank Fleischer
- Boehringer Ingelheim Pharma GmbH & Co. KG, Birkendorfer Straße 65, 88397, Biberach, Germany
| | - Meinhard Kieser
- Institute of Medical Biometry and Informatics, University of Heidelberg, Im Neuenheimer Feld 130.3, 69120, Heidelberg, Germany.
| |
Collapse
|
7
|
Molins E, Cobo E, Ocaña J. Two-stage designs versus European scaled average designs in bioequivalence studies for highly variable drugs: Which to choose? Stat Med 2017; 36:4777-4788. [PMID: 28853164 DOI: 10.1002/sim.7452] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2016] [Revised: 07/22/2017] [Accepted: 08/07/2017] [Indexed: 11/06/2022]
Abstract
The usual approach to determine bioequivalence for highly variable drugs is scaled average bioequivalence, which is based on expanding the limits as a function of the within-subject variability in the reference formulation. This requires separately estimating this variability and thus using replicated or semireplicated crossover designs. On the other hand, regulations also allow using common 2 × 2 crossover designs based on two-stage adaptive approaches with sample size reestimation at an interim analysis. The choice between scaled or two-stage designs is crucial and must be fully described in the protocol. Using Monte Carlo simulations, we show that both methodologies achieve comparable statistical power, though the scaled method usually requires less sample size, but at the expense of each subject being exposed more times to the treatments. With an adequate initial sample size (not too low, eg, 24 subjects), two-stage methods are a flexible and efficient option to consider: They have enough power (eg, 80%) at the first stage for non-highly variable drugs, and, if otherwise, they provide the opportunity to step up to a second stage that includes additional subjects.
Collapse
Affiliation(s)
- Eduard Molins
- Department of Statistics and Operations Research, Universitat Politècnica de Catalunya, Barcelona, Spain
| | - Erik Cobo
- Department of Statistics and Operations Research, Universitat Politècnica de Catalunya, Barcelona, Spain
| | - Jordi Ocaña
- Department of Genetics, Microbiology and Statistics, Universitat de Barcelona, Barcelona, Spain
| |
Collapse
|
8
|
Inflation of Type I Error in the Evaluation of Scaled Average Bioequivalence, and a Method for its Control. Pharm Res 2016; 33:2805-14. [DOI: 10.1007/s11095-016-2006-1] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2016] [Accepted: 07/19/2016] [Indexed: 11/26/2022]
|
9
|
Queckenberg C, Erlinghagen V, Baken BCM, Van Os SHG, Wargenau M, Kubeš V, Peroutka R, Novotný V, Fuhr U. Pharmacokinetics and pharmacogenetics of capecitabine and its metabolites following replicate administration of two 500 mg tablet formulations. Cancer Chemother Pharmacol 2015; 76:1081-91. [PMID: 26242222 DOI: 10.1007/s00280-015-2840-6] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2015] [Accepted: 07/27/2015] [Indexed: 12/27/2022]
Abstract
PURPOSE To describe concentration versus time profiles of capecitabine and its metabolites 5'-DFUR, 5'-DFCR and 5-FU, depending on tablet formulation and on frequent and/or relevant genetic polymorphisms of cytidine deaminase, dihydropyrimidine dehydrogenase, thymidylate synthase and methylenetetrahydrofolate reductase (MTHFR). METHODS In 46 cancer patients on chronic capecitabine treatment, who voluntarily participated in the study, individual therapeutic doses were replaced on four consecutive mornings by the study medication. The appropriate number of 500 mg test (T) or reference (R) capecitabine tablets was given in randomly allocated sequences TRTR or RTRT (replicate design). Average bioavailability was assessed by ANOVA. RESULTS Thirty female and 16 male patients suffering from gastrointestinal or breast cancer (mean age 53.4 years; mean dose 1739 mg) were included. The T/R ratios for AUC0-t(last) and C max were 96.7 % (98 % CI 90.7-103.2 %) and 87.2 % (98 % CI 74.9-101.5 %), respectively. Within-subject variability for AUC0-t(last) and C max (coefficient of variation for R) was 16.5 and 30.2 %, respectively. Similar results were seen for all metabolites. No serious adverse events occurred. For the MTHFR C677T (rs1801133) genotype, an increasing number of 677C alleles showed borderline correlation with an increasing elimination half-life of capecitabine (p = 0.043). CONCLUSIONS The extent of absorption was similar for T and R, but the rate of absorption was slightly lower for T. While such differences are not considered as clinically relevant, formal bioequivalence criteria were missed. A possible, probably indirect role of the MTHFR genotype in pharmacokinetics of capecitabine and/or 5-FU should be investigated in further studies.
Collapse
Affiliation(s)
- Christian Queckenberg
- Department of Pharmacology, University of Cologne, Cologne, Germany. .,Clinical Trials Centre Cologne, Medical Faculty, University of Cologne, Gleueler Str. 269, 50935, Cologne, Germany.
| | - V Erlinghagen
- Department of Pharmacology, University of Cologne, Cologne, Germany
| | | | | | - M Wargenau
- M.A.R.C.O. GmbH & Co. KG, Düsseldorf, Germany
| | - V Kubeš
- Quinta-Analytica S.r.o., Prague, Czech Republic
| | - R Peroutka
- Quinta-Analytica S.r.o., Prague, Czech Republic
| | - V Novotný
- Quinta-Analytica S.r.o., Prague, Czech Republic
| | - U Fuhr
- Department of Pharmacology, University of Cologne, Cologne, Germany
| |
Collapse
|