Multispectral and Molecular Docking Studies Reveal Potential Effectiveness of Antidepressant Fluoxetine by Forming π-Acceptor Complexes

Poor mood, lack of pleasure, reduced focus, remorse, unpleasant thoughts, and sleep difficulties are all symptoms of depression. The only approved treatment for children and adolescents with major depressive disorder (MDD) is fluoxetine hydrochloride (FXN), a serotonin selective reuptake inhibitor antidepressant. MDD is the most common cause of disability worldwide. In the present research, picric acid (PA); dinitrobenzene; p-nitro benzoic acid; 2,6-dichloroquinone-4-chloroimide; 2,6-dibromoquinone-4-chloroimide; and 7,7′,8,8′-tetracyanoquinodimethane were used to make 1:1 FXN charge-transfer compounds in solid and liquid forms. The isolated complexes were then characterized by elemental analysis, conductivity, infrared, Raman, and 1H-NMR spectra, thermogravimetric analysis, scanning electron microscopy, and X-ray powder diffraction. Additionally, a molecular docking investigation was conducted on the donor moiety using FXN alone and the resulting charge transfer complex [(FXN)(PA)] as an acceptor to examine the interactions against two protein receptors (serotonin or dopamine). Interestingly, the [(FXN)(PA)] complex binds to both serotonin and dopamine more effectively than the FXN drug alone. Furthermore, [(FXN)(PA)]–serotonin had a greater binding energy than [FXN]–serotonin. Theoretical data were also generated by density functional theory simulations, which aided the molecular geometry investigation and could be beneficial to researchers in the future.


Introduction
Depression is a common mental illness that affects almost 300 million people of all ages around the world [1]. Depression's effects can last a long time or recur, and they can have a substantial influence on a person's ability to function and live a happy life. Major depressive disorder (MDD) is the primary cause of disability worldwide. At least one distinct depressive episode lasting at least two weeks, as well as evident changes in mood, interests, and pleasure, cognitive changes, and vegetative symptoms, characterize this severe condition. Girls are more likely than boys to be affected by MDD, with a global incidence of 3.0% and 1.8%, respectively [2].
MDD is the second leading cause of chronic illness among all medical diseases according to the years lived with disability [3]. Additionally, MDD is associated with a higher risk of contracting diseases, such as diabetes, heart disease, and stroke, which ultimately enhances the severity of the disease [4]. Notably, death by suicide can be induced by MDD. It is estimated that up to 50% of around 800,000 suicides committed annually around the world take place during a depressed episode [5], and MDD patients are approximately 20-fold more likely than the general population to commit suicide [6].
Fluoxetine hydrochloride (FXN) was first described as a selective serotonin (5-HT)uptake inhibitor in 1974 [7] and was the initial representative of a new group of antidepressant medications identified as selective serotonin reuptake inhibitors. In fact, FXN was the first of this class to be released in the United States. Notably, inhibitors of this class are still among the most often prescribed antidepressants today [8]. Moreover, FXN is the only approved treatment for pediatric and adolescent patients with MDD [9].
Charge transfer (CT), which involves an electron donor interacting with several acceptors, is now a significant factor in the study of drug-receptor binding mechanisms [10] and several biological fields [11]. As a result, CT interactions of some acceptors have been successfully used in pharmacokinetic analysis [12], with extensive studies conducted on CT-receptor compounds for broad applications [13].
Chemists, biologists, and pharmacists are all interested in the chemistry of CT interactions and complex building, particularly when it comes to pharmaceuticals. Their interest in CT complexes is related to their high value and numerous uses in a variety of industries and technological, chemical, biological, and pharmaceutical fields. For instance, CT complexes have been employed in biological, electrical, optical, and magnetic experiments, as well as in investigations of pharmaceutical receptor-binding mechanisms [14][15][16]. The development of quick, dependable, and straightforward techniques for the qualitative and quantitative detection of drugs in bulk and/or pharmaceutical dose forms was also based on the production of CT complexes [17,18].
In the current study, the CT complexes of FXN with six different types of π-acceptors were synthesized and thoroughly characterized at different levels in both solid and liquid forms. The interactions between two protein receptors (serotonin or dopamine) and ligands were also examined by molecular docking using Autodock Vina software. To understand how receptor-ligand interactions work, several molecular dynamic simulations were run. Multiple degrees of comparisons were made between CT complexes of FXN and the FXN drug alone against the serotonin receptor.

Multispectroscopic Studies
Recently, there has been a lot of focus on the development of stable CT complexes that come from the interaction of electron acceptors with chemical or biological molecules. The important physical and chemical characteristics of these complexes are the source of this interest. Compared with the methods of drug determination described previously in the literature, the CT complexation is an important approach that is less complicated, more affordable, and more effective [29]. Understanding the interactions between medications and receptors as well as the mechanisms behind drug action may be possible through the study of drug CT complexes [30].
In the present research, a molar ratio of 1:1 was used for the elemental analyses of all synthesized FXN-CT complexes, which are easily soluble in dimethylformamide and dimethylsulfoxide. The molar conductance values of the CT complexes of FXN revealed a slightly electrolytic nature (37-56 Ω −1 cm −1 mol −1 ) upon the creation of positive and negative datives anions owing to hydrogen bonding.

Multispectroscopic Studies
Recently, there has been a lot of focus on the development of stable CT complexes that come from the interaction of electron acceptors with chemical or biological molecules. The important physical and chemical characteristics of these complexes are the source of this interest. Compared with the methods of drug determination described previously in the literature, the CT complexation is an important approach that is less complicated, more affordable, and more effective [29]. Understanding the interactions between medications and receptors as well as the mechanisms behind drug action may be possible through the study of drug CT complexes [30].
In the present research, a molar ratio of 1:1 was used for the elemental analyses of all synthesized FXN-CT complexes, which are easily soluble in dimethylformamide and dimethylsulfoxide. The molar conductance values of the CT complexes of FXN revealed a slightly electrolytic nature (37-56 Ω −1 cm −1 mol −1 ) upon the creation of positive and negative datives anions owing to hydrogen bonding.
The infrared spectrum of FXN had different stretching vibration bands at 3500 cm −1 (N-H stretching), 3200-3000 cm −1 (C−H stretching), 1600-1400 cm −1 (C=C stretching), 1300 cm −1 (C-O stretching), and 1200-1000 cm −1 (C-F stretching) [31]. After CT complexation, the υ(N-H) band was absent, and the δ(NH) band shifted to lower wavenumbers (1611-1629 cm −1 ). These alterations could result from a single pair of -NH electrons interacting with the six acceptors. The n-π* CT transitions caused the C=C stretching frequency to decrease from 1590 cm −1 to 1569-1537 cm −1 . In addition, the stretching vibration of υ(C-N) shifted to 1349-1324 cm −1 in all six acceptors with the FXN drug complexes according to the interaction between -NH and the acceptors (Figure 1). Few bands were observed from 2373-2452 cm −1 due to intermolecular hydrogen bonding between the basic nitrogen of FXN as an n-donor (D) and the π-acceptor (A), resulting in CT complexes of the n-π type (D-A) ( Figure 2) [32]. Few bands were observed from 2373-2452 cm −1 due to intermolecular hydrogen bonding between the basic nitrogen of FXN as an n-donor (D) and the π-acceptor (A), resulting in CT complexes of the n-π type (D-A) ( Figure 2) [32]. 1 HNMR spectra of the FXN with the six π-acceptors complexes are shown in Figure  3. The -NH peak of the secondary amine of FXN shifted to a lower field (9.11-9.58 ppm) due to the involvement of the -NH proton on the receptor fragment (FXN). Aromatic protons also shifted downfield in the spectra of the FXN complexes. Such a result strengthens the transfer of electrons from the lone electron pair on the nitrogen atom of the secondary amine group -NH to the acceptors and aromatic rings of the FXN donor via intermolecular hydrogen bonding. The form of the different (FXN)(π-acceptors) is dependent on the current acceptor according to electron micrographs due to the varying chemical compositions. EDX spectra confirm the existence of carbon, oxygen, and halogens (chlorine and fluorine) in the FXN 1 HNMR spectra of the FXN with the six π-acceptors complexes are shown in Figure 3. The -NH peak of the secondary amine of FXN shifted to a lower field (9.11-9.58 ppm) due to the involvement of the -NH proton on the receptor fragment (FXN). Aromatic protons also shifted downfield in the spectra of the FXN complexes. Such a result strengthens the transfer of electrons from the lone electron pair on the nitrogen atom of the secondary amine group -NH to the acceptors and aromatic rings of the FXN donor via intermolecular hydrogen bonding. 1 HNMR spectra of the FXN with the six π-acceptors complexes are shown in Figure  3. The -NH peak of the secondary amine of FXN shifted to a lower field (9.11-9.58 ppm) due to the involvement of the -NH proton on the receptor fragment (FXN). Aromatic protons also shifted downfield in the spectra of the FXN complexes. Such a result strengthens the transfer of electrons from the lone electron pair on the nitrogen atom of the secondary amine group -NH to the acceptors and aromatic rings of the FXN donor via intermolecular hydrogen bonding. The form of the different (FXN)(π-acceptors) is dependent on the current acceptor according to electron micrographs due to the varying chemical compositions. EDX spectra confirm the existence of carbon, oxygen, and halogens (chlorine and fluorine) in the FXN  (Figure 4). The homogeneity and consistency of particle morphologies in synthesized FXN CT complexes imply that the morphological phases of the [(FXN)(πacceptor)] complexes have a consistent matrix, with particle sizes ranging from 50 to 500 µm ( Figure 4).
Molecules 2022, 27, x FOR PEER REVIEW 5 of 20 complexes ( Figure 4). The homogeneity and consistency of particle morphologies in synthesized FXN CT complexes imply that the morphological phases of the [(FXN)(π-acceptor)] complexes have a consistent matrix, with particle sizes ranging from 50 to 500 μm ( Figure 4).

Molecular Docking Studies
Serotonin 5-HT2A (PDB ID: 6A94) and dopamine (PDB ID: 6CM4) protein receptors were docked against all synthesized CT complexes, resulting in the optimum docking pose. For comparison reasons, FXN was utilized as control. Interestingly, five of six CT complexes showed higher docking scores than that of FXN alone at both receptors (Table  1). The interaction data and best docking pose of [(FXN)(PA)]-serotonin and (FXN)-serotonin are displayed in Figure 5 and Table 2.

Molecular Docking Studies
Serotonin 5-HT2A (PDB ID: 6A94) and dopamine (PDB ID: 6CM4) protein receptors were docked against all synthesized CT complexes, resulting in the optimum docking pose. For comparison reasons, FXN was utilized as control. Interestingly, five of six CT complexes showed higher docking scores than that of FXN alone at both receptors (Table 1). Table 1. Docking score of all [(FXN)(π-acceptor)] complexes and FXN drug alone into two protein receptors [serotonin (PDB ID: 6A94) and dopamine (PDB ID: 6CM4)].

Hydrogen Bonds, Ionizability, Hydrophobicity, and Aromatic Surfaces
The ligand-receptor complexes were assessed using DiscoveryStudio (DS) software. Mutual interactions were explored, and different surfaces were created around the ligands [35].

Hydrogen Bonds, Ionizability, Hydrophobicity, and Aromatic Surfaces
The ligand-receptor complexes were assessed using DiscoveryStudio (DS) software. Mutual interactions were explored, and different surfaces were created around the ligands [35].  Herein, fluorine atoms play a significant role in hydrogen-bond formation. The hydrogen-atom acceptor area is indicated in green, and the donor area is indicated in pink for the amino acid residues at the hydrogen-bond surface in Figure 8a and Figure S1a. The presence of the hydrophilicity features of the receptor around the ligand is confirmed by the hydrophobicity surface (Figure 8b and Figure S1b). Further, the aromatic face/edge surface (Figure 8c and Figure S1c, orange/blue = face/edge) was revealed using the docking outputs. The ionization surface reflected the acidic and basic propensity (Figure 8d and Figure S1d, blue color = basic and red color = acidic) [36].

MD Simulation and Structural Stability Analysis upon Ligand Binding
The highest docking scores from CTpS and FXNS were used as the starting structures to achieve the MD simulation at a time of 100 ns ( Figure 10). The best-docked result was used to investigate the binding process at the receptor's active region in well-defined water environments. MD simulation data were processed for structural stability studies by determining the root mean square deviation (RMSD). As demonstrated in the RMSD graphic, both CTpS and FXNS achieved steady conformation after~55 ns, with suitable RMSDs of 2.01 Å and 2.21 Å, respectively ( Figure 10). Molecules 2022, 27, x FOR PEER REVIEW 9 of 20 Herein, fluorine atoms play a significant role in hydrogen-bond formation. The hydrogen-atom acceptor area is indicated in green, and the donor area is indicated in pink for the amino acid residues at the hydrogen-bond surface in Figures 8a and S1a. The presence of the hydrophilicity features of the receptor around the ligand is confirmed by the hydrophobicity surface (Figures 8b and S1b). Further, the aromatic face/edge surface (Figures 8c and S1c, orange/blue = face/edge) was revealed using the docking outputs. The ionization surface reflected the acidic and basic propensity (Figures 8d and S1d, blue color = basic and red color = acidic) [36].

MD Simulation and Structural Stability Analysis upon Ligand Binding
The highest docking scores from CTpS and FXNS were used as the starting structures to achieve the MD simulation at a time of 100 ns ( Figure 10). The best-docked result was used to investigate the binding process at the receptor's active region in well-defined water environments. MD simulation data were processed for structural stability studies by determining the root mean square deviation (RMSD). As demonstrated in the RMSD graphic, both CTpS and FXNS achieved steady conformation after ~55 ns, with suitable RMSDs of 2.01 Å and 2.21 Å , respectively ( Figure 10). An RMSD value of <3.0 Å is widely considered to be the most acceptable and implies that the complex is suitably stable [37]. Lower RMSD values for CTpS and FXNS suggest that the ligand reduces protein flexibility, which blocks a conformational change, indicating that CTpS develops a more stable combination. The results are consistent with the theory that ligand-receptor interactions reduce the distance between protein chains and bring them closer together (Figure 11) [38]. An RMSD value of <3.0 Å is widely considered to be the most acceptable and implies that the complex is suitably stable [37]. Lower RMSD values for CTpS and FXNS suggest that the ligand reduces protein flexibility, which blocks a conformational change, indicating that CTpS develops a more stable combination. The results are consistent with the theory that ligand-receptor interactions reduce the distance between protein chains and bring them closer together ( Figure 11) [38]. An RMSD value of <3.0 Å is widely considered to be the most acceptable and implies that the complex is suitably stable [37]. Lower RMSD values for CTpS and FXNS suggest that the ligand reduces protein flexibility, which blocks a conformational change, indicating that CTpS develops a more stable combination. The results are consistent with the theory that ligand-receptor interactions reduce the distance between protein chains and bring them closer together ( Figure 11) [38]. The standard deviation and average distance between each pair of amino acids for each conformation were displayed using RR distance maps ( Figure 12) [39,40]. The standard deviation and average distance between each pair of amino acids for each conformation were displayed using RR distance maps ( Figure 12) [39,40]. The red and blue parts depict residue pairings through the highest distance differences between two amino acids, while the white diagonal represents residue pairings with a distance of zero ( Figure 12). The radius of gyration (Rg) values for CTpS and FXNS were 25.62 and 26.12, respectively. The Rg for CTpS and FXNS declined during the simulation, indicating that the structures were becoming further compressed ( Figure 13). The red and blue parts depict residue pairings through the highest distance differences between two amino acids, while the white diagonal represents residue pairings with a distance of zero ( Figure 12). The radius of gyration (Rg) values for CTpS and FXNS were 25.62 and 26.12, respectively. The Rg for CTpS and FXNS declined during the simulation, indicating that the structures were becoming further compressed ( Figure 13). The red and blue parts depict residue pairings through the highest distance differences between two amino acids, while the white diagonal represents residue pairings with a distance of zero ( Figure 12). The radius of gyration (Rg) values for CTpS and FXNS were 25.62 and 26.12, respectively. The Rg for CTpS and FXNS declined during the simulation, indicating that the structures were becoming further compressed ( Figure 13).

Hydrogen-Bond Analysis
The amount of hydrogen-bond interactions that have arisen in ligand-receptor combinations (FXNS and CTpS) were plotted against time using a grid-search at grid = 25 × 11 × 14 and rcut = 0.35 ( Figure 14).

Hydrogen-Bond Analysis
The amount of hydrogen-bond interactions that have arisen in ligand-receptor combinations (FXNS and CTpS) were plotted against time using a grid-search at grid = 25 × 11 × 14 and rcut = 0.35 ( Figure 14).

Solvent Accessibility Surface Area Analysis
The solvent accessibility surface area (SASA) was found to change because of the binding of ligand to receptor (Figure 15).

Solvent Accessibility Surface Area Analysis
The solvent accessibility surface area (SASA) was found to change because of the binding of ligand to receptor ( Figure 15).
The hydrogen bonds between ligands (FXN) or [(FXN)(PA)] with serotonin receptors were calculated (Figure 14). The average number of H-bonds per timeframe was 0.656 of 252,486 for FXN-serotonin complex and was 0.216 of 252,466 for [(FXN)(PA)]-serotonin complex. Overall, the receptor-ligand interaction was found to markedly enhance the Hbonds, with more bonds in [(FXN)(PA)]-serotonin than FXN-serotonin.

Solvent Accessibility Surface Area Analysis
The solvent accessibility surface area (SASA) was found to change because of the binding of ligand to receptor (Figure 15).    Blue represents the electropositive region, while red represents electronegative areas. These findings point to the molecule's preferred binding sites for electrophilic and nucleophilic charges [42]. The MEP surface is mapped with the color scale from deep red (−9.028 Blue represents the electropositive region, while red represents electronegative areas. These findings point to the molecule's preferred binding sites for electrophilic and nucleophilic charges [42]. The MEP surface is mapped with the color scale from deep red (−9.028 × 10 −2 ) to deep blue (+9.028 × 10 −2 ) ( Figure 16) [43].
TD-DFT was used to examine the nature of electronic transitions in water. The TD-DFT yielded two electronic absorption bands at 436 and 451 nm. HOMO → LUMO and HOMO-1 → LUMO were allocated at 451 and 436 nm, respectively. The electron acceptors were mainly LUMO and electron donors were HOMO, as observed in the FXN moiety of the CT complex [(FXN)(PA)]. Figure 17 represents the HOMO and LUMO spatial arrangements and their gaps and associated energies. Blue represents the electropositive region, while red represents electronegative areas. These findings point to the molecule's preferred binding sites for electrophilic and nucleophilic charges [42]. The MEP surface is mapped with the color scale from deep red (−9.028 × 10 −2 ) to deep blue (+9.028 × 10 −2 ) ( Figure 16) [43].
TD-DFT was used to examine the nature of electronic transitions in water. The TD-DFT yielded two electronic absorption bands at 436 and 451 nm. HOMO → LUMO and HOMO-1 → LUMO were allocated at 451 and 436 nm, respectively. The electron acceptors were mainly LUMO and electron donors were HOMO, as observed in the FXN moiety of the CT complex [(FXN)(PA)]. Figure 17 represents the HOMO and LUMO spatial arrangements and their gaps and associated energies. The molecular orbital energy level diagram of the CT complex [(FXN)(PA)] is presented in Figure 18.   [44]. Based on these DFT results, some theoretical molecular parameters related to chemical reactivity in water as solvent are presented in Table 3 [45,46].  The HOMO-LUMO and HOMO-1-LUMO gaps (∆E) for [(FXN)(PA)] were 3.646 and 3.768 eV, respectively [44]. Based on these DFT results, some theoretical molecular parameters related to chemical reactivity in water as solvent are presented in Table 3 [45,46].   Figure 19). The energy difference between HOMO and LUMO has an impact on how chemically stable molecules are. As previously reported, ligands with a smaller band gap are soft in nature, have low kinetic stability, and higher chemical reactivity [47]. The opposite is true for ligands with a wide energy gap, which are harder by nature, more unstable, and more chemically reactive [45,47].  [42,44,48,49].
At room temperature, we mixed all of the combinations for 1 h. We used filtering to remove the solid products, subsequently washed the sample with minimal volumes of chloroform, then dehydrated under vacuum with CaCl2.

Molecular Docking
We used OpenBabelIGUI software version 2.4.1 [51] to obtain the structure of FXN and the six CT complexes in PDBQT format. We used MMFF94 force field and conjugate gradient optimization algorithm to diminish the energy of the structures [52]. We used the conjugate gradient optimization algorithm using PyRx-Python prescription 0.8 for 500 steps. We employed the RCSB protein data repository to obtain the 3D crystal structures of serotonin 5-HT2A (PDB ID: 6A94) and dopamine (PDB ID: 6CM4) [53]. To add polar hydrogen atoms to the receptors and determine their Kollman charges, we utilized the Autodock Tool [54]. To assign partial charges, we employed the Geistenger approach. We used Autodock Vina to carry out the docking calculations [55]. We used a grid of 54 × 40 × 42 and center x = 13.562, center y = 0.255, center z = 61.097 with an exhaustiveness value of 8. We selected the docked positions on the basis of good docking scores and examined the interactions using DS Visualizer.

MD Simulation Study
We accomplished the MD simulation using the best-docked complexes of the receptor and ligand with a high docking score for FXN alone and [(FXN)(PA)]. We attained the topologies and parameters of the ligands by CGenFF with CHARMM-GUI [56,57]. We used online server CHARMM-gui to insert the DPPC (dipalmitoylphosphatidylcholine) In order to make the solid CT complexes, we mixed 0.309 gm (1 mmol) fluoxetin hydrochloride with 1 mmol of each π-acceptor in 20 mL of chloroform solvent [50].
At room temperature, we mixed all of the combinations for 1 h. We used filtering to remove the solid products, subsequently washed the sample with minimal volumes of chloroform, then dehydrated under vacuum with CaCl 2 .

Molecular Docking
We used OpenBabelIGUI software version 2.4.1 [51] to obtain the structure of FXN and the six CT complexes in PDBQT format. We used MMFF94 force field and conjugate gradient optimization algorithm to diminish the energy of the structures [52]. We used the conjugate gradient optimization algorithm using PyRx-Python prescription 0.8 for 500 steps. We employed the RCSB protein data repository to obtain the 3D crystal structures of serotonin 5-HT2A (PDB ID: 6A94) and dopamine (PDB ID: 6CM4) [53]. To add polar hydrogen atoms to the receptors and determine their Kollman charges, we utilized the Autodock Tool [54]. To assign partial charges, we employed the Geistenger approach. We used Autodock Vina to carry out the docking calculations [55]. We used a grid of 54 × 40 × 42 and center x = 13.562, center y = 0.255, center z = 61.097 with an exhaustiveness value of 8. We selected the docked positions on the basis of good docking scores and examined the interactions using DS Visualizer.

MD Simulation Study
We accomplished the MD simulation using the best-docked complexes of the receptor and ligand with a high docking score for FXN alone and [(FXN)(PA)]. We attained the topologies and parameters of the ligands by CGenFF with CHARMM-GUI [56,57]. We used online server CHARMM-gui to insert the DPPC (dipalmitoylphosphatidylcholine) membrane. We added lower-and upper leaflet with 72 DPPC molecules. We used SPC water models that extended 10 Å from the receptor to examine the receptor-ligand configurations in a rectangular box [58]. We administered 52 K + and 57 Cl − ions (0.15 M salt) to neutralize the systems and reproduce physiological salt concentrations ( Figure 21). In the NPT/NVT equilibration run, we subjected both systems to periodic boundary conditions at constant temperature (300 K) and pressure (1.0 bar) using Leap-frog MD integrator for 100 ns simulation time [59].

Density Functional Theory
We performed DFT analysis by Gaussian 09RevD.01 package [64]. We used the calculations to obtain a stable molecular geometry and study the electronic transitions in the CT complex [(FXN)(PA)] from a theoretical perspective. Furthermore, we applied B3LYP/6-311G++ basic to obtain the optimized structure of the CT complexes [65]. Then, we assessed the electrostatic potential map (MEP), LUMO, and HOMO spatial plots of the [(FXN)(PA)] complex [66]. The system's chemical stability is partly determined by the boundary molecular orbitals. We also investigated several structure-based molecular properties. We also used this theory to calculate structure-based molecular parameters, such as bond lengths, angles, atomic charges, total energy, electronic properties, and protein stability. We used ChemCraft 1.5 software for visualization [67].

Conclusions
We evaluated the 1:1 complexes of FXN with six π-acceptor complexes in both solid and liquid phases. Several spectroscopic analyses were used to characterize the isolated complexes. Additionally, molecular docking demonstrated that [(FXN)(PA)] has the maximum immunomolecular docking score towards the serotonin [5-HT2A (PDB ID: 6A94)] To remove inappropriate links inside the system, we reduced energy using the steepest descent approach with 5000 steps [60]. We employed the GROMACS program to accomplish the trajectory investigation [61].
We assessed the root mean square deviation by gmx rms tool. We studied hydrogen bonding by Gmxhbond instrument. We used the gmx gyrate and gmxsasa tools to calculate the gyration radius and solvent-accessible surface area. We created plots with Grace Software and created pictures with Polo/VMD [62,63].

Density Functional Theory
We performed DFT analysis by Gaussian 09RevD.01 package [64]. We used the calculations to obtain a stable molecular geometry and study the electronic transitions in the CT complex [(FXN)(PA)] from a theoretical perspective. Furthermore, we applied B3LYP/6-311G++ basic to obtain the optimized structure of the CT complexes [65]. Then, we assessed the electrostatic potential map (MEP), LUMO, and HOMO spatial plots of the [(FXN)(PA)] complex [66]. The system's chemical stability is partly determined by the boundary molecular orbitals. We also investigated several structure-based molecular properties. We also used this theory to calculate structure-based molecular parameters, such as bond lengths, angles, atomic charges, total energy, electronic properties, and protein stability. We used ChemCraft 1.5 software for visualization [67].

Conclusions
We evaluated the 1:1 complexes of FXN with six π-acceptor complexes in both solid and liquid phases. Several spectroscopic analyses were used to characterize the isolated complexes. Additionally, molecular docking demonstrated that [(FXN)(PA)] has the maximum immunomolecular docking score towards the serotonin [5-HT2A (PDB ID: 6A94)] receptor, and was more effective than the FXN reactant donor alone. It was also found that at 100 ns MD simulation, [(FXN)(PA)]-serotonin was more stable than (FXN)-serotonin. Furthermore, the theoretical data from the DFT calculations helped to explore the molecular geometry of the CT complexes along with their important physical parameters. The stability of the CT complex [(FXN)(PA)] was also studied by DFT calculations.