{'fraction_non_alphanumeric': 0.035561640176305594, 'fraction_numerical': 0.018999599305462803, 'mean_word_length': 4.727481353987378, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 20, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Optimization of S-Nitrosocaptopril Monohydrate Storage Conditions Based on Response Surface Method

From unstable crystals to relatively stable monohydrate crystals, many researchers have been working on S-nitrosocaptopril for more than two decades. S-nitrosocaptopril monohydrate (Cap-NO·H2O) is a novel crystal form of S-nitrosocaptopril (Cap-NO), and is not only a nitric oxide (NO) donor, but also an angiotensin-converting enzyme inhibitor (ACEI). Yet, a method for long-term storage has never been reported. In order to determine the optimal storage conditions, Plackett–Burman (PB) design was performed to confirm the critical factors. Response surface methodology (RSM) was employed to determine the optimal Cap-NO·H2O storage condition, based on the rough interval determined by the path of steepest ascent experiment. The optimized storage condition was denoted as nitrogen purity of 97%, temperature of −10 °C and 1.20 g deoxidizer. In this case, a final preservation rate of 97.91 ± 0.59% could be obtained. In specific storage conditions, Cap-NO·H2O was found to be stable for at least 6 months in individual PE package, procreating a potentially applicable avenue.


Introduction
S-nitrosocaptopril is a S-nitrosylated captopril developed more than 20 years ago [1], but its technical bottleneck of large-scale synthesis was recently overcome [2]. Moreover, S-nitrosocaptopril monohydrate (Cap-NO·H 2 O) is a novel and stable compound based on S-nitrosocaptopril (Cap-NO), which is not only a nitric oxide (NO) donor but also an angiotensin-converting enzyme inhibitor (ACEI). Cap-NO·H 2 O has great potential for the study and possible treatment of numerous diseases such as pulmonary hyperattention [2][3][4], hypertension regardless of renin status, angina pectoris [5] and congestive heart failure [6], acute respiratory distress syndrome [7] and preventing postoperative circulation tumor cells metastasis [8].
One of the most common ways of degradation of S-nitrosocaptopril Monohydrate (Cap-NO·H 2 O) during storage is oxidation reaction [9]. When Cap-NO·H 2  Although the content of -S-/N≡O+ is only 6-10%, it will attack the -S-N bond of Cap-NO, resulting in a knock-on effect of Cap-NO degradation [2]. At the same time, if there are oxygen molecules in the environment, oxygen free radicals will be generated and interact with Active Pharmaceutical Ingredients (API), impurities and other auxiliary materials to generate more free radicals and accelerate the degradation of Cap-NO·H 2 O. Although the Molecules 2021, 26, 7533 2 of 13 precise mechanism of the interaction of the oxidation process is not very clear, the reaction process is generally an auto-oxidation process [10], including three processes of initiation, circulation and termination.

Initiation:
In· + R-H→R· + In-H Circulation: R· + O 2 →ROO· R· + R-H→R-H· R-OO·+ R-H→ROOH + R· Termination: 2ROO→Product Among them is an unknown free radical initiator, and R-H represents pharmaceutical substances, excipients, or other pollutants. Whether directly reacting with the drug or indirectly acting through other excipients, oxygen molecules participated in the propagation step and the catalytic cycle of the oxidative degradation of the pharmaceutical substance in the drug formulation. The oxidative degradation of API leads to a decrease in drug efficacy and an increase in oxidative degradation products.
The oxidative degradation of API leads to the decrease of efficacy and the increase of oxidative degradation products (impurities). The rate of oxidative degradation depends on the nature of API and the sensitivity of storage conditions (such as temperature, humidity, oxygen concentration, light and time, etc.). It is known that the harmful effects of the oxidation process of Cap-NO·H 2 O include product discoloration (from red to white), changes in dissolution rate, red-brown gas generation and bad odors, etc. [2]. The key point is that the generated oxidative degradation products may lead to the adverse reactions of API, or the effects on the efficacy of the medicine.
To solve this problem, the common solutions include: (1) adding non-mixed antioxidants or chelating agents and (2) placing it in an inert gas atmosphere [11]. On this basis, the API is stored in a constant and low temperature environment. Although the form of adding antioxidants or chelating agents to drugs is more common, there are fewer studies on the control and application of multi-influencing factors under conditions of nitrogen content. This may be because it is difficult to quantitatively evaluate the oxygen content of a packaging system. However, inert gas protection is a direct method to effectively isolate the API and its auxiliary materials from oxygen. By removing one of the two key reactants in the oxidation cycle from the packaging, the degradation of susceptible formulations can be minimized, and the shelf life of oxygen-sensitive commodities will be prolonged.
PBD is an experimental program developed by R.L. Plackett and J.P. Burman in 1946 [12]. In order to improve the quality control process, the design could be used to study the influence of design parameters on the system state, so as to make reasonable decisions. The orthogonal array designed by Plackett and Burman (PB) can produce unbiased estimates of all major influences in the least possible design. However, the initial estimate of the optimum values derived from the fractional factorial design could be far from the true optimum values. Thus, the method of steepest ascent was usually introduced to move rapidly to the general vicinity of the actual optimum [13]. Then, response surface methodology (RSM) should be conducted to obtain an optimal response as a statistical multi-response optimization method [14]. RSM is an effective statistic tool with minimized experimental trials and a 3D description of the interactions of parameters [15].
In the present paper, the Plackett-Burman method was performed to screen significant impacts on API storage out of five factors based on previous studies on Cap-NO·H 2 O stability. Then, the response surface method was used to obtain the optimal storage conditions of the small test level that provided guidance for the storage conditions of the industrial scale-up stage.

Plackett-Burman Experimental Results and Analysis
As indicated in Figure 1, Cap-NO·H 2 O (10 µg/mL) was observed as eluted at the wavelength of 215 nm, whose retention time was 9.65 min.

Plackett-Burman Experimental Results and Analysis
As indicated in Figure 1, Cap-NO·H2O (10 μg/mL) was observed as eluted at the wavelength of 215 nm, whose retention time was 9.65 min. The experimental results based on the Plackett-Burman experimental grouping design were shown in Table 1.
The results of ANOVA analysis are shown in Table S1 (Supplementary Material). The significance of the model (P value) was below 0.0001, which indicated that the experimental model was significant; namely, the model had a high degree of fit and was suitable for predicting the importance of each factor. The order of the influence of each influencing factor on the preservation rate of Cap-NO·H2O was "X4 > X3 > X2 > X5 > X1"; namely, the The experimental results based on the Plackett-Burman experimental grouping design were shown in Table 1.
The results of ANOVA analysis are shown in Table S1 (Supplementary Material). The significance of the model (p value) was below 0.0001, which indicated that the experimental model was significant; namely, the model had a high degree of fit and was suitable for predicting the importance of each factor. The order of the influence of each influencing factor on the preservation rate of Cap-NO·H 2 O was "X 4 > X 3 > X 2 > X 5 > X 1 "; namely, the influence of each factor on the preservation rate of Cap-NO·H 2 O was as follows: Temperature > Nitrogen purity > Deoxidizer > Desiccant > Light. The significance of the temperature was below 0.0001, which indicated that the temperature had a very significant effect on the preservation rate of Cap-NO·H 2 O. The significance of the nitrogen purity of the storage atmosphere was below 0.05, which indicated that the nitrogen purity had a significant impact on the preservation rate of Cap-NO·H 2 O. While the significance of deoxidizer was 0.0596 (larger than 0.05), which was not statistically significant, its value was near the critical value. It may be that the number of repetitions was not enough, so it should be considered as one of the key factors.
According to the residual shown in Figure 2a, the residual value was calculated based on the difference between the predicted value and the actual value. In this figure, the actual residual value was plotted along the horizontal x-axis, while the predicted expected value was plotted along the y-axis. It could be seen that basically all the values fell near the straight line, indicating that the residuals obeyed a normal distribution that could intuitively support the results of the above model analysis. As shown in Figure 2b, the Pareto chart ranked the absolute values of each factor's standardized effects in descending order. Among them, only the temperature was above the significant line, the nitrogen purity was above the significant line, and the other factors and virtual factors were not significant, which could directly support the model analysis results.
influence of each factor on the preservation rate of Cap-NO·H2O was as follows: Temp ature > Nitrogen purity > Deoxidizer > Desiccant > Light. The significance of the temp ature was below 0.0001, which indicated that the temperature had a very significant ef on the preservation rate of Cap-NO·H2O. The significance of the nitrogen purity of storage atmosphere was below 0.05, which indicated that the nitrogen purity had a sig icant impact on the preservation rate of Cap-NO·H2O. While the significance of deoxid was 0.0596 (larger than 0.05), which was not statistically significant, its value was near critical value. It may be that the number of repetitions was not enough, so it should considered as one of the key factors.
According to the residual shown in Figure 2a, the residual value was calculated ba on the difference between the predicted value and the actual value. In this figure, the tual residual value was plotted along the horizontal x-axis, while the predicted expec value was plotted along the y-axis. It could be seen that basically all the values fell n the straight line, indicating that the residuals obeyed a normal distribution that could tuitively support the results of the above model analysis. As shown in Figure 2b, the reto chart ranked the absolute values of each factor's standardized effects in descend order. Among them, only the temperature was above the significant line, the nitrogen rity was above the significant line, and the other factors and virtual factors were not nificant, which could directly support the model analysis results.

Results and Analysis of the Path of Steepest Ascent
According to the results of the PB experiment, the insignificant factors (desiccant light) were excluded. A climbing experiment was performed on the remaining fact temperature, nitrogen purity, and deoxidizer. The experimental results were shown Table 2. According to the maximum condition of the experimental results, the stor temperature of the center point of the response surface experiment was −8 °C, the nitro purity was 95%, and the deoxidizer was 1.13 g.

Results and Analysis of the Path of Steepest Ascent
According to the results of the PB experiment, the insignificant factors (desiccant and light) were excluded. A climbing experiment was performed on the remaining factors: temperature, nitrogen purity, and deoxidizer. The experimental results were shown in Table 2. According to the maximum condition of the experimental results, the storage temperature of the center point of the response surface experiment was −8 • C, the nitrogen purity was 95%, and the deoxidizer was 1.13 g.

Results and Analysis of Response Surface Experiments
According to the response surface experimental design grouping, the results of the preservation rate of Cap-NO·H 2 O were shown in Table 3. Table 3. Results of Box-Behnken response surface design.

Number
Temperature X 1 Nitrogen Purity X 2 Deoxidizer X 3 Preservation Rate %Y Based on the above experimental results, ANOVA was used for statistical analysis, and the results were shown in Table S2. According to the statistical results, a ternary quadratic regression equation was fitted and shown as follows: The R 2 of the model was 0.9843, the adjusted R 2 was 0.9531 and the expected R 2 was 0.8162. The difference between the expected R 2 and adjusted R 2 was less than 0.2, indicating that the model was reasonable. From the analysis results in the table, it could be seen that the significance of the model was 0.0001, which was extremely significant, and the lack-of-fit term was 0.148145, which was not significant, indicating that the model could better fit and predict the experimental results. Under the interaction of multiple factors, the P value of temperature, nitrogen purity, and the deoxidizer were all less than 0.01, indicating that those three factors have a very significant impact on the storage of Cap-NO·H 2 O. The P values of nitrogen purity and deoxidizer were less than 0.05 in the two interaction factors, indicating that the interaction between nitrogen purity and deoxidizer had a significant influence in this storage model. However, the interaction between temperature and the other two factors was not significant in this model. As shown in Figure 3, the contour plot and the 3D surface of the center point, which was set by the climbing experiment, deviated from the actual experimental results. However, the equation had a stable point, when X 1 was 10.0417, X 2 was 96.4796 and X 3 was 1.17574, then the output Y was 96.4796. Namely, when the stable point conditions were storage temperature of −10.0417 • C, nitrogen purity of 96.4796% and deoxidizer of 1.17574 g, then the predicted preservation rate (Y) was 96.4796. Meanwhile, the maximum value of Y was 98.296, when X 1 was 10.042, X 2 was 96.480 and X 3 was 1.176, then the output Y was 98.296. When storage temperature was −10.042 • C, nitrogen purity was 96.480% and deoxidizer was 1.176 g, the predicted preservation rate (Y) was 98.296. From the perspective of the steepness of the contour, the temperature had the greatest influence on the response of the system, followed by the nitrogen purity and the deoxidizer. X1 was 10.042, X2 was 96.480 and X3 was 1.176, then the output Y was 98.296. When storage temperature was −10.042 °C, nitrogen purity was 96.480% and deoxidizer was 1.176 g, the predicted preservation rate (Y) was 98.296. From the perspective of the steepness of the contour, the temperature had the greatest influence on the response of the system, followed by the nitrogen purity and the deoxidizer.  According to the prediction of the model, the actual experimental conditions were adjusted to storage temperature of −10 • C, nitrogen purity of 97%, and deoxidizer of 1.20 g. Under this condition, three sets of repetitions were performed, and the final true value was 97.91 ± 0.59%, as shown in Table 4, which was basically consistent with the predicted value. Those results show that the model had certain guiding significance for the practical preservation of Cap-NO·H 2 O.

Discussion
Through the PB experiment, it was known that the key influencing factors for the preservation of Cap-NO·H 2 O were temperature, nitrogen purity of the storage environment, and the deoxidizer. However, the desiccant and light were not the key factors, which seems to differ from the known properties of the functional group -SNO [16]. The possible reasons for this result were as follows: (1) the experimental result may come from the ultraviolet light, which is the key factor for photosensitivity of Cap-NO·H 2 O [17], and the PA/PE composite sealing packaging bag may block ultraviolet light; and (2) the Cap-NO·H 2 O used in the experiment had low moisture content after drying, so it was difficult to obtain moisture from the environment under the sealed conditions of small space and relatively stable ambient temperature. During the experiment, the use of desiccant had no significant influence on the stable storage of Cap-NO·H 2 O crystals.
Three-dimensional (3D) response surface and two-dimensional (2D) contour plots were provided as graphical representations of the regression equation (Figure 3). As shown in Figure 3a,b, the effect of temperature on storage changed more rapidly than that of deoxidant and nitrogen purity, since the curved surface of temperature was much steeper than the curved surface of the latter. The factor of nitrogen purity only had a relatively large effect at lower levels of purity with steep and sharply dissected slopes (Figure 3a,c), while the change of deoxidant in intervals tends to be gentler with gradual slopes (Figure 3b,c). The elliptical shape of all the contour plots showed that three pairwise interactions were significant, particularly when the factor of temperature was included (Figure 3b,d).
The results of the response surface showed that under the interaction of temperature, nitrogen purity and deoxidizer, all three factors had a significant effect on the stable storage of Cap-NO·H 2 O crystals, but only the interaction of nitrogen purity and deoxidizer were statistically significant (Table S2). The reasons were as follows: (1) temperature was a key factor. When the ambient temperature was high, the impurity coming from the synthesis of Cap-NO·H 2 O or the -SNO group was relatively active, and the form content of -S-/N≡O+ would increase and attack the functional groups of other molecules, causing chain degradation reactions. At this time, the control of the oxygen content and humidity of the environment would no longer suppress the occurrence of this process. Moreover, (2) the interactive relationship between the nitrogen purity and the deoxidizer were remarkable because when the deoxidizer absorbed the oxygen in the environment, the nitrogen purity in the circumstance increased relatively.
However, when the use of deoxidizer increased, the preservation rate of Cap-NO·H 2 O crystals would decrease. This may be because the selected mixed deoxidizer was a mixture of iron (Fe) powder. Although there was packaging to isolate Cap-NO·H 2 O crystals and deoxidizer, when the amount of deoxidizer is too high, some of it may exist as aerosol ions in the packaging environment, and metal ions could degrade the functional group of Cap-NO·H 2 O-SNO.
In summary, in the subsequent pilot trials and expanded production of drugs, due to the full consideration of the packaging factors of Cap-NO·H 2 O and the temperature control in the storage and transportation links, the common materials PA/PE can block ultraviolet radiation. The composite material can play a better role in blocking UV degradation of Cap-NO·H 2 O. The drying process of Cap-NO·H 2 O was an important step in the production process. This step can greatly save the cost of subsequent drug storage. At the same time, we could further consider replacing the deoxidizer or improving the packaging of the deoxidizer, so that the package would not only isolate the impact of deoxidizer on the API, but also deoxidize the environment, or replace the non-metal ion type deoxidizer. On the other hand, the experimental results of the response surface also gave us an inspiration of the application level, namely, in this trial with the deoxidizer, the high purity of nitrogen in the nitrogen-filled atmosphere was not a necessary condition.

The Synthesis of Cap-NO·H 2 O
The synthesis of Cap-NO·H 2 O refers to the experimental procedures of published paper [2].

The Preparation and Use of Nitrogen Environment
A specific ratio of nitrogen was prepared by the improved drainage method [18]. As shown in Figure 4a, we filled the graduated gas cylinder with deionized water and buckled it upside down into a vessel containing deionized water. According to the ratio requirement of the experimental group, the corresponding gas with a certain volume ratio was filled from the gas charging port, and the corresponding volume of water was discharged to obtain the nitrogen-containing gas with a corresponding volume ratio. was filled from the gas charging port, and the corresponding volume of water was discharged to obtain the nitrogen-containing gas with a corresponding volume ratio. The calculation method of charging gas was shown in the following formula: where: x-Response volume of air; x ≥ 0. y-Response volume of nitrogen; y > 0. Q-The percentage of nitrogen in the gas mixture; Q > 78. After a simple transform, the equation can be denoted as follows: After the gas collection was complete, we placed the bottle body upright. We put glass wool on the bottom of the drying tower to prevent the desiccant above the waist from falling, filled the drying tower with CaCl2 desiccant above the waist, then applied petroleum jelly to the upper bottle and rotated it to prevent air leakage, and aligned the air outlet with the tower body. The connecting port was used to prevent the moisture from affecting the API in the package. Deionized water was injected from the water injection port, and the nitrogen-containing gas was discharged and injected into the API package. The specific link method is shown in Figure 4b.

Nitrogen-Filled Packaging Method
Nitrogen filling and sealing adopted a small multifunctional continuous vacuum pumping/inflating packaging machine (LF1080B). The operating parameters were as follows: power was 1.15 kW, power source was 220 V 50 Hz and sealing speed was 0-12 m/min. A total of 2 g Cap-NO·H2O was weighed and placed into the packaging bag. The procedure was as follows: (1) Nitrogen filling: connected the pressure control valve to the nitrogen cylinder (tighten), kept the pressure regulating valve closed, turned on the cylinder switch, and adjusted the pressure to 0.5-1 kPa. Turned on the packaging machine and the power switch, then turned on the heating switch to set the heating temperature to 220 °C. When the temperature reached the set temperature, turned on the fan, air extraction and nitrogen flushing. Adjusted the nitrogen inlet pressure valve to make the pressure appropriate (about 0.1-0.3 kPa).
(2) Mixed gas filling: connected the drying tower and the ultrapure water machine system to the water injection port for continuous water injection.
(3) Charging air: did not connect the air port. Packing: Inserted the suction tube into the side of the packaging bag containing the sample (away from the conveyor belt), then stepped on the foot switch to start pumping. The calculation method of charging gas was shown in the following formula: where: x-Response volume of air; x ≥ 0. y-Response volume of nitrogen; y > 0. Q-The percentage of nitrogen in the gas mixture; Q > 78. After a simple transform, the equation can be denoted as follows: After the gas collection was complete, we placed the bottle body upright. We put glass wool on the bottom of the drying tower to prevent the desiccant above the waist from falling, filled the drying tower with CaCl 2 desiccant above the waist, then applied petroleum jelly to the upper bottle and rotated it to prevent air leakage, and aligned the air outlet with the tower body. The connecting port was used to prevent the moisture from affecting the API in the package. Deionized water was injected from the water injection port, and the nitrogen-containing gas was discharged and injected into the API package. The specific link method is shown in Figure 4b.

Nitrogen-Filled Packaging Method
Nitrogen filling and sealing adopted a small multifunctional continuous vacuum pumping/inflating packaging machine (LF1080B). The operating parameters were as follows: power was 1.15 kW, power source was 220 V 50 Hz and sealing speed was 0-12 m/min. A total of 2 g Cap-NO·H 2 O was weighed and placed into the packaging bag. The procedure was as follows: (1) Nitrogen filling: connected the pressure control valve to the nitrogen cylinder (tighten), kept the pressure regulating valve closed, turned on the cylinder switch, and adjusted the pressure to 0.5-1 kPa. Turned on the packaging machine and the power switch, then turned on the heating switch to set the heating temperature to 220 • C. When the temperature reached the set temperature, turned on the fan, air extraction and nitrogen flushing. Adjusted the nitrogen inlet pressure valve to make the pressure appropriate (about 0.1-0.3 kPa). (2) Mixed gas filling: connected the drying tower and the ultrapure water machine system to the water injection port for continuous water injection. (3) Charging air: did not connect the air port.
Packing: Inserted the suction tube into the side of the packaging bag containing the sample (away from the conveyor belt), then stepped on the foot switch to start pumping.
After the suction was completed, moved the suction tube to the conveyor belt and kept for 1-3 s until it was full of gas, so that it could enter the conveyor belt for sealing.
After packaging, we adjusted the set temperature to room temperature (without turning off the fan), and turned off the air extraction and inflation switches. When the temperature dropped below 80 • C, turned off the heating power and fan, we then turned off the power and the cylinder switches, and removed the pressure control valve.

Chromatographic Conditions
We weighed precisely 10mg of dried Cap-NO·H 2 O, added 1ml of methanol to dissolve it, diluted it by 10,000 times (10 µg/mL), and then filtered it with a 25 µm filter into an automatic sample bottle. The chromatographic separation was performed at an HPLC instrument (2695-2489, Waters, USA) equipped with an ODS column (250 mm × 4.6 mm × 5 µm) at the wavelength of 215 nm. Column temperature was set at 40 • C, and the injection volume was 20 µL. Isocratic elution was run for 20 min, and the mobile phase consisted of 0.1% phosphoric acid aqueous solution and methanol (50/50, v/v) at a flow rate of 1 mL/min. The PB experiment in this paper adopted a 3-level experimental design, therefore n = 3. It adopted an experimental design table with n = 12 times, and the number of dummy variables was 6. It measured the percentage of remaining amount of Cap-NO·H 2 O after 180 days of storage under the control of other variables, such as light, deoxidizer, nitrogen purity, temperature and desiccant. The factor design was shown in Table 5. After screening several significant factors, the key factors were selected in the follow-up experiments. Exposure experiment: one set of packaging bags was made of brown packaging bags and wrapped in light-proof paper, and the other set was made of transparent packaging bags; deoxidizer experiment: 1.5 g deoxidizer was put into one group before the package was sealed, and no deoxidizer was placed in the other group; nitrogen filling experiment: one group was filled with nitrogen in the package, the other group was filled with air; temperature experiment: one group was placed in a refrigerator at −20 • C, and the other group was exposed to room temperature; desiccant experiment: 1.2 g desiccant was put into one group, and no desiccant was placed in the other group. Each experimental group was repeated 3 times.

The Experiment Design of the Path of Steepest Ascent
According to the PB experiment, after removing the non-critical factors, in order to quickly find the best state point interval of the remaining key factors, the path of steepest ascent was adopted. The first-order linear fitting formula of the steepest climbing experiment was shown as follows:ŷ where: y-Optimal response value; b i -Variation coefficient of influence factor; x i -Impact factors. The step length formula was shown as follows: where: ∆x j -The variable step with the largest absolute value of the regression coefficient; b j -Absolute maximum regression coefficient; b i -Calculation coefficient.

The Experiment Design of Response Surface
The response surface method (RSM) is a classical statistical modeling technique dedicated to evaluating the interaction between a single experimental group factor and one or more experimental group factors. Therefore, in order to obtain a reliable model, each factor needs to first be tested. In this paper a central composite experimental design [19] was adopted to optimize the storage conditions of Cap-NO·H 2 O with nitrogen purity, storage temperature, deoxidizer quality, desiccant quality, and light as the influencing factors. Preliminary studies have shown that nitrogen purity, storage temperature and desiccant were important factors affecting the stable storage of Cap-NO·H 2 O. Therefore, these variables were set as X 1 , X 2 , and X 3 , respectively. The variables followed the formula below: where: x i -Independent variable; X i -The real value of the independent variable; X i -The central value of the independent variable; ∆X i -Step length. The independent variable for Box-Behnken factor level design in this paper was shown in Table 6. The maximum storage capacity of Cap-NO·H 2 O was adopted as the response valuê Y i (U). Regression analysis was carried out on the acquired data by using Design Expert 10. The predicted response was calculated with the following second-order polynomial: where: Y i -Predicted response value; x i , x j -Independent variables; β i -Linear coefficient; β ij -Interaction term coefficient; β ii -Square coefficient.

Statistical Methods
Statistical calculations in this paper used the statistical software Design Expert 10 and IBM SPSS Statistics 23.
Supplementary Materials: The following are available online. Table S1: Significance analysis of Plackett-Burman experimental design results; Table S2