Integration of Fuzzy Matter-Element Method and 3D-QSAR Model for Generation of Environmentally Friendly Quinolone Derivatives

The environmental pollution of quinolone antibiotics (QAs) has caused rising public concern due to their widespread usage. In this study, Gaussian 09 software was used to obtain the infrared spectral intensity (IRI) and ultraviolet spectral intensity (UVI) of 24 QAs based on the Density Functional Theory (DFT). Rather than using two single-factor inputs, a fuzzy matter-element method was selected to calculate the combined effects of infrared and ultraviolet spectra (CI). The Comparative Molecular Field Analysis (CoMFA) was then used to construct a three-dimensional quantitative structure–activity relationship (3D-QSAR) with QAs’ molecular structure as the independent variable and CI as the dependent variable. Using marbofloxacin and levofloxacin as target molecules, the molecular design of 87 QA derivatives was carried out. The developed models were further used to determine the stability, functionality (genetic toxicity), and the environmental effects (bioaccumulation, biodegradability) of these designed QA derivatives. Results indicated that all QA derivatives are stable in the environment with their IRI, UVI, and CI enhanced. Meanwhile, the genetic toxicity of the 87 QA derivatives increased by varying degrees (0.24%–29.01%), among which the bioaccumulation and biodegradability of 43 QA derivatives were within the acceptable range. Through integration of fuzzy matter-element method and 3D-QSAR, this study advanced the QAs research with the enhanced CI and helped to generate the proposed environmentally friendly quinolone derivatives so as to aid the management of this class of antibiotics.


Introduction
Quinolone antibiotics (QAs) are a new class of artificially synthesized anti-infective drugs. They are widely used in the treatment of various infectious diseases due to their broad bactericidal spectrum, strong antibacterial properties, low toxicity, few side effects, and low price [1]. They are a class of antibiotics widely applied in veterinary medicine. After usage, QAs will remain in the animal as their original molecular structures or pass through feces and urine in the form of metabolites. These QAs and their metabolites will excrete into the external environment [2]. According to Xander et al. [3], 47 research papers have reported the distribution of QAs in the environment and found that their residues are concentrated worldwide including in Pearl River, Huangpu River, Xiaoqing River, Hai River, Baiyangdian Lake, and Dongting Lake in China; densely populated watersheds such as Lake Ontario in North America, Taft River and in Europe, and Mekong River in Asia. Zhang et al. [4] used the III fugacity model to study the distribution characteristics of QAs in 58 watersheds in China. The results showed that after metabolism in humans and animals, the total amount of QAs excreted from feces and urine is 5.4 × 10 7 kg. Only 0.02 × 10 7 kg are removed in wastewater treatment, while the remaining 5.38 × 10 7 kg still directly enter the receiving environment. Chen et al. [5] used gas chromatography and mass spectrometry to detect QA concentrations in 19 key groundwater monitoring wells in Beijing. The average concentrations of the QAs ciprofloxacin and norfloxacin were found to be 4.9 and 0.2 ng, respectively. In addition, studies have shown that the levels of norfloxacin, ciprofloxacin, and enoxacin in groundwater in Spain, were 64. 13, 38.93, 44.47 ng, respectively [6]. Norfloxacin, ciprofloxacin, and enoxacin in the groundwater of a Swedish pharmaceutical company reached 31, 14,000, and 1900 ng, respectively [7]. These studies have shown that the high residual concentrations of QAs are a threat to the environment. QAs accumulate in various aquatic organisms, however, the concentration of different QAs are quite different. QAs, including ofloxacin (2.72 µg/kg), enrofloxacin (762.34 µg/kg), and ciprofloxacin (3.08 µg/kg) were found in bivalve samples [8]. A mean value of ofloxacin (5.58 µg/kg), enrofloxacin (3.08 µg/kg), ciprofloxacin (4.17 µg/kg), and norfloxacin (23.8 µg/kg) were found in the muscles of fish, turtles, and birds from Baiyangdian Lake, China [9]. Li et al. [10] found that among 22 antibiotics, QAs had the highest concentration in 190 mollusks samples with a mean concentration of 86.76 µg/kg dry weight. The investigation of antibiotics in mollusks from Bohai Sea, China, indicated that the mean concentrations of QAs in the mollusks were in the order of "NOR > OFL > CIP > FLE > SAR > LOM > ENR > DIF". The bioaccumulation capacity of QAs is highly relevant to their molecular structures [10]. Therefore, it is of critical importance that the relationship between QAs molecular structures and their characteristics can be identified, so that environmentally friendly QA derivatives can be generated.
3D quantitative structure-activity relationship (3D-QSAR) models have been utilized to design and develop potent drugs by correlating 3D-structural features of the chemicals with properties of interest [11]. Among the existing 3D-QSAR methods, the Comparative Molecular Field Analysis (CoMFA) and the Comparative Molecular Similarity Indices Analysis (CoMSIA) are extensively used in the current practice of rational drug design [12]. 3D-QSAR has been applied to design environmentally friendly molecules. Wang et al. [13] designed a pentachlorophenol molecule with lower bioaccumulation by 3D-QSAR. Tong et al. [14] studied the design of HIV protease inhibitors by using a certain 3D-QSAR method (CoMFA or CoMSIA). Gu et al. [15] predicted the octanol-water partition coefficient for polychlorinated naphthalenes through 3D-QSAR models. However, there are still very limited studies about the use of 3D-QSAR to predict the environmentally friendly properties of QAs. Zhao et al. [16] combined QSAR/QSPR with molecular docking to examine the biodegradability of C20-carbonyl and the C21-carboxyl groups of fluoroquinolones. Zhao et al. [17] designed new fluoroquinolones using SYBYL-X 2.0. After changing the molecular groups, the photodegradability of new fluoroquinolones increased from 15.04% to 40.92% [17].
To construct a 3D-QSAR model for QA analysis, effective methods for detection of QAs in environmental samples are needed so as to obtain infrared and ultraviolet spectral characteristics of QAs to build the database. Infrared and ultraviolet spectral based analytical methods have been well developed for QA detection. Claine [18], Efthimiadou et al. [19], Skyrianou et al. [20], and Zampakou et al. [21] conducted infrared spectroscopy analysis to detect QAs. Bailac et al. [22] detected quinolones in chicken tissues by liquid chromatography based on ultraviolet spectral analysis. Zampakou et al. [21] studied the interaction of quinolone antimicrobial complexes with calf-thymus deoxyribonucleic acid investigated by ultraviolet spectroscopy. However, 3D-QSAR model development needs a large amount of infrared and ultraviolet spectral characteristics of QAs as inputs. Obtaining such information is expensive and time-consuming. Rather than using experimental analysis to generate infrared and ultraviolet spectral data, theoretical calculation methods [23,24] can assist QAs quantification to speed up the database generation. Among them, the Density Functional Theory (DFT) based calculation has been widely adopted. Qiu [25] et al. used DFT to calculate the infrared and ultraviolet spectra of phthalate esters and their derivatives, and used 3D-QSAR to analyze the spectral changes of phthalates before and after substitution. Yang [26] et al. used DFT to calculate the infrared spectra of 36 polybrominated biphenyl molecules and recorded the highest infrared vibration intensity of each molecule.
DFT can help to obtain either infrared or ultraviolet spectra of QAs. These infrared or ultraviolet spectra information can then feed a 3D-QSAR model as a single-factor input. However, instead of using two single-factor inputs (i.e., infrared spectra intensities and ultraviolet spectra intensities), an infrared-ultraviolet integrated-factor is desired to increase the efficiency and accuracy of a 3D-QSAR model. It is thus essential to introduce an effective method for integrating infrared and ultraviolet spectral characteristics accurately. The fuzzy matter-element analysis theory can serve this purpose. It can solve the multiple parameter evaluation problem of incompatibility through establishing the corresponding matter-element [27]. Yang et al. [28] used the fuzzy matter-element method to assess the water resources carrying capacity of six regions in Huai River, China. Liu and Zou [29] used the fuzzy matter-element evaluation method to assess water quality, obtaining very similar results to the official reported water quality. Compared with the traditional method for superscale calculation, the fuzzy matter-element analysis method has lower workload and would overcome the adverse effects from abnormal values. The integration of fuzzy matter-element method and 3D-QSAR models has a great potential for generating more accurate results, though the topic has not been tackled in previous studies.
In this study, 24 QAs were selected to build the 3D-QSAR models, and among them, marbofloxacin and levofloxacin were selected as the target molecules to be modified. DFT was used obtain infrared and ultraviolet spectra intensities of 24 QAs. The fuzzy matter-element method was adopted to deal with the magnitude of single spectrum (infrared spectral intensity (IRI) and ultraviolet spectral intensity (UVI)) intensities, and the combined effect (CI) values that were obtained. In addition, a 3D-QSAR model based on CI was constructed, and the substituted sites of QAs (improving the CI values) were precisely obtained based on the contour map. The associated CI, IRI, and UVI were investigated and compared. Eventually, the stability, functionality (genetic toxicity), and environmental friendliness (bioaccumulation and biodegradability) of all QA derivatives were evaluated in order to generate more environmentally friendly QA derivatives. This study considers both functionality and environmental friendliness of QA derivatives which leads to a tool for a generation of new antibiotics in the future.

Data Sources
The theoretical calculations of the infrared and ultraviolet spectra of the 24 QAs studied in this paper are based on Gaussian-09 software (Gaussian Inc., Wallingford, CT, USA). Gaussian is the most widely used software package for computational and quantum chemistry [30]. Gaussian can be used to process larger molecular weight structures such as anthraquinone, polyhydroxylated anthroquinones, and polyfluorinated dibenzo-p-dioxins for research and calculation [24,31]. DFT is a classical computational method of electronic structure theory which has been frequently used [31]. In most DFT calculations today, the combination of Gaussian orbits is commonly used to represent atomic orbits [32]. The infrared spectrum and ultraviolet spectrum are optimized for the structure of QAs in a gaseous environment at the B3LYP/6-311G (d) level. The infrared spectrum is the optimal structure and its value is calculated and corrected by a correction factor of 0.9614. The ultraviolet spectrum is based on the optimized ground-state and excited-state geometry, calculated using the DFT method [32].

Construction of Composite Fuzzy Matter-Element Matrix about QAs' CI
The elementary elements of an ordered triad R = (M, C, q) are described by the three matter-elements: quinolones (M), detection type (C), and spectral intensity (q) [33]. There are two types of detection methods of quinolone M: infrared spectrum (C 1 ) and ultraviolet spectrum (C 2 ), corresponding spectral intensities q 1 , q 2 , . . . , q 24 , and R is called m-dimensional fuzzy matter-element. When all the quantities in the composite matter-element (R m ) have ambiguity, then R m is called a composite fuzzy matter-element and can be recorded as: where R mn is m dimensional matter-element for 24 QAs, C j is the detection type of j-th characteristic (j = 1, 2), M i the i-th molecule (i = 1, 2, . . . , 24), q ij is the i-th spectral intensity corresponding to the j-th detection type of a molecule.

Determination of Subordination Membership Degree of QAs' CI
The subordination membership degree indicates the fuzzy value corresponding to the single effect parameter of the infrared and ultraviolet spectrums belonging to the standard value of each single effect parameter. The degree of membership is represented by η ij . The fuzzy matter-element matrix R' mn of the preferential membership can be written as: where R' mn is the priority membership fuzzy matter-element; η ij is the priority membership of the j-th detection type of the i-th molecule belonging to the standard sample, i = 1, 2, . . . , 24; j = 1, 2. The priority membership can be recorded as: where Maxq ij represents the maximum value of the spectral intensity in each detection type corresponding to each molecule; Minq ij represents the minimum value of the spectral intensity in each detection type corresponding to each molecule. It is easier to detect the effect of infrared and ultraviolet spectra, when the value of the effects is larger. Therefore, the largest index will be used for further calculation.

Determination of Standard Fuzzy Matter-Element and Variance Compound Fuzzy Matter-Element Matrix for the CI of QAs' IRI and UVI
The m-dimensional standard fuzzy matter-element R 0m is the maximum or minimum value of the spectral strength of the superiority of the spectral intensity in each detection type in the fuzzy matter-of-priority membership. The form is as follows: where M 0 represents the standard sample; η 0j represents the maximum value of the membership degree of the j-th detection type.
R ∆ of variance compound fuzzy matter element can be written as: where ∆ i j = (η 0j -η ij ) 2 ; because the largest subordination membership degree of this paper is the most optimum, η 0j equals to 1.

Determination of Weighting for Evaluation Index of the CI of QAs' IRI and UVI
The weight determination methods include a subjective weighting approach and objective weighting approach [34]. The subjective weighting approach in this paper uses the empirical method and the objective weighting approach uses the entropy method. The entropy method is calculated according to the following steps: (1) Construct a judgment matrix R for the infrared spectrum and ultraviolet spectrum of 24 QAs, and normalize R to obtain a normalized judgment matrix A. Determine the entropy of the j-th detection type according to the definition of entropy: where S j is the entropy of the j-th detection type; a ij is the fuzzy value of each item in the normalized judgment matrix.
(2) Calculate the entropy weighting of each detection type, it can be recorded as follows: where w j is the entropy of the j-th detection type, j = 1, 2.
(3) Determine the comprehensive weights of the comprehensive effect parameters of QAs' infrared and ultraviolet spectra: The objective weighting w j determined by the entropy method is combined with the subjective weighting θ j determined by the empirical method. The comprehensive weighting w' j of each detection type is finally determined as: where ξ is the preference coefficient of the subjective weighting ξ∈[0,1]. A larger ξ value indicates that the research focuses more on subjective weighting; conversely, a lower ξ indicates that decision makers should pay more attention to objective weighting.

Euclidean Closeness Calculation for the CI of QAs' IRI and UVI
The Euclidean closeness was used to comprehensively analyze the infrared and ultraviolet spectra of 24 QAs. The calculation formula of Euclidean closeness e i can be determined as follows:

Construction of 3D-QSAR Model to Generate Environmentally Friendly QAs
Based on the database generated from Section 2.2, a 3D-QSAR model was constructed in this paper. The 3D-QSAR analysis was performed using SYBYL-X2.0 software (Tripos Company, St. Louis, USA) [35]. The comprehensive effects of infrared and ultraviolet spectra of 24 QAs were calculated based on the fuzzy matter-element method. When calculating the parameters of the CoMFA field, partial least squares (PLS) analysis was used to establish the relationship between the structure of the target compound and the biological activity [36]. When using PLS analysis, the Leave-One-Out method was used to cross-validate the training set compounds and calculate the cross-validation coefficient q 2 and the number of optimal principal components n [37]. Non-cross-validation regression (No Validation) was then used to perform regression analysis and calculate the non-cross-validation coefficient r 2 , standard deviation SEE, and test value F, in order to complete the establishment of the CoMFA model [38]. In addition, single-effect models of infrared and ultraviolet spectra also follow the above operations to construct CoMFA models of single-effects of infrared and ultraviolet spectra of QAs.

Calculation of QAs' IRI and UVI based on DFT
There are many types of QAs, many of which have different behavioral characteristics and hazards in the environment. Therefore, determining the types of QAs that exist in the environment and extracting spectral information assists prediction and evaluation of the magnitude of the hazards caused by such substances to the environment. In this study, the UVI of 24 QAs were obtained using the DFT method based on the optimized ground-state and excited-state geometry. The IRI of 24 QAs were calculated using the Gaussian 09 software with a correction factor of 0.9614. The intensities of both the infrared and ultraviolet spectrum were optimized in a gaseous environment at the B3LYP/6-311G (d) level. Table 1 lists the calculation results of infrared and ultraviolet spectral intensities of 24 QAs. Two characteristics of QAs, infrared and ultraviolet spectrums, were applied and the associated fuzzy membership degree, standard fuzzy matter-element and variance composite fuzzy matter-element matrix of the two spectral intensities were calculated using the fuzzy matter-element method according to Equations (1)- (9). The calculation results are summarized in Table 2, which were used to further calculate the CI of QAs.

Comprehensive Weightings for CI of QAs' IRI and UVI
The data of the infrared and ultraviolet spectra of QAs are derived from the results of DFT calculations in quantum chemical calculations. Both detection methods have shown similar physical significance when detecting quinolones [39,40]. Therefore, the subjective weighting for both infrared and ultraviolet spectra was set as 0.5. In addition, the obtained data was normalized before calculating the objective weightings. The specific results are shown in Table 3. Table 3. The results of objective weightings and comprehensive weightings of QAs' infrared and ultraviolet spectra.

Entropy
Entropy Weights Comprehensive Weights From the results, it can be seen that the objective weights calculated by the entropy method are 0.5002 and 0.4998, and the comprehensive weights are 0.5001 and 0.4999, which are close to 0.5.

CI of QAs' IRI and UVI
In this study, the Euclidean closeness was calculated by using Equation (9) with the processed single-effect parameters of 24 QAs representing the comprehensive effects of infrared and ultraviolet spectra of these QAs. As presented in Table 4, the CI values obtained based on the fuzzy matter-element method varied from 0.1613 to 0.8958. The dynamic range of the CI values is 5.6, greater than that stated in previous studies [41][42][43][44][45]. It indicated that the CI data set of QAs' IRI and UVI could be employed to build the associated 3D-QSAR model based on the comprehensive effect.

Construction of the 3D-QSAR Model
The CI of 16 QAs were randomly selected as data sources. Temafloxacin with the highest CI (e 24 = 0.7296) was chosen as the target molecule, 13 QAs as the training set, and the remaining 4 QAs as the testing set (Temafloxacin exists in both the training set and the testing set). Based on these comprehensive effects, a 3D-QSAR model was constructed. SYBYL-X2.0 (Tripos Company, USA) was used in this study to select the lowest energy conformation of the molecule as the dominant stable conformation, and G-H electrical charges (Gasteiger-Huckel) were loaded. The maximum number of optimizations was 10,000 using the Powell energy gradient method, and the energy convergence was limited to 0.005 kJ/mol [46,47]. The optimized molecules were stored in the database for alignment. All molecules were aligned with the pharmacophore characteristic elements in the marked area shown in Figure 1 as the basic skeleton. The structures of selected target molecules, marbofloxacin and levofloxacin are also shown in Figure 1. The CoMFA module will be used to establish the infrared and ultraviolet spectra suitable for QAs by using the basic skeleton for comprehensive effects prediction.

Construction of the 3D-QSAR Model
The CI of 16 QAs were randomly selected as data sources. Temafloxacin with the highest CI (e24 = 0.7296) was chosen as the target molecule, 13 QAs as the training set, and the remaining 4 QAs as the testing set (Temafloxacin exists in both the training set and the testing set). Based on these comprehensive effects, a 3D-QSAR model was constructed. SYBYL-X2.0 (Tripos Company, USA) was used in this study to select the lowest energy conformation of the molecule as the dominant stable conformation, and G-H electrical charges (Gasteiger-Huckel) were loaded. The maximum number of optimizations was 10,000 using the Powell energy gradient method, and the energy convergence was limited to 0.005 kJ/mol [46,47]. The optimized molecules were stored in the database for alignment. All molecules were aligned with the pharmacophore characteristic elements in the marked area shown in Figure 1 as the basic skeleton. The structures of selected target molecules, marbofloxacin and levofloxacin are also shown in Figure 1. The CoMFA module will be used to establish the infrared and ultraviolet spectra suitable for QAs by using the basic skeleton for comprehensive effects prediction.

Performance Evaluation of the 3D-QSAR model
The evaluation results of 3D-QSAR model show that in the CoMFA model, the cross-validation coefficient q 2 is 0.67 (q 2 > 0.5), and the best principal component n is 3, indicating that the model has a good prediction ability. The cross-validation coefficient R 2 is 0.984 (>0.9), the standard error of estimate (SEE) is 0.023, and the F-test value is 179.271, indicating that the model has a good ability to fit and predict [48].
Golbraikh and Ropsha [49] confirmed that the strict QSAR model verification procedures should include internal and external verification. The use of internal verification parameters such as q 2 cannot assess the quality of the model. Therefore, external verification methods need to be applied. External verification methods are one of the most valuable verification methods, and were applied to evaluate the predictive ability of the obtained model. The overall predictive ability of the CoMFA model was externally verified by predicting the activity of the compounds in the independent testing set [49].

Performance Evaluation of the 3D-QSAR Model
The evaluation results of 3D-QSAR model show that in the CoMFA model, the cross-validation coefficient q 2 is 0.67 (q 2 > 0.5), and the best principal component n is 3, indicating that the model has a good prediction ability. The cross-validation coefficient R 2 is 0.984 (>0.9), the standard error of estimate (SEE) is 0.023, and the F-test value is 179.271, indicating that the model has a good ability to fit and predict [48].
Golbraikh and Ropsha [49] confirmed that the strict QSAR model verification procedures should include internal and external verification. The use of internal verification parameters such as q 2 cannot assess the quality of the model. Therefore, external verification methods need to be applied. External verification methods are one of the most valuable verification methods, and were applied to evaluate the predictive ability of the obtained model. The overall predictive ability of the CoMFA model was externally verified by predicting the activity of the compounds in the independent testing set [49]. The predictive ability of the model is represented by the predicted r 2 pred , and its calculation formula is as follows: In the formula, PRESS refers to the sum of the squared deviations between the calculated and predicted values in the test set, and SD refers to the sum of the squared deviations between the calculated values of the compounds in the test set and the average values of the calculated values of the compounds in the training set.
The constructed CoMFA model was used to predict the activity of the test set molecules, and external verification was performed based on the prediction results ( Table 5). The results show that the interaction test coefficient r 2 pred of the external prediction set is 0.9695 (>0.6) [50]. Tropsha method was also used to evaluate the external prediction ability of this model. The calculation results indicated that r was 0.9932, r 2 was 0.9864 (>0.6), k was 1.0269 (0.85 < k < 1.15), k' was 0.9712 (0.85 < k' < 1.15), r 0 2 was 0.893, and r 0 2 was 0.932. Furthermore, (r 2 −r 0 2 )/r 2 = −0.0060 (<0.1) and (r 2 −r 0 2 )/r 2 = 0.0115 (<0.1) The above parameters all met the external verification requirements [51,52], therefore, the constructed 3D-QSAR model has shown a satisfactory external prediction capability.

Determination of Substitution Characteristics Based on the Contour Maps
In the CoMFA model, the contributions of the steric and electrostatic fields are 68.10% and 31.90%, respectively, suggesting that the steric effect and electrical distribution of the groups will affect CI of QAs. In the steric field, the green area indicates that the introduction of bigger molecular groups in this area can improve the CI of QAs, while the yellow area indicates that the introduction of bigger molecular groups in this area can reduce the comprehensive effect of QAs. In the electrostatic field, the blue area indicates that the addition of positive groups is beneficial to improve the CI of QAs, and the red area indicates that the addition of negative groups is beneficial to the comprehensive effect of QAs.
This study uses marbofloxacin and levofloxacin as examples to modify the substituents at positions 1 (CH 3 ) and 2 (CH 3 ), respectively. From the contour maps of marbofloxacin, it can be seen that the area near the 1-position substituent is mainly blue, indicating that the introduction of a positive group at the 1-position substituent is beneficial to improve the QAs' CI; the green and red regions near the substituent at position 2 indicate that the introduction of bigger molecular groups and negative groups at the substituent in this position is conducive to improving the comprehensive effect value of the infrared and ultraviolet spectra of QAs. From the contour maps of levofloxacin, it can be seen that the vicinity of the substituent at position 1 is mainly blue, indicating that the introduction of a positively charged group at the substituent in this position is beneficial to improve the CI of QAs. The regions near the 2-position substituent are mainly green and red, indicating that the introduction of bigger molecular groups and negative groups at the 2-position substituent is conducive to improving the CI of QAs (Figure 2).  Table S1.

Molecular Spectral Characteristics
In this paper, the 3D-QSAR model was used to evaluate the CI of 87 QA derivatives. In addition, a single-effect model of infrared spectrum and ultraviolet spectrum was established in order to further verify the feasibility and rationality of the fuzzy matter-element method. Therefore, the single-effect model was used to evaluate the infrared spectrum intensity and ultraviolet spectrum single-effect of 87 QA derivatives (Table 6).
By analyzing the data in Table 6, it can be seen that when the comprehensive effects of 87 QA derivatives were enhanced, the single effects of the IRI and UVI were also enhanced. The average value of the comprehensive effect enhancement range was 186.46%, the average value of the singleeffect enhancement amplitude of the infrared spectrum was 69.45%, and the average value of the single-effect enhancement amplitude of the ultraviolet spectrum was 1398.39%. Figure S1 was further generated to reflect the relationship between each QA structure and the associated activities (i.e., IRI, UVI, and CI) with marbofloxacin and levofloxacin used as the target molecules. Results indicated that the molecular design of QA derivatives showed positive effect on all activities with obvious enhancement of IRI, UVI, and CI values.  Table S1. In this paper, the 3D-QSAR model was used to evaluate the CI of 87 QA derivatives. In addition, a single-effect model of infrared spectrum and ultraviolet spectrum was established in order to further verify the feasibility and rationality of the fuzzy matter-element method. Therefore, the single-effect model was used to evaluate the infrared spectrum intensity and ultraviolet spectrum single-effect of 87 QA derivatives (Table 6).
By analyzing the data in Table 6, it can be seen that when the comprehensive effects of 87 QA derivatives were enhanced, the single effects of the IRI and UVI were also enhanced. The average value of the comprehensive effect enhancement range was 186.46%, the average value of the single-effect enhancement amplitude of the infrared spectrum was 69.45%, and the average value of the single-effect enhancement amplitude of the ultraviolet spectrum was 1398.39%. Figure S1 was further generated to reflect the relationship between each QA structure and the associated activities (i.e., IRI, UVI, and CI) with marbofloxacin and levofloxacin used as the target molecules. Results indicated that the molecular design of QA derivatives showed positive effect on all activities with obvious enhancement of IRI, UVI, and CI values. The contour map is used as the basis for obtaining accurate modification information in the 3D-QSAR model. Its spatial distribution and the size of each colored region are closely related to the activity of the designed derivative molecule. Therefore, this paper used marbofloxacin's position 2 as an example to compare the combined effect model of the infrared and ultraviolet spectra with the contour map of the single-effect model. A qualitative mechanism analysis of the spectral enhancement of QA derivatives was conducted in order to reveal the effect of the modified information presented in the contour map on the degree of spectral enhancement of the derivative molecule (Figure 3).

The Enhancement of CI, IRI, and UVI of QA Derivatives' The Quantitative Mechanism Analysis of the Spectral Enhancement of QA Derivatives Based on the Contour Map
The contour map is used as the basis for obtaining accurate modification information in the 3D-QSAR model. Its spatial distribution and the size of each colored region are closely related to the activity of the designed derivative molecule. Therefore, this paper used marbofloxacin's position 2 as an example to compare the combined effect model of the infrared and ultraviolet spectra with the contour map of the single-effect model. A qualitative mechanism analysis of the spectral enhancement of QA derivatives was conducted in order to reveal the effect of the modified information presented in the contour map on the degree of spectral enhancement of the derivative molecule ( Figure 3). According to the contour map of the steric field in the comprehensive effect model, the substituent at position 2 is surrounded by a green region, however, the other positions are also surrounded by green regions. Therefore, it can be considered that the steric field is not significant to the contribution of the 2-position substituent. In the infrared-effect model, there is no colored area near the 2-position substituent, thus, the steric field can be considered to have no contribution to the 2-position substituent in the model. In the ultraviolet-effect model, the green region is near the 2position substituent, and it is the only green region distribution in the contour map. Therefore, the steric field is considered to have the most significant contribution to the 2-position substitution in the model. According to this, it can be qualitatively considered that the order of the contribution rate of the steric field to the modified position is UVI > CI > IRI.
According to the contour map of the steric field, the red area is distributed near the 2-position substituent; however, because the other positions are also close to red areas, the electrostatic field can be considered not significant to the contribution of the 2-position substituent. In the infrared-effect model, there is no colored area near the 2-position substituent, thus, the electrostatic field can be considered to have no contribution to the 2-position substituent in the model; in the ultraviolet-effect model, the red region is near the 2-position substituent, and it is the only red region distribution in the contour map. Therefore, the electrostatic field is considered to have the most significant contribution to the 2-position substitution in the ultraviolet-effect model. According to this, it can be qualitatively considered that the order of the contribution rate of the electrostatic field to the modified position is the UV single effect model > combined effect model > infrared single effect model.
In summary, according to the modified information provided by the contour map, it can be seen that the order of the contribution rate of the steric and electrostatic fields in the three groups of models According to the contour map of the steric field in the comprehensive effect model, the substituent at position 2 is surrounded by a green region, however, the other positions are also surrounded by green regions. Therefore, it can be considered that the steric field is not significant to the contribution of the 2-position substituent. In the infrared-effect model, there is no colored area near the 2-position substituent, thus, the steric field can be considered to have no contribution to the 2-position substituent in the model. In the ultraviolet-effect model, the green region is near the 2-position substituent, and it is the only green region distribution in the contour map. Therefore, the steric field is considered to have the most significant contribution to the 2-position substitution in the model. According to this, it can be qualitatively considered that the order of the contribution rate of the steric field to the modified position is UVI > CI > IRI.
According to the contour map of the steric field, the red area is distributed near the 2-position substituent; however, because the other positions are also close to red areas, the electrostatic field can be considered not significant to the contribution of the 2-position substituent. In the infrared-effect model, there is no colored area near the 2-position substituent, thus, the electrostatic field can be considered to have no contribution to the 2-position substituent in the model; in the ultraviolet-effect model, the red region is near the 2-position substituent, and it is the only red region distribution in the contour map. Therefore, the electrostatic field is considered to have the most significant contribution to the 2-position substitution in the ultraviolet-effect model. According to this, it can be qualitatively considered that the order of the contribution rate of the electrostatic field to the modified position is the UV single effect model > combined effect model > infrared single effect model.
In summary, according to the modified information provided by the contour map, it can be seen that the order of the contribution rate of the steric and electrostatic fields in the three groups of models are the UV single factor model > combined effect model > infrared single effect model. The results are consistent with spectral enhancement (the average value of the combined effect enhancement amplitude is 186.46%, the average value of the infrared spectrum single effect enhancement amplitude is 69.45%, and the ultraviolet spectrum single effect enhancement amplitude is 1398.39%). This shows that the spectral enhancement amplitude of QA derivatives has a certain internal relationship with the contribution of steric and electrostatic fields to its modified positions. The larger the contribution rate, the larger the corresponding spectral enhancement amplitude.
The Quantitative Mechanism Analysis of the Spectral Enhancement of QA Derivatives Based on the Modified Positions and the Properties of Substituted Groups In addition to the distribution of colored regions in the contour map, the spectral enhancement amplitude is also affected by the properties of modified positions and the properties of substituted groups. Therefore, the number of modified positions, modified substitutes, and the properties of substituted groups were studied in this paper in order to reveal its intrinsic relationship with the spectral enhancement amplitude (Table 7).
According to the analysis of the data in Table 7, it was found that among the single substitution modification features, seven QA derivatives modified at the 1-position, such as Derivative 1, have an average enhancement range of the CI of 72.06%. The effect enhancement range is 60.02%, and the UV spectrum single effect enhancement range is 339.31%. Derivative 7 and the other 10 QA derivatives modified at the 2-position substituent have an average enhancement range of 113.08%. The single-effect enhancement of the infrared spectrum is 50.80%, and the single-effect enhancement of the ultraviolet spectrum is 421.09%. In the double-substitution modification feature, 70 QA derivatives, such as Derivative 14, are simultaneously modified at the 1-and 2-positions. The average enhancement of the CI is 208.39%, the IRI is 73.05%, and UVI is 1643.92%. From the results, it can be observed that under each substitution feature, the magnitude of the spectral enhancement amplitude in the three groups of models is still the UV single effect model > composite effect model > infrared single effect model, and the results are similar to the size distribution of the color regions in the contour map consistency.
In the comprehensive effect model, the spectral enhancement of QA derivatives obtained by modification at the 1-position, 2-position, and double substitution at 1-and 2-positions were 72.06%, 113.08%, and 208.39%, respectively. In the infrared single-effect model, the spectral enhancement of QA derivatives obtained by modification at the 1-position, 2-position, and double substitution at 1-and 2-positions were 60.02%, 50.80%, and 73.05%, respectively. In the effect model, the spectral enhancement of the QA derivatives obtained by modification at the 1-position, 2-position, and double substitution at 1-and 2-positions were 339.31%, 421.09%, and 1643.92%, respectively. By comparison, it can be found that, with the exception of the infrared single-effect model, the magnitudes of the spectral enhancement amplitudes in the remaining two models are both double-substitution features > 2-position substituent single-substitution features > 1-position substituent single-substitution features. In addition, the nature of the modified group selected for the 1-position substituent increased by 2.92%, the 2-position substituent was 61.93%, and the 1-and 2-position substituents were modified at the same time to 79.43%. The spectral enhancement amplitude has universal consistency, which indicates that the spectral enhancement amplitude of QA derivatives has a certain correlation with the properties of modification sites and groups. As the number of modification sites increases, the nature of the modified molecular groups becomes more prominent, and the spectrum will have a greater magnitude of enhancement.

Molecular Stability Evaluation of QA Derivatives
The 87 QA derivatives designed based on the CoMFA model are used to characterize the stability of 87 QA derivatives in order to further derivatize QAs.
The positive frequency value of the molecule can directly reflect whether the molecule can stably exist in the environment. When the positive frequency value is greater than zero, it indicates that the molecule can exist stably in the environment; otherwise, it cannot [24,53]. In this paper, DFT was used to calculate the positive frequency values of 87 QAs derivatives to verify their stability (Table 8). From the above results, the positive frequency values of the 87 QA derivatives are all greater than zero, indicating that all 87 QA derivatives designed in this study can stably exist in the environment.

Functional and Environmental Friendliness Evaluation of QA derivatives
Zhao et al. [16,17] selected -lgLOEC (LOEC is the lowest observed effective concentration) to construct a hologram quantitative structure-activity relationship (HQSAR) model of the QAs genetic toxicity. QSAR models were designed independently to predict the genotoxicity, bioaccumulation, and biodegradability of the 87 QA derivatives (Table 9).  Genotoxicity refers to how QAs can selectively inhibit two enzymes that play a role in DNA synthesis in bacteria. Topoisomerase II and IV are two enzymes which interfere with the replication, transcription, repair, and recombination of bacterial DNA. These interferences make it impossible for bacteria to pass down genetic information, which lead to the increase in genotoxicity of QAs, making them conducive to improving the medicinal effect of this class of drugs [54]. According to the prediction results of genotoxicity, compared with marbofloxacin and levofloxacin, the genetic toxicity of 87 QA derivatives increased to varying degrees, and the increase range was 0.24%-29.01%.
The n-octanol-water partition coefficient (K ow ) is one of the most important property parameters for studying the environmental behavior of organic matter. This parameter can simulate the distribution behavior of organic matter in lipid and water, and characterize the bioaccumulative ability of organic matter in environmental media. Organic substances with logK ow greater than 5.0 are banned worldwide [55]. According to the bioaccumulation prediction results, compared to marbofloxacin and levofloxacin, the bioaccumulation of 49 QA-derived molecules has been reduced to varying degrees, with a decrease of 0.25% to 60.47%, and the remaining 38 QA-derived molecules increased. The bioaccumulation of the molecule has increased to varying degrees, but its logK ow value is much lower than 5.0, which is below the banned threshold. Therefore, the bioaccumulation of all 87 QA derivatives are considered to be within the acceptable range.
The biodegradability is expressed by the LibDock Score (LDS) of the oxidoreductase enzyme of Phanerochaete chrysosporium and the QAs in the aerobic process of urban sewage treatment plants. The value of the LDS can represent the biodegradability of QAs [45]. Based on the prediction of biodegradability, it can be seen that compared to marbofloxacin and levofloxacin, the biodegradability of 34 QA derivatives has increased to varying degrees, with an increase of 0.12%-2.52%; for the remaining 53 QA derivatives, the biodegradability of biomolecules has decreased, with a decrease range of 0.12%-16.87%. Among them, the decline of nine QA derivatives is less than 1%, and their biodegradability can be considered unchanged. A total of 43 QA derivatives with increased or unchanged biodegradability can be selected.
In summary, among the 87 QA derivatives designed in this study, the properties of function (genetic toxicity) and environmentally friendliness (bioaccumulation, biodegradability) of 43 QAs are acceptable.

Conclusions
Based on DFT and fuzzy matter-element method, the CI of 24 quinolone antibiotics (QAs) were calculated in this study. A 3D-QSAR model was then constructed by using CI as the dependent variable. A total of 87 QA derivatives with enhanced CI and single effects were designed. All of the designed QA derivatives can stably exist in the environment, and the functionality and environmental friendliness of 43 QA derivatives are within the acceptable range. Derivative 3 may be considered as an alternative of marbofloxacin with increased genotoxicity (10.83%) and biodegradability (0.47%), and decreased bioaccumulation (24.83%). Derivative 74 with the improved genotoxicity (17.49%) and biodegradability (0.29%), as well as the decreased bioaccumulation (6.58%), shows a better performance than levofloxacin. In addition, molecules combined with the modification information provided by the contour map were analyzed to identify why the CI of QA derivatives was enhanced. The research output provides theoretical support for obtaining novel antibiotic drug molecules that are easy to detect and less harmful to the environment and human health.