Next Article in Journal
Ampholytic and Polyelectrolytic Starch as Matrices for Controlled Drug Delivery
Previous Article in Journal
Physiologically Relevant In Vitro-In Vivo Correlation (IVIVC) Approach for Sildenafil with Site-Dependent Dissolution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comprehensive Study of Intermediate and Critical Quality Attributes for Process Control of High-Shear Wet Granulation Using Multivariate Analysis and the Quality by Design Approach

1
School of Pharmacy, Sungkyunkwan University, Suwon 440-746, Korea
2
Department of Pharmaceutical Engineering, Inje University, Gyeongnam 621-749, Korea
*
Authors to whom correspondence should be addressed.
Pharmaceutics 2019, 11(6), 252; https://doi.org/10.3390/pharmaceutics11060252
Submission received: 22 April 2019 / Revised: 19 May 2019 / Accepted: 30 May 2019 / Published: 1 June 2019

Abstract

:
A robust manufacturing process and the relationship between intermediate quality attributes (IQAs), critical quality attributes (CQAs), and critical process parameters (CPPs) for high-shear wet granulation was determined in this study. Based on quality by the design (QbD) approach, IQAs, CQAs, and CPPs of a telmisartan tablet prepared by high-shear wet granulation were determined and then analyzed with multivariate analysis (MVA) to evaluate mutual interactions between IQAs, CQAs, and CPPs. The effects of the CPPs on the IQAs and CQAs were quantitatively predicted with empirical models of best fit. The models were used to define operating space, and an evaluation of the risk of uncertainty in model prediction was performed using Monte Carlo simulation. MVA showed that granule size and granule hardness were significantly related to % dissolution. In addition, granule FE (Flow Energy) and Carr’s index had effects on tablet tensile strength. Using the manufacture of a clinical batch and robustness testing, a scale-up from lab to pilot scale was performed using geometric similarity, agitator torque profile, and agitator tip speed. The absolute biases and relative bias percentages of the IQAs and CQAs generated by the lab and pilot scale process exhibited small differences. Therefore, the results suggest that a risk reduction in the manufacturing process can be obtained with integrated process parameters as a result of the QbD approach, and the relationship between IQAs, CQAs, and CPPs can be used to predict CQAs for a control strategy and SUPAC (Scale-Up and Post-Approval Guidance).

1. Introduction

Pharmaceutical manufacturing processes are complex and controlled by multi-factorial relationships between materials and process parameters. Therefore, it is difficult to investigate the pharmaceutical manufacturing process with the OFAT (one-factor-at-a-time) method. Combination approaches, such as experimental design and MVA, could present a comprehensive and intensive understanding of the relevant multi-factorial relationships between process parameters and product quality attributes [1,2]. MVA is the statistical method that can reduce or classify variables using correlation between variables. It also provides visual effects for multivariate data to enable an understanding of the relationship between distribution features, patterns, natural groupings, and relationships variables [1,3,4,5]. In pharmaceutical experimental design, there are numerous variables, such as material properties, manufacturing processes, quality attributes, and intermediate quality attributes, which occur during the manufacturing process. It is difficult to study the relationships among such a large number of variables because only a limited number of variables can be processed through experimental design analysis [1]. Generally, the number of control factor is less than 3–4, and the level of control factors is more than 3 levels in the optimization stage. If the number of control factors increases or the number of level increases, the number of experiment increases exponentially [6]. In experimental design, the relationship between control factors and response factors is statistically analyzed with analysis of variance. However, it is difficult to define the relationship between response factors using experimental design. MVA can be regarded as a complementary method to experimental design, which presents additional information about the product and processes, such as the effect of intermediate quality attributes on product quality attributes. Therefore, when MVA and experimental design are combined, the integrated MVA is a more powerful tool to elucidate the complex relationships in the pharmaceutical manufacturing process.
High-shear wet granulation is the most commonly used technique in the pharmaceutical industry because of various advantages, such as lower use of binding agents, shortened processing times, the formation of reproducible granules with a uniform size distribution, and formation of dense granules suitable for tableting [7]. Agitator speed, chopper speed, agitator shape, chopper shape, liquid to solid ratio, liquid addition rate, massing time, drying temperature, and drying time could be considered process parameters of high-shear wet granulation. The process parameters have an effect on product quality attributes (e.g., dissolution profile, content uniformity, assays, and tablet hardness, among others), as well as intermediate quality attributes (e.g., granule size, granule hardness, granule flow energy, and granule density, among others) [8,9,10,11,12,13]. Improper agitator and chopper speeds can result in uncontrolled granule growth [8,13,14,15,16]. This results in undesired granule porosity and granule strength. In addition, a high agitator speed presents a high shearing force to powder bed, which might produce smaller granules with narrow size distribution [17,18]. Wet massing time might be one of the most important parameters in the granulation process because the parameter allows the coalescence and consolidation of granules by kneading the wet mass [8,19,20]. The wet massing time should be optimized to obtain desired granule properties. The failure of massing time control can produce weak granules or undesired density granules, which makes them difficult to handle for next unit operation. Moreover, the wet massing time may have an effect on granule size and size distribution. Improper wet massing time may produce a bimodal granule size distribution or monomodal granule size distribution [12,21,22]. Additionally, the liquid addition rate is related to the onset of the nucleation stage and operation time [9]. The liquid is uniformly distributed on a powder bed at a slow addition rate. This might decrease the powder agglomeration time. On the contrary to this, a fast addition rate may not present uniform distribution of the liquid on a powder bed, which leads to lump formation with a wide granule size distribution [23,24,25,26]. The improper process can result in undesirable intermediate quality attributes, which could be significantly related to poor product quality attributes.
Based on quality by the design (QbD) implementation and various studies of the QbD approach, in the early stage of development, the formulation and the processes are, in general, optimized using lab-scale equipment. Subsequently, the formulation is fixed during the scale-up [27,28,29,30]. Thus, the formulation should be robust and has to lead to the same quality of product using the lab-scale and pilot or commercial-scale equipment. Therefore, the scale-up process should be significantly controlled to present a target quality for the product. Although there are various scale-up strategies, general considerations in the scale-up process are dimensional analysis and the principle of similarity [28,31,32]. Dimensional analysis is based on the assumption that a mathematical formulation of a physical process should be valid in any system of dimensions. It is a tool to produce dimensionless numbers that completely describe the process. The dimensionless numbers most commonly used to describe the high-shear wet granulation are Newton, Froude, and Reynolds. The Newton number is related to the drag force acting on the agitator and quantified power required to overcome friction in fluid flow in a stirred vessel. The Froude number is the interplay of the centrifugal force and the centripetal force produced by the wall. The Reynolds number is related to the inertial force and the viscous force [28,31,32]. The principle of similarity can be explained with geometrical similarity, kinematic similarity, and dynamic similarity. The dimensionless numbers might have the same numerical value under the condition of similarity. However, the scale-up strategy of high-shear wet granulation is particularly difficult and often problematic because of the inherent heterogeneous nature of the materials used and the various process parameters are complexly interrelated. When a scale-up process of high-shear wet granulation is not well performed, the IQAs might change. A change in the granule size, granule size distribution, granule density, compressibility, or granule flow energy might significantly affect product quality, such as tablet hardness, dissolution, content uniformity, disintegration time, tablet friability, and assays [33]. If the improper scale-up of high-shear wet granulation is not detected at the end product test stage, it could result in huge economic losses and wasted time.
Although the QbD approach is essentially applied for the development of pharmaceutical products, most studies of QbD approach for the development of pharmaceutical products were focused on the relationships between control factors (critical material attributes or CPPs) and response factors (drug product CQAs). In addition, IQAs after unit operation might not be significantly investigated to define the relationship with drug product CQAs in the QbD approach [33,34,35,36,37,38,39]. This study was conducted to investigate the robust manufacturing process with respect to the relationships among IQAs, critical quality attributes (CQAs), and critical process parameters (CPPs) for high-shear wet granulation. Based on the QbD approach, the control and response factors were determined by risk assessment, which is based on our knowledge and experience with related processes and the target profiles of the control tablet. The model drug was telmisartan, and a response surface design was used to obtain the optimal process parameters of high-shear wet granulation to prepare telmisartan granules using three control factors: agitator speed, massing time, and binding suspension addition rate. The quantitative effects of the different control factors on the response factors were accurately predicted using empirical models of best fit. The models and Monte Carlo simulation were used to identify the operating space. The relationships between IQAs and CQAs were analyzed with principal component analysis (PCA), where partial least squares (PLS) was used to evaluate the relationships among IQAs, CQAs, and CPPs. To manufacture a clinical batch and conduct robustness testing, a scale-up from lab to pilot scale was performed using geometric similarity, agitator torque profile, and agitator tip speed. Moreover, a pharmacokinetics study was performed with the clinical batch.

2. Materials and Methods

2.1. Materials

Telmisartan, D-mannitol, meglumine, and magnesium stearate were received from Jeil Pharmaceutical Co., Ltd. (Seoul, Korea). Polyvinylpyrrolidone (PVP K25) and crospovidone (Kollidon®® CL-F) were obtained from BASF AG (Ludwigshafen, Germany). Dicalcium phosphate and sodium hydroxide were obtained from Sigma-Aldrich Co. (St. Louis, MO, USA). Microcrystalline cellulose (MCC 101) was obtained from DFE Pharma (Dusseldorf, Germany). F-MELT®® (Type C) was obtained from Fuji Chemical Industries Co., Ltd. (Toyama, Japan). All other reagents were of analytical or HPLC grade and were used as received.

2.2. Risk Assessment

IQAs, CQAs, and CPPs were established by risk assessment, which was based on our knowledge and experience with the related process and the target profiles of the control tablet. Failure mode and effects analysis (FMEA) was applied to quantify the degree of risk connected with these process parameters [34,36]. Based on the total risk priority number (RPN), a score matrix presented the priority and/or level of attention required to be applied to each step. Risk quantification was performed by ranking severity (S), probability (P), and detectability (D), which were rated individually from 1 to 5. Severity was related to the results of the failure mode and considered the worst result of the failure, determined by the degree of effect to repair the failure. The likelihood of failure was represented as the probability and detectability, which indicated the ability to determine and detect the existence of a failure mode, respectively. The RPN was calculated by multiplying the three scoring columns [34]. RPN was classified as low risk (1–45), medium risk (46–90), and high risk (91–125). Moreover, a pareto chart was combined with the RPN to identify the control factors, which were evaluated using the experimental design.

2.3. Experimental Design to Investigate High-Shear Wet Granulation in Lab Scale

To generate the experimental design, Design-Expert® software (version 10; Stat-Ease Inc., Minneapolis, MN, USA) was used. Response surface design was used to identify the optimal process parameters of high-shear wet granulation to prepare telmisartan granules. The three control factors were agitator speed (p1), massing time (p2), and binding suspension addition rate (p3). IQAs, such as granule hardness (q1), Carr’s index of granules (q2), granule size (q3), granule true density (q4), granule bulk density (q5), and granule FE (q6), and CQAs, such as tablet tensile strength (q7) and % dissolution (q8q11) were evaluated as responses of experimental design, which were determined during the risk assessment. The best-fit empirical model was used to predict the relationship between control factors and response factors by comparing statistical parameters, such as the probability value (p value), the multiple correlation coefficient (R2), and adjusted multiple correlation coefficient (adjusted R2). Contour plots for control factors versus response factors were also presented using the software. To determine the optimized process parameters, the empirical models, which interpreted the relationship between control and response factors, were used. The target response values were the following: maximum tablet tensile strength within the lower (≥ 125 N/cm2) limit; minimum Carr’s index within the upper (≤ 15) limit; optimal intermediate ranges of dissolved drug at 5 min (72.29–83.29%), 10 min (78.68–88.58%), 15 min (85.00–92.10%), and 30 min (87.21–94.47%). In addition, 95% confidence intervals for the optimal condition were used for the control strategy.
To prepare the test tablet, telmisartan granules were generated by high-shear wet granulation, and Table 1 presents the granule components. Solving agents were sodium hydroxide and meglumine. PVP K25 was used as the binding agent. The telmisartan suspension was prepared by the following: i) dissolving solving agents in 50% (v/v) ethanol at ambient temperature, ii) introducing telmisartan and binding agent into the i) solution at ambient temperature, and stirring at 50 rpm until it became transparent, iii) the suspension was then stored until use at ambient temperature. To remove any aggregates, F-MELT®, D-mannitol dicalcium phosphate, crospovidone, and MCC were filtered with a #30-mesh sieve and then blended for 3 min at 200 rpm; the powder mixture from each experiment was 475 g. A high-shear granulator (Mixer-granulator P 1, DIOSNA Dierks & Söhne GmbH, Germany) with a capacity of 1.23 L was used to prepare the granulation. The agitator speed, massing time, and binding suspension addition rate were set to the predetermined values of the experimental design, and the chopper speed was 2500 rpm. The end point of granulation was determined by monitoring agitator torque. The agitator torque profile provides wet granulation process stages, such as stages 1–4 [40]. Stage 1 represents the initiation of wetting and hence the agitator torque remains relatively constant. Stage 2 represents that agglomeration begins and agitator torque increases. Stage 3 involves the generation of useful granules and a plateau in agitator torque. Stage 4 corresponds to over-wetting with fluctuations in agitator torque. The granules were dried in an oven at 60 °C for 1 h. IQAs, such as granule hardness, Carr’s index, granule true density, granule bulk density, granule FE, and granule particle size were evaluated using the dried granules. The dried granules were introduced into a stainless steel double cone blender (Kopamtec, Anyang, Korea) with a capacity of 20,000 cm3 and a fill-weight of approximately 30% of the mixer volume. Magnesium stearate was introduced into the blender and then blended at 50 rpm for 10 min. The resulting granules were compressed on a single-punch hydraulic laboratory press (Ichihashi Seiki Co., Ltd., Kyoto, Japan) at 15 MPa using a rectangle-shaped punch (15 × 8 mm). CQAs were evaluated using the tablets.

2.4. Measurement of Quality Attributes

2.4.1. Measurement of Granule Hardness and Tablet Tensile Strength

The granule strength was measured with a texture analyzer (TA.XT plusC from Stable MicroSystems, Godalming, UK), which has a 500 N force capacity, 0.1 g force resolution, and is equipped with a 50 kg·f loadcell. The dry granules were individually compressed with a 10 mm steel probe. The trigger force was set as 0.001 N. The force versus distance were recorded to calculate the total work (area under the curve). More than 30 granules were evaluated to obtain statistically reliable results. The means of 30 determinations were used to calculate the granule strength for each experiment.
The tablet tensile strength was evaluated by compressing diametrically using a texture analyzer (TA.XT plusC, Stable Micro Systems Ltd., Godalming, UK), and tensile strength was calculated using Fall and Newton’s method [41] as shown in Equation (1):
T e n s i l e   s t r e n g t h = F π R T
where F is the force required to break the tablet (N), R is the radial length of the tablet (mm), and T is the thickness of the tablet (mm). The tablet was compressed with a 20 mm steel probe. The trigger force was set as 0.001 N. The force versus distance were recorded to calculate the maximum force. The means of five determinations were used to calculate the tablet tensile strength for each experiment.

2.4.2. Measurement of Granule Density and Carr’s Index

An MT-1000 instrument (Seishin Enterprise Co., Tokyo, Japan) was used to evaluate the bulk and tap densities of the granules. To determine the bulk density, the granules were passed through a #18 sieve to break up agglomerates that might have formed during storage. An excess of granules was introduced into a 100 mL cylinder until it overflowed and then the top was carefully scraped to remove the excess granules from the top of the cylinder. The granule weight was accurately measured. The tapped density was determined with a 100 mL cylinder. The granules were introduced into the cylinder and then the cylinder was tapped with 250 drops per min with a fixed drop of 14 ± 2 mm. Each analysis was repeated three times. Carr’s index was determined using Equation (2).
Car   r     s   Index = tapped   density bulk   density tapped   density × 100 %
The means of three determinations was used to calculate the granule density and Carr’s index for each experiment.
The true density of the granules was evaluated with a helium pycnometer (AccuPyc 1330; Micromeritics Instrument Co., Norcross, GA, USA). Before measurement, the instrument calibration was performed with a standard steel ball to increase the accuracy of the pycnometer. The granule weight was accurately measured, and it was introduced into the sample cell. The sample cell was filled with helium gas and the gas was discharged into a second empty cell. The granule volume was calculated by measuring pressure by filling the sample cell with helium gas. The means of three determinations were used to calculate the granule true density for each experiment.

2.4.3. Measurement of Particle Size

The laser diffraction method was used to evaluate the particle size distribution of the granules. Approximately, a 5 g granule size was introduced into a HELOS particle size analyzer (Clausthal-Zellerfeld, Germany) ranging from 0.1 µm to 8750 µm. The particle size distribution diameters (d10, d50, and d90) were determined. The measurements were repeated three times and the mean value of d50 was used as particle size in this study.

2.4.4. Powder Property Analysis Using a Rheometer

A Freeman FT4 powder rheometer (Freeman Technology, Micrometitics Instrument Corporation, Norcross UK) was used to analyze granule flow energy (FE). A conditioning cycle was completed prior to every test to remove the variability introduced by the operator during loading of the sample and any residual compaction from previous tests. A blade (stainless steel blade of diameter Db = 23.5 mm) moved downward into and through a powder bed contained within a cylindrical vessel with an inside diameter Dv = 25 mm and volume of 25 mL. The blade speed was 100 mm/s with a −10° helix angle. As the blade moved through the sample, the FT4 measured both rotational and vertical resistances, in the form of torque and force, respectively. It was the composite of these two values that quantified the granule FE. It was calculated with the Equation (3):
granule   FE = 0 L ( T R tan α +   F b a s e ) d H
where R is the agitator blade radius, T is the torque acting on the agitator blade, α is the helix angle, and L is the penetration depth.

2.4.5. In Vitro Dissolution Test

According to the USP Apparatus 2 guidelines (paddle method), the in vitro dissolution test was conducted (DT 720, ERWEKA GmbH, Germany) with 900 mL of dissolution buffer (50 mM phosphate buffer, pH 7.5) maintained at 37 ± 0.5 °C. The paddle speed was maintained with 75 rpm. At predetermined time intervals, sample aliquots (5 mL) were withdrawn with 15 mL syringes. These were filtered through a 0.45 µm membrane syringe filter. Quantitative analysis for drug content was performed using an HPLC system (Agilent Technologies, Santa Clara, CA, USA) with UV detection at a wavelength of 237 nm. A Supelco Discovery C18 column (4.6 × 250 mm, 5 µM) (Sigma-Aldrich, St. Louis, MO, USA) was used and maintained at 40 °C. The mobile phase was a 45:55 volume mixture of aqueous buffer (pH 2.8, 20 mM phosphate buffer prepared with ammonium dihydrogen phosphate and phosphoric acid) and methanol. The injection volume was 20 µL and the flow rate was 1.0 mL/min. The cumulative percentage of the released drug was calculated, and the results are presented as the mean value of at four tablets for each experiment.

2.5. Multivariate Analysis between IQAs, CQAs, and CPPs

The combination of experimental design and MVA has grown to be a more powerful system. The commonly used methods include PCA and PLS. PCA is a statistical method that visualizes the relationship between independent variables (X) by shortening multidimensional data to be analyzed into two or three dimensional data. The PLS is a statistical method for establishing the relationship between the independent variable (X) and the dependent variable (Y) based on the principle of PCA described above. That is, by applying PCA to each of the X and Y variables, T and U, consisting of the score values of X and Y, are projected into a new space to obtain a linear regression model [2,3,4,42,43].
In this study, PCA and PLS were performed with SIMCA© software (Sartorius Stedim Biotech., version 15, Umeå, Sweden) to evaluate the mutual relation between both IQAs of processes and drug product CQAs, and variable relationships in the experimental design of process parameters for high-shear wet granulation to prepare telmisartan granules. Moreover, the best fitted regression models between IQAs and CQAs were constructed with a prediction interval (95% PI), a confidence interval (95% CI), and a p value to be applicable at the pilot scale.

3. Results and Discussion

3.1. Risk Assessment

A risk assessment of the manufacturing process was performed to identify the high risk process parameters that might have significant effects on the CQAs of the drug product. The IQAs of the output material from the manufacturing process steps that impacted the final drug product CQAs were simultaneously identified. The quality target product profile (QTPP) was defined with the analysis of the control tablet and pharmacopeia [33]. The dosage form, dosage design, route of administration, and stability were based on the analysis of the control tablet because generic drug products are asked to have the same or better quality for pharmaceutical equivalence requirements. Therefore, the above product profiles were not established as CQAs in this study.
Generally, the IQAs (Carr’s index, granule size, granule hardness, granule true density, granule bulk density, and granule FE) could potentially be influenced by the process parameters and these might be significantly related to the CQAs. Carr’s index could indicate the flowability of the granule. A poorly flowable granule has a Carr index greater than 25, whereas a Carr’s index below 15 is considered to be an indication of good flowability [44]. It is calculated with bulk density and tapped density, which could be significantly influenced by the process parameters [8,19,34,45]. Based on our experience related with the development of similar processes, the particle size distribution of granules and granule hardness manufactured by high-shear wet granulation could significantly influence drug product quality, such as dissolution, tablet tensile strength, content uniformity, and bioavailability [8,19,46,47]. However, the target values of the IQAs, except Carr’s index, could not be determined before the development of the mutual effect with the CQAs. Therefore, the target values were not defined in the risk assessment at early development stages.
The CQAs (friability, tablet tensile strength, content uniformity (CU), assay, and dissolution) could potentially be influenced by process parameters. Friability might be critically related to the safety and efficacy of the drug product. More than 1.0% (w/w) of weight loss can cause a significant impact on patient safety and efficacy [33]. However, friability of the test tablet is controlled by stricter standards as seen in the target value; minimized to less than 0.25% (w/w). Tablet tensile strength could affect the dissolution profile, pharmacokinetic parameters, and disintegration, which is directly related to the safety and efficacy of the drug product. A strong or weak tablet indicates excessive bonding between API and excipients, which can lead to an improper dissolution profile. The bonding could be influenced by process parameters. Variability in content uniformity or assay is related to the safety and efficacy of the drug product. Process variables could affect the assay or content uniformity of the drug product. The dissolution profile could be directly related to the safety and efficacy of the drug product, which significantly influences bioavailability. Therefore, the above mentioned the IQAs and CQAs should be evaluated, and the mutual interactions between the IQAs and CQAs were analyzed in detail during process parameter development studies.
As shown in Table 2, FMEA was performed to identify the effects of CPPs on CQAs and IQAs. Previous experience with these process steps and information about the control tablet from the published literature [8,9,14,19,40,41,48] or prior analysis were used to determine the degree of risk associated with process parameters and its potential to impact the CQAs and IQAs. For the risk assessment, three categories (high, medium, and low) were mapped based on the degree of RPN. Moreover, a pareto chart was constructed with RPN as shown in Figure 1. The RPN and the cumulative total number were presented in descending order by the bars and line, respectively. The RPN of agitator speed, massing time, and binding suspension addition rate were 125, 125, and 100, respectively. Moreover, the cumulative percent of the three process parameters was 83.7%. Based on the risk assessment, the agitator speed, binding suspension addition rate, and massing time were identified as potentially having a risk impact on CQAs and IQAs and these were evaluated in detail with experimental design.

3.2. Effect of High-Shear Wet Granulation Process Parameters on IQAs and CQAs

A response surface design was used to evaluate the effect of the three control factors on the three response factors. Table 3 presents the experimental matrix for the 15 experiments and experiment results. The mathematical empirical models were generated with coded values for factor levels. A positive coefficient indicated a synergistic effect on the response factor, whereas an antagonistic effect on the response factor was a negative coefficient. In addition, a larger coefficient indicated a stronger effect of the independent variable of the response [49,50,51].
However, CQAs, such as friability, content uniformity, and assay, were excluded to determine operating space because all run orders demonstrated acceptable values, such as content uniformity (ranging from 1.2–2.5% RSD), friability (<0.2% weight loss when compressed using 15 MPa of force), and assay (ranging from 98.6–100.5%). Moreover, none of the process parameters showed a significant impact on these CQAs.

3.2.1. Significant Factors for Granules Hardness (q1)

Granule hardness is one of the most important parameters that can affect the disintegration time, tensile strength, and dissolution rate of tablets. Moreover, granules should have sufficient strength to resist compression. For this reason, granule hardness was evaluated in this study with a texture analyzer. As shown in Table 3, the granule strength manufactured by different process parameters ranged from 0.90–5.00 N. The effect of process parameters on granule hardness was presented as the empirical model described by Equation (4).
q 1 = 2.16 + 1.49 p 1 0.64 p 2 + O .29 p 3 + 0.55 p 1 2
All of the factors had significant effects on granule hardness (p < 0.05). Agitator speed (p1) and massing time (p2) had positive effects, whereas binding suspension addition rate (p3) had a negative effect on granule hardness. This is presented in the contour plot in Figure 2a. The plot shows that higher agitator speeds and longer wet massing times increased granule hardness; granule hardness was also increased at low water addition rates. Higher agitator speeds and longer massing times could increase granule density and hardness because high shearing forces consistently exist between the particles [8,13,14,52]. Thereore, the bonding strength between particles strengthens the resistance to the separating forces.

3.2.2. Significant Factors for Carr’s Index (q2)

The flowability of the granules is significantly influenced by the granulation process [53]. Carr’s index is frequently used in pharmaceutics as an indication of the flowability of a powder. In high flowability granules, the bulk density and tapped density would be close in value; therefore, Carr’s index would be small. In contrast, in poor flowability powders with greater interparticle interactions, the difference between the bulk and tapped densities would be larger; therefore, Carr’s index would be larger. Carr’s index greater than 25 is considered an indication of poor flowability, whereas below 15 indicates good flowability [17,33]. As shown in Table 3, the Carr’s indexes manufactured by different process parameters ranged from 9.00–16.39, which could be interpreted as an indication of good flowability [54]. Regression analysis of the response factor, with significance at p < 0.05, was used to generate the empirical model described by Equation (5).
q 2 = 13.05 2.25 p 1 1.45 p 2
Agitator speed (p1) and massing time (p2) had significantly negative effects on granule flowability. The effect of agitator speed and massing time on Carr’s index of the granules is presented in the contour plot in Figure 2b. An increase in agitator speed and massing time could generate high shearing forces between the particles. Consequently, the growth and densification of the granules would progress. Therefore, Carr’s index decreases with increased agitator speed and massing time.

3.2.3. Significant Factors for Granule Size (q3)

The process parameters in the high-shear wet granulation process might significantly influence granule size. Granule size distribution might be significantly related to desired qualities in key areas, such as product uniformity, dissolution, flow, hardness, and bioavailability. Granule size is also used to optimize downstream processing as it relates to blending and compression. As shown in Table 3, the granule sizes manufactured by different process parameters ranged from 170.09–393.59 µm, representing the surface mean diameter. Regression analysis of the response factor, with a significance value of p < 0.05, was used to generate the empirical model described by Equation (6).
q 3   = 264.61 + 87.66 p 1 + 19.05 p 2
Equation (6) and the contour plot in Figure 2c show that the agitator speed (p1) had the largest effect on the surface mean diameter, followed by the massing time (p2). Both variables have a synergistic effect on the surface mean diameter. Long wet massing times and high agitator speeds may result in granule growth by coalescence, but large granules may undergo breakage until a steady state of granule size is reached. This implies that agitator speed and massing time have greater potential for controlling the granulation process [19,55,56].

3.2.4. Significant Factors for Granule True Density (q4) and Bulk Density (q5)

The granule true density after high-shear wet granulation was determined using a pycnometer. As shown in Table 3, the granule true density manufactured by different process parameters ranged from 1.409–1.491 g/cm3. Regression analysis of the response factor, with a significance value of p < 0.05, was used to generate the empirical model described by Equation (7).
q 4   = 1.45 + 0.0025 p 1 0.0173 p 2 0.0024 p 3 + 0.0092 p 1 p 2 + 0.0285 p 2 p 3
The equation shows that the significant factors affecting granule true density were agitator speed (p1), binding suspension addition rate (p3), massing time (p2), the mutual interactions between agitator speed and massing time, and mutual interactions between binding suspension, addition rate, and massing time. The contour plot presented in Figure 2d shows the effect of binding suspension addition rate and massing time on granule true density. It is evident that the granule true density increased with decreasing binding suspension addition rate and massing time (negative effect). Meanwhile, these two parameters exhibited a strong positive interaction (i.e., massing time showed a larger impact on the granule true density when using a larger binding suspension addition rate).
As shown in Table 1, the granule bulk density manufactured by different process parameters ranged from 0.467–0.479 g/mL. Regression analysis of the response factor, with a significance value of p < 0.05, was used to generate the empirical model described by Equation (8).
q 5   = 0.4715 + 0.0015 p 1 + 0.0024 p 2 + 0.002 p 3 0.0030 p 1 p 2 + 0.0034 p 2 p 3
The equation shows that agitator speed (p1), binding suspension addition rate (p3), massing time (p2), the mutual interactions between them were the significant factors affecting granule bulk density. The effect of binding suspension addition rate and massing time on the granule bulk density is presented in Figure 2e as a contour plot. The granule bulk density increased with increasing binding suspension addition rate and massing time (positive effect). These two parameters also exhibited a strong positive interaction. Meanwhile, the mutual interaction between agitator speed and massing time had a negative effect on the granule bulk density.

3.2.5. Significant Factors for Granule FE (q6)

Recently, the flow properties of pharmaceutical powders have been characterized with the rheological measurements [57,58,59,60]. As the blade passed through a predetermined distance into a powder bed, the torque and force on the blade were measured [60,61]. As shown in Table 3, the granule FE values are ranged from 121–220. Regression analysis of the response factor, with p < 0.05, was used to generate the empirical model described by Equation (9).
q 6 = 175.44 30.22 p 1 19.54   p 2 + 3.58 p 3 + 8.97 p 2 p 3
The equation shows that the significant factors affecting granule FE were agitator speed (p1), massing time (p2), binding suspension addition rate (p3), and the mutual interactions between binding suspension addition rate and massing time. The contour plot presented in Figure 2f shows the effect of blade speed and massing time on the granule FE. Based on the coefficients of the empirical model, the granule flow energy decreased with increasing blade speed and massing time (negative effect).

3.2.6. Significant Factors for Tablet Tensile Strength (q7)

Tablet dosage forms should have a certain strength to resist the mechanical stress of handling during manufacture, packaging, and shipping [31]. In addition, previous studies reported that tablet strength was related to tablet disintegration time, dissolution profiles, and bioavailability [62]. Table 3 shows that the tablet tensile strengths of different formulations ranged from 125.88–145.51 N/cm2. Regression analysis of the response factor, with a significance of p < 0.05, was used to generate the empirical model described by Equation (10).
q 7   = 137.54 + 6.13 p 1 + 3.25 p 2 0.70 p 3 1.59 p 1 p 2 1.56 p 1 p 3
This equation and the contour plot in Figure 3a show that agitator speed (p1) had the highest effect on tablet tensile strength, followed by massing time (p2). The positive coefficients of agitator speed and massing time indicate that higher values for these parameters result in higher tensile strength. binding suspension addition rate (p3) had a negative effect on tensile strength. The mutual interactions between the agitator speed and binding suspension addition rate, and between the agitator speed and massing time had negative effects on the tensile strength of the tablet.

3.2.7. Significant Factors for Dissolution (q8q11)

Table 3 shows the dissolution results of the telmisartan tablets manufactured with different process parameters. Regression analysis of the response factors, with a significance value of p values < 0.05, was used to generate the empirical models described by Equations (11)–(14).
q 8 = 78.45 10.86 p 1 1.35 p 2 7.24 p 1 2
q 9 = 82.51 5.97 p 1 1.89 p 2
q 10 = 88.04 4.07 p 1 1.64 p 2 3.05 p 1 2
q 11 = 89.50 2.69 p 1 1.10 p 2
As shown in Equations (11)–(14) and the contour plots in Figure 3b–e, the agitator speed (p1) and massing time (p2) could have significant, negative effects on dissolution, whereas the binding suspension addition rate (p3) may have no effect on dissolution. From 5–30 min, the larger p1 value could indicate that dissolution over the initial 30 min is significantly controlled by agitator speed. The blade presents the mechanical energy to generate granules, whereas the mechanical energy generated by the chopper breaks up granules that are generated. Thus, the response of the granules, such as granule size and growth rate could be significantly influenced by the energy resulting from the agitator and chopper speed [63,64]. Generally, larger granules have a lower dissolution profile than smaller granules because smaller granules have a larger surface area in contact with the dissolution medium. Therefore, increased agitator speed could result in increased granule size, negatively influencing the dissolution profile. The massing time also had a negative effect on the dissolution profile. A high massing time could present the mechanical energy required for mixing the feed powder. This mechanical energy could increase the granule size; thus, the dissolution profile could be negatively influenced by the massing time.

3.3. Optimal Process Parameters and Monte Carlo Simulations

The optimization of the process parameters was conducted with statistically validated empirical equations, which were used to interpret the relationship between the control and response factors. The optimization condition for the high-shear granulation process presented the following: maximum tablet tensile strength within the lower (≥ 125 N/cm2) limit; minimum Carr’s index within the upper (≤ 15) limit; optimal intermediate ranges of dissolved drug at 5 min (72.29–83.29%), 10 min (78.68–88.58%), 15 min (85.00–92.10%), and 30 min (87.21–94.47%). In addition, a 95% confidence interval of the optimal condition was used for the control strategy. Based on these conditions, the responses were then combined to reveal the overall optimum region. This optimum region of the high-shear granulation process is shown in Figure 4a, in which the yellow area indicates the values of the control factors that achieve the desired responses simultaneously. To investigate the robustness of the optimum region, Monte Carlo simulations were performed with MODDE® software (Sartorius Stedim Biotech., version 12.0.1, Umeå, Sweden). To find the most robust process parameters in the operating space, the setting point should have a DPMO (Defects Per Million Opportunities) around or less than the specification. A DMPO result of 1000 means 0.1% risk of failure. The result is presented in Figure 4b with the probability of failure. The green area shows a quality certainty of 99.9%. Furthermore, as the probability of failure increased, the area was represented with red color. This suggested that it was not safe to obtain the target response at a higher agitator speed and higher massing time. However, agitator speeds ranging 700–900 rpm and massing times ranging 2.5–3.5 min in high-shear wet granulation had a low probability of failure less than 1%. To validate the optimal process parameters, these settings were performed at the lab scale. The spray rate was fixed with 5.20 mL/min. The agitator speed was selected 700, 800, and 900 rpm in the operating space. The massing time was also selected in the operating space as 2.5, 3.0, 3.5 min. Table 4 presents the experimental results with the target values of the response factors. Additionally, differences in biases and relative bias percentages between the response values generated by the optimal process parameters and the target values were small. The difference in biases relative bias percentages were lower than 2.20 and 5.08. This suggests that a risk reduction for the high-shear wet granulation at the lab scale can be obtained using QbD approach that controlled the development process from the initial stage with RA to robust operating space.

3.4. Multivariate Analysis between CQAs of Drug Product and IQAs of Process

The purpose of PCA is to identify relationships between IQAs and CQAs in the experimental design. All 15 experimental design samples were included in the PCA. 11 variables were analyzed, including IQAs, such as granule size, granule hardness, Carr’s index, granule FE, bulk density, and true density, and CQAs, such as tensile strength and % dissolution at each time point. Two PC models were fitted to the variables. As shown in Figure 5, the first and second PCs (P[1] and P[2], respectively) explained 76.0% and 11.3% of the overall variability in the dataset, respectively. R2 is used as a measure of how well the variance of the data is explained in the model, and Q2 is used as a measure of how well the variance of the data is predicted in the model. The closer each value is to 1.0, the better the fit of the model. A score scatter plot in Figure 6a suggests there are no groups and no outliers outside the ellipse (95% confidence interval). In addition, the center point samples (run order 2, 5, and 15) were located near each center, indicating good reproducibility. Run order 1 and 3 were located further away from the other run orders. Score contribution plots of run order 1 and 3 are shown in Figure 7. It identifies why these differ more and which variables are the greatest contributing factors to the differences. The variables with larger positive and lesser negative values had significant effects on differentiating run order 1 and 3. The variables were shown to have completely opposite effects. This suggests that agitator speeds ranging 400–1200 rpm and massing time ranging 1–5 min in high-shear wet granulation have significant effects on IQAs and CQAs.
The dataset can be displayed on a loading plot. The loading value reflects the weight of the variable relative to the principal component. In the score and loading plot, if the variable is highly correlated with the principal component, a positive loading value or a negative loading value is assigned. The range is −1 to +1. If the value is near 0, it indicates that the variable does not provide information to this main component. The loading plot can be used to identify the covariance between variables and to interpret the pattern observed in the score plot. As shown in Figure 6b, the loading scatter plot can elucidate the relationship between the IQAs and CQAs. In the case of P[1], which fully explains the degree of variance in the overall data, tensile strength (CQA), with a strong positive loading value, and Carr’s index and granule Fe, with strong negative loading values, can be regarded as having a negative relationship, because they are located on opposite sides. In addition, it can be seen that granule size and granule hardness, with strong positive loading values for P[1], and are located on the opposite side of % dissolution with a strong negative loading value; thus, they also have a negative relationship.
The purpose of PLS is to identify how all available process parameters and IQAs impact CQAs. PLS was used to define the relationship between X variables (process parameters and IQAs) and Y variables (CQAs). Figure 8 presents a predicted versus observed plot for CQAs from calibration models, which is presented by a solid black line. The R2 had % dissolution at 5 min, 10 min, and 15 min and tensile strength at 0.8805, 0.9367, 0.8364 and 0.8959, respectively. This might indicate a good model fit. As shown in Figure 9, the VIP (variable importance) presents the importance variables in the model both with respect to Y (CQAs) and X (process parameters and IQAs). VIP is normalized and the average squared VIP value is 1. Variables in the model with a VIP >1 are important to the model [1]. It appears that granule hardness, granule size, granule FE, Carr’s index, and agitator speed were important variables with respect to CQAs. Furthermore, massing time, true density, bulk density, and binding suspension addition rate did not have a significant effect on CQAs.
Based on the MVA, the relationship between IQAs and CQAs (granule size versus % dissolution, granule hardness versus % dissolution, granule FE versus tensile strength, and Carr’s index versus tensile strength) were fitted to establish prediction models with 95% confidence intervals and 95% prediction intervals. Figure 10, Figure 11 and Figure 12 show the fitted line plots. The models contained an adequate number of observations throughout the entire range of all the predictor values. The scatter points generally followed the regression line. There did not appear to be any curvature in the data. The p values for the regression models were less than the significance level of 0.05. The S value represents the standard deviation of the distance between the data values and the fitted values. Therefore, the S value can present how well the model describes the response. The lower the value of S, the better the model describes the response. The S values of the prediction models were lower than 3.719, which indicated that the models could accurately predict the response. Moreover, the higher R-Sq (> 80.7%) also indicated a good model fit. Therefore, using these models, CQAs could be accurately predicted with IQAs. In addition, the target values for IQAs could be established with these models.

3.5. Scale-Up to Validate the Mutual Effects between CPPs, IQAs, and CQAs

Process development using the QbD approach was studied at the lab scale (granulator capacity is 1.0 L). To validate the mutual effects between CPPS, IQAs, and CQAs, and manufacture clinical batch, a scale-up from lab to pilot scale was performed using geometric similarity, agitator torque profile, and agitator blade tip speed. The optimal process obtained from the lab-scale experimental design was used. Based on the experimental design and MVA, IQAs, and CQAs were significantly affected by agitator speed and massing time, whereas binding suspension addition rate did not have a significant effect. Therefore, the agitator blade tip speeds were kept constant regardless of granulator size. The agitator speed used in the scale-up is shown in Table 5. This was calculated from the agitator blade tip speed and the diameters of the agitator blade. Massing time was determined by the agitator torque profile. In various literature, agitator torque profile has been used for a scale-up strategy and granules might have similar properties when agitator torque profiles were similar [30,40,65,66]. The agitator torque profiles at the lab scale and pilot scale are shown in Figure 13. Based on the profiles, it was indicated that granule properties at the lab scale and pilot scale might be similar even with the same mixing time, which was 3 min. A geometric similarity was obtained using a similar granulator (P 6). The Diosna granulator allows bowls of different capacities in the range 1-6 L to be easily interchanged, which means that different batch sizes can be manufactured using the same basic equipment [24,28,29,67]. In addition, the fill levels were maintained at a similar ratio, such as 38.6% at the lab scale and 39.2% at the pilot scale.
Table 6 presents the validated results of IQAs and CQAs at the lab and pilot scale process. The absolute biases and relative bias percentages between the IQAs and CQAs generated by lab and pilot scale process had small differences, which were less than 9.09% (relative biases of friability). The IQA and CQA values might demonstrate similar target values, such as maximum tablet tensile strength within the lower (≥ 135 N/cm2) and upper (≤ 145 N/cm2) limits; minimum Carr’s index within the lower (≥ 9) and upper (≤ 15) limits; optimal intermediate ranges of dissolved drug at 5 min (72.29–83.29%), 10 min (78.68–88.58%), 15 min (85.00–92.10%), and 30 min (84.47–94.47%); minimum friability less than 0.25%; optimal intermediate ranges of assay (95–105%); minimum content uniformity less than 5.0% RSD. This suggests that the integrated process parameters obtained from experimental design at the lab scale could be applied at the pilot scale with the scale-up strategy. Additionally, the similar relationship between IQAs, CQAs, and CPPs was shown regardless of batch size, it can be used to predict CQAs for control strategy.

4. Conclusions

Based on the control strategy, the development of process parameters for high-shear wet granulation was studied at the lab scale with the QbD approach. Risk assessment was conducted to determine the control and response factors and it was based on a primary knowledge and the target values of the control tablet. A response surface design was performed to obtain the optimal process parameters. The predictable empirical equations of best fit were generated with the quantitative effects of the different control factors on the response factors. The operating space at the lab scale was obtained with the empirical equations. The robustness of the operating space was identified with a Monte Carlo simulation. The relationship between IQAs, CQAs, and CPPs was investigated with MVA. The validation results at the pilot scale suggested that the integrated process parameters obtained from experimental design at the lab scale could be applied at the pilot scale with the scale-up strategy. This study suggests that the risk reduction for the manufacturing process could be obtained with integrated process parameters as a result of the QbD approach, and the relationship among IQAs, CQAs, and CPPs could be used to predict CQAs for control strategy and SUPAC (Scale-Up and Post-Approval Guidance).

Author Contributions

Conceptualization: D.H.C.; Methodology: D.H.C., Software: D.H.C.; Validation: J.K.H.; Formal analysis: J.K.H.; Investigation: J.K.H.; Resources: D.H.C.; Data curation: D.H.C.; Writing-original draft preparation: J.K.H.; Writing-review and editing: B.S.S., D.H.C.; Visualization: B.S.S.; Supervision: B.S.S., D.H.C.; Project administration: D.H.C.; Funding acquisition: D.H.C.

Funding

This work was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2016R1D1A1B03931513).

Conflicts of Interest

The authors report no conflicts of interest.

References

  1. Huang, J.; Kaul, G.; Cai, C.; Chatlapalli, R.; Hernandez-Abad, P.; Ghosh, K.; Nagi, A. Quality by design case study: An integrated multivariate approach to drug product and process development. Int. J. Pharm. 2009, 382, 23–32. [Google Scholar] [CrossRef] [PubMed]
  2. Wu, H.; Tawakkul, M.; White, M.; Khan, M.A. Quality-by-Design (QbD): An integrated multivariate approach for the component quantification in powder blends. Int. J. Pharm. 2009, 372, 39–48. [Google Scholar] [CrossRef] [PubMed]
  3. Rajalahti, T.; Kvalheim, O.M. Multivariate data analysis in pharmaceutics: A tutorial review. Int. J. Pharm. 2011, 417, 280–290. [Google Scholar] [CrossRef] [PubMed]
  4. Gabrielsson, J.; Lindberg, N.O.; Lundstedt, T. Multivariate methods in pharmaceutical applications. J. Chemom. 2002, 16, 141–160. [Google Scholar] [CrossRef]
  5. Naelapää, K.; Allesø, M.; Kristensen, H.G.; Bro, R.; Rantanen, J.; Bertelsen, P. Increasing process understanding by analyzing complex interactions in experimental data. J. Pharm. Sci. 2009, 98, 1852–1861. [Google Scholar] [CrossRef] [PubMed]
  6. Choi, D.H.; Shin, S.; Viet Truong, N.K.; Jeong, S.H. A new experimental design method to optimize formulations focusing on a lubricant for hydrophilic matrix tablets. Drug Dev. Ind. Pharm. 2012, 38, 1117–1127. [Google Scholar] [CrossRef] [PubMed]
  7. Hansuld, E.M.; Briens, L. A review of monitoring methods for pharmaceutical wet granulation. Int. J. Pharm. 2014, 472, 192–201. [Google Scholar] [CrossRef]
  8. Badawy, S.I.F.; Menning, M.M.; Gorko, M.A.; Gilbert, D.L. Effect of process parameters on compressibility of granulation manufactured in a high-shear mixer. Int. J. Pharm. 2000, 198, 51–61. [Google Scholar] [CrossRef]
  9. Benali, M.; Gerbaud, V.; Hemati, M. Effect of operating conditions and physico–chemical properties on the wet granulation kinetics in high shear mixer. Powder Technol. 2009, 190, 160–169. [Google Scholar] [CrossRef]
  10. Butensky, M.; Hyman, D. Rotary Drum Granulation. An Experimental Study of the Factors Affecting Granule Size. Ind. Eng. Chem. Fundam. 1971, 10, 212–219. [Google Scholar] [CrossRef]
  11. Liu, H.; Wang, K.; Schlindwein, W.; Li, M. Using the Box-Behnken experimental design to optimise operating parameters in pulsed spray fluidised bed granulation. Int. J. Pharm. 2013, 448, 329–338. [Google Scholar] [CrossRef] [PubMed]
  12. Mangwandi, C.; Albadarin, A.B.; Al-Muhtaseb, A.a.H.; Allen, S.J.; Walker, G.M. Optimisation of high shear granulation of multicomponent fertiliser using response surface methodology. Powder Technol. 2013, 238, 142–150. [Google Scholar] [CrossRef]
  13. Rahmanian, N.; Ghadiri, M.; Jia, X.; Stepanek, F. Characterisation of granule structure and strength made in a high shear granulator. Powder Technol. 2009, 192, 184–194. [Google Scholar] [CrossRef]
  14. Cavinato, M.; Bresciani, M.; Machin, M.; Bellazzi, G.; Canu, P.; Santomaso, A. The development of a novel formulation map for the optimization of high shear wet granulation. Chem. Eng. J. 2010, 164, 350–358. [Google Scholar] [CrossRef]
  15. Hegedűs, Á.; Pintye-Hódi, K. Influence of the type of the high-shear granulator on the physico-chemical properties of granules. Chem. Eng. Process. 2007, 46, 1012–1019. [Google Scholar] [CrossRef]
  16. Pathare, P.B.; Baş, N.; Fitzpatrick, J.J.; Cronin, K.; Byrne, E.P. Effect of high shear granulation process parameters on the production of granola cereal aggregates. Biosyst. Eng. 2011, 110, 473–481. [Google Scholar] [CrossRef]
  17. Kumar, A.; Gernaey, K.V.; Beer, T.D.; Nopens, I. Model-based analysis of high shear wet granulation from batch to continuous processes in pharmaceutical production—A critical review. Eur. J. Pharm. Biopharm. 2013, 85, 814–832. [Google Scholar] [CrossRef]
  18. Kumaresan, T.; Joshi, J.B. Effect of impeller design on the flow pattern and mixing in stirred tanks. Chem. Eng. J. 2006, 115, 173–193. [Google Scholar] [CrossRef]
  19. Badawy, S.I.F.; Narang, A.S.; LaMarche, K.; Subramanian, G.; Varia, S.A. Mechanistic basis for the effects of process parameters on quality attributes in high shear wet granulation. Int. J. Pharm. 2012, 439, 324–333. [Google Scholar] [CrossRef]
  20. Shi, L.; Feng, Y.; Sun, C.C. Massing in high shear wet granulation can simultaneously improve powder flow and deteriorate powder compaction: A double-edged sword. Eur. J. Pharm. Sci. 2011, 43, 50–56. [Google Scholar] [CrossRef]
  21. Mangwandi, C.; Adams, M.; Hounslow, M.; Salman, A. Effect of impeller speed on mechanical and dissolution properties of high-shear granules. Chem. Eng. J. 2010, 164, 305–315. [Google Scholar] [CrossRef]
  22. Mirza, Z.; Liu, J.; Glocheux, Y.; Albadarin, A.B.; Walker, G.M.; Mangwandi, C. Effect of impeller design on homogeneity, size and strength of pharmaceutical granules produced by high-shear wet granulation. Particuology 2015, 23, 31–39. [Google Scholar] [CrossRef] [Green Version]
  23. Bouwman, A.; Henstra, M.; Westerman, D.; Chung, J.; Zhang, Z.; Ingram, A.; Seville, J.; Frijlink, H. The effect of the amount of binder liquid on the granulation mechanisms and structure of microcrystalline cellulose granules prepared by high shear granulation. Int. J. Pharm. 2005, 290, 129–136. [Google Scholar] [CrossRef] [PubMed]
  24. Chitu, T.M.; Oulahna, D.; Hemati, M. Wet granulation in laboratory-scale high shear mixers: Effect of chopper presence, design and impeller speed. Powder Technol. 2011, 206, 34–43. [Google Scholar] [CrossRef]
  25. Iveson, S.M.; Litster, J.D.; Hapgood, K.; Ennis, B.J. Nucleation, growth and breakage phenomena in agitated wet granulation processes: A review. Powder Technol. 2001, 117, 3–39. [Google Scholar] [CrossRef]
  26. Iveson, S.M.; Wauters, P.A.; Forrest, S.; Litster, J.D.; Meesters, G.M.; Scarlett, B. Growth regime map for liquid-bound granules: Further development and experimental validation. Powder Technol. 2001, 117, 83–97. [Google Scholar] [CrossRef]
  27. Faure, A.; York, P.; Rowe, R. Process control and scale-up of pharmaceutical wet granulation processes: A review. Eur. J. Pharm. Biopharm. 2001, 52, 269–277. [Google Scholar] [CrossRef]
  28. Levin, M. How to Scale-Up a Wet Granulation End Point Scientifically; Academic Press: Cambridge, MA, USA, 2015; ISBN 978-0-12-803522-1. [Google Scholar]
  29. Landin, M.; York, P.; Cliff, M.; Rowe, R.; Wigmore, A. Scale-up of a pharmaceutical granulation in fixed bowl mixer-granulators. Int. J. Pharm. 1996, 133, 127–131. [Google Scholar] [CrossRef]
  30. Leuenberger, H. New trends in the production of pharmaceutical granules: The classical batch concept and the problem of scale-up. Eur. J. Pharm. Biopharm. 2001, 52, 279–288. [Google Scholar] [CrossRef]
  31. Parikh, D.M. Handbook of Pharmaceutical Granulation Technology, 2nd ed.; Taylor and Francis Group: Oxfordshire, UK, 2005; pp. 7–77. ISBN 978-0-8247-2647-8. [Google Scholar]
  32. Parikh, D.M. Handbook of Pharmaceutical Granulation Technology, 3rd ed.; informa: London, UK, 2016; ISBN 978-1-4398-0789-7. [Google Scholar]
  33. Lee, A.R.; Kwon, S.Y.; Choi, D.H.; Park, E.S. Quality by Design (QbD) approach to optimize the formulation of a bilayer combination tablet (Telmiduo®) manufactured via high shear wet granulation. Int. J. Pharm. 2017, 534, 144–158. [Google Scholar] [CrossRef]
  34. Aksu, B.; Paradkar, A.; de Matas, M.; Özer, Ö.; Güneri, T.; York, P. A quality by design approach using artificial intelligence techniques to control the critical quality attributes of ramipril tablets manufactured by wet granulation. Pharm. Dev. Technol. 2013, 18, 236–245. [Google Scholar] [CrossRef] [PubMed]
  35. Basalious, E.B.; El-Sebaie, W.; El-Gazayerly, O. Application of pharmaceutical QbD for enhancement of the solubility and dissolution of a class II BCS drug using polymeric surfactants and crystallization inhibitors: Development of controlled-release tablets. AAPS PharmSciTech 2011, 12, 799–810. [Google Scholar] [CrossRef] [PubMed]
  36. Charoo, N.A.; Ali, A.A. Quality risk management in pharmaceutical development. Drug Dev. Ind. Pharm. 2013, 39, 947–960. [Google Scholar] [CrossRef] [PubMed]
  37. Cogdill, R.P.; Drennen, J.K. Risk-based quality by design (QbD): A Taguchi perspective on the assessment of product quality, and the quantitative linkage of drug product parameters and clinical performance. J. Pharm. Innov. 2008, 3, 23–29. [Google Scholar] [CrossRef]
  38. Xie, L.; Wu, H.; Shen, M.; Augsburger, L.L.; Lyon, R.C.; Khan, M.A.; Hussain, A.S.; Hoag, S.W. Quality-by-design (QbD): Effects of testing parameters and formulation variables on the segregation tendency of pharmaceutical powder measured by the ASTM D 6940-04 segregation tester. J. Pharm. Sci. 2008, 97, 4485–4497. [Google Scholar] [CrossRef] [PubMed]
  39. Xu, X.; Khan, M.A.; Burgess, D.J. A quality by design (QbD) case study on liposomes containing hydrophilic API: I. Formulation, processing design and risk assessment. Int. J. Pharm. 2011, 419, 52–59. [Google Scholar] [CrossRef] [PubMed]
  40. Cavinato, M.; Bresciani, M.; Machin, M.; Bellazzi, G.; Canu, P.; Santomaso, A.C. Formulation design for optimal high-shear wet granulation using on-line torque measurements. Int. J. Pharm. 2010, 387, 48–55. [Google Scholar] [CrossRef]
  41. Staniforth, J.N.; Baichwal, A.R.; Hart, J.P. Interpretation of creep behaviour of microcrystalline cellulose powders and granules during compaction. Int. J. Pharm. 1987, 40, 267–269. [Google Scholar] [CrossRef]
  42. Andersson, M.; Ringberg, A.; Gustafsson, C. Multivariate methods in tablet formulation suitable for early drug development: Predictive models from a screening design of several linked responses. Chemom. Intell. Lab. Syst. 2007, 87, 125–130. [Google Scholar] [CrossRef]
  43. Dumarey, M.; Goodwin, D.J.; Davison, C. Multivariate modelling to study the effect of the manufacturing process on the complete tablet dissolution profile. Int. J. Pharm. 2015, 486, 112–120. [Google Scholar] [CrossRef]
  44. Choi, D.H.; Kim, N.A.; Chu, K.R.; Jung, Y.J.; Yoon, J.H.; Jeong, S.H. Material properties and compressibility using Heckel and Kawakita equation with commonly used pharmaceutical excipients. J. Pharm. Investig. 2010, 40, 237–244. [Google Scholar] [CrossRef]
  45. Thapa, P.; Choi, D.H.; Kim, M.S.; Jeong, S.H. Effects of granulation process variables on the physical properties of dosage forms by combination of experimental design and principal component analysis. Asian J. Pharm. Sci. 2018. [Google Scholar] [CrossRef]
  46. Oka, S.; Kašpar, O.; Tokárová, V.; Sowrirajan, K.; Wu, H.; Khan, M.; Muzzio, F.; Štěpánek, F.; Ramachandran, R. A quantitative study of the effect of process parameters on key granule characteristics in a high shear wet granulation process involving a two component pharmaceutical blend. Adv. Powder Technol. 2015, 26, 315–322. [Google Scholar] [CrossRef]
  47. Rahmanian, N.; Naji, A.; Ghadiri, M. Effects of process parameters on granules properties produced in a high shear granulator. Chem. Eng. Res. Des. 2011, 89, 512–518. [Google Scholar] [CrossRef]
  48. Cavinato, M.; Andreato, E.; Bresciani, M.; Pignatone, I.; Bellazzi, G.; Franceschinis, E.; Realdon, N.; Canu, P.; Santomaso, A.C. Combining formulation and process aspects for optimizing the high-shear wet granulation of common drugs. Int. J. Pharm. 2011, 416, 229–241. [Google Scholar] [CrossRef]
  49. Huang, Y.B.; Tsai, Y.H.; Yang, W.C.; Chang, J.S.; Wu, P.C.; Takayama, K. Once-daily propranolol extended-release tablet dosage form: Formulation design and in vitro/in vivo investigation. Eur. J. Pharm. Biopharm. 2004, 58, 607–614. [Google Scholar] [CrossRef] [PubMed]
  50. Huang, Y.B.; Tsai, Y.H.; Lee, S.H.; Chang, J.S.; Wu, P.C. Optimization of pH-independent release of nicardipine hydrochloride extended-release matrix tablets using response surface methodology. Int. J. Pharm. 2005, 289, 87–95. [Google Scholar] [CrossRef]
  51. Fouad, S.A.; Basalious, E.B.; El-Nabarawi, M.A.; Tayel, S.A. Microemulsion and poloxamer microemulsion-based gel for sustained transdermal delivery of diclofenac epolamine using in-skin drug depot: In vitro/in vivo evaluation. Int. J. Pharm. 2013, 453, 569–578. [Google Scholar] [CrossRef]
  52. Rahmanian, N.; Ghadiri, M.; Ding, Y. Effect of scale of operation on granule strength in high shear granulators. Chem. Eng. Sci. 2008, 63, 915–923. [Google Scholar] [CrossRef]
  53. Javaheri, H.; Carter, P.; Elkordy, A. Wet granulation to overcome liquisolid technique issues of poor flowability and compactibility: A study to enhance glibenclamide dissolution. J. Pharm. Drug Dev. 2014, 1, 501–512. [Google Scholar] [CrossRef]
  54. Kumar, V.; de la Luz Reus-Medina, M.; Yang, D. Preparation, characterization, and tabletting properties of a new cellulose-based pharmaceutical aid. Int. J. Pharm. 2002, 235, 129–140. [Google Scholar] [CrossRef]
  55. Knight, P.C.; Johansen, A.; Kristensen, H.G.; Schæfer, T.; Seville, J.P.K. An investigation of the effects on agglomeration of changing the speed of a mechanical mixer. Powder Technol. 2000, 110, 204–209. [Google Scholar] [CrossRef]
  56. Ohno, I.; Hasegawa, S.; Yada, S.; Kusai, A.; Moribe, K.; Yamamoto, K. Importance of evaluating the consolidation of granules manufactured by high shear mixer. Int. J. Pharm. 2007, 338, 79–86. [Google Scholar] [CrossRef] [PubMed]
  57. Chitu, T.M.; Oulahna, D.; Hemati, M. Rheology, granule growth and granule strength: Application to the wet granulation of lactose–MCC mixtures. Powder Technol. 2011, 208, 441–453. [Google Scholar] [CrossRef]
  58. Klausner, J.F.; Chen, D.; Mei, R. Experimental investigation of cohesive powder rheology. Powder Technol. 2000, 112, 94–101. [Google Scholar] [CrossRef]
  59. Bharadwaj, R.; Ketterhagen, W.R.; Hancock, B.C. Discrete element simulation study of a Freeman powder rheometer. Chem. Eng. Sci. 2010, 65, 5747–5756. [Google Scholar] [CrossRef]
  60. Freeman, R. Measuring the flow properties of consolidated, conditioned and aerated powders—A comparative study using a powder rheometer and a rotational shear cell. Powder Technol. 2007, 174, 25–33. [Google Scholar] [CrossRef]
  61. Hare, C.; Zafar, U.; Ghadiri, M.; Freeman, T.; Clayton, J.; Murtagh, M. Analysis of the dynamics of the FT4 powder rheometer. Powder Technol. 2015, 285, 123–127. [Google Scholar] [CrossRef] [Green Version]
  62. Liebermann, H.; Lachman, L.; Kanig, J. The Theory and Practice of Industrial Pharmacy, 4th ed.; CBS Publishers and Distributors: New Delhi, India, 2013; ISBN 978-81-239-2289-8. [Google Scholar]
  63. Knight, P. An investigation of the kinetics of granulation using a high shear mixer. Powder Technol. 1993, 77, 159–169. [Google Scholar] [CrossRef]
  64. Sunada, H.; Hasegawa, M.; Makino, T.; Sakamoto, H.; Fujita, K.; Tanino, T.; Kokubo, H.; Kawaguchi, T. Study of standard tablet formulation based on fluidized-bed granulation. Drug Dev. Ind. Pharm. 1998, 24, 225–233. [Google Scholar] [CrossRef]
  65. Koide, T.; Nagato, T.; Kanou, Y.; Matsui, K.; Natsuyama, S.; Kawanishi, T.; Hiyama, Y. Detection of component segregation in granules manufactured by high shear granulation with over-granulation conditions using near-infrared chemical imaging. Int. J. Pharm. 2013, 441, 135–145. [Google Scholar] [CrossRef] [PubMed]
  66. Bhargava, A.; Jelen, P. Lactose solubility and crystal growth as affected by mineral impurities. J. Food Sci. 1996, 61, 180–184. [Google Scholar] [CrossRef]
  67. Chitu, T.M.; Oulahna, D.; Hemati, M. Wet granulation in laboratory scale high shear mixers: Effect of binder properties. Powder Technol. 2011, 206, 25–33. [Google Scholar] [CrossRef]
Figure 1. Pareto chart of the risk priority number (RPN). Individual values are represented in descending order by bars, and the cumulative total is represented by the line. The left vertical axis is the RPN. The right vertical axis is the cumulative percentage of the RPN.
Figure 1. Pareto chart of the risk priority number (RPN). Individual values are represented in descending order by bars, and the cumulative total is represented by the line. The left vertical axis is the RPN. The right vertical axis is the cumulative percentage of the RPN.
Pharmaceutics 11 00252 g001
Figure 2. Contour plots correlating response factors (IQAs) to optimize process parameters of high-shear wet granulation for the telmisartan test tablet. (a) Granule hardness, (b) Carr’s index, (c) Granule size, (d) True density, (e) Bulk density, and (f) Granule flow energy (FE).
Figure 2. Contour plots correlating response factors (IQAs) to optimize process parameters of high-shear wet granulation for the telmisartan test tablet. (a) Granule hardness, (b) Carr’s index, (c) Granule size, (d) True density, (e) Bulk density, and (f) Granule flow energy (FE).
Pharmaceutics 11 00252 g002
Figure 3. Contour plots correlating response factors (CQAs) to optimize process parameters of high-shear wet granulation for the telmisartan test tablet. (a) Tensile strength, (b) % dissolution at 5 min, (c) % dissolution at 10 min, (d) % dissolution at 15 min, and (e) % dissolution at 30 min.
Figure 3. Contour plots correlating response factors (CQAs) to optimize process parameters of high-shear wet granulation for the telmisartan test tablet. (a) Tensile strength, (b) % dissolution at 5 min, (c) % dissolution at 10 min, (d) % dissolution at 15 min, and (e) % dissolution at 30 min.
Pharmaceutics 11 00252 g003
Figure 4. (a) Operating spaces (overlay plot) to optimize process parameters of high-shear wet granulation for the telmisartan test tablet. The yellow area is the operating space with 95% confidence intervals. (b) Monte Carlo simulation results for the operating space. The green area shows a quality certainty of 99.9.
Figure 4. (a) Operating spaces (overlay plot) to optimize process parameters of high-shear wet granulation for the telmisartan test tablet. The yellow area is the operating space with 95% confidence intervals. (b) Monte Carlo simulation results for the operating space. The green area shows a quality certainty of 99.9.
Pharmaceutics 11 00252 g004aPharmaceutics 11 00252 g004b
Figure 5. Cumulative R2 (R2X(cum)) and Q2 (Q2(cum)) of response factors (IQAs and CQAs) for two principal components.
Figure 5. Cumulative R2 (R2X(cum)) and Q2 (Q2(cum)) of response factors (IQAs and CQAs) for two principal components.
Pharmaceutics 11 00252 g005
Figure 6. (a) Score scatter plot with two principal components (t[1] and t[2]) and (b) loading plot (p[1] and p[2]) presenting variable relationships.
Figure 6. (a) Score scatter plot with two principal components (t[1] and t[2]) and (b) loading plot (p[1] and p[2]) presenting variable relationships.
Pharmaceutics 11 00252 g006
Figure 7. Score contribution plots showing variables contributing to the difference between run order 1 (a) or run order 3 (b) and the average of all run orders.
Figure 7. Score contribution plots showing variables contributing to the difference between run order 1 (a) or run order 3 (b) and the average of all run orders.
Pharmaceutics 11 00252 g007
Figure 8. Predicted versus measured plots for tensile strengths (CQAs). (a) % dissolution at 5 min, (b) % dissolution at 10 min, (c) % dissolution at 15 min, and (d) tensile strength.
Figure 8. Predicted versus measured plots for tensile strengths (CQAs). (a) % dissolution at 5 min, (b) % dissolution at 10 min, (c) % dissolution at 15 min, and (d) tensile strength.
Pharmaceutics 11 00252 g008
Figure 9. VIP (variable importance) plot presenting very important variables in terms of CQAs.
Figure 9. VIP (variable importance) plot presenting very important variables in terms of CQAs.
Pharmaceutics 11 00252 g009
Figure 10. The fitted line plots of % dissolution ((a) 5 min, (b) 10 min, (c) 15 min, and (d) 30 min) versus granule hardness. The black dots, red line, blue dotted line, and black dotted line represent the experimental results, best fitted regression, 95% confidence intervals, and 95% prediction intervals, respectively.
Figure 10. The fitted line plots of % dissolution ((a) 5 min, (b) 10 min, (c) 15 min, and (d) 30 min) versus granule hardness. The black dots, red line, blue dotted line, and black dotted line represent the experimental results, best fitted regression, 95% confidence intervals, and 95% prediction intervals, respectively.
Pharmaceutics 11 00252 g010
Figure 11. The fitted line plots of % dissolution ((a) 5 min, (b) 10 min, (c) 15 min, and (d) 30 min) versus granule size. The black dots, red line, blue dotted line, and black dotted line represent the experimental results, best fitted regression, 95% confidence intervals, and 95% prediction intervals, respectively.
Figure 11. The fitted line plots of % dissolution ((a) 5 min, (b) 10 min, (c) 15 min, and (d) 30 min) versus granule size. The black dots, red line, blue dotted line, and black dotted line represent the experimental results, best fitted regression, 95% confidence intervals, and 95% prediction intervals, respectively.
Pharmaceutics 11 00252 g011
Figure 12. The fitted line plots of (a) tensile strength versus granule FE and (b) tensile strength versus Carr’s index. The black dots, red line, blue dotted line, and black dotted line represent the experimental results, best fitted regression, 95% confidence intervals, and 95% prediction intervals, respectively.
Figure 12. The fitted line plots of (a) tensile strength versus granule FE and (b) tensile strength versus Carr’s index. The black dots, red line, blue dotted line, and black dotted line represent the experimental results, best fitted regression, 95% confidence intervals, and 95% prediction intervals, respectively.
Pharmaceutics 11 00252 g012
Figure 13. Agitator torque profiles at the lab scale (black line) and pilot scale (red line). The agitator powder consumption is plotted versus liquid/solid (L/S) ratio (%).
Figure 13. Agitator torque profiles at the lab scale (black line) and pilot scale (red line). The agitator powder consumption is plotted versus liquid/solid (L/S) ratio (%).
Pharmaceutics 11 00252 g013
Table 1. Formulation of telmisartan granules.
Table 1. Formulation of telmisartan granules.
PurposeExcipientAmount (mg/tablet)
Active pharmaceutical ingredientTelmisartan80
Solubilizing agentNaOH6.7
Meglumine24
Binding agentPVP K2512
Thicking agentD-mannitol166
MCC 10139
Dicalcium phosphate9.3
Disintegrant agentF-melt typeC20
Crospovidone118
Lubricant agentSt-Mg5
Total480
Table 2. Risk assessment based on prior knowledge and experience, and information about the control tablet from published literature or analyses. (RPN: risk priority number; IQAs: intermediate quality attributes; CQA: critical quality attribute; CU: content uniformity; FE: flow energy)
Table 2. Risk assessment based on prior knowledge and experience, and information about the control tablet from published literature or analyses. (RPN: risk priority number; IQAs: intermediate quality attributes; CQA: critical quality attribute; CU: content uniformity; FE: flow energy)
Unit ProcessFailure ModeSPDRPNRisk DegreeRelated IQAs and CQAs
High shear granulation for test tabletAgitator speed555125Highdissolution, CU, tensile strength, Friability, granule size, granule hardness, Carr’s index, granule density, granule FE
Chopper speed52550Moderatedissolution, CU, tensile strength, Friability, granule size, granule hardness, Carr’s index, granule density, granule FE
Solvent spray rate455100Highdissolution, CU, tensile strength, Friability, granule size, granule hardness, Carr’s index, granule density, granule FE
Massing time555125Highdissolution, CU, tensile strength, Friability, granule size, granule hardness, Carr’s index, granule density, granule FE
Drying temp.3139Low tensile strength, Friability, granule hardness, Carr’s index, granule density
Drying time3139Lowtensile strength, Friability, granule hardness, Carr’s index, granule density
Table 3. Experimental design of high-shear wet granulation process parameters with three input control factors and values of the response factors for different process parameters
Table 3. Experimental design of high-shear wet granulation process parameters with three input control factors and values of the response factors for different process parameters
Run OrderControl FactorsResponse Factors
CPPIntermediate QAsDrug Product CQAs
Agitator Speed (rpm)Massing Time (min)Spray Rate (mL/min)Granule Hardness (N)Carr’s IndexGranule Size (μm)Granule True Density (g/cm3)Granule Bulk Density (g/mL)Granule FETensile Strength (N/cm2)% Dissolution
5 min10 min15 min30 min
p1p2p3q1q2q3q4q5q6q9q10q11q12q13
1120055.24.969.00393.591.4380.469121145.1983.2589.7791.2593.54
280035.22.0613.40256.001.4570.471180138.0080.2184.5687.1289.94
340015.20.9016.39186.941.4780.468220125.8882.5688.1289.2291.45
480013.71.4615.27232.851.4920.471205134.1582.2587.5888.6791.21
580035.22.1413.29265.951.4550.472179138.2578.1183.5587.4589.79
680056.72.5312.33299.631.4650.479166140.1279.8784.9789.9791.22
7120015.23.2811.74328.041.4630.470158142.7477.4582.5886.6489.12
840055.21.5513.65197.121.4160.479184134.5577.9582.1286.2589.84
940033.71.3414.97170.091.4560.471201131.5879.4584.2189.3190.55
10120036.73.5611.56359.091.4540.473155141.0077.7582.9887.1290.05
1140036.71.0715.54172.721.4430.475209133.4278.5884.9789.5491.21
12120033.75.0010.27347.431.4580.468138145.5158.4373.1279.5984.58
1380016.71.5414.07245.791.4340.468189132.8762.5578.9582.6587.89
1480053.73.2910.86141.801.4090.469146255.6660.5475.5581.2386.57
1580035.22.0813.40138.071.4490.471180258.2359.8774.6280.2285.56
Table 4. Optimal process parameters, target values, and validated results at the lab scale for the response factors
Table 4. Optimal process parameters, target values, and validated results at the lab scale for the response factors
Optimal Process ParametersResponse Factors
x1x2x3q1q2q7q8q9q10
Agitator Speed (rpm)Massing Time (min)Spray RRate (mL/min)Carr’s indexTensile Strength (N/cm2)Dissolution at 5 min (%)Dissolution at 10 min (%)Dissolution at 15 min (%)Dissolution at 30 min (%)
7002.55.2013.87133.179.582.3988.5590.26
Target values13.97135.181.0584.4789.2791.45
Absolute viases0.102.001.552.080.721.19
Relative biases (%)4.501.472.072.540.821.33
7003.55.2013.12134.7279.2582.2387.1387.9
Target values13.24136.9280.3783.5388.4589.89
Absolute viases0.122.201.121.301.321.99
Relative biases (%)5.081.611.501.621.542.29
9002.55.2012.79137.7175.4680.7286.1687.78
Target values12.85138.3575.6281.4887.2489.10
Absolute viases0.060.640.160.761.081.32
Relative biases (%)2.330.460.220.951.261.51
9003.55.2012.08139.6674.1879.6185.9387.64
Target values12.12139.7974.9480.5486.4288.55
Absolute viases0.040.130.760.930.490.91
Relative biases (%)1.430.091.041.160.571.03
80035.2013137.4177.6981.5887.5588.59
Target values13.04137.5478.4582.5188.0489.50
Absolute viases0.040.130.760.930.490.91
Relative biases (%)1.430.091.041.160.571.03
Table 5. Process parameter settings for the scale-up strategy of high-shear wet granulation.
Table 5. Process parameter settings for the scale-up strategy of high-shear wet granulation.
Process ParametersScale
LabPilot
Used granulatorP1P6
Granulator capacity
(L)
1.236.05
Batch size
(g)
4752375
Fill level
(%)
38.639.2
Agitator speed corresponding
to agitator tip speed of 6.0 m/s
800 rpm480 rpm
Massing time
(min)
33
Spray rate
(ml/min)
5.25.2
Table 6. Validated results for IQAs and CQAs in the scale up.
Table 6. Validated results for IQAs and CQAs in the scale up.
Quality AttributesScaleAbsolute BiasesRelative Biases (%)
LabPilot
IQAsGranule hardness (N)2.172.250.083.69
Carr’s index13.1213.450.332.52
Granule size (μm) 264.23275.4511.224.25
Granule true density (g/cm3)1.4531.4570.0040.28
Granule bulk density (g/mL)0.4760.4790.0030.63
Granule FE17517321.16
CQAsTensile strength (N/cm2)137.55138.250.700.51
% dissolution at 5 min78.5479.150.610.78
% dissolution at 10 min82.5582.960.410.50
% dissolution at 15 min88.1488.350.210.24
% dissolution at 30 min89.5290.080.560.63
Friability (%)0.120.110.019.09
Content uniformity (%)0.780.740.045.41
Assay (%)101.43100.780.650.64

Share and Cite

MDPI and ACS Style

Han, J.K.; Shin, B.S.; Choi, D.H. Comprehensive Study of Intermediate and Critical Quality Attributes for Process Control of High-Shear Wet Granulation Using Multivariate Analysis and the Quality by Design Approach. Pharmaceutics 2019, 11, 252. https://doi.org/10.3390/pharmaceutics11060252

AMA Style

Han JK, Shin BS, Choi DH. Comprehensive Study of Intermediate and Critical Quality Attributes for Process Control of High-Shear Wet Granulation Using Multivariate Analysis and the Quality by Design Approach. Pharmaceutics. 2019; 11(6):252. https://doi.org/10.3390/pharmaceutics11060252

Chicago/Turabian Style

Han, Jong Kwon, Beom Soo Shin, and Du Hyung Choi. 2019. "Comprehensive Study of Intermediate and Critical Quality Attributes for Process Control of High-Shear Wet Granulation Using Multivariate Analysis and the Quality by Design Approach" Pharmaceutics 11, no. 6: 252. https://doi.org/10.3390/pharmaceutics11060252

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop