Synthesis, Crystal Structure, DFT Calculations, Hirshfeld Surface Analysis and In Silico Drug-Target Profiling of (R)-2-(2-(1,3-Dioxoisoindolin-2-yl)propanamido)benzoic Acid Methyl Ester

The work here reflects synthesis, DFT studies, Hirshfeld charge analysis and crystal data exploration of pharmacologically important (R)-2-(2-(1,3-dioxoisoindolin-2-yl)propanamido)benzoic acid methyl ester (5) to understand its properties for further chemical transformations. The methyl anthranilate (2) was produced by the esterification of anthranilic acid in an acidic medium. The phthaloyl-protected alanine (4) was rendered by the fusion of alanine with phthalic anhydride at 150 °C, followed by coupling with (2) furnished isoindole (5). The characterization of products was performed using IR, UV-Vis, NMR and MS. Single-crystal XRD also verified the structure of (5) in which N-H⋯O bonding stabilizes the molecular configuration of (5), resulting in the formation of S(6) hydrogen-bonded loop. The molecules of isoindole (5) are connected in the form of dimers, and the π⋯π stacking interaction between aromatic rings further stabilizes the crystal packing. DFT studies suggest that HOMO is over the substituted aromatic ring, the LUMO is present mainly over the indole side, and nucleophilic and electrophilic corners point out the reactivity of the product (5). In vitro and in silico analysis of (5) shows its potential as an antibacterial agent targeting DNA gyrase and Dihydroorotase from E. coli and tyrosyl-tRNA synthetase and DNA gyrase from Staphylococcus aureus.


Introduction
Isoindoles belong to an important class of heterocyclic compounds of pharmaceutical importance. The substituted isoindoles have significant therapeutic potential and remain a long-standing focus of dynamic scientific interest. It is a large group of salutary agents showing anti-HIV [1], anticancer [2], antipsychotic [3], antianxiety [4,5], antiinflammatory [6], antimicrobial [7,8], antiarrhythmic [9] and anticonvulsant [10] activities. The molecular modification and derivatization of available drugs have led to enormously improved medicines.
We have already successfully synthesized chiral seven-membered rings of benzoxazepines [11]. During our current studies, attempts were made to synthesize chiral 1,4benzodiazepines. We have tried to develop a simple method of preparation of asymmetric benzodiazepines by coupling methyl ester of anthranilic acid with phthaloyl-protected (S)-alanine. The reaction yielded a crystalline product as an exclusive product with a significant yield. Spectroscopic characterization indicated that the desired product was not formed; instead, an isoindole (5) was formed, which failed to be further converted into benzodiazepine 7 (Scheme 1). The resulting substituted isoindole 5 is an attractive pharmacophore for different activities and is important in medicines and drug chemistry. Considering the structure of the isoindole with anthranilic acid moiety, isoindole moiety and amino acid residue, we decided to apply DFT and Hirshfeld studies on the synthesized product to have a deeper insight into the chemical behavior and reactivity of this interesting product to further derivatize it as a potential drug. Considering the pharmacological importance mentioned above, we also decided to check the pharmacological screening of the achieved isoindole. a long-standing focus of dynamic scientific interest. It is a large group of salutary agents showing anti-HIV [1], anticancer [2], antipsychotic [3], antianxiety [4,5], anti-inflammatory [6], antimicrobial [7,8], antiarrhythmic [9] and anticonvulsant [10] activities. The molecular modification and derivatization of available drugs have led to enormously improved medicines.
We have already successfully synthesized chiral seven-membered rings of benzoxazepines [11]. During our current studies, attempts were made to synthesize chiral 1,4-benzodiazepines. We have tried to develop a simple method of preparation of asymmetric benzodiazepines by coupling methyl ester of anthranilic acid with phthaloyl-protected (S)-alanine. The reaction yielded a crystalline product as an exclusive product with a significant yield. Spectroscopic characterization indicated that the desired product was not formed; instead, an isoindole (5) was formed, which failed to be further converted into benzodiazepine 7 (Scheme 1). The resulting substituted isoindole 5 is an attractive pharmacophore for different activities and is important in medicines and drug chemistry. Considering the structure of the isoindole with anthranilic acid moiety, isoindole moiety and amino acid residue, we decided to apply DFT and Hirshfeld studies on the synthesized product to have a deeper insight into the chemical behavior and reactivity of this interesting product to further derivatize it as a potential drug. Considering the pharmacological importance mentioned above, we also decided to check the pharmacological screening of the achieved isoindole.

Results and Discussions
Anthranilic acid was esterified to form compound 2 to protect its carboxy-terminal (C-terminal). Further, the coupling of ester 2 with phthalamide 4 formed an expected coupled product 5. The deprotection (ester and phthaloyl moieties) followed by conversion of resulting carboxylic acid 6 to COCl2 and the intramolecular cyclization was expected to construct diazepine ring 7, which was abandoned due to the conversion of 5 back into 4 instead of the desired benzodiazepine 7 (Scheme 1) under applied conditions for cyclization reaction. The substituted isoindole 5 was prepared during the proposed route adopted for synthesizing benzodiazepine 7, which was characterized by using different spectroscopic and mass spectrometric techniques followed by DFT studies.

Characterization
The esterification of anthranilic acid employing acid catalysis formed 2 as purplecolored oil with a 50% yield by procedure mentioned in the literature [12]. The absence of the O-H signal and the presence of a broad C=O signal in the IR spectrum indicated successful esterification. Although the 1 H-NMR spectrum of ester 2 showed an aromatic re-

Results and Discussions
Anthranilic acid was esterified to form compound 2 to protect its carboxy-terminal (C-terminal). Further, the coupling of ester 2 with phthalamide 4 formed an expected coupled product 5. The deprotection (ester and phthaloyl moieties) followed by conversion of resulting carboxylic acid 6 to COCl 2 and the intramolecular cyclization was expected to construct diazepine ring 7, which was abandoned due to the conversion of 5 back into 4 instead of the desired benzodiazepine 7 (Scheme 1) under applied conditions for cyclization reaction. The substituted isoindole 5 was prepared during the proposed route adopted for synthesizing benzodiazepine 7, which was characterized by using different spectroscopic and mass spectrometric techniques followed by DFT studies.

Characterization
The esterification of anthranilic acid employing acid catalysis formed 2 as purplecolored oil with a 50% yield by procedure mentioned in the literature [12]. The absence of the O-H signal and the presence of a broad C=O signal in the IR spectrum indicated successful esterification. Although the 1 H-NMR spectrum of ester 2 showed an aromatic region similar to that of the reactant, three additional proton singlet peaks appeared, supporting the formation of methyl anthranilate 2. The molecular ion peak of ester appeared at 151 amu, having a 14 amu increase on the molecular ion peak of anthranilic acid, confirming successful esterification. The MeO radical and CO were lost successively, supporting a logical fragmentation pattern. The rearranged fragment (β-keto acid) generated via the MeOH removal from [M +• ] is attributed as the base peak (119 amu) (Scheme 2) (Supplementary Materials; Compound characterization data). gion similar to that of the reactant, three additional proton singlet peaks appeared, supporting the formation of methyl anthranilate 2. The molecular ion peak of ester appeared at 151 amu, having a 14 amu increase on the molecular ion peak of anthranilic acid, confirming successful esterification. The MeO radical and CO were lost successively, supporting a logical fragmentation pattern. The rearranged fragment (β-keto acid) generated via the MeOH removal from [M +• ] is attributed as the base peak (119 amu) (Scheme 2) (Supplementary Materials; Compound characterization data). It was essential to protect the NH2 terminal of (S)-alanine 2 to avoid interference during the coupling reaction. The N-protection of (S)-alanine 3 was carried out by the fusion of amino acid with phthalic anhydride [13]. The removal of unreacted phthalic anhydride was easy after the completion of the reaction due to sublimation. The reaction is highly attractive for being environmentally friendly and efficient with a high yield (70-80%). This condensation resulted in the formation of crystalline 4 with significant yields (Table 1). Compound 4 exhibited molecular ion peaks that were in accordance with the expected cyclized product for successful N-protection. The loss of OH radical and CO or CO2 gave daughter fragments, confirming the formation of the required product (Scheme 3). The signals of four aromatic region protons and the presence of peaks in the aliphatic region of 1 H-NMR also suggested successful N-protection. In product 4, H 2 appeared at 7.80-7.90 ppm as a multiplet (Table 2). Concrete evidence for N-protection of 3 was shown by single-crystal XRD ( Figure 1).
It was essential to protect the NH 2 terminal of (S)-alanine 2 to avoid interference during the coupling reaction. The N-protection of (S)-alanine 3 was carried out by the fusion of amino acid with phthalic anhydride [13]. The removal of unreacted phthalic anhydride was easy after the completion of the reaction due to sublimation. The reaction is highly attractive for being environmentally friendly and efficient with a high yield (70-80%). This condensation resulted in the formation of crystalline 4 with significant yields (Table 1). Compound 4 exhibited molecular ion peaks that were in accordance with the expected cyclized product for successful N-protection. The loss of OH radical and CO or CO 2 gave daughter fragments, confirming the formation of the required product (Scheme 3).
porting the formation of methyl anthranilate 2. The molecular ion peak of ester appeared at 151 amu, having a 14 amu increase on the molecular ion peak of anthranilic acid, confirming successful esterification. The MeO radical and CO were lost successively, supporting a logical fragmentation pattern. The rearranged fragment (β-keto acid) generated via the MeOH removal from [M +• ] is attributed as the base peak (119 amu) (Scheme 2) (Supplementary Materials; Compound characterization data). It was essential to protect the NH2 terminal of (S)-alanine 2 to avoid interference during the coupling reaction. The N-protection of (S)-alanine 3 was carried out by the fusion of amino acid with phthalic anhydride [13]. The removal of unreacted phthalic anhydride was easy after the completion of the reaction due to sublimation. The reaction is highly attractive for being environmentally friendly and efficient with a high yield (70-80%). This condensation resulted in the formation of crystalline 4 with significant yields (Table 1). Compound 4 exhibited molecular ion peaks that were in accordance with the expected cyclized product for successful N-protection. The loss of OH radical and CO or CO2 gave daughter fragments, confirming the formation of the required product (Scheme 3). The signals of four aromatic region protons and the presence of peaks in the aliphatic region of 1 H-NMR also suggested successful N-protection. In product 4, H 2 appeared at 7.80-7.90 ppm as a multiplet (Table 2). Concrete evidence for N-protection of 3 was shown by single-crystal XRD ( Figure 1). The signals of four aromatic region protons and the presence of peaks in the aliphatic region of 1 H-NMR also suggested successful N-protection. In product 4, H 2 appeared at 7.80-7.90 ppm as a multiplet (Table 2). Concrete evidence for N-protection of 3 was shown by single-crystal XRD ( Figure 1). The acid chloride of N-protected alanine 4 was synthesized by refluxing it (1 h) with SOCl 2 . This acid chloride was then coupled with ester 2 by heating for 4 h at 40 • C, yielding colorless crystals of isoindole 5.
The IR spectrum of the isoindole 5 showed distinctive absorbance analogous to the three C=O groups at 1770, 1707, 1595 cm −1 and one N-H at 3269 cm −1 . The 1 H-NMR identified eight aromatic protons. Phthaloyl 4H, in the same electronic environment, appeared as two multiplets, each with 2H. The aliphatic proton signals also agreed with the expected values, indicating a doublet of three CH 3 protons, a quartet of H 1 and a singlet relative to 3H of OMe. The anthranilate portion showed the presence of four aromatic protons. The splitting pattern of this portion was similar to that of methyl anthranilate 2; however, all these protons were relatively deshielded post-coupling, such as the H 3 signal that jumped most downfield from 6.75 ppm to 8.70 ppm.
Molecules 2023, 28, x FOR PEER REVIEW The acid chloride of N-protected alanine 4 was synthesized by refluxing it (1 SOCl2. This acid chloride was then coupled with ester 2 by heating for 4 h at 40 °C, y colorless crystals of isoindole 5.
The IR spectrum of the isoindole 5 showed distinctive absorbance analogou three C=O groups at 1770, 1707, 1595 cm -1 and one N-H at 3269 cm -1 . The 1 H-NM tified eight aromatic protons. Phthaloyl 4H, in the same electronic environment, ap as two multiplets, each with 2H. The aliphatic proton signals also agreed with pected values, indicating a doublet of three CH3 protons, a quartet of H 1΄ and a relative to 3H of OMe. The anthranilate portion showed the presence of four ar protons. The splitting pattern of this portion was similar to that of methyl anthran however, all these protons were relatively deshielded post-coupling, such as the H 3 that jumped most downfield from 6.75 ppm to 8.70 ppm.
The low field shift is justified due to the functional group interconversion of to the amido group; the amino group has +R while the amido group shows the −R The field effect that deshields the proton H 3″ is also a reason for the ppm shift ( Fig   Figure 2. The electronic environment within the isoindole 5.
The 13 C-NMR analysis also supported the coupled product structure. The fo The low field shift is justified due to the functional group interconversion of amino to the amido group; the amino group has +R while the amido group shows the −R effect. The field effect that deshields the proton H 3 is also a reason for the ppm shift ( Figure 2). The acid chloride of N-protected alanine 4 was synthesized by refluxing it (1 SOCl2. This acid chloride was then coupled with ester 2 by heating for 4 h at 40 °C, y colorless crystals of isoindole 5.
The IR spectrum of the isoindole 5 showed distinctive absorbance analogou three C=O groups at 1770, 1707, 1595 cm -1 and one N-H at 3269 cm -1 . The 1 H-NM tified eight aromatic protons. Phthaloyl 4H, in the same electronic environment, ap as two multiplets, each with 2H. The aliphatic proton signals also agreed with pected values, indicating a doublet of three CH3 protons, a quartet of H 1΄ and a relative to 3H of OMe. The anthranilate portion showed the presence of four a protons. The splitting pattern of this portion was similar to that of methyl anthran however, all these protons were relatively deshielded post-coupling, such as the H that jumped most downfield from 6.75 ppm to 8.70 ppm.
The low field shift is justified due to the functional group interconversion o to the amido group; the amino group has +R while the amido group shows the −R The field effect that deshields the proton H 3″ is also a reason for the ppm shift ( Fig   Figure 2. The electronic environment within the isoindole 5. The 13 C-NMR analysis also supported the coupled product structure. The fo (CO2Me, C 3 , C 1 and C 2′ ) appeared at 172.0, 168.8, 168.0 and 167.9 ppm, respe whereas eight CHs of the aromatic region and four q Cs appeared in the range 120. The 13 C-NMR analysis also supported the coupled product structure. The four C=O (CO 2 Me, C 3 , C 1 and C 2 ) appeared at 172.0, 168.8, 168.0 and 167.9 ppm, respectively, whereas eight CHs of the aromatic region and four q Cs appeared in the range 120.5-134.2 ppm. The The high-resolution mass spectrum ultimately showed a peak of [M + Na] + at 375.0943 amu, which agrees with the targeted product's calculated mass. In the LR EIMS, the [M] +• appeared at 352 amu. The further fragmentation mode and base peak also strengthened the proof for production of 5. The MeOH removal gave a daughter fragment at 320 amu (Scheme 4). Further reduction of 5 to yield 6 was a complete failure, and further conversion into 7 could not be carried out.

DFT Exploration
The three-dimensional structure of compound (5), as shown in Figure 3 was modeled using Gauss View, and the geometry was then optimized at the PBE0-D3BJ/def2-TZVP level of density functional theory. The structure of compound (5) is a true minimum, as indicated by the absence of any imaginary frequency. Comparing the calculated structure with the crystallographic information shows that both results correlate well. For instance, the N8-C21 bond length is 1.460 Å, while N6-C28 is 1.396 Å in the crystal structure. The calculated bond lengths are 1.439 Å and 1.385 Å, with negligible deviation from the experimental bond lengths ( Figure 3).  (5) at PBE0-D3BJ/def2-TZVP level of theory. The grey color shows carbon in 3D models; red is for oxygen, and blue shows nitrogen atoms. Hydrogen atoms are omitted for clarity.

Frontier Molecular Orbital Analysis (FMO) and Polarizability
FMO is one of the important tools for analyzing reactivity and describing many other physical parameters, including the energies of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbitals (LUMO) [14]. The calculated HOMO-LUMO energy difference in the present case was 4.59 eV. Figure 4 shows the isosurface plots of HOMO and LUMO of the compound (5). It can be seen that HOMO is localized over the substituted aromatic ring (right side), and the LUMO exists mainly over the indole side. This further suggests that the nucleophilic side of the molecule is the right side (according to Figure 4), while the indole side is electrophilic. The hyperpolarizability (β) of compound (1) is 223.05 a.u., which shows that the molecule has a moderate nonlinear optical response ( Figure 4). Further reduction of 5 to yield 6 was a complete failure, and further conversion into 7 could not be carried out.

DFT Exploration
The three-dimensional structure of compound (5), as shown in Figure 3 was modeled using Gauss View, and the geometry was then optimized at the PBE0-D3BJ/def2-TZVP level of density functional theory. The structure of compound (5) is a true minimum, as indicated by the absence of any imaginary frequency. Comparing the calculated structure with the crystallographic information shows that both results correlate well. For instance, the N8-C21 bond length is 1.460 Å, while N6-C28 is 1.396 Å in the crystal structure. The calculated bond lengths are 1.439 Å and 1.385 Å, with negligible deviation from the experimental bond lengths (Figure 3). strengthened the proof for production of 5. The MeOH removal gave a daughter fragment at 320 amu (Scheme 4). Further reduction of 5 to yield 6 was a complete failure, and further conversion into 7 could not be carried out.

DFT Exploration
The three-dimensional structure of compound (5), as shown in Figure 3 was modeled using Gauss View, and the geometry was then optimized at the PBE0-D3BJ/def2-TZVP level of density functional theory. The structure of compound (5) is a true minimum, as indicated by the absence of any imaginary frequency. Comparing the calculated structure with the crystallographic information shows that both results correlate well. For instance, the N8-C21 bond length is 1.460 Å, while N6-C28 is 1.396 Å in the crystal structure. The calculated bond lengths are 1.439 Å and 1.385 Å, with negligible deviation from the experimental bond lengths ( Figure 3).  (5) at PBE0-D3BJ/def2-TZVP level of theory. The grey color shows carbon in 3D models; red is for oxygen, and blue shows nitrogen atoms. Hydrogen atoms are omitted for clarity.

Frontier Molecular Orbital Analysis (FMO) and Polarizability
FMO is one of the important tools for analyzing reactivity and describing many other physical parameters, including the energies of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbitals (LUMO) [14]. The calculated HOMO-LUMO energy difference in the present case was 4.59 eV. Figure 4 shows the isosurface plots of HOMO and LUMO of the compound (5). It can be seen that HOMO is localized over the substituted aromatic ring (right side), and the LUMO exists mainly over the indole side. This further suggests that the nucleophilic side of the molecule is the right side (according to Figure 4), while the indole side is electrophilic. The hyperpolarizability (β) of compound (1) is 223.05 a.u., which shows that the molecule has a moderate nonlinear optical response ( Figure 4).  (5) at PBE0-D3BJ/def2-TZVP level of theory. The grey color shows carbon in 3D models; red is for oxygen, and blue shows nitrogen atoms. Hydrogen atoms are omitted for clarity.

Frontier Molecular Orbital Analysis (FMO) and Polarizability
FMO is one of the important tools for analyzing reactivity and describing many other physical parameters, including the energies of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbitals (LUMO) [14]. The calculated HOMO-LUMO energy difference in the present case was 4.59 eV. Figure 4 shows the isosurface plots of HOMO and LUMO of the compound (5). It can be seen that HOMO is localized over the substituted aromatic ring (right side), and the LUMO exists mainly over the indole side. This further suggests that the nucleophilic side of the molecule is the right side (according to Figure 4), while the indole side is electrophilic. The hyperpolarizability (β) of compound (1) is 223.05 a.u., which shows that the molecule has a moderate non-linear optical response (Figure 4).

Hirshfeld Charge Analysis
Hirshfeld charge analysis provides an insight into the spatial distribution of electron density in space throughout the molecule. This aids in understanding the intermolecular interactions within crystal packing [15]. Table 3 shows Hirshfeld charges of a few key atoms of the compound (5). It may be deduced that the oxygen (O1, O2, O3) and nitrogen (N8) atoms on the isoindole side are comparatively more negatively charged than the atoms on the other side. This, in turn, causes a positive charge on the adjacent carbon (C9, C20) atoms in the former case than those in the latter. However, C38 on the substituted aromatic ringside is the most positively charged due to the attachment of two oxygen atoms (Table 3).  Figure 3; * H.C. = Hirshfeld charge.

Crystallographic Data
The ultimate evidence for coupled product 5 came from the single-crystal XRD structure. In 5 ( Figure 5, Table 4

Hirshfeld Charge Analysis
Hirshfeld charge analysis provides an insight into the spatial distribution of electron density in space throughout the molecule. This aids in understanding the intermolecular interactions within crystal packing [15]. Table 3 shows Hirshfeld charges of a few key atoms of the compound (5). It may be deduced that the oxygen (O1, O2, O3) and nitrogen (N8) atoms on the isoindole side are comparatively more negatively charged than the atoms on the other side. This, in turn, causes a positive charge on the adjacent carbon (C9, C20) atoms in the former case than those in the latter. However, C38 on the substituted aromatic ringside is the most positively charged due to the attachment of two oxygen atoms (Table 3). Atom numbering (in the format symbol-number) is from Figure 3; * H.C. = Hirshfeld charge.

Crystallographic Data
The ultimate evidence for coupled product 5 came from the single-crystal XRD structure. In 5 ( Figure 5, Table 4     The ring comprising groups A and C has a dihedral angle of 73.27 (3) • . The C-atom (C10) is at the distance of 0.8382 (3) Å from the plane of moiety B. The important bond angles and bond lengths are specified in Table 5. Table 5. Selected bond lengths in angstrom and bond angles in degrees in 5.
In addition to H-bonding, the crystal packing is further stabilized by offset π· · · π stacking interaction as shown in Figure 7. A strong offset stacking π· · · π interaction is formed between the phenyl ring Cg2 (C2-C7) of a molecule located in an asymmetric location and the phenyl ring Cg2 (C2-C7) located at (1-x, 1-y, 1-z) with an inter-centroid distance of 3.928 Å and ring offset of 0.424 Å, where Cg is the ring centroid. In general, the offset π· · · π stacking interaction is found between aromatic rings with inter-centroid distances ranging from 3.928 to 4.252 Å. The ring offset ranges from 0.776 to 2.060 Å. In addition to H-bonding, the crystal packing is further stabilized by offset π⋯π stacking interaction as shown in Figure 7. A strong offset stacking π⋯π interaction is formed between the phenyl ring Cg2 (C2-C7) of a molecule located in an asymmetric location and the phenyl ring Cg2 (C2-C7) located at (1-x, 1-y, 1-z) with an inter-centroid distance of 3.928 Å and ring offset of 0.424 Å, where Cg is the ring centroid. In general, the offset π⋯π stacking interaction is found between aromatic rings with inter-centroid distances ranging from 3.928 to 4.252 Å. The ring offset ranges from 0.776 to 2.060 Å.  In single crystals, the crystal packing is stabilized by the non-covalent interactions. For the exploration of these interactions, the Hirshfeld surface (HS) is performed on Crystal Explorer, version 17.5 [16]. Color coding shows the non-covalent interactions on a Hirshfeld surface displayed on dnorm [17]. Red spots on HS indicate contacts for which the distance between atoms is smaller than the sum of the Van der Waal radii. White and blue regions represent weak and negligibly small interactions, respectively, as shown in Figure   Figure 7. Graphical presentations of offset π· · · π interaction for (5) (hydrogen atoms are not shown to avoid ambiguity, and Å is used for distance calculation).
In single crystals, the crystal packing is stabilized by the non-covalent interactions. For the exploration of these interactions, the Hirshfeld surface (HS) is performed on Crystal Explorer, version 17.5 [16]. Color coding shows the non-covalent interactions on a Hirshfeld surface displayed on d norm [17]. Red spots on HS indicate contacts for which the distance between atoms is smaller than the sum of the Van der Waal radii. White and blue regions represent weak and negligibly small interactions, respectively, as shown in Figure 8a. No red spot was found over the HS near the imine NH group, indicating the absence of N-H· · · N and N-H· · · O types of intermolecular H-bonding. Red spots on the HS are found near carbonyl O-atoms (O1/O2) of isoindoline-1,3-dione group A, CH of methoxy group, CH of carbon directly attached with group A and CH of the methyl group (C10), indicating that these atoms are engaged in H-bonding. Hirshfeld surface analysis can be used to explore the π· · · π stacking interaction between aromatic rings. For this, HS is plotted over the shape index. On this HS, the aromatic rings are seen to be involved in π· · · π stacking interaction due to the existence of consecutive triangular regions of red and blue colour (Figure 8b). Intermolecular interactions in terms of individual interatomic contacts can be inspected by plotting 2D plots in which de represents the separation of HS to the closest nucleus outside the HS, whereas this distance from HS to the closest nucleus inside HS is represented by di. Figure 9a shows a 2D plot for overall interactions. The central sky-blue zones of nearly triangular shape show that π⋯π stacking interaction is present in the crystal packing. Each interatomic contact is determined by considering its reciprocal contact. The interatomic contacts of type H⋯H, O⋯H, C⋯C and C⋯H are found to be the significant participants in crystal packing with %age contribution of 44.7% (Figure 9b  Intermolecular interactions in terms of individual interatomic contacts can be inspected by plotting 2D plots in which d e represents the separation of HS to the closest nucleus outside the HS, whereas this distance from HS to the closest nucleus inside HS is represented by d i . Figure 9a shows a 2D plot for overall interactions. The central sky-blue zones of nearly triangular shape show that π· · · π stacking interaction is present in the crystal packing. Each interatomic contact is determined by considering its reciprocal contact. The interatomic contacts of type H· · · H, O· · · H, C· · · C and C· · · H are found to be the significant participants in crystal packing with %age contribution of 44.7% (Figure 9b), 29.9% (Figure 9c), 10.3% (Figure 9d) and 10.1% (Figure 9e), respectively. The interatomic contacts that have minor contribution to crystal packing are O· · · C (2.1%, Figure 9f), N· · · C (2%, Figure 9g), O· · · O (0.4%, Figure 9h), N· · · H (0.3%, Figure 9i) and O· · · N (0.1%, Figure 9j).
To learn more about crystal packing, computations were performed to determine how all the atoms present inside the HS interacts with the atom present outside the HS. Such types of interactions are quantized, as shown in Figure 10a, in which ALL represent all the atoms present inside HS. ALL-H interaction is the strongest interface, with a percentage impact of 66.5%. Other such interactions are ALL-N, ALL-O and ALL-C, with respective percentage contributions of 1.3%, 15.4% and 16.8% (Figure 10a). Likewise, the interaction of an atom within HS with atoms of molecules in its surroundings was investigated. In such interactions, ALL represents all the atoms of molecules present in the HS surrounding. H-ALL interaction is found to be the strongest interaction, with a percentage impact of  centage impact of 66.5%. Other such interactions are ALL-N, ALL-O and ALL-C, with respective percentage contributions of 1.3%, 15.4% and 16.8% (Figure 10a). Likewise, the interaction of an atom within HS with atoms of molecules in its surroundings was investigated. In such interactions, ALL represents all the atoms of molecules present in the HS surrounding. H-ALL interaction is found to be the strongest interaction, with a percentage impact of 63.3%. Other such interactions are N-ALL, O-ALL and C-ALL, with respective percentage impacts of 1.3%, 17.5% and 18% (Figure 10b). In the unit cell, the percentage space occupied by the molecules define the crystal packing strength. If most of the space in the unit cell of a single crystal is occupied by voids, then it means crystal packing is weak. The void analysis was performed by assuming that the electron density around every atom is spherical. Graphical view of voids in compound 5 is shown in Figure 11. The volume of voids in the unit cell was 254.51 Å. In the unit cell, 14.6% of space was occupied by voids; in other words, 85.4% was occupied by atoms indicating that crystal packing is strong and has no large cavity. In the unit cell, the percentage space occupied by the molecules define the crystal packing strength. If most of the space in the unit cell of a single crystal is occupied by voids, then it means crystal packing is weak. The void analysis was performed by assuming that the electron density around every atom is spherical. Graphical view of voids in compound 5 is shown in Figure 11. The volume of voids in the unit cell was 254.51 Å. In the unit cell, 14.6% of space was occupied by voids; in other words, 85.4% was occupied by atoms indicating that crystal packing is strong and has no large cavity.

General Experimental Procedure
Freshly distilled analytical-grade organic solvents were used. The EtOAc, CHCl 3 , MeOH and EtOH were obtained from commercial sources. For crystallization and chro-matography, a fraction of distilled hexane at 60-80 • C was employed. All solid chemicals were used as such.
Gallenkamp was used to record melting points (England, Serial MFB 595 010M) and reported uncorrected. For UV spectra, a Thermo Spectronic (UV-1700) spectrometer was used, and IR spectra were taken on IR Prestige 21 FTIR spectrophotometer, Shimadzu (Japan) at High Tech Laboratory, UoS, Sargodha. Molar absorptivity ε and absorption maxima are both reported. Equation ε = A/cl, was used to find ε values.
LR EIMS (low-resolution electron impact mass spectra) were recorded on Fisons Instrument VG Autospec mass spectrometer at ICCBS (International Center for Chemical and Biological Sciences), HEJRIC, University of Karachi, Karachi, while electrospray ionization mass spectra (ESI MS) were recorded using Q-TOF Ultima API (Micromass) at BMSF (Biomedical Mass Spectrometry Facility), UNSW (University of New Southwales), Sydney, Australia.
The NMR spectra of our products were recorded at ICCBS, HEJRIC, University of Karachi, Karachi, Pakistan, and at the School of Chemistry, University of New South Wales (UNSW), Sydney, Australia. The 1 H NMR spectral data are shown as follows: chemical shifts in ppm downfield with reference to TMS (δ); proton number/count; multiplicity due to neighboring protons; observed J (coupling constant) in Hertz (Hz). Multiplicities are recorded in the form of s (singlet), bs (broad singlet), d (doublet), t (triplet), q (quartet), p (quintet), m (multiplet), dd (doublet of doublet), dt (doublet of triplet) and their combinations. Chemical shifts δ in the 13 C NMR are reported downfield from TMS in ppm and show all identifiable carbons.
For recording XRD data to elucidate the structure of crystals, Bruker Kappa APEX 11 CCD diffractometer at the Department of Physics, UoS, was used. For the radiation source, fine-focus sealed tube monochromator graphite (Mo Kα radiation, λ = 0.71073 Å) was used.
We employed density functional tools under the realm of computational studies using the software Gaussian 09 [18]. We used PBE0 hybrid functional [19] and triple ζ basis set def2-TZVP [20] for our calculations. Dispersion interactions were included during optimization by adding Grimme's empirical dispersion correction with Becke-Johnston damping (D3BJ) [17,21,22]. Calculations for frequency were performed to confirm the structure as a minimum on the potential energy surface. (2) The anthranilic acid 1 (1.00 g, 7.3 mmol, 1 eq) was dissolved in MeOH (7 mL), and 5.1 mL conc. H 2 SO 4 (9.3 g, 0.095 mol, 13 eq) was added. The resulting solution was refluxed for 72 h while the reaction was monitored by TLC to ensure the completion of the reaction. A 1N aq. NaHCO 3 solution was used to neutralize the solution and was extracted with n-hexane (2 × 20 mL). The n-hexane layer was filtered, dried over anhydrous Na 2 SO 4 and rotary concentrated to obtain a purple-colored ester 2 (50%) in the liquid phase. R f ; 0.6 (n-hexane: EtOAc, 1:1); density, 1.17 g/cm 3  The (S)-alanine 3 (1 g, 0.02 mol, 1 eq) was fused with phthalic anhydride (0.03 mol, 1.1 eq) at about 155 • C for 60 min. The reaction completion was confirmed by the disappearance of reactants on TLC. The ninhydrin dip confirmed the absence of unreacted amino acids. The reaction flask was cooled to ambient temperature, and the unreacted crystallized phthalic anhydride was scratched off from the flask wall. The rest of the solid product was recrystallized from EtOH, and H 2 O yielded colourless prisms such as crystals of phthalamide 4 (750).  (5) The SOCl 2 (0.7 mL, 7 eq) was added to phthalamide 4 (300 mg, 0.85 mmol, 1 eq) and heated at 100 • C until a clear solution was obtained. The addition of ester 2 (1 eq, 0.18 mL) was carried out to the reaction solution and was further heated for 5 h at 50 • C. The cooled solution was partitioned between H 2 O and CHCl 3 (3×15 mL). The organic layer was filtered, dried over anhydrous Na 2 SO 4 and concentrated under reduced pressure, which resulted in a pale, yellow-colored solid. It was crystallized from C 2 H 5 OH to form white-colored crystals (40%).

In Vitro and In Silico Analysis
Disk diffusion assay was performed as previously described [23]. Clinical isolates of E. coli and S. aureus were used to assess antibacterial activities, while Ciprofloxacin (1 mg/mL) was used as a standard drug. In silico studies were performed as described previously by our group. The protein structures were downloaded from the Protein Data Bank (PDB, https://www.rcsb.org/, accessed on 12 December 2022), and YASARA Structure ver. 22.10 was used for molecular docking with parameters used before [24]. The ligand-protein docking was visualized with PyMol ver. 2.5.4 [25] and LigPlot+ ver. 2.2.5 [26].

Conclusions
The structure of isoindole 5 is verified by the SC-XRD technique. The molecular configuration of isoindole 5 is stabilized due to N-H· · · O bonds that form a S (6) hydrogenbonded loop. The molecules are connected in the form of dimers by C-H· · · O bonds to make a R 2 2 (11) loop. The π· · · π stacking interaction between aromatic rings with intercentroid distance range from 3.602 to 4.252 Å helps with additional stability of crystal packing. DFT calculations were performed for compound (5). These provided us with an understanding of the geometrical features. FMO analysis was performed to observe the reactivity pattern of the compound and also to calculate the energy difference between the HOMO and the LUMO. Hyperpolarizability studies showed that the compound showed a moderate non-linear optical response. DFT studies showed HOMO's localization is localized over the substituted aromatic ring (right side), and the LUMO exists mainly over benzazole. This further suggests that the nucleophilic side of the molecule is the right side (according to Figure 6), while the indole side is electrophilic. A study of Hirshfeld charges helped with the analysis of charge distribution over the molecule. All the data points out the chemical reactivity of the product (5). Preliminary studies regarding in vitro and in silico antimicrobial activities of (5) against E. coli and S. aureus shows this compound can be a good nucleus for developing a further antibacterial compound.