Abstract

Background. Tablets are still the most preferred means of drug delivery. The search for new and improved direct compression tablet excipients is an area of research focus. This is because the direct compression method overcomes the drawbacks of granulation methods of tablet production. It exempts several treatment steps associated with the granulation methods. The requirements for the powders to be directly compressible include flowability, low friction tendency, compressibility, and fast disintegration capacity. Taro Boloso-I is a new variety of Colocasia esculenta (L. Schott) yielding 67% more than a previously reported variety (Godare) in Ethiopia. This study is aimed at enhancing the flowability while keeping the compressibility and compactibility of the pregelatinized Taro Boloso-I starch. Methods. Central composite design was used for the optimization of two factors which were the temperature and duration of pregelatinization against 4 responses. The responses were angle of repose, Hausner’s ratio, Heckel’s yield pressure, and tablet breaking force. Results and Discussions. An increase in the temperature resulted in decrease in both the angle of repose and the Hausner ratio and that of time decreased angle of repose as well. The Heckel yield pressure was observed to increase with increasing levels of both temperature and time. The pregelatinized starch prepared by heating 15% slurry of Taro Boloso-I starch at the pregelatinization temperature of 66.22°C for 20 min showed desired flow property and compressibility. Conclusions. Pregelatinized Taro Boloso-I starch could be regarded as a potential direct compression excipient in terms of flowability, compressibility, and compactibility. The PGTBIS could perform better as filler and binder in direct compression tablets than the Starch 1500® in terms of compactibility.

1. Introduction

Solid dosage forms, especially tablets, are still the most preferred means of drug delivery [1, 2]. The search for new and improved direct compression tablet excipients is an area of research focus. This is because the direct compression method overcomes the drawbacks of granulation methods of tablet production. It exempts several treatment steps associated with the granulation methods. The requirements for the powders to be directly compressible include flowability, low friction tendency, compressibility, and fast disintegration capacity [37]. While starch is one of the leading polymers of preference for use as a pharmaceutical excipient, it has notable limitations including poor flowability and compressibility. However, native Taro Boloso-I starch is unique in that it has appreciable compressibility properties [8] while having very poor flow properties [9]. Unless some physical modifications such as pregelatinization or chemical modifications are employed to enhance the flow property, many of the potential functionalities such as direct compression properties of the native starch could not be utilized [10, 11]. Pregelatinization is claimed to improve the flowability of starches [12, 13]. It has been reported that the direct compression properties of starches would be enhanced following changes in density, particle size and shape, and moisture content [14, 15] by pregelatinization better than acetylation [16]. The native Taro Boloso-I starch has low amylose to amylopectin ratio ( to ). Additionally, it has a high onset, peak, and endset temperatures of gelatinization, namely, 68.40°C, 75.46°C, and 84.40°C, respectively. It is comprised of polyhedral/angular granules existing in A-type polymorphism [9]. With the assumption that the process of pore reduction during compression of powders follows linear kinetics [17, 18], the Heckel model is the most accepted method for analyzing volume reduction with a pressure exerted on powders [1921]. It involves plotting tablet densification versus the applied pressure. The linear portion of the graph has slope and intercept both of which are material constants. It represents the plastic deformation of the powder. This region is obtained by considering the line of best fit. Most of the time, the line of best fit is in the range of the compression pressure 50-150 MPa. The reciprocal of its slope () is equal to the Heckel yield pressure (), the measure of the ability of the powder to deform plastically [19]. In the current study, the flow property and compressibility of pregelatinized Taro Boloso-I starch (PGTBIS) at different values of temperature and duration of pregelatinization were evaluated. The Heckel yield pressure () (using out-die-tablet Heckel model), angle of repose, compressibility index, and Hausner ratio were used as responses. The factors were optimized using the circumscribed central composite design (CCD) response surface methodology (RSM) [2224].

2. Materials and Methods

2.1. Materials

Taro Boloso-I was obtained from Areka Agricultural Research Institute, located at Areka (300 km South of Addis Ababa), Wolaita, Ethiopia; sodium hydroxide (BDH Poole Co, UK) and sodium chloride (Sörensen, Lauren, Denmark) were used as obtained.

2.2. Preparation of Pregelatinized Starch

Taro Boloso-I starch (NTBIS) was extracted as per the methods described by Balla et al. [9] and pregelatinized using a method used by Odeku et al. [12]. Accordingly, a total of 13 aqueous slurries (15%, ) of NTBIS were heated in a water bath at specified temperatures (61.86–89.14°C) with continuous heating and uniform stirring for specified periods of time (11.72–68.28 min). The pregelatinized starches were then dried at 40 for 48 h and powdered in a laboratory grinder (Pulverisette 2, Fritsch, Germany) so that all the powders passed through a 224 μm aperture sieve. Finally, the samples were stored separately in tightly sealed glass containers.

2.3. Determination of Size Distribution, Density, Compressibility, and Flow Properties

The size distribution of PGTBIS particles was determined by using leather light diffractometer (Malvern Mastersizer 2000, Germany). For determination of bulk and tapped densities, the standard cylinder and tapped densitometer (ERWEKA, type SVM, Germany), respectively, were used. Bulk and tapped densities were used for the calculation of Carr’s index and Hausner’s ratio [2526]. For the determination of true density (), xylene displacement method was used with equation (1). Tablet true density () was computed as the weighted average of the component true densities. For the study of flow properties, the angle of repose and flow rate were determined by using the funnel method. where , , , and stand for the weight of starch, the specific gravity of xylene (0.865), the weight of the flask filled with xylene, and the weight of flask with starch and xylene filled.

Heckel’s plots: a single punch machine (Korsch AG XP1 K0010288, Germany) equipped with a round flat-faced stainless steel die cavity with a diameter of 10 cm was used for the preparation of the compacts. In advance of the compression, the punch faces and the die wall were lubricated with magnesium stearate suspension in 95% ethanol. Eight tablets of  mg weight were prepared at each of 8 different compression pressures 38.22, 76.43, 114.65, 152.87, 191.08, 229.30, 267.52, and 305.58 MPa. Compact weight, tablet breaking force (TBF), diameter, and thickness of out-of-die tablets were measured 24 h after compression [22]. Compact density (tablet density, ) was calculated using its weight (), diameter (), and thickness (). Out-of-die tablet relative density () at any compression pressure was calculated as the ratio of compact density to true density and used for the calculation of the tablet porosity. The negative natural logarithm of tablet porosity was taken as its densification (equation (2)). Similarly, compression pressure was calculated from compression force in kN and tablet diameter (). Finally, graph of versus pressure (in MPa) was plotted to analyze compression properties of the powder compacts (equation (3)). The linear portion of the plot was obtained from the line of best fit using the Origin Software 7 (Origin LabTM Corporation, USA) such that the correlation coefficient was reasonably close to 1. Constant “” (extrapolated intercept of the straight portion of the Heckel plot) was used to calculate the relative density just before deformation () (equation (4)). per se, along with relative density at a point applied pressure equals zero (), was used to calculate the phase of densification at low pressure, (equation (5)) [22]. had been determined from the ratio of bulk density to true density.

2.4. Compactibility

The compactibility was assessed by using the method used elsewhere by Mitrevej et al. [18] and also by Mužíková and Kubíčková [27]. The mean tablet breaking force values of 300 mg of the pregelatinized Taro Boloso-I starch according to the study design were tested 24 hr after compression at 12 kN. The compression force 12 kN was selected for the optimization formulations because of two reasons. First, the densification of most of the compacts of the investigational starches including that of the starches modified at the central replica was observed to be plastic at or beyond 12 kN according to the Heckel plots besides it resulted in the strong compacts.

2.5. Optimization

Experimental designs: using circumscribed CCD, the effects of temperature and time of pregelatinization as factors on angle of repose, Hausner’s ratio, Heckel’s number, and TBF of the compacts prepared were evaluated. There were 13 experimental runs as shown in Table 1.

Polynomial regression algorithms were created using ANOVA output according to Design-Expert 8.0.7.1 (Stat-Ease, Corp. Australia) and used to predict the responses. Moreover, the factors were optimized by using numerical and graphical methods of optimization [23].

2.6. Validation and Comparison of Optimum Formulation

For validation of the models selected, tablet formulations at three points different from the design points were prepared, evaluated, and compared with the predicted values. For validation of the optimization result, the optimized product was prepared and compared as powder and in tablet formulations with the native starch.

2.7. Statistical Analyses

Statistical analyses including one-way ANOVA were applied by using Design-Expert 8.0.7.1 software (Stat-Ease, Corp. Australia) for the optimizations. The responses were expressed in polynomial models in terms of temperature and time of pregelatinization. All the results of the direct measurements were presented as arithmetic deviation. Origin version 7 (Origin LabTM Corporation, USA) was used to determine the lines of best fit. The target limit of the significance of statistical data was 95% CI.

3. Results and Discussion

3.1. Densities and Flow Properties

The density and flow properties of PGTBIS are presented in Table 2. The modifications increased bulk and tapped densities. However, it decreased true density most likely due to the hydrothermal disruption-related loosening of molecular packing and diffusion of amylose from the starch granules [28]. These effects were the highest at 89.12°C for 40 min as this temperature was higher than the endset temperature (84.40°C) [9] accompanied by complete disruption of the granular structure of the starch. The least decrease in true density was observed at 60.86°C for 40 min which shows that the temperature of gelatinization is associated directly with a decrease in its effect on the true density. The effects of true density reduction were lower also at 65°C for 20 min and 65°C for 60 min which are the temperature points below the onset temperature (68.40°C) where structural disruption of the starch granules is very insignificant. It gives the impression that the effect of gelatinization on the PGTBIS has some association with onset, peak, and endset temperatures of gelatinization which were 68.40°C, 75.46°C, and 84.40°C, respectively [9]. Generally, except for the points of pregelatinization at 60.86°C for 40 min and 65°C for 20 min, pregelatinization was shown to enhance flow property as it decreased the Carr index, Hausner ratio, and angle of repose and enabled the powder to flow. This was consistent with the published literature elsewhere [12].

3.2. Heckel’s Plots

Compressibility leads to compactibility which per se leads to tabletability [29, 30]. The compressibility of pregelatinized starches was studied by Heckel’s models. The results of Heckel’s plots and parameters were as shown in Figure 1 and Table 3, respectively.

Figure 1 shows that all the Heckel plots start with nonlinear portions indicating that the initial densification is predominantly due to the rearrangement of particles. The next segment for all of the plots is linear suggesting a region where the dominant mechanism of densification is plastic deformation. The relevant parameters of the initial nonlinear and linear regions including spans of pressure are described in Table 3. The Heckel plot of the powder of the design point 85°C for 60 min was observed to have ups and downs. This is a sign of brittle fracture as a mechanism of deformation explained by secondary rearrangement phases of particles [18].

The results in Table 3 imply that the starch modified at 65°C for 20 min, 85°C for 20 min, 89.12 °C for 40 min, 75°C for 11.72 min, and 75°C for 68.28 min attained maximum densification before compression (relative density ) within 0-38.22 MPa span of compression pressure. At 38.22 MPa onwards, the particles of these powders start to deform and possibly form bonding. All the rest of the points including the 5 replicates at 75°C for 40 min start particulate deformation beyond this pressure observed at the point 76.43 MPa onwards. Likewise, plastic deformation of the powder extends highest (229.30 MPa) for the modification at 65°C for 60 min and 61°C for 20 min followed by all the 5 central replicate points and 89.12°C for 40 min (191.08 MPa). The maximum pressure for the plastic deformation was the least (114.65MPa) for the starches modified at 85°C for 20 min, 75°C for 12 min, and 75°C for 68 min. The largest region (152.87 MPa) of plastic deformation was observed at 65°C for 60 min. The Heckel yield pressure () was calculated from the inverse of the slope () of the linear portions. Accordingly, the modification at 65°C for 20 min had the smallest Heckel yield pressure ( MPa) followed by 141.66 MPa at 60.12°C for 40 min. The explanation might be that the temperature was too low to change the native nature of the starch granules which are small enough to be compressible [31, 32]. This indicates that the pregelatinization of the starch at 65°C for 20 min results in the most plastic deforming PGTBIS as evidenced by literature [21]. On the other hand, the least compressible powder ( MPa) was observed at 85°C for 20 min followed by 250.31 MPa observed at 65°C for 60 min and 249 MPa observed at 75°C for 68.28 min. Constant is related to the particle rearrangement and die filling before deformation [22]. Constant “” is related to particle rearrangement and die filling just before deformation and the bonding of discrete particles. It was used to calculate the total densification hitherto (). According to the results, the powder at 60.86°C for 40 min showed the highest value for (1.4844), the least being 0.7205 observed at 89.12°C for 40 min. The corresponding relative density () values were observed to be 0.7734 and 0.5135 indicating the powders to be the most and the least compressible, respectively. The phase of densification at low pressure () was calculated from the value of the and the value of densification due to die filling (). The powder at 60.86°C for 40 min showed the largest (0.4368). Importantly, exhibits a phase of densification at low pressure and shows further densification due to particle fragmentation [18].

3.3. Compactibilities

The TBF and friability of compacts are described in Table 4.

The highest value of TBF was observed at 65°C for 20 min followed by the central replica indicating the compactibility. The modifications at 65°C for 20 min and 75°C for 12 min and all the 5 replicates at 75°C for 40 min had revealing appreciable cohesion of the compact particles [31, 33].

3.4. Analysis of Trends of Responses

Adequate mathematical models suitable for demonstrating the trends of responses graphically and mathematically were selected using the Design-Expert software. The values of angle of repose, Hausner’s ratio, Heckel’s yield pressure, and TBF values of the tablets were fed into the Design-Expert software. Then, for each of the responses, comparison of , adjusted , predicted , and predicted residual sum of square (PRESS) values of linear, two-factor interaction (2FI), quadratic, and cubic models were undertaken. All the responses were found to best satisfy the quadratic polynomial model except for Heckel’s yield pressure. The model of best fit for Heckel’s yield pressure was shown to be the two-factor interaction model. These respective models which were suggested to satisfy the best had greater values (0.9846, 0.9470, 0.9394, 0.9585, and 0.9949), adjusted values (0.9736, 0.9091, 0.8961, 0.9447, and 0.9913) closer to 1 and in more reasonable agreement with the predicted values (0.9064, 0.7707, 0.7707, 0.8780, and 0.9778), and smaller predicted residual sum of square (PRESS) values (29.51, 0.0034, 0.0033, 3517.29, and 4.09) than any other nonaliased models [23]. The adequacy of these models for predicting the influences of the factors on the responses was verified by employing the ANOVA test. The model significance tests were such that and lack-of-fit tests were insignificant () for all the five responses. Moreover, the adequate precision values, 36.86, 20.032, 12.564, 23.306, and 54.469 of the angle of repose, Hausner ratio, Heckel yield pressure, and TBF, respectively, showed that the signals are adequate and the respective models are quite valid. Two extra tests including the normal probability plots versus predicted values and the internally studentized residuals were also considered to prove the validity of the selected regression models. The normal probability plots of residuals versus predicted values reasonably approximated the normal plot of the predicted values. Similarly, the internally studentized residual values were <3 units away from zero and also randomly scattered for all the cases. These conditions further confirmed that the selected regression models were valid for the intended presentation of the trend of changes in responses in terms of the factors.

Having selected the models of best fits with statistically valid adequacy, equivalent mathematical equations were used to reveal the individual and interaction effects along with response surfaces. With this concern, equations (6) and (7) are the equations of angle of repose in terms of actual and coded factors, respectively, and used to predict the actual and relative impacts of each factor on angle of repose. Figure 2 shows the trend of angle of repose with changes in the temperature and time of pregelatinization diagrammatically. where and are the actual values of temperature and time, respectively. where and are the coded values of temperature and time, respectively.

Obviously, increase in both temperature and time decreases the angle of repose () and (), respectively. This indicates enhancement of flow property supported by published report elsewhere about the effect of the modification by Adedokun and Itiola [24]. Temperature had positive quadratic effect on the angle of repose with coefficients of 1.93 ().

For making prediction of the Hausner ratio for each of the levels of factors, the equation of the Hausner ratio in terms of actual values of factors is shown in equation (8). The same equation in terms of coded levels of the factors is given in equation (9). It is preferred for the estimation of the relative impact of each of the factors on the Hausner ratio. The trend of the Hausner ratio with changing levels of factors is diagrammatically depicted in Figure 3. where and are the actual values of temperature and time, respectively. where and are the coded values of temperature and time, respectively.

The increase in temperature decreases the Hausner ratio as opposed to the time of pregelatinization. With coefficients of 0.0234 () and -0.0241 () of temperature and time of pregelatinization, respectively, positive and negative quadratic effects on the Hausner ratio were observed.

The prediction of Heckel’s yield pressure could be made by using its equation in terms of the actual values of each of the factors (equation (10)). The relative impacts of each of the variables on the Heckel yield pressure could also be determined by equation (11) in which it is expressed in terms of coded levels of each of the variables. The 3D graphical demonstration of changes in the Heckel yield pressure with changing levels of the factors is presented in Figure 4. where , , and are the actual values of Heckel’s yield pressure, temperature, and time, respectively. where , , and are the coded values of Heckel’s yield pressure, temperature, and time, respectively.

The Heckel yield pressure was observed to increase with increasing levels of both temperature (coefficient of 25.94, ) and time (coefficient of 25.87, ). The interaction of the two factors antagonized the Heckel yield pressure with a factor of 55.44 ().

In the same fashion, the compactibility property with changing levels of actual values of the factors in the design space was predictable using equation (12). The relative impacts of each of the factors on the TBF were evident from its equation in terms of the coded levels of factors (equation (13)). The trend TBF with changing levels of the factors is diagrammatically presented in 3D (Figure 5). where , , and stand for tablet breaking force and actual values of temperature and time, respectively. where , , and stand for tablet breaking force and coded values of temperature and time, respectively.

Both temperature and time of pregelatinization exhibited negative effect on the square root of TBF with coefficients of -3.13 () and -2.71 (), respectively. The time of pregelatinization was observed to have a negative quadratic effect on the square root of TBF with a coefficient of -2.53 ().

Finally, the models selected were validated by pregelatinization of the same starch at three points other than the design points and comparing experimentally determined levels of responses with the predicted values (Table 5).

None of the responses in the selected data points showed as high as ±5% error with respect to the predicted value verifying that the models selected could demonstrate the trend of responses with changing levels of the factors within the design space.

3.5. Optimization

To optimize the factors towards the goals of the responses, the numerical (Figures 6 and 7(a)) and graphical (Figure 7(b)) methods of optimization were used [23].

The region of design space protruding upwards in Figure 7(a) and shaded with yellow color in Figure 7(b) represents the set of points fulfilling the criteria limited in the constraints. In these diagrams, the desirability flag and selected point show the point suggested as optimum with respect to all the responses considered. As a result, the optimum condition of pregelatinization was 66.22°C for 20 min resulting in the angle of repose of 28.95°, HR of 1.20, Heckel yield pressure of 101.35 MPa, and TBF of the tablet of 135.45 N with the desirability of 0.910 indicating that this point is the most desirable with respect to the goals of each of the responses.

3.6. Validation of the Optimization and Comparative Evaluation of Optimized Starch

In order to validate the optimization result, the optimized pregelatinized Taro Boloso-I starch (PGTBIS) was prepared by heating 15% slurry of the NTBIS at 66.22°C for 20 minutes. Then, the angle of repose, Hausner ratio, Heckel Yield pressure, and tablet breaking strength of the PGTBIS and the Starch 1500® were determined. The PGTBIS was compared with the NTBIS and Starch 1500® in terms of the bulk density, tapped density, true density, Hausner ratio, and Carr index which are presented in Table 6.

The particle size, [3, 4], of PGTBIS determined using laser light diffractometer (176.65 μm) was by far higher than that of NTBIS reported in our previous publication (μm) [9, 34]. Both bulk density ( g/ml) and tapped density ( g/ml) of PGTBIS were higher than NTBIS () may be due to changes in particle sizes and shapes [35]. Density values are related to total, interparticle, and intraparticle porosities which in turn indicate the properties of excipients like compressibility and tablet disintegrating potential and are also related to flow properties [31]. The higher bulk density is advantageous for it reduces fill volume during tableting. The angle of repose value of PGTBIS () was higher than that of the Starch 1500® () () indicating that it might have lower flowability. However, both the PGTBIS and Starch 1500® were within the range of excellent flowability index (25-30°) in USP/NF [36]. The Carr index of PGTBIS (%) was comparable to that of the Starch 1500® (%) () but noticeably lower than that of NTBIS (%) (). These values of Carr’s index gave a good promise regarding the compressibility and hence the tabletability of the powder [36]. The reason for low Carr’s index could be low interparticular interaction which favors compressibility. It was in the range of values claimed to be fair in flow according to the USP/NF classification of the flowability [37]. The reason for enhancement of the flow property of PGTBIS compared to that of the NTBIS might be due to increment of the particle size more than 100-folds as the experimental findings clearly show. This is consistent with literature [38]. Observation of the true density values determined for use in Heckel’s parameters has shown that the pregelatinization slightly decreased the true density from  g/ml to  g/ml (). This could possibly be owing to hydrothermal disruption of crystal structure and diffusion of amylose molecules out of the granules [28]. Accordingly, PGTBIS has statistically comparable density to that of Starch 1500® () (). The Heckel plot of the PGTBIS was constructed and is shown in Figure 8. The shape of the Heckel plot and parameters from its linear segment (slope, , and intercept, ) were used for the exploration of mechanisms of deformation, the extent of plastic deformation, and the influence of rearrangement or fragmentation during pressure-induced densification of powders. Initial densification of PGTBIS gives the impression that it is due to rearrangement or fragmentation of the powder particles revealed by the nonlinearity. The linear segment () of the plot suggests that the granules deform plastically in the range compression pressure 76.43-229.30 MPa. The Heckel yield pressure () was calculated from the inverse of the slope ( MPa-1) and equals 98.81 MPa with error of only 2.6% from the predicted value 101.45 MPa. The Heckel yield pressure reported in this study was lower than the Heckel yield pressure values of MCC, mannitol, lactose, and dicalcium phosphate reported elsewhere [39]. The implication is the PGTBIS powder could have better plasticity and compressibility. Moreover, the Heckel constant “” (1.1039) shows the total densification of the powder bed before deformation and corresponds to . The phase of densification by rearrangement before the particles of powder start individual deformation and possible bonding was found to be .

For the insight of compactibility profiles, TBF and friability of tablets compressed at 12 kN of pure PGTBIS, NTBIS, and the Starch 1500® were investigated (Table 7). The TBF of the pure PGTBIS tablet was  N which shows that the mean value was almost similar to the predicted value (135.4 N). The TBF value of the native starch compressed at the same pressure was  N which shows that the decrease in the compactibility of the Taro Boloso-I starch following pregelatinization was not statistically significant (). The tablets of Starch 1500® compressed at the same compression pressure had lower TBF value () than that of PGTBIS (). Besides, the friability of tablets of PGTBIS was lower than that of the Starch 1500® () and statistically comparable to the friability of compacts of the NTBIS (). The TBF and friability data observed in this study suggest that PGTBIS would have better compactibility than the Starch 1500®. The confirmation results verify that the analysis and optimization processes were well represented and valid in the mathematical models and the optimization. Since the pregelatinized starch has ensured to have good flowability and compactibility, it can be considered as a promising direct compression excipient for pharmaceutical tablet dosage forms. In comparison to the Starch 1500® as a filler and binder in direct compression tablets, the PGTBIS has comparable flow properties but better compactibility. The better compactibility might be because of formation of interparticular bridges to better extent with all other necessary preconditions such as presence of optimum moisture [40]. The PGTBIS could therefore be considered in the direct compression tablet formulations instead of the Starch 1500®.

4. Conclusions and Recommendations

An increase in both temperature and time of pregelatinization of Taro Boloso-I starch enhanced its flowability. The Heckel yield pressure was observed to increase with the increase in levels of both temperature and time. Both temperature and time of pregelatinization were found to have decreasing effects on compressibility/compactibility. Pregelatinization had a trade-off effect on compressibility and flowability. Heating at 66.22°C for 20 min was the optimum point of pregelatinization for balanced and acceptable flowability and compressibility/compactibility. At this point, sufficient flowability and compressibility/compactibility were achieved indicated by the angle of repose of , Hausner ratio of , Heckel yield pressure of 104.4 MPa, and TBF of 138.0 N at a compression force of 12 kN. Accordingly, it was discovered that the PGTBIS starch could be taken as a potential direct compression excipient. Its novelty is that the PGTBIS could perform better as filler and binder in direct compression tablets than the Starch 1500® in terms of compactibility. Further investigations including brittle fracture index, Young’s modulus, toughness, lubricant sensitivity, and dilution potential studies by using specific drugs should be performed towards its application in pharmaceutical industries.

Data Availability

The authors confirm that all the data associated with this paper are available upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Acknowledgments

The authors would like to thank Areka Town Administration and School of Pharmacy, Addis Ababa University, for this is part of the work presented to its institutional repository as a thesis [41] (although there are significant improvements). The authors also appreciate Areka Agricultural Research Center for the material donations and Ethiopian Pharmaceutical Manufacturing Share Company (EPHARM) for the access to laboratory facilities. The abstract of this article has also been shortlisted as the one to be presented in proceedings of the 2nd Global Conference on Pharmaceuticals and Clinical Research to be held in 2023, Dubai, UAE [42].