Optimization of Spray Drying Process Parameters for the Preparation of Inhalable Mannitol-Based Microparticles Using a Box-Behnken Experimental Design

Inhalation is used for local therapy of the lungs and as an alternative route for systemic drug delivery. Modern powder inhalation systems try to target the required site of action/absorption in the respiratory tract. Large porous particles (LPPs) with a size >5 μm and a low mass density (usually measured as bulk or tapped) of <0.4 g/cm3 can avoid protective lung mechanisms. Their suitable aerodynamic properties make them perspective formulations for deep lung deposition. This experiment studied the effect of spray-drying process parameters on LPP properties. An experimental design of twelve experiments with a central point was realized using the Box–Behnken method. Three process parameters (drying temperature, pump speed, and air speed) were combined on three levels. Particles were formed from a D-mannitol solution, representing a perspective material for lung microparticles. The microparticles were characterized in terms of physical size (laser diffraction), aerodynamic diameter (aerodynamic particle sizer), morphology (SEM), and densities. The novelty and main goal of this research were to describe how the complex parameters of the spray-drying process affect the properties of mannitol LPPs. New findings can provide valuable data to other researchers, leading to the easy tuning of the properties of spray-dried particles by changing the process setup.


Introduction
Pulmonary administration mediated by inhalation is an important non-invasive route of drug delivery [1,2]. This administration brings the possibility of both local and systemic effects using a wide range of different drugs [3,4]-from small molecules (beta-agonists or corticosteroids [5,6]) to large proteins (insulin [7]) or nucleic acids [1]. The lungs represent a target organ with several advantages: direct drug delivery to the site of action, bypassing the liver's first-pass metabolism, large absorption surface, rich blood supply, reduced enzymatic activity, rapid onset of action, high bioavailability, and lower effective dose needed, etc. [1,8,9]. Delivery of the drug in the form of dry powder is more beneficial than liquid formulations. Dry powder inhalers are propellant-free, low-cost, portable, and easy to operate (no need for hand-and-breath coordination), and the drug's dry state improves the formulation's stability [2,10,11]. Dry particle usage brings more versatility of their properties control via particle engineering (particle size, surface morphology, mass density, stability, aerodynamic properties, charge, etc.) [12][13][14].
The particles have to possess the required properties for lung deposition. Their physical and, more importantly, aerodynamic diameter (mass median aerodynamic diameter-MMAD), are crucial among the others. MMAD is characterized as the diameter of a spherical particle of unit density with the same settling velocity as the particle of interest, diabetes mellitus. Moreover, it may be incompatible with specific APIs (e.g., formoterol or proteins) as it interacts with their functional groups due to its reducing sugar function. For these reasons, the pharmaceutical industry is looking for alternative materials for dry powder formulation. Different carbohydrates (trehalose, raffinose), polyols (mannitol, sorbitol, erythritol), or polymers (chitosan, PLGA) have been investigated [10]. For this experimental work, mannitol was chosen for its non-animal origin and lower hygroscopicity than lactose. It is a non-reducing polyalcohol with a sweet taste, which could indicate proper dosing. It is also approved as a GRAS substance by the FDA and the European Regulatory Committee [32]. It is currently marketed for diagnostic (Aridol TM ) and therapeutic (Bronchitol TM ) purposes as a dry powder inhaler [10].
Nowadays, the optimization of process/formulation parameters represents an important step in research, and the design of experiment (DoE) approach is getting great attention in pharmaceutical technology [33]. Response surface methodology (RSM) is widely used to understand the relationship between input variables and response(s) or to find optimal process settings while simultaneously achieving high efficiency and low costs. A suitable tool for RSM application when having 2-4 independent factors is the Box-Behnken design (BBD) [34], which requires fewer factor combinations than a central composite design and is rotatable or nearly rotatable, still sufficient to fit a quadratic model. BBD is also beneficial when avoiding the extreme factor level combinations is desirable.
The influence of process parameters on mannitol particles properties was investigated using different DoE approaches and statistical methods. For example, Littringer et al. used a full-factorial design and stated that higher temperatures lead to a smoother surface and lower feed concentration. High atomizer rotation speeds created smaller particles [35].
Kramek-Ramanowska et al. concluded that the minimum particle size is mainly determined by a low drying temperature and a high gas flow rate through the nozzle using the same DoE type [36]. According to Li et al., aerosol performance is better with increasing pump speed rate based on DoE analysis [37]. With the application of BBD, Guimarães et al. showed the effect of outlet temperature and feed concentration on particle size and shape [38].
The present study investigated the influences of three less-discussed process parameters of spray drying (drying temperature, pump speed, and air speed) on the various responses (properties of LPPs, especially particle size characteristics). A BBD and analysis of DoE data were employed to provide the best mathematical models of the mentioned dependencies using multiple linear regression (MLR) and to visualize results via perspective plots. This study's main goals and novelty were to describe the influences of studied process variables on mannitol LPP features and to provide usable new findings to other research groups trying to tune the properties of spray-dried particles by changing the process setup.

Spray Drying Process
A LabPlant spray dryer (SD-06, LabPlant, Filey, UK) with a drying chamber (215 mm in diameter and 500 mm long) was used for particle formulation. A spray solution of 10% (w/w) D-mannitol (Penta, CZ; Mr = 182.18) in distilled water was prepared in laboratory conditions at temperatures up to 25 • C. For atomization, a two-fluid, 2 mm nozzle was used. The drying conditions for the experiment were combined into three levels ( Table 1). The inlet temperature of the drying gas was set to 100 • C, 120 • C, or 140 • C, respectively. The feed rate of the primary solution was determined by pump speeds of 5 mL/min, 10 mL/min, and 15 mL/min, respectively. The last process parameter air speed was 3.0 m/s, 3.6 m/s, and 4.2 m/s. Constant atomizing pressure was set to 300 kPa. The final product was separated through a cyclone and stored in glass bottles in an exicator for a maximum of 10 days before evaluation.

Particle Size and Aerodynamic Diameter
Each sample was dispersed in 96% ethanol. The suspension was analyzed by laser diffraction (LA-960, Horiba, Kyoto, Japan). Before measuring, a small amount (<1 g) of each sample was dispersed in ethanol, sonicated for 30 s, and stirred for maximal homogeneity of the suspension. Refractive index settings were: ethanol RI = 1.360, mannitol RI = 1.333, and imaginary index iRI = 0, respectively. Measurement was performed according to the recommended settings [39]. Each measurement of the laser diffraction median size (S, µm) was performed in three repetitions and presented as a mean with geometric standard deviation.
Aerodynamic properties were measured by an APS (aerodynamic particle sizer-TSITM model 3321, TSI Incorporated, St. Paul, MN, USA). The method is based on measuring the time of flight of a particle over a well-defined distance between two lasers with detectors. Firstly, the sample is dispersed by the compressed air using the TSI small-scale powder disperser 3433 USA (set parameters: air speed = 5 L/min, air pressure = 200 kPa). Then, the sample goes through the nozzle, where the flow of the particles is accelerated. The velocity of the particles is measured in time in the detection area. Measurement was performed according to the recommended settings [40]. Mass median aerodynamic diameter (MMAD, µm) and numeric median aerodynamic diameter (NMAD, µm) were detected. Each measurement was performed in three repetitions and presented as a mean with geometric standard deviation.

Bulk and Tapped Density, Hausner Ratio
For density measurements, 5 mL of each sample was used. The bulk density was calculated from the weight and volume according to the European Pharmacopeia. Tapped density was measured by Erweka ® model SVM102 (DE). The volumes of samples were reported after 10, 500, 1250, and 2500 taps. The Hausner ratio was calculated using tapped (after 2500 taps) and bulk density values [41].

Particle Morphology
The morphology and surface structure of mannitol LPPs were characterized by scanning electron microscopy (SEM). Samples were placed onto aluminum stubs with doublesided adhesive carbon tape, coated with a 10 nm gold layer (sputter Q150R, Quantum Technologies, London, UK), and visualized using an SEM (MIRA3, Tescan, Brno, Czechia). Obtained signals of the samples were produced by secondary electrons (SE) at 5 kV voltage and 5 kx magnification. Particles were classified on a three-point scale as G-good, A-acceptable, and B-bad ( Table 2). The integrity of the surface, spherical shape, uniformity of size distribution, and presence of fragments were monitored. Spherical, integral, uniformly wrinkled, and size-uniform particles are marked as good-G. Oval, integral with damages, non-homogeneously wrinkled, or partially wrinkled particles with larger differences in size are marked as acceptable-A. Particles with an irregular shape, a high proportion of damaged particles, wide size distribution, and a high fragment content are marked as bad-B. To build a mathematical model, the morphology was expressed on the coded scale 1-3 (3 = G-good, 2 = A-acceptable, 1 = B-bad). Despite the discrete nature of the obtained response, the assumption is that the particle morphology grade for various process parameter settings within the DoE region can take on values approximately in the range of 1-3.

Experimental Design
The three-factorial and a three-level BBD with three replications of the center point were applied for DoE, resulting in 15 experimental runs to assess the effect of input parameters on various product quality attributes. The process parameters: drying temperature (X 1 ; in the range of 100-140 • C), pump speed (X 2 ; in the range of 5-15 mL/min), air speed (X 3 ; in the range of 3.0-4.2 m/s) were selected as studied input variables. Each factor was divided into three levels coded as low (−1), medium (0), and high (+1), as listed in Table 1, which were chosen based on previous experiments with the production of LPPs used in experiments on an artificial lung model to simulate particle deposition during different breath cycles [42].
Utilized BBD is a cubic design, where the 12 experimental runs are represented by the mid-point of each edge of the 3D cube, and the 13th point is the center point coded as a set of factors: 0, 0, 0. The experimental runs were carried out randomly to minimize the effect of unexplained variability on observed responses. Replications of the center point were added to ensure process stability and check the model curvature. They were dispersed as evenly as possible throughout the BBD matrix (the 1st, the 8th, and the 15th experimental run). The measured characteristics of microparticles in the present study were S (Y 1 , µm), MMAD (Y 2 , µm), NMAD (Y 3 , µm), relative standard deviation (RSD) of MMAD (Y 4 , %), particle morphology characterized by SEM (Y 5 ; in the coded scale 1-3; higher grade is better), and flow properties (bulk density, tapped density, and Hausner ratio).

Data Analysis
DoE analysis in this study included the following steps (generalized procedure):

2.
Choosing a suitable mathematical model for the obtained data.

3.
Testing the assumptions required for ANOVA and regression analysis (MLR), primarily using visual assessment of graphical outputs (histograms, residual graphs); eventually, exclusion of the outliers from the analysis and subsequent building of a new model.

4.
Simplification of the model by gradual backward elimination of insignificant terms (assessment of their p-values in ANOVA table) while monitoring and comparing values of coefficient of determination R 2 , adjusted R 2 , Akaike information criterion, and p-value of the original and simplified model to achieve the best quality of the model fit. Percentage prediction error (PPE) was assessed to ensure the validity of the generated regression equation. PPE for each experimental run was calculated as: (observed value-predicted value)/predicted value × 100 (%).

6.
Interpretation of obtained regression equation and graphical outputs: investigation of the effects of independent factors (process conditions) and their possible interactions on the response (quality attributes of microparticles) and the determination of optimal process conditions for the desired response.
The results of MLR were fitted to a quadratic polynomial model given by equation (Equation (1)): where Y stands for the estimated response; X 1 , X 2 , and X 3 are independent variables; b 0 is an intercept/constant; b 1 , b 2 , and b 3 are regression coefficients for linear terms; b 12 , b 13 , and b 23 are regression coefficients for interaction/cross product terms; b 11 , b 22 and b 33 are regression coefficients for quadratic terms. Standardized regression coefficients (β) were also determined to compare the influence of individual terms on the response. Data analysis was performed with a significance level of α = 0.05, so the effects at p < 0.05 were considered statistically significant. The design of the experiment and subsequent analysis of the obtained data were carried out in R software, version 4.2.1 [43].

Results and Discussion
All results obtained are presented in Table 2. Measured were: laser diffraction median size value (Y 1 ; S, µm) and mass median aerodynamic diameter (Y 2 ; MMAD, µm) as the most pronounced parameters in the research. Numerical median aerodynamic diameter (Y 3 ; NMAD, µm) provides important information about the numerical distribution of particles. Y 4 represents the RSD of MMAD. Finally, scanning electron microscopy was performed as an important indicator of surface and morphological properties (Y 5 ; SEM, values: G-good, A-acceptable, B-bad). Other measured values are bulk density (g/cm 3 ), tapped density (g/cm 3 ), and Hausner ratio.
To assess the combined effects of three factors (drying temperature, pump speed, and air speed) on the various responses (particle size parameters and morphology), a BBD of 15 experimental runs in conjunction with subsequent data analysis was used. By performing multiple regression analysis on the experimental data, a second-order polynomial equation for each response variable was obtained (expressed in tabular arrangement). The resulting MLR models, including regression coefficients for each term, corresponding p-values, and selected model characteristics (R 2 , adjusted R 2 , p-value), are summarized in Table 3. The mean values of each response variable for each experimental run, as well as the values predicted by MLR models (from Table 3) and errors of prediction, are reported in Table 4. This evaluation is presented for quantities for which a significant influence of the process parameters on the final value was revealed using statistical testing. The raw data, with indicated outliers excluded from the subsequent analysis, are presented in Table S1. Table 3. MLR models for selected variables: estimated regression coefficients (b) with p-values, standardized regression coefficients (β) and the basic model characteristics.

Laser Diffraction Median Size-S (Y 1 )
The MLR model of S was built using 15 experimental runs with three repeated measurements. Three outliers were excluded from the analysis based on the residual plots. The quadratic model, including all terms, was chosen as the best-fitting model. The R 2 value of 0.813 means that the model can clarify 81.3% of the variability, and the associated p-value of less than 0.001 indicates the model is highly significant (Table 3). PPE for most runs is in the range of several units up to tens of percent, except for run 11-in this DoE region, the predictive ability of the model sharply decreases (extremely high PPE value, 748.9%). The mean PPE value is relatively high (67.6%). If run 11 is not considered, the mean PPE value is 18.9%, which is acceptable. Based on the mentioned characteristics, the model fitting can be considered satisfactory for response determination, especially in identifying independent variables' effects on the S value. Based on the standardized regression coefficient values (Table 3, β values), drying temperature and air speed (linear and quadratic terms in both cases) have the most considerable effect on the response. Pump speed and interactions between variables also contribute to the model (Figure 1). model fitting can be considered satisfactory for response determination, especially in identifying independent variables' effects on the S value. Based on the standardized regression coefficient values (Table 3, β values), drying temperature and air speed (linear and quadratic terms in both cases) have the most considerable effect on the response. Pump speed and interactions between variables also contribute to the model (Figure 1). The interrelationships between process parameters and S, closely defined by the quadratic equation, are depicted via perspective plots, where the effects of two input variables and their interaction at the middle level of the third variable are visualized. The middle level of the third input variable is marked as "slice at" in the graph, and it is calculated as the mean for all measurements after excluding outliers; therefore, it may not exactly correspond to the medium level (0) from Table 1. It is further described, like other monitored microparticle properties, in the Results and Discussion section. It could be seen that there was a considerable variation in the S value with different process parameter settings (Figure 1). The perspective plot in Figure 1A shows that S rapidly increased up to the values above 30 μm with increasing drying temperature (X1) or pump speed (X2), as well as when combining the high values of these variables.
This observation follows the Broadhead et al. study, which suggested that the material's rapid agglomeration can cause the formation of larger particles at higher temperatures [44]. In the study published by Stahl et al., it was observed that increasing temperature significantly influences physical particle size. High inlet temperature leads to the formation of the largest particles. Another finding of this study was that particle size decreased when air speed increased. This effect is caused because higher air speed produces more energy for aerosolization (conversion of the primary dispersion into the droplets), resulting in smaller droplets [45]. This finding was confirmed partly in our study due to the particles' physical size affection by other process parameters (not set on the same level in the compared study).
As shown in Figure 1, using the low to medium values of X1 and X2 seems to be the most appropriate for achieving particles of around 10 μm. The saddle shape of the response surface in Figure 1B,C suggests a more complex dependency, where low or conversely high air speed (X3) in combination with low drying temperature (X1) or low pump speed (X2) appears to be the experimental settings leading to the S less than 10 μm. The experimental region around the center point results in a response slightly above 10 μm and thus represents the optimal conditions for preparing LPPs suitable for lung deposition.
From the measured sample set, microparticles with the optimal S value are provided by central point-experimental run 13 (mean and SD for 13a, 13b and 13c − S: 13.34 ± 4.26 The interrelationships between process parameters and S, closely defined by the quadratic equation, are depicted via perspective plots, where the effects of two input variables and their interaction at the middle level of the third variable are visualized. The middle level of the third input variable is marked as "slice at" in the graph, and it is calculated as the mean for all measurements after excluding outliers; therefore, it may not exactly correspond to the medium level (0) from Table 1. It is further described, like other monitored microparticle properties, in Section 3. It could be seen that there was a considerable variation in the S value with different process parameter settings (Figure 1). The perspective plot in Figure 1A shows that S rapidly increased up to the values above 30 µm with increasing drying temperature (X 1 ) or pump speed (X 2 ), as well as when combining the high values of these variables.
This observation follows the Broadhead et al. study, which suggested that the material's rapid agglomeration can cause the formation of larger particles at higher temperatures [44]. In the study published by Stahl et al., it was observed that increasing temperature significantly influences physical particle size. High inlet temperature leads to the formation of the largest particles. Another finding of this study was that particle size decreased when air speed increased. This effect is caused because higher air speed produces more energy for aerosolization (conversion of the primary dispersion into the droplets), resulting in smaller droplets [45]. This finding was confirmed partly in our study due to the particles' physical size affection by other process parameters (not set on the same level in the compared study).
As shown in Figure 1, using the low to medium values of X 1 and X 2 seems to be the most appropriate for achieving particles of around 10 µm. The saddle shape of the response surface in Figure 1B,C suggests a more complex dependency, where low or conversely high air speed (X 3 ) in combination with low drying temperature (X 1 ) or low pump speed (X 2 ) appears to be the experimental settings leading to the S less than 10 µm. The experimental region around the center point results in a response slightly above 10 µm and thus represents the optimal conditions for preparing LPPs suitable for lung deposition.

Mass Median Aerodynamic Diameter-MMAD (Y 2 )
BBD, including 15 experimental runs with 3 repeated measurements, was used to obtain an MLR model of MMAD. Based on the residual plots, three outliers were excluded from the data analysis. MLR proposes the quadratic model, including all terms, as the best-fitting model. The R 2 value, which is found to be 0.810, suggests that the model can explain 81.0% of the total variability. The detected p-value of less than 0.001 implies that  (Table 3). The low mean PPE (1.9%) indicates the high predictive ability of the MLR model. Considering these facts, it can be concluded that the resulting model is adequate for MMAD response prediction under various experimental conditions within the range of the analyzed BBD. It can be seen from Table 3 that the importance of the effect of individual process parameters on MMAD decreases in the following order: drying temperature > air speed > pump speed (applies to both linear and quadratic terms of the mentioned variables).
The interrelationships between the dependent variable and independent parameters were further elucidated using perspective plots. Response surface as the combination of two process parameters by considering the third parameter at a middle level is presented in Figure 2. MMAD dropped to a value around 7 µm when the drying temperature (X 1 ) and air speed (X 3 ) decreased and pump speed (X 2 ) increased, while under other conditions within the experimental design it reached values up to 8.2-8.4 µm. Especially high X 1 values (Figure 2A,B) and a combination of X 2 low level and X 3 high level ( Figure 2C) resulted in a quite high MMAD value. On the other hand, the center point and the region around (0, 0, 0) correspond to the appropriate process conditions leading to the small particles. best-fitting model. The R value, which is found to be 0.810, suggests that the model can explain 81.0% of the total variability. The detected p-value of less than 0.001 implies that the regression equation is significant ( Table 3). The low mean PPE (1.9%) indicates the high predictive ability of the MLR model. Considering these facts, it can be concluded that the resulting model is adequate for MMAD response prediction under various experimental conditions within the range of the analyzed BBD. It can be seen from Table 3 that the importance of the effect of individual process parameters on MMAD decreases in the following order: drying temperature > air speed > pump speed (applies to both linear and quadratic terms of the mentioned variables).
The interrelationships between the dependent variable and independent parameters were further elucidated using perspective plots. Response surface as the combination of two process parameters by considering the third parameter at a middle level is presented in Figure 2. MMAD dropped to a value around 7 μm when the drying temperature (X1) and air speed (X3) decreased and pump speed (X2) increased, while under other conditions within the experimental design it reached values up to 8.2-8.4 μm. Especially high X1 values (Figure 2A,B) and a combination of X2 low level and X3 high level ( Figure 2C) resulted in a quite high MMAD value. On the other hand, the center point and the region around (0, 0, 0) correspond to the appropriate process conditions leading to the small particles.  Table 2). Values were achieved at medium and low process parameter values for air speed and drying temperature, while pump speed was at a middle and high level, except in run 9, where the pump speed was at level −1.
Maltesen et al. used a similar experimental setup in their study aimed at inhalation insulin particle preparation. Process parameters (drying temperature, pump speed, and air speed) were set at three levels. Like our experiment, only small differences in aerodynamic diameters between samples were observed, a range of 2.57-3.71 µm (compared against one level of insulin concentration. It was selected according to the largest number of samples for a more accurate comparison). However, the lowest values were also achieved with medium settings of the process parameters [46].

Numeric Median Aerodynamic Diameter-NMAD (Y 3 )
For NMAD, the MLR model using 15 experimental runs of BBD with 3 repetitions was calculated. The exclusion of two outlier measurements was based on the residual plot indication. The MLR method provided a prediction equation including a quadratic term only for pump speed as a statistically significant coefficient; the quadratic effect of other variables (drying temperature and air speed) and their mutual interaction have not been assessed as significant by ANOVA. The quality of the resulting model was evaluated by R 2 value, p-value, and PPE. The coefficient of determination (0.679) does not indicate a very strong goodness of fit. In contrast, the mean value of PPE is low enough (4.2%), which, on the contrary, indicates a good prediction ability ( Table 3). The overall model p-value of less than 0.001 demonstrates that the regression equation is significant. Therefore, the resulting model can still be used to determine the dependencies occurring in the data matrix. Based on the regression coefficient values listed in Table 3, it can be deduced that the quadratic term of the pump speed represents the greatest influence on the NMAD. The interaction terms and linear terms of drying temperature and air speed have a lower effect on the response, but they are still significant.
Response surfaces in Figure 3 were plotted based on the quadratic model, with all terms as the combination of two independent variables at a middle level of the third one. The shape of the response surface in Figure 3A indicates that NMAD decreases with increasing drying temperature (X 1 ) at any value of pump speed (X 2 ). Still, this effect is most pronounced at X 2 medium level. A similar dependence can be observed in Figure 3C, where there was a descending response with the increment in air speed (X 3 ) at various X 2 . Still, again, the steepest descent manifested at X 2 medium level. As shown in Figure 3B, a visual analysis of the depicted response surface plot indicates that the NMAD decreases as the X 1 and X 3 increase. For all response surfaces, the NMAD values decreased to approximately 3 µm while the upper parts of the response surface were close to 4 µm. Although the dependences determined by MLR, in this case, differ from the response surfaces for the variable Y 2 (compare Figures 2 and 3), the center point and the region around (0, 0, 0) can again be evaluated as a suitable combination of process parameters resulting in an optimal response. variables (drying temperature and air speed) and their mutual interaction have not been assessed as significant by ANOVA. The quality of the resulting model was evaluated by R 2 value, p-value, and PPE. The coefficient of determination (0.679) does not indicate a very strong goodness of fit. In contrast, the mean value of PPE is low enough (4.2%), which, on the contrary, indicates a good prediction ability ( Table 3). The overall model pvalue of less than 0.001 demonstrates that the regression equation is significant. Therefore, the resulting model can still be used to determine the dependencies occurring in the data matrix. Based on the regression coefficient values listed in Table 3, it can be deduced that the quadratic term of the pump speed represents the greatest influence on the NMAD. The interaction terms and linear terms of drying temperature and air speed have a lower effect on the response, but they are still significant.
Response surfaces in Figure 3 were plotted based on the quadratic model, with all terms as the combination of two independent variables at a middle level of the third one. The shape of the response surface in Figure 3A indicates that NMAD decreases with increasing drying temperature (X1) at any value of pump speed (X2). Still, this effect is most pronounced at X2 medium level. A similar dependence can be observed in Figure 3C, where there was a descending response with the increment in air speed (X3) at various X2. Still, again, the steepest descent manifested at X2 medium level. As shown in Figure 3B, a visual analysis of the depicted response surface plot indicates that the NMAD decreases as the X1 and X3 increase. For all response surfaces, the NMAD values decreased to approximately 3 μm while the upper parts of the response surface were close to 4 μm. Although the dependences determined by MLR, in this case, differ from the response surfaces for the variable Y2 (compare Figures 2 and 3), the center point and the region around (0, 0, 0) can again be evaluated as a suitable combination of process parameters resulting in an optimal response. The experimental NMAD results ranged from 2.68 to 4.14 μm (Table 2), which can be considered satisfying. However, RSDs of NMAD reached enormous values (over 100%) in most cases, making these results inaccurate. However, the spray drying method generates polydisperse particles, which correspond to higher values of NMAD RSD [47]. The experimental NMAD results ranged from 2.68 to 4.14 µm (Table 2), which can be considered satisfying. However, RSDs of NMAD reached enormous values (over 100%) in most cases, making these results inaccurate. However, the spray drying method generates polydisperse particles, which correspond to higher values of NMAD RSD [47].

Relative Standard Deviation of Mass Median Aerodynamic Diameter-MMAD RSD (Y 4 )
An MLR model of MMAD RSD was built on BBD with 15 experimental runs and 3 repeated measurements. Based on the residual plots, two outliers were excluded from the data analysis, and another three measurements that were excluded during MMAD model development were not included in the calculations either. A quadratic model including all terms was generated and assessed as the most accurate and meaningful model. The R 2 value of 0.775 implies that 77.5% of the variability can be interpreted by the model, which is acceptable. The statistical significance of the model was confirmed by the p-value (p < 0.001) ( Table 3). Predicted errors of all runs within 5% (the mean value of 2.2%) demonstrate the excellent fit of the MLR model. All the results mentioned above clearly show that the quality of the model is satisfactory and reliable, especially for examining data dependencies in the DoE region. As presented in Table 3, the linear and quadratic coefficients of drying temperature were the terms that most affected the response. Other quadratic terms and interactions between pump speed and air speed contribute less to the model but are still statistically significant, as can also be inferred from

Relative Standard Deviation of Mass Median Aerodynamic Diameter-MMAD RSD (Y4)
An MLR model of MMAD RSD was built on BBD with 15 experimental runs and 3 repeated measurements. Based on the residual plots, two outliers were excluded from the data analysis, and another three measurements that were excluded during MMAD model development were not included in the calculations either. A quadratic model including all terms was generated and assessed as the most accurate and meaningful model. The R 2 value of 0.775 implies that 77.5% of the variability can be interpreted by the model, which is acceptable. The statistical significance of the model was confirmed by the p-value (p < 0.001) ( Table 3). Predicted errors of all runs within 5% (the mean value of 2.2%) demonstrate the excellent fit of the MLR model. All the results mentioned above clearly show that the quality of the model is satisfactory and reliable, especially for examining data dependencies in the DoE region. As presented in Table 3, the linear and quadratic coefficients of drying temperature were the terms that most affected the response. Other quadratic terms and interactions between pump speed and air speed contribute less to the model but are still statistically significant, as can also be inferred from Figure 4. Since the regression model had three independent variables, one variable was fixed at a constant at a middle level for each perspective plot ( Figure 4). Values of MMAD RSD range in the tens of percent and increase to a maximum of about 40%, especially in the area around the central point, as can be seen from the shape of the response surface in Figures 4A,B. With a gradual increase in drying temperature (X1) at any pump speed (X2) or any air speed (X3), a decrease in MMAD RSD to values approaching 30% can be observed. From Figure 4C, an interaction between X2 and X3 can be deduced. High values of MMAD RSD are achieved, especially at high levels of both independent variables. A decrease in MMAD RSD can be obtained with a combination of experimental conditions: high X2 at the low level of X3 or, conversely, low X2 at the high level of X3. Therefore, the region around the central point (0, 0, 0) appears less suitable for particle size variability.
Experimentally, MMAD RSD is in the range of 30.0-42.5%, representing corresponding values for the spray drying method. Spray drying, like most methods, leads to the formation of polydisperse particles with a variety of distributions in their sizes. This parameter is, therefore, important just in the overall spectrum of evaluation. The most pronounced parameters for LPPs are laser diffraction median size, mass median aerodynamic size, and particle morphology [48]. Since the regression model had three independent variables, one variable was fixed at a constant at a middle level for each perspective plot ( Figure 4). Values of MMAD RSD range in the tens of percent and increase to a maximum of about 40%, especially in the area around the central point, as can be seen from the shape of the response surface in Figure 4A,B. With a gradual increase in drying temperature (X 1 ) at any pump speed (X 2 ) or any air speed (X 3 ), a decrease in MMAD RSD to values approaching 30% can be observed. From Figure 4C, an interaction between X 2 and X 3 can be deduced. High values of MMAD RSD are achieved, especially at high levels of both independent variables. A decrease in MMAD RSD can be obtained with a combination of experimental conditions: high X 2 at the low level of X 3 or, conversely, low X 2 at the high level of X 3 . Therefore, the region around the central point (0, 0, 0) appears less suitable for particle size variability.
Experimentally, MMAD RSD is in the range of 30.0-42.5%, representing corresponding values for the spray drying method. Spray drying, like most methods, leads to the formation of polydisperse particles with a variety of distributions in their sizes. This parameter is, therefore, important just in the overall spectrum of evaluation. The most pronounced parameters for LPPs are laser diffraction median size, mass median aerodynamic size, and particle morphology [48].

Particle Morphology Characterized by SEM (Y 5 )
For particle morphology, an MLR model based on 15 experimental runs of BBD was established. In terms of quadratic coefficients, the MLR method provided a prediction equation that included a statistically significant quadratic term for only air speed. In contrast, the quadratic effect of pump speed and drying temperature was insignificant according to the ANOVA table. Model competency was validated using the following parameters: R 2 value, p-value, and PPE. The value of the determination coefficient (0.866) shows that the model fitting is relatively good. The mean PPE (13.7%) indicates acceptable magnitudes of differences between experimental and calculated data ( Table 3). The results may also be affected by the relatively low number of experimental values (absence of repeated measurements for each run), given the nature of the discussed microparticles' properties. The model p-value of 0.012 proves the statistical significance of the identified dependencies. Therefore, the found equation can be used to display the regression equation and to explain these relationships. From Table 3, it can be seen that the air speed has a considerable influence on the response. The lower magnitude of other terms indicates their lower effects; however, some are still statistically significant. Figure 5 displays the surface profiles as the response of varying values of two independent factors at a fixed level (0) of the third one, whereas the quadratic model with all terms was considered. Better particle morphology can be achieved at the middle level of pump speed (X 2 ) and air speed (X 3 ), irrespective of the drying temperature (X 1 ), as can be seen in Figure 5A,B. From Figure 5C, it can be interpreted that the combination of medium-to-high levels of X 2 and X 3 leads to the maximum particle morphology grade, shown in the plot as a peak value of 3. Thus, it can be concluded that the experimental conditions corresponding to the region around the center point (0, 0, 0) result in the microparticles with the best particle morphology. equation that included a statistically significant quadratic term for only air speed. In contrast, the quadratic effect of pump speed and drying temperature was insignificant according to the ANOVA table. Model competency was validated using the following parameters: R 2 value, p-value, and PPE. The value of the determination coefficient (0.866) shows that the model fitting is relatively good. The mean PPE (13.7%) indicates acceptable magnitudes of differences between experimental and calculated data ( Table 3). The results may also be affected by the relatively low number of experimental values (absence of repeated measurements for each run), given the nature of the discussed microparticles' properties. The model p-value of 0.012 proves the statistical significance of the identified dependencies. Therefore, the found equation can be used to display the regression equation and to explain these relationships. From Table 3, it can be seen that the air speed has a considerable influence on the response. The lower magnitude of other terms indicates their lower effects; however, some are still statistically significant. Figure 5 displays the surface profiles as the response of varying values of two independent factors at a fixed level (0) of the third one, whereas the quadratic model with all terms was considered. Better particle morphology can be achieved at the middle level of pump speed (X2) and air speed (X3), irrespective of the drying temperature (X1), as can be seen in Figure 5A,B. From Figure 5C, it can be interpreted that the combination of medium-to-high levels of X2 and X3 leads to the maximum particle morphology grade, shown in the plot as a peak value of 3. Thus, it can be concluded that the experimental conditions corresponding to the region around the center point (0, 0, 0) result in the microparticles with the best particle morphology. The experimental finding showed that samples 7 (X1: −1, X2: 0, X3: +1), 9 (X1:0, X2: −1, X3: −1), 12 (X1: 0, X2: +1, X3: 01), 13 a, b, c (X1: 0, X2: 0, X3: 0) have process parameters set to provide particle sizes S (parameter Y1, μm) in the range 11.08 ± 3.65 μm-13.34 ± 4.26 μm, which is the closest to ten targeted micrometers. Particle morphology provided by SEM pointed out spherical, polydisperse particles with minimal damage. These particles are marked as good (G) except for sample 9. In Figure 6, sample 9 (X1: 0, X2: −1, X3: −1) spherical particles with a larger size distribution, locally wrinkled with few holes present, and very few fragments can be seen. Similar observations were observed in samples 2 (X1: +1, X2: −1, X3: 0), 4 (X1: +1, X2: +1, X3: 0), and 8 (X1: +1, X2: 0, X3: +1). These particles are marked as acceptable-A. Spherical polydisperse particles, with a small presence of holes, with suitable shapes were observed in samples 1 (X1: −1, X2: −1, X3: 0), 3 (X1: −1, X2: +1, X3: 0), 6 (X1: Figure 5. Three-dimensional (3D) perspective plots: response surface of particle morphology (Y 5 , the scale 1-3, the higher grade is better) as a function of drying temperature (X 1 , • C), pump speed (X 2 , mL/min), and air speed (X 3 , m/s).
Pharmaceutics 2023, 15, x FOR PEER REVIEW 13 of 16 irregular shape, a high proportion of damage, wide size distribution, a high fragment content, and a tendency to agglomerate can be observed. These are marked as bad (B). The influence of drying temperature, pump speed, and air speed on bulk, tapped density, and Hausner ratio (and mutual interactions of process factors) have not been confirmed as statistically significant using ANOVA (all p > 0.05).
All samples showed densities < 0.4 g/cm 3 , which fulfills the requirement for LPPs according to current knowledge ( Table 2). Bulk density values ranged between 0.25 and 0.31 g/cm 3 and tapped densities were between 0.34 and 0.41 g/cm 3 . These small differences The influence of drying temperature, pump speed, and air speed on bulk, tapped density, and Hausner ratio (and mutual interactions of process factors) have not been confirmed as statistically significant using ANOVA (all p > 0.05).
All samples showed densities < 0.4 g/cm 3 , which fulfills the requirement for LPPs according to current knowledge (Table 2). Bulk density values ranged between 0.25 and 0.31 g/cm 3 and tapped densities were between 0.34 and 0.41 g/cm 3 . These small differences are not considered statistically significant. The Hausner ratio values as the parameter of powder flow (important for processing powders into capsules for DPI filling) ranged from 1.25 to 1.47, which corresponds to acceptable to bad flow behavior [1,41,49].
In summary, as follows from the discussion of the individual responses and the significance of the regression model coefficients (Table 3), the crucial factor influencing the resulting particle properties is drying temperature. Increasing drying temperature leads to the formation of larger particles in physical diameter (S) and aerodynamic size (MMAD). Higher feed rate (pump speed) led to the formation of larger particles in physical diameter. Increasing the air speed tends to produce particles with a larger MMAD.
In the promising samples, the MMAD was around 7 µm, which is still greater than the optimal range (1-5 µm). However, sample 13 (central point) represents particles with the lowest MMAD and a passable S. The central point appears to have an optimal setting of process parameters for expanding and incorporating other parameters (e.g., concentration and atomization pressure) or for fine-tuning the properties of the spray-dried LPPs to reach the deep lung region.

Conclusions
The main purpose of this study was to determine the spray-drying process parameters affecting mainly the physical size (S) of the prepared LPPs and their aerodynamic properties (MMAD, NMAD). Other parameters, such as particle morphology, bulk, tapped density, and HR, were less important. The particles' size and aerodynamic diameter were significantly affected by the variation of all three considered process parameters (drying temperature, pump speed, and air speed). In general, the effect of spray-drying process parameters on the properties of mannitol LPPs was described and evaluated. The region around the DoE central point (the drying temperature of 120 • C, the pump speed of 10 mL/min, and the air speed of 3.6 m/s) was found optimal for achieving the desired response of the monitored particles' characteristics. The innovative presented data can provide other research groups with a description of how to tune spray-dried particles properties by changing the process setup.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/pharmaceutics15020496/s1, Table S1: Results of experiment for variables measured in 3 replications: raw data with indicated outliers.