Application of Plackett–Burman Design for Spectrochemical Determination of the Last-Resort Antibiotic, Tigecycline, in Pure Form and in Pharmaceuticals: Investigation of Thermodynamics and Kinetics

Tigecycline (TIGC) reacts with 7,7,8,8-tetracyanoquinodimethane (TCNQ) to form a bright green charge transfer complex (CTC). The spectrum of the CTC showed multiple charge transfer bands with a major peak at 843 nm. The Plackett–Burman design (PBD) was used to investigate the process variables with the objective being set to obtaining the maximum absorbance and thus sensitivity. Four variables, three of which were numerical (temperature—Temp; reagent volume—RV; reaction time—RT) and one non-numerical (diluting solvent—DS), were studied. The maximum absorbance was achieved using a factorial blend of Temp: 25 °C, RV: 0.50 mL, RT: 60 min, and acetonitrile (ACN) as a DS. The molecular composition that was investigated using Job’s method showed a 1:1 CTC. The method’s validation was performed following the International Conference of Harmonization (ICH) guidelines. The linearity was achieved over a range of 0.5–10 µg mL−1 with the limits of detection (LOD) and quantification (LOQ) of 166 and 504 ng mL−1, respectively. The method was applicable to TIGC per se and in formulations without interferences from common additives. The application of the Benesi–Hildebrand equation revealed the formation of a stable complex with a standard Gibbs free energy change (∆G°) value of −26.42 to −27.95 kJ/mol. A study of the reaction kinetics revealed that the CTC formation could be best described using a pseudo-first-order reaction.


Introduction
Tetracyclines (TCs) are a large family of antibiotics that are used both within veterinary and therapeutic rehearsals [1]. Tigecycline (TIGC), as shown in Scheme 1, is the newest member of the tetracycline family that belongs to its third generation [2,3]. TIGC is the first antibacterial drug to be classified under glycylcyclines [4]. Until 2019, TIGC was listed as an 'essential medicine' by the World Health Organization (WHO) [5,6]. As a (9-t-butylglycylamido) derivative of the parent minocycline, TIGC possesses an enhanced ability to overcome the two resistance mechanisms that the TCs could encounter: ribosomal protection and the TC-specific efflux pump acquisition [7]. low possibilities for renal and other organ toxicity, the use of TIGC could help reduce the hospital burden on the other broad-spectrum antimicrobials [7]. Nevertheless, the administration of TIGC was related to some incidents of unknown deaths. As of 2010, a black box warning for TIGC was consequently released [10]. As a result, the use of TIGC as a last-resort drug was reserved for conditions in which an alternate treatment is not suitable. Moreover, with the release of TIGC into aquatic bodies that likely contain the Tet(X) gene, the degradation of this antibiotic possibly causes the further emergence of microbial resistance [2,3,11]. Therefore, finding a simple and sensitive approach for sensing even ultra-low concentrations of TIGC is crucial. A survey of the literature shows that the majority of the efforts in this regard were dedicated to studying TIGC in terms of its mechanism of action, pharmacokinetics, and pharmacodynamics, with fewer efforts being dedicated to finding In this itinerary, TIGC is a broad-spectrum antibiotic that kills bacterial cells by inhibiting protein synthesis. TIGC is commonly utilized to cure several diseases such as those caused by skin and intra-abdominal pathogens [7]. Yet, first and foremost, TIGC is used to cure infections instigated by multiple-drug resilient pathogens [8]. Recently, TIGC has been considered as a promising candidate for treating acute myeloid leukemia [9]. With an unknown influence on the emergence of resistance in hospitalized patients and low possibilities for renal and other organ toxicity, the use of TIGC could help reduce the hospital burden on the other broad-spectrum antimicrobials [7]. Nevertheless, the administration of TIGC was related to some incidents of unknown deaths. As of 2010, a black box warning for TIGC was consequently released [10]. As a result, the use of TIGC as a last-resort drug was reserved for conditions in which an alternate treatment is not suitable. Moreover, with the release of TIGC into aquatic bodies that likely contain the Tet(X) gene, the degradation of this antibiotic possibly causes the further emergence of microbial resistance [2,3,11].
Therefore, finding a simple and sensitive approach for sensing even ultra-low concentrations of TIGC is crucial. A survey of the literature shows that the majority of the efforts in this regard were dedicated to studying TIGC in terms of its mechanism of action, pharmacokinetics, and pharmacodynamics, with fewer efforts being dedicated to finding a simple analytical technique for the determination of TIGC. The reported techniques mainly included chromatographic-based approaches [12][13][14][15][16] and to a lesser extent spectrophotometric [17] and spectrofluorimetric [18,19] approaches. By and large, the chromatographic approaches are sophisticated and require well-trained staff. Moreover, and to the best of our information, all the reported techniques for the analysis of TIGC were univariate-based, where a single variable is investigated at a time while the rest are kept constant. Employing this practice for optimizing process variables implies the absence of the overall vision of the process with a greater consumption of chemicals and resources.
Charge transfer complexes (CTC) are usually formed via molecular interactions of electron donors and electron acceptors with the subsequent development of a colored complex that exhibits different behaviors and properties [20][21][22]. 7,7,8,8-tetracyanoquinodimethane (TCNQ) is a strong electron acceptor with four cyano-groups and a π-conjugation system, as shown in Scheme 2 [21].
low possibilities for renal and other organ toxicity, the use of TIGC could help reduce the hospital burden on the other broad-spectrum antimicrobials [7]. Nevertheless, the administration of TIGC was related to some incidents of unknown deaths. As of 2010, a black box warning for TIGC was consequently released [10]. As a result, the use of TIGC as a last-resort drug was reserved for conditions in which an alternate treatment is not suitable. Moreover, with the release of TIGC into aquatic bodies that likely contain the Tet(X) gene, the degradation of this antibiotic possibly causes the further emergence of microbial resistance [2,3,11]. Scheme 1. Chemical structure of TIGC. Therefore, finding a simple and sensitive approach for sensing even ultra-low concentrations of TIGC is crucial. A survey of the literature shows that the majority of the efforts in this regard were dedicated to studying TIGC in terms of its mechanism of action, pharmacokinetics, and pharmacodynamics, with fewer efforts being dedicated to finding a simple analytical technique for the determination of TIGC. The reported techniques mainly included chromatographic-based approaches [12][13][14][15][16] and to a lesser extent spectrophotometric [17] and spectrofluorimetric [18,19] approaches. By and large, the chromatographic approaches are sophisticated and require well-trained staff. Moreover, and to the best of our information, all the reported techniques for the analysis of TIGC were univariate-based, where a single variable is investigated at a time while the rest are kept constant. Employing this practice for optimizing process variables implies the absence of the overall vision of the process with a greater consumption of chemicals and resources.
Charge transfer complexes (CTC) are usually formed via molecular interactions of electron donors and electron acceptors with the subsequent development of a colored complex that exhibits different behaviors and properties [20][21][22]. 7,7,8,8-tetracyanoquinodimethane (TCNQ) is a strong electron acceptor with four cyano-groups and a πconjugation system, as shown in Scheme 2 [21]. In the current study, TCNQ will be used as the π-acceptor to react with TIGC as an electron donor to produce an intensely colored CTC. The approach we are offering herein is unique in terms of utilizing CTC for the sensitive and selective determination of TIGC, as well as being the first multivariate-based approach for the determination of TIGC with all the benefits this platform offers. In addition to saving time, efforts, and lowering the possible number of experimentations, a huge amount of data will be generated, an issue that helps to illustrate the reaction under investigation. Therefore, the conclusions can be Scheme 2. Chemical structure of TCNQ.
In the current study, TCNQ will be used as the π-acceptor to react with TIGC as an electron donor to produce an intensely colored CTC. The approach we are offering herein is unique in terms of utilizing CTC for the sensitive and selective determination of TIGC, as well as being the first multivariate-based approach for the determination of TIGC with all the benefits this platform offers. In addition to saving time, efforts, and lowering the possible number of experimentations, a huge amount of data will be generated, an issue that helps to illustrate the reaction under investigation. Therefore, the conclusions can be considered highly irrevocable. It is noteworthy to mention that the Plackett-Burman design (PBD) is a two-level fraction factorial screening design that is used when only the main variables are concerned [23][24][25][26].
The aim of this work is to determine TIGC in its pure form and in the pharmaceutical dosage forms via a sensitive and selective spectrochemical approach. Moreover, the optimum conditions that maximize the reaction response, absorbance in this case, can be realized using PBD as a strategy. In this itinerary, four variables (temperature, reagent vol- ume 'RV', reaction time 'RT', and the diluting solvent 'DS') will be assessed. One response variable will be measured: absorbance at 843 nm (Y ctc ). Therefore, the optimized approach can be applied in quality control and quality assurance laboratories especially in developing countries where the chance of using sophisticated approaches such as chromatography is not always plausible. Job's method of continuous variation will be applied to determine the best stoichiometry for the TIGC-TCNQ interaction [27]. The limits of detection (LOD) and quantification (LOQ) will be determined and the whole process will be assessed in terms of accuracy, precision, repeatability, ruggedness, and will follow the International Conference of Harmonization (ICH) recommendations [28].

Absorption Spectra
The interaction of TCNQ as a π-acceptor with TIGC as an n-electron donor has produced an intensely colored green complex that shows three major peaks at 843, 744, and 680 nm, as shown in Figure 1. These peaks can be ascribed to the development of TCNQ radical anions as a result of the transfer of n-electrons from TIGC to TCNQ in the organic solvent that was used (acetonitrile, ACN). The absorption band at 843 nm was chosen as the analytical wavelength considering the sensitivity and blank absorbance. sign (PBD) is a two-level fraction factorial screening design that is used when only the main variables are concerned [23][24][25][26].
The aim of this work is to determine TIGC in its pure form and in the pharmaceutical dosage forms via a sensitive and selective spectrochemical approach. Moreover, the optimum conditions that maximize the reaction response, absorbance in this case, can be realized using PBD as a strategy. In this itinerary, four variables (temperature, reagent volume 'RV', reaction time 'RT', and the diluting solvent 'DS') will be assessed. One response variable will be measured: absorbance at 843 nm (Yctc). Therefore, the optimized approach can be applied in quality control and quality assurance laboratories especially in developing countries where the chance of using sophisticated approaches such as chromatography is not always plausible. Job's method of continuous variation will be applied to determine the best stoichiometry for the TIGC-TCNQ interaction [27]. The limits of detection (LOD) and quantification (LOQ) will be determined and the whole process will be assessed in terms of accuracy, precision, repeatability, ruggedness, and will follow the International Conference of Harmonization (ICH) recommendations [28].

Absorption Spectra
The interaction of TCNQ as a π-acceptor with TIGC as an n-electron donor has produced an intensely colored green complex that shows three major peaks at 843, 744, and 680 nm, as shown in Figure 1. These peaks can be ascribed to the development of TCNQ radical anions as a result of the transfer of n-electrons from TIGC to TCNQ in the organic solvent that was used (acetonitrile, ACN). The absorption band at 843 nm was chosen as the analytical wavelength considering the sensitivity and blank absorbance.

Charge Transfer Reaction: Proposed Pathway
The interaction of TIGC (D) and TCNQ (A) is probably based on the reaction of the nitrogen atom of TIGC as a donor with the π-acceptor system of TCNQ to produce a CTC. The resulting product will consequently dissociate into radical anions as determined by

Charge Transfer Reaction: Proposed Pathway
The interaction of TIGC (D) and TCNQ (A) is probably based on the reaction of the nitrogen atom of TIGC as a donor with the π-acceptor system of TCNQ to produce a CTC. The resulting product will consequently dissociate into radical anions as determined by the solvent polarity. In the current investigation, can the polar solvent, will assist in a complete electron transfer from a TIGC to TCNQ moiety with the formation of intensely colored radical anions, as shown in Scheme 3 [29][30][31][32]. the solvent polarity. In the current investigation, can the polar solvent, will assist in a complete electron transfer from a TIGC to TCNQ moiety with the formation of intensely colored radical anions, as shown in Scheme 3 [29][30][31][32]. The structure of TIGC shows five ionizable groups: three basic nitrogens and two acidic hydroxyl groups [2,33,34]. In TIGC, the n-electron donation could have occurred via the aliphatic amine group (carboxamide side chain of Ring A), or the secondary amine group in the side chain (pKa = 8.9), or the tertiary amino group at ring A (pKa = 9.5). However, steric hinderance and the unavailability or delocalization of electrons to the benzene ring should be taken into consideration. The structure of TIGC shows five ionizable groups: three basic nitrogens and two acidic hydroxyl groups [2,33,34]. In TIGC, the n-electron donation could have occurred via the aliphatic amine group (carboxamide side chain of Ring A), or the secondary amine group in the side chain (pK a = 8.9), or the tertiary amino group at ring A (pK a = 9.5). However, steric hinderance and the unavailability or delocalization of electrons to the benzene ring should be taken into consideration.

Assessment of Reaction Conditions
As previously revealed, the PBD was the design used to optimize the reaction conditions. Four variables (three numerical and one non-numerical) and a single response variable were considered. The target was set to maximize the response. The PBD generates a number of runs (N) that is a multiple of four.
The PBD is the design of choice when merely the key factors are concerned. The PBD is one of the ordinarily used designs in assessing the robustness while validating an analytical method. The key reason for its use is that the PBD only focuses on the main variables [24,35,36]. The screened factors are shown in Table 1, while the design matrix is shown in Table 2 along with the measured and the predicted responses.

Response Transformation and Modelling
The assessment of the substantial variables was performed using different approaches. The Pareto chart of standardized effects was one of the approaches that were used accompanied by an analysis of variance (ANOVA) testing at a 95.0 confidence interval (CI). As shown in Figure 2, variables that exceed the reference line are seen as statistically significant. The reaction time (RT, C) was the most influencing variable while the reagent volume (RV, B) was statistically insignificant. Similar conclusions can be derived from the ANOVA testing shown in Table 3 where only the RV has a probability value (p-value) that exceeds the significance level (α = 0.05), implying that the variable is statistically insignificant. Table 3 also shows that the lack-of-fit values were close to 1.000, thus implying the goodness-of-fit of the obtained data.

Response Transformation and Modelling
The assessment of the substantial variables was performed using different approaches. The Pareto chart of standardized effects was one of the approaches that were used accompanied by an analysis of variance (ANOVA) testing at a 95.0 confidence interval (CI). As shown in Figure 2, variables that exceed the reference line are seen as statistically significant. The reaction time (RT, C) was the most influencing variable while the reagent volume (RV, B) was statistically insignificant. Similar conclusions can be derived from the ANOVA testing shown in Table 3 where only the RV has a probability value (pvalue) that exceeds the significance level (α = 0.05), implying that the variable is statistically insignificant. Table 3 also shows that the lack-of-fit values were close to 1.000, thus implying the goodness-of-fit of the obtained data.   It is important to point out that these data were collected after response transformation, applying the Box-Cox response transformation where the value of λ (transformation factor) was 0.50 as per the following equation, Equation (1) [37,38]: The output of this processing is the mathematical model shown in Equation (2) where Ct. Pt. is the central point. The model summary shows the coefficient of determination, R 2 = 0.9665, R 2 (adjusted) = 0.9572, and R 2 (predicted) = 0.9378. The high value of R 2 reflects the fact that the proposed model well-fit the obtained data. The high value of R 2 (predicted) infers that the model has a good prediction capability for new observations. This is well-reflected and can be concluded from Table 2, where the values of the predicted absorbance were very close to the observed ones as inferred from the small values of RE. The difference between R 2 (adjusted) and R 2 (predicted) was less than 10%, implying the absence of model over-fitting. The normal probability plot of residuals was used to verify the postulation that the residuals are normally distributed. Figure 3 shows that Anderson-Darling (AD) indicator value, which is used to determine how favorably data obey a certain distribution, is less than the p-value, signifying the good distribution of the data [39].
Pure Error 7 0.03792 0.005417 Total 23 1.60803 * DF is degrees of freedom; SS is sum of squares, and MS is mean of squares. Adj SS is the adjusted sum of squares, and Adj MS is the adjusted mean of squares.
It is important to point out that these data were collected after response transformation, applying the Box-Cox response transformation where the value of λ (transformation factor) was 0.50 as per the following equation, Equation (1) [37,38]: The output of this processing is the mathematical model shown in Equation (2): where Ct. Pt. is the central point. The model summary shows the coefficient of determination, R 2 = 0.9665, R 2 (adjusted) = 0.9572, and R 2 (predicted) = 0.9378. The high value of R 2 reflects the fact that the proposed model well-fit the obtained data. The high value of R 2 (predicted) infers that the model has a good prediction capability for new observations. This is well-reflected and can be concluded from Table 2, where the values of the predicted absorbance were very close to the observed ones as inferred from the small values of RE. The difference between R 2 (adjusted) and R 2 (predicted) was less than 10%, implying the absence of model over-fitting. The normal probability plot of residuals was used to verify the postulation that the residuals are normally distributed. Figure 3 shows that Anderson-Darling (AD) indicator value, which is used to determine how favorably data obey a certain distribution, is less than the p-value, signifying the good distribution of the data [39].

Impact of the Diluting Solvent
As a part of the investigation, the impact of the polarity of the diluting solvent on the measured response was investigated. The individual value plot was used to check the presence of outliers and the distribution coverage, as shown in Figure 4. A visual inspection of the represented data shows that the use of ACN has resulted in a slightly higher absorbance compared to methanol (MeOH). As a part of the investigation, the impact of the polarity of the diluting solvent on the measured response was investigated. The individual value plot was used to check the presence of outliers and the distribution coverage, as shown in Figure 4. A visual inspection of the represented data shows that the use of ACN has resulted in a slightly higher absorbance compared to methanol (MeOH). A further investigation was conducted using a 2-sample t-test. The purpose of conducting this test was to measure the effect of the used solvent on the absorbance of the CTC. The obtained results implied that the difference between MeOH and ACN on the CTC that was formed was statistically insignificant, where the p-value was greater than the significance level (α = 0.05). In addition, the mean of the measured response in the case A further investigation was conducted using a 2-sample t-test. The purpose of conducting this test was to measure the effect of the used solvent on the absorbance of the CTC. The obtained results implied that the difference between MeOH and ACN on the CTC that was formed was statistically insignificant, where the p-value was greater than the significance level (α = 0.05). In addition, the mean of the measured response in the case of ACN was not markedly higher in comparison to that obtained using MeOH, as shown in Table 4.

Optimization Phase: Contour Plots
The contour plots, Figure 5, display the relationship between two factors and the response surface, which appears as contours. Figure 5 shows a sample contour plot for the relationship between the reaction time (RT) and the temperature (Temp), and the absorbance at 843 nm. The legend on the right-hand side shows that dark green zones are the zones in which absorbance exceeds 1.5. This absorbance value was attained using a factorial combination of RT: 53-60 min and by heating samples the at 47-50 • C.

Optimization Phase: Individual Desirability Function
An optimization plot (desirability function plot; Figure is not shown) was utilized to find the optimum conditions that maximize the absorbance. The obtained individual desirability (d) was 1.0000, indicating that the optimal conditions are favorable to obtain the maximum absorbance of the CTC. The optimum conditions (which will be used for construction of calibration curve) to obtain an absorbance value of 1.000 are 25 °C, 0.50 mL of TCNQ, 60 min as RT, and ACN as DS.

Linear Range and Sensitivity
The calibration curve of the CTC formation reaction shows a linear relationship between the concentration (0.5-10 µg mL −1 ) and absorbance recorded at λmax = 843 nm. The linear relationship is presented by Equation (3): The value of R 2 was relatively high, reflecting the linearity of obtained data. The an-

Optimization Phase: Individual Desirability Function
An optimization plot (desirability function plot; Figure is not shown) was utilized to find the optimum conditions that maximize the absorbance. The obtained individual desirability (d) was 1.0000, indicating that the optimal conditions are favorable to obtain the maximum absorbance of the CTC. The optimum conditions (which will be used for construction of calibration curve) to obtain an absorbance value of 1.000 are 25 • C, 0.50 mL of TCNQ, 60 min as RT, and ACN as DS.

Linear Range and Sensitivity
The calibration curve of the CTC formation reaction shows a linear relationship between the concentration (0.5-10 µg mL −1 ) and absorbance recorded at λ max = 843 nm. The linear relationship is presented by Equation (3): The value of R 2 was relatively high, reflecting the linearity of obtained data. The analytical data are summarized in Table 5. The sensitivity was assessed using the limit of detection (LOD) and limit of quantification (LOQ), determined using Equations (4) and (5): where SD represents the standard deviation of the blank and a is the slope of the straight line of the calibration curve of TIGC. The LOD and LOQ values are listed in Table 5. The values of LOD and LOQ indicate that the proposed method is extremely sensitive, with a nanogram level detection limit. Table 5. Analytical parameters for the determination of TIGC-TCNQ CTC at λ max = 843 nm.

Accuracy and Precision
The accuracy and precision were measured at 95.0% CI for three concentrations of TIGC as bulk powder and in the pharmaceutical dosage form. The intra-day precision was evaluated by measuring each sample three times in the same day, while the inter-day precision was assessed over three different days. The results showed that the % relative standard deviation (%RSD) was <3%, indicating good precision. The accuracy was expressed in terms of %RE. Table 6 summarizes the obtained data and shows the accuracy and precision of the current method. The accuracy was additionally evaluated by applying the standard addition technique, as shown in Table 6. Calibration and standard addition methods were applied to the pharmaceutical formulation and the %mean recovery values for both methods were high (100.50 ± 2.91% for the calibration method and 100.80 ± 2.46% for the standard addition method), as shown in Table 6, indicating that the matrix (commonly added adjuvants and diluents) has a minor effect.

Comparison to the Reference Method
The proposed method was compared to a reference spectrophotometric method in terms of the student t-test and an F-statistic, as shown in Table 7. The results show that the obtained tand Fvalues were less than the tabulated values signifying that the proposed technique is comparable to the reference approach [17].

Determination of the Reaction Stoichiometry
The molecular composition of the resultant CTC was determined by conducting Job's method of continuous variation. A plot of the relationship between the absorbance of the resulting complex versus the mole fraction (V r /V r + V d ) was drawn, where V r is the volume of TCNQ and V d is the volume of TIGC-the Figure is not shown. The maximum amount of the complex formed was observed at a mole fraction of 0.50; therefore, the best stoichiometric ratio for the reaction of TIGC and TCNQ was 1:1.

Investigation of the Reaction Thermodynamics
The formation constant for the charge transfer reaction of TIGC and TCNQ was calculated using the Benesi-Hildebrand equation [40][41][42][43] as follows (Equation (6)): where  Table 8. and A AD , ε AD , and are the absorbance, the molar absorptivity, and the formation constant of the CTC, respectively. The Benesi-Hildebrand equation can be applied for the 1:1 CTC by keeping [Ao] lower than [Do]. In other words, [Do] must be 5-10 times higher than [Ao]. The thermodynamic parameters were obtained by studying the reaction at four different temperatures (25,40,55, and 70 °C). Plots of [ ] ⁄ versus 1 [ ] ⁄ were sketched, as shown in Figure 6. Straight-line equations were acquired for the CTC at different temperatures. The slope of the straight-line is 1/ and the intercept is 1 ⁄ . The obtained figures for and ε AD , determined at 25, 40, 55, and 70 °C, are displayed in Table 8.    The obtained data show that K AD c and ε AD at room temperature (298 K) have reached high values of 7.90 × 10 4 L mol −1 and 4.96 × 10 4 L mol −1 cm −1 , thus corroborating the stability of the formed CTC. These high values could be ascribed to the high donating ability of TIGC, prominent electron affinity of TCNQ, and the high electric permittivity of ACN. Moreover, these high values demonstrate the dissociation of the outer sphere CTC (TIGC-TCNQ) resulting in the formation of radical ions of the two species with high electrostatic attraction, as shown in Scheme 3 [44][45][46][47].
According to the obtained data at the different temperatures, the complex formation was thermodynamically favored as reflected by the immense negative value of the standard Gibbs free energy change (∆G • ), ranging between −26.42 and −27.95 kJ/mol. The value of ∆G • decreases (less negative) as the temperature increases, implying that a complex formation might not be favored at a higher temperature. The formation constant values, K AD c , confirm these findings. The values of K AD c , calculated at 298, 313, 328, and 343 K, were employed to compute the standard enthalpy of formation (∆H • ) employing Van't Hoff plots as depicted in Figure 7 (summary is shown in Table 8) and utilizing Equation (7) [48,49]: (7) [48,49]: In this equation, R denotes the gas constant (8.314 J/mol. K), and T stands for the temperature (Kelvin). Plotting log against the reciprocal of the absolute temperature, as shown in Figure 7, has resulted in a linear relationship with a slope of −∆H°⁄(2.303 R). The determined ∆H° is a high negative value, signifying that the formation of the TIGC-TCNQ CTC is exothermic in nature. The standard entropy change (∆S°) was obtained using Equation (8): The value of ∆S° was low (−31.93 ± 0.37 J/K. mol), indicating the simple composition of the CTC without the interference of the solvent. In this equation, R denotes the gas constant (8.314 J/mol. K), and T stands for the temperature (Kelvin). Plotting log K AD c against the reciprocal of the absolute temperature, as shown in Figure 7, has resulted in a linear relationship with a slope of −∆H • ⁄(2.303 R).
The determined ∆H • is a high negative value, signifying that the formation of the TIGC-TCNQ CTC is exothermic in nature. The standard entropy change (∆S • ) was obtained using Equation (8): The value of ∆S • was low (−31.93 ± 0.37 J/K. mol), indicating the simple composition of the CTC without the interference of the solvent. Table 8 also shows that the value of ε AD (molecular extinction coefficient) increases as the temperature increases, as shown in Figure 8. Conversely, the value of ε AD is inversely proportional to the formation constant, K AD c .
Pharmaceuticals 2022, 15, x FOR PEER REVIEW 13 of 19 Table 8 also shows that the value of ε AD (molecular extinction coefficient) increases as the temperature increases, as shown in Figure 8. Conversely, the value of ε AD is inversely proportional to the formation constant, . Figure 8. The relationship between the temperature (K) and the molecular extinction coefficient (ε).

Evaluation of the Reaction Kinetics
The kinetic studies were conducted by measuring the absorbance of the resultant CTC after reaction times of 5, 15, 30, 45, 60, 75, and 90 min. The obtained data show that the absorbance of the CTC increases as the reaction time increases. This behavior was utilized to assess the reaction kinetics. Employing the formerly outlined optimum settings,

Evaluation of the Reaction Kinetics
The kinetic studies were conducted by measuring the absorbance of the resultant CTC after reaction times of 5,15,30,45,60,75, and 90 min. The obtained data show that the absorbance of the CTC increases as the reaction time increases. This behavior was utilized to assess the reaction kinetics. Employing the formerly outlined optimum settings, the initial rates were obtained from the slopes of absorbance-time graphs. The reaction order was then established corresponding to the donor reactant (TIGC) in the presence of a fixed concentration of acceptor (TCNQ). In this study, several kinetic models were used, including the initial rate, rate constant, fixed concentration, and fixed time approaches [50][51][52][53][54]. Moreover, these techniques were evaluated, and the choice of the best technique depended merely on the relevance, LOD, LOQ, and linear range of the acquired data.

Initial Rate Approach
Following the previously outlined optimal conditions, the reaction's initial rates were established by utilizing the slopes of the absorbance-time plots. The reaction order can be determined by plotting reaction rates versus the initial absorbance. Consequently, the rate of reaction can be represented by Equation (9): In this equation, K symbolizes the rate constant, and m + n signifies the overall reaction order with respect to TCNQ (m) and TIGC (n). The initial rate of the CTC formation reaction follows a pseudo-first-order reaction rate and can be expressed as follows, in Equation (10): Log (rate) = log K = log∆A/∆t = log k + n log [C] (10) In this equation, A represents the absorbance and t is the reaction time (s). A regression of the log (rate) versus log [TIGC] is provided Equation (11): Log (rate) = log K = log ∆A/∆t = 1.1742 + 0.9977 log C, R 2 = 0.9526 (11) Therefore, K = 14.93 s −1 and the slope = n = 0.9977 1, substantiating that the CTC formation follows the pseudo-first-order. Based on these data, the reaction rate was [TIGC]-dependent and would be best portrayed using Equation (12): In Equation (12), K represents the pseudo-first-order rate constant, and n denotes the reaction order related to [TIGC].

Fixed-Time Approach
In this approach, calibration curves were drawn in the range of 1-8 µg mL −1 of TIGC at reaction times of 5,15,30,45,60,75, and 90 min. The regression equations for each reaction time were calculated and the results are shown in Table 9. From the results, the best linearity was obtained at 60 min, confirming the findings of the optimization phase of the PBD.  Graphs of absorbance versus time (min) were plotted and a straight line was depicted to intersect with the largest possible number of curves. Another graph of 1/t (s) versus concentration (M) was thenceforth obtained, Equation (14):

Assessment of the Kinetics Procedures
A comparison between the four procedures proposed for the determination of the reaction kinetics is shown in Table 10. It is evident that the initial rate and fixed time methods showed the widest linear range with the lowest LOD and LOQ. Yet, the rate constant method shows less linearity compared to the ones shown for the fixed-time and fixed-concentration methods. S b = SD of slope, ±tS b = confidence limit for slope, S a = SD of intercept, ±tS a = confidence limit for intercept, S y/x = SD of the regression.

Instrumentation and Software
The software employed to constitute the PBD matrix was Minitab ® 19 (Minitab ® Inc., State College, PA, USA). Deionized water utilized in the current investigation was obtained from a Millipore-Q water system (Burlington, MA, USA). Absorbance was determined by means of a UV-Vis spectrophotometer (Agilent diode-array, Agilent, Santa Clara, CA, USA) equipped with 10 mm quartz cuvettes (matched). Samples were heated to the required temperatures when needed using a thermostatically controlled water bath.

Reagents and Standards
Stock solutions of TIGC (0.05%) and TCNQ (0.1%) were freshly prepared daily by weighing specific amounts of TIGC and TCNQ powders and diluting them with acetonitrile to the mark in 100 mL volumetric flasks. Working solutions were obtained via serial dilutions of the stock solutions using the same solvent.

Materials
Tigecycline (TIGC, Batch# AT108181901) was the product of Biosynth ® Carbosynth Ltd. (Compton, Berkshire, UK). The electron acceptor 7,7,8,8-tetracyanoquinodimethane (TCNQ, purity 98%) was purchased from Sigma-Aldrich (St. Louis, MO, USA). Acetonitrile procured from BDH Chemicals (BDH Laboratory Supplies, Poole, UK) and methanol from Merck KGaA (Darmstadt, Germany) were used as received with no further purification. Tygacil ® vials were purchased from local pharmacy stores in Cairo, Egypt. The studied variables and their proposed levels are listed in Table 1. Twenty-four experimental runs took place (with 4 center points). Aliquots of 0.05% stock solution (500 µg mL −1 of TIGC) were used in all experimental runs. A volume of 0.1% TCNQ, as shown in the scenario revealed in Table 2, was added to TIGC solution and volume was completed to the mark using the suitable diluting solvent. Type of diluting solvent, temperature at which reaction took place, and the reagent volume are all shown for each run in the design matrix table, as shown in Table 2. Reagent blanks were prepared and measured similarly. The absorbance of the resulting green colored solutions was measured at λ max = 843 nm. To create the calibration curve, different concentrations of the stock solution were prepared from the TIGC stock solutions by serial dilution to obtain a final concentration in the range of 0.5-10 µg mL −1 , and the calibration curve was constructed using the same procedure and implementing the optimal conditions.

Procedure for the Formulation
Tygacil ® vials (product of Patheon Italia S.P.A., Italy), labelled to contain 50 mg TIGC per vial, were the formulation of choice. The lyophilized content of the vial was further crushed and an amount of 30.8 mg of the powdered material (equivalent to 10.0 mg TIGC) was accurately weighed, dissolved in ACN, filtered, and transferred into a 50 mL volumetric flask. The volume was completed using ACN and formulation stock solution of 200 µg mL −1 TIGC was then ready for further analysis.

Standard Addition Method
Ten samples were made by inserting a fixed amount of 30 µL (200 µg mL −1 ) of the formulated drug solution to ten volumetric flasks (labelled as S0-S10). Different volumes (10-180 µL) of TIGC (0.05%) were then added followed by applying the optimum conditions and measuring the absorbance at λ max = 843 nm.

Procedure for Job's Method
Job's method [27] was applied to ascertain the molar ratio for the interaction of TIGC with TCNQ. Equimolar solutions (0.854 mM) of TIGC and TCNQ were prepared. A total of ten samples (10 mL) were prepared where the total volume was kept at 2.0 mL of both drug and reagent for each solution, with the aid of acetonitrile as a solvent. Absorbance of the prepared set was recorded at λ max = 843 nm contrasted to a reagent blank.

Investigation of the Reaction Kinetics
Seven sets were prepared with a [TIGC] of 1-8 µg mL −1 . Absorbance of every group was recorded at 5, 15, 30, 45, 60, 75, and 90 min as described under the general procedure.

Conclusions
A simple, sensitive, and selective methodology was developed to generate a CTC involving tigecycline as the electron donor and TCNQ as a πacceptor. The Plackett-Burman Design (PBD) was implemented to attain the optimum process variables and the maximum response-absorbance of the CTC in this case. A Pareto chart and ANOVA testing were used to determine the statistical significance of the studied variables. The results showed that the RV was not statistically significant. The desirability function plot showed that the optimum conditions that could be used to maximize the response were Temp: 25 • C, RV: 0.50 mL, RT: 60 min, and ACN as DS. Job's method of continuous variation showed that the complex has a molecular composition of 1:1. ICH recommendations were used to assess the developed technique. The obtained data show excellent accuracy and precision with no significant differences compared to the reference method. Interferences from the common excipients and additives were not observed. Therefore, the developed approach could be utilized for the routine analysis of tigecycline in its pure form and in formulations. The calibration curve obtained under the optimum conditions was rectilinear in the range 0.5-10 µg mL −1 . The stability of the CTC complex was determined using thermodynamic studies; thus, the Benesi-Hildebrand equation and Van't Hoff plots showed the formation of a stable complex with a formation constant of 1.05 × 10 4 − 7.90 × 10 4 L mol −1 and a molar absorptivity of 4.96 × 10 4 − 8.89 × 10 4 L mol −1 cm −1 . The Reaction kinetics were studied and revealed that the interaction between TIGC and TCNQ follows a pseudo-first-order reaction, and the best linearity was obtained at RT: 60 min.