Green Biosynthesis of Silver Nanoparticles Using Annona glabra and Annona squamosa Extracts with Antimicrobial, Anticancer, Apoptosis Potentials, Assisted by In Silico Modeling, and Metabolic Profiling

Annona glabra L. (AngTE) and Annona squamosa L. (AnsTE) fruits have been widely used in cancer treatment. Accordingly, their extracts were used to synthesize silver nanoparticles via a biogenic route (Ang-AgNPs) and (Ans-AgNPs), respectively. Chemical profiling was established using UPLC-QTOF-MS/MS. All species were tested for anticancer activity against human cervical cancer cells (HeLa), prostate adenocarcinoma metastatic (PC3), and ovary adenocarcinoma (SKOV3) using sulphorhodamine B assay. Apoptosis was determined using Annexin flow cytometry along with cell cycle analysis and supported by a molecular docking. The antibacterial and synergistic effect when combined with gentamicin were evaluated. A total of 114 compounds were tentatively identified, mainly acetogenins and ent-kaurane diterpenes. AnsTE and Ans-AgNPs had the most potent cytotoxicity on HeLa and SKOV3 cells, inducing a significant apoptotic effect against all tumor cells. The AnsTE and Ans-AgNPs significantly arrested PC3, SKOV3, and HeLa cells in the S phase. The nanoparticles demonstrated greater antibacterial and antifungal activities, as well as a synergistic effect with gentamicin against P. aeruginosa and E. coli. Finally, a molecular docking was attempted to investigate the binding mode of the identified compounds in Bcl-2 proteins’ receptor, implying that the fruits and their nanoparticles are excellent candidates for treating skin infections in patients with ovarian or prostatic cancer.


Introduction
Annona squamosa L. and A. glabra L., two species in the family Annonaceae, have been widely investigated for their therapeutic potential, considering that the major components of Annonaceous plants are acetogenins. Both species contain various phytochemicals molecular formulas, major fragment ions, and the identities of the identified compounds are recorded in Table S1. The UPLC chromatograms of A. squamosa and A. glabra ethanolic extracts are shown in Figure S1.

Acetogenins
Natural substances known as annonaceous acetogenins (ACGs) are only present in the Annonaceae species. They are fatty acid derivatives that have an alkanol moiety attached to the 2,4-disubstituted -lactone ring at position 2, where a long hydrocarbon chain that frequently contains a number of oxygenated moieties such as hydroxyl, acetoxy, carbonyl, and/or a double bond typically contains variable numbers (1)(2)(3)(4) of tetrahydrofuran rings (THF). The ACGs can also be divided into subgroups based on how many and how the THF rings are arranged within each molecule. These subgroups include mono-THF, adjacent bis-THF, nonadjacent bis-THF, nonadjacent THF or non-THF, and adjacent tris-THF. Based on the accurate mass measurements, fragmentation patterns, and available literature, 77 different ACGs were tentatively identified in the extracts of A. squamosa and A. glabra using UPLC-QTOF technique.
According to the amount of THFs, the detected ACGs were divided into four subclasses: adjacent bis-THF, non-adjacent bis-THF, mono-THF with one adjacent-OH group, and mono THF with OH groups on both sides. In their mass spectra, ACGs are mainly identified as sodium adducts [M+Na] + . The removal of the lactone moieties results in the ACGs' MS 2 spectra exhibiting distinctive fragment ions (−112 amu). The hydrocarbon chain's length, the positions of the hydroxyl groups, and the positions of the THF rings may all be calculated thanks to the fragmentation around the hydroxyl groups. The amount of hydroxyl groups contained in the molecule is shown by the subsequent loss of water (−18 amu).
Squadiolin C is one of the mono-THFs with an adjacent OH group that can serve as a model. It contained a protonated molecular ion with the formula C 35 H 65 O 7 + at m/z 597.4731 (peak 61, Table S1). The fragmentation behavior of Squadiolin C showed characteristic ions due to the successive loss of four water molecules from the parent ion. This was obvious from the ions formed at m/z 579, 561, 543, and 525, indicating the presence of four OH groups in its structure. The remaining acetogenins (Table S1) were identified following the same dissociation pathways and were compared with the literature [18,22,23].
The mono-THF with two adjacent hydroxyl groups were also identified in the ethanolic extracts of the two fruits. Annonacinone (peak 79, Table S1) is an example. It exhibited a molecular ion at m/z 617.4393 [M+Na] + , corresponding to the molecular formula C 35 H 62 O 7 . Its mass spectra revealed fragment ions at m/z 577 [(M+H) + -H 2 O] + , 559 [(M+H) + -2H 2 O] + , and 541 [M+H-3H 2 O] + , which suggests the existence of three OH groups, caused by the sequential loss of three water groups.
Bullatacin is an example of the nearby bis-THF ACGS that have been found (peak 87, Table S1). It had a molecular ion with the chemical formula C 37 H 67 O 7 + at m/z 623.4863. The Na adduct is what causes the peak at m/z 645.4683. The consecutive losses of water molecules from the parent ion produced the mass fragments at m/z 587 [(M+H) + -2H 2 O] + and 569 [(M+H) + -3H 2 O] + , indicating the existence of three hydroxyl groups.
Squamostatin B (peak 30, Table S1) was one of the discovered non-adjacent bis-THF ACGS. It displayed a sodium adduct peak at m/z 661. 4632

Diterpenes
Diterpenoids called ent-kauranes are typically found in Annonaceous plants and play a supporting role in the production of gibberellin plant growth hormones. They exhibit a variety of biological properties, including anti-inflammatory, antibacterial, and  /MS fragmentation and published literature,  twenty ent-kauranes were inferred to be present in the ethanolic extracts. Their MS/MS  spectra were dominated by the diagnostic fragment ions created by the loss of water and  carbon monoxide molecules. A molecular ion at m/z 363.2499 (C 22 H 35 O 4 + ) and a base peak at m/z 289 due to the loss of the (CH 2 OCOCH 3 ) + ion were both present in Annoglabasin B (compound 70, Table S1). At m/z 271 and 234, respectively, where a water molecule and subsequently a CO group were successively lost, other distinguishing ions were also visible. Annoglabasin D (compound 17, Table S1) had a protonated molecular ion at m/z 391.2489 corresponding to the formula C 23 H 35 O 5 + . In its mass spectrum it gives the fragment ions at m/z 331, 300, and 272 due to successive loss of acetoxy, methoxy, and CO groups, respectively.
Annoglabasin E (compound 48 in Table S1) included a protonated molecular ion with the chemical formula C 20 H 33 O 3 + at m/z 321.2398. The molecular ion's removal of the CH 2 OH group is indicated by the base peak at m/z 289 that was recorded. At m/z 271 and 243, respectively, further fragment ions were discovered because of the loss of an H 2 O and a CO group.
Annoglabasin F (compound 21, Table S1) showed a molecular ion at m/z 379.2423 corresponding to the formula C 22 H 35 O 5 + and a base peak at m/z 361 due to the loss of hydroxyl group. Additionally, the peaks that appeared at m/z 302, 271, and 243 were formed from the successive loss of COOCH 3 , OCH 3 , and CO groups, respectively. Following the same fragmentation pathway, the remaining ent-kauranes were assigned and their identities were recorded in Table S1.

UV Spectroscopy
UV-vis spectroscopy is the first analysis parameter used as an indicator for the production of stable silver nanoparticles. Due to nanoparticles' selective UV absorbance, it is thought to be the quickest and easiest way to achieve accurate evaluation for nanoparticle synthesis. AgNPs exhibit a distinctive optical reflectivity that causes them to interact strongly with certain light wavelengths. Their collective electron oscillation results in a surface Plasmon resonance (SPR) absorption band that is produced by free electrons. AgNP absorption is influenced by the dielectric medium as well as chemical surroundings, particle shape, and particle size. The produced Ang-AgNPs and Ans-AgNPs were subjected to UV measurements, which revealed absorbance at 392 nm and 414 nm, respectively ( Figure 1). By using FTIR spectroscopy measurements, it was possible to estimate which functional chemical groups of the AngTE and AnsTE oversaw the reduction reaction that con-

FTIR
By using FTIR spectroscopy measurements, it was possible to estimate which functional chemical groups of the AngTE and AnsTE oversaw the reduction reaction that contributed to the creation of the silver nanoparticles. Peaks at 3480, 2920, and 1610 cm −1 indicate the functional groups in AnsTE, which are OH, C aliphatic, and C=O of phenolic acids and flavonoids. Peaks at 3470, 2920, and 1610 cm −1 represent the functional groups in AngTE. While the peaks at 1450 cm −1 and 1440 cm −1 in AngTE and AnsTE, respectively, demonstrated the presence of polyphenols and aromatic compounds, the peaks at 1050 cm −1 and 1100 cm −1 in AngTE and AnsTE, respectively, confirmed the presence of secondary OH groups ( Figure 2).

High-Resolution Transmission Electron Microscope (HR-TEM)
By employing A. glabra and A. squamosa extracts as a reducing agent, naturally occurring silver nanoparticles were produced. The spherical shape of the AgNPs was confirmed by transmission electron microscopy. In the same sequence, the two Annona species produce AgNPs with average particle sizes of 7.11 and 6.63 nm for Ang-AgNPs and Ans-AgNPs [8] (Figure 3). Selected Area Electron Diffraction Pattern (SAED) evidenced the crystalline nature of both Ang-AgNPs and Ans-AgNPs.

Zeta Potential and DLS
Nanoparticles stability, nanocharge, and surrounding media size were indicated to the Ang-AgNPs and Ans-AgNPs by Zeta potential technique. The application of the zeta potential technique was employed to evaluate the surface nano charges. The zeta potential gave a value of −17.7.0 ± 0.7 mV and −6.0 ± 0.09 mV to the nanoparticles Ang-AgNPs and Ans-AgNPs, respectively, formed through A. glabra and A. squamosa. The stabilities of the formed nanosilver particles are indicated by the negative charges of the Ang-AgNPs and Ans-AgNPs ( Figure 4).

Zeta Potential and DLS
Nanoparticles stability, nanocharge, and surrounding media size were indicated to the Ang-AgNPs and Ans-AgNPs by Zeta potential technique. The application of the zeta potential technique was employed to evaluate the surface nano charges. The zeta potential gave a value of −17.7.0 ± 0.7 mV and −6.0 ± 0.09 mV to the nanoparticles Ang-AgNPs and Ans-AgNPs, respectively, formed through A. glabra and A. squamosa. The stabilities of the formed nanosilver particles are indicated by the negative charges of the Ang-AgNPs and Ans-AgNPs ( Figure 4).   Diffraction silver nano peaks obtained for the formed AgNPs by A.glabra, were represented by peaks on the 2 Figure 4. Zeta potential measurements of the biosynthesized silver nanoparticles using the extracts, AngTE and AnsTE; (A,B): particles size distribution of Ang-AgNPs and Ans-AgNPs, respectively, (C,D): particles charges of Ang-AgNPs and Ans-AgNPs, respectively.

X-ray Diffraction (XRD)
Diffraction silver nano peaks obtained for the formed AgNPs by A.glabra, were represented by peaks on the 2 Ꝋ scale at 38.45°, 44.31°, 64.80°, and 77.81°, respectively ascribed to (111), (200), (220), and (311) planes ( Figure 5), while silver nano peaks obtained for the formed AgNPs by A. squamosa, were represented by peaks on the 2 Ꝋ scale at 38.24°, 45.06°, 64.55°, and 77.73°, respectively ascribed to (111), (200), (220), and (311) planes ( Figure 5). These results are evidence of the formation of pure silver nanoparticles in crystalline nature for both Ang-AgNPs and Ans-AgNPs. The findings are consistent with previous reports on synthesis of silver nanoparticles [24].  These results are evidence of the formation of pure silver nanoparticles in crystalline nature for both Ang-AgNPs and Ans-AgNPs. The findings are consistent with previous reports on synthesis of silver nanoparticles [24].

Scanning Electron Microscope (SEM)
SEM can precisely demonstrate the surface appearance of the biosynthesized silver nanoparticles, the particle size, shape, and distribution of the nanoparticles, in addition to determining of the micro or nanoscale formation. The SEM analysis of the biosynthesized silver nanoparticles from A. glabra and A. squamosa revealed the spherical nature and aggregate tendency of Ang-AgNPs and Ans-AgNPs ( Figure 6).

Antimicrobial Evaluation of the Extracts and Their Biogenic Nanoparticles
In the present study, the results show that the extracts possess antimicrobial activities against the tested microorganisms. The Ans-AgNPs and Ang-AgNPs showed more activ-

Scanning Electron Microscope (SEM)
SEM can precisely demonstrate the surface appearance of the biosynthesized silver nanoparticles, the particle size, shape, and distribution of the nanoparticles, in addition to determining of the micro or nanoscale formation. The SEM analysis of the biosynthesized silver nanoparticles from A. glabra and A. squamosa revealed the spherical nature and aggregate tendency of Ang-AgNPs and Ans-AgNPs ( Figure 6).

Scanning Electron Microscope (SEM)
SEM can precisely demonstrate the surface appearance of the biosynthesized silver nanoparticles, the particle size, shape, and distribution of the nanoparticles, in addition to determining of the micro or nanoscale formation. The SEM analysis of the biosynthesized silver nanoparticles from A. glabra and A. squamosa revealed the spherical nature and aggregate tendency of Ang-AgNPs and Ans-AgNPs ( Figure 6).

Antimicrobial Evaluation of the Extracts and Their Biogenic Nanoparticles
In the present study, the results show that the extracts possess antimicrobial activities

. Antimicrobial Evaluation of the Extracts and Their Biogenic Nanoparticles
In the present study, the results show that the extracts possess antimicrobial activities against the tested microorganisms. The Ans-AgNPs and Ang-AgNPs showed more activity than that observed by AngTE and AnsTE activities. The MIC of AngTE and AnsTE was shown different values ranging from 16.5 until 8.25 mg/mL in (Figure 7), however Ans-AgNPs and Ang-AgNPs had MIC values varying between 0.75 mg/mL and 0.093 mg/mL. The results indicated that the Ans-AgNPs and Ang-AgNPs have stronger activity than that of AngTE and AnsTE. The Ans-AgNPs has the highest activity with a lower MIC value (0.093 mg/mL) against C. albicans while Ang-AgNPs has recorded the highest activity against P. aeruginosa with 0.18 mg/mL (Table S2).
Pharmaceuticals 2021, 14, x FOR PEER REVIEW 9 of 24 value (0.093 mg/mL) against C. albicans while Ang-AgNPs has recorded the highest activity against P. aeruginosa with 0.18 mg/mL (Table S2). The green synthesized nanoparticles were more active materials towards the P. aeruginosa and E. coli, so we went deeper to evaluate the synergistic potential of the nanoparticles with gentamicin and vancomycin by evaluating the zone of inhibition (Tables 1-3).
In a combination assay of Ans-AgNPs and Ang-AgNPs with gentamicin, they have shown a higher synergistic effect towards P. aeruginosa and E. coli, respectively. However, the assay of Ans-AgNPs with vancomycin did not show a higher synergistic effect towards S. aureus than standard antibiotic vancomycin. The difference in the synergistic effect between the different bacterial strains might be due to the thin peptidoglycan layer in case of Gram-negative strains, which facilitates the penetration of the nanoparticle extracts (Ans-AgNPs and Ang-AgNPs).   The green synthesized nanoparticles were more active materials towards the P. aeruginosa and E. coli, so we went deeper to evaluate the synergistic potential of the nanoparticles with gentamicin and vancomycin by evaluating the zone of inhibition (Tables 1-3). Table 1. Synergistic activity of Ans-AgNPs (4B/FM) and gentamicin against the bacterial species.

Reference Zone of Inhibition
Pseudomonas aeruginosa 10 20 14 15 E. coli 10 19 14 15 In a combination assay of Ans-AgNPs and Ang-AgNPs with gentamicin, they have shown a higher synergistic effect towards P. aeruginosa and E. coli, respectively. However, the assay of Ans-AgNPs with vancomycin did not show a higher synergistic effect towards S. aureus than standard antibiotic vancomycin. The difference in the synergistic effect between the different bacterial strains might be due to the thin peptidoglycan layer in case of Gram-negative strains, which facilitates the penetration of the nanoparticle extracts (Ans-AgNPs and Ang-AgNPs). Table 2. Synergistic activity of Ans-AgNPs (4D/FM) and Gentamicin against the bacterial species.  Table 3. Synergistic activity of Ans-AgNPs (4D/FM) and vancomycin against S. aureus.

Ans-AgNPs Ans-AgNPs + Vancomycin Vancomycin
Staphylococcus aureus 15 12 17 The synergistic effect towards P. aeruginosa and E. coli were proven through the measuring of the combination's (formed silver nanoparticles+ gentamicin) inhibition zone and compared to the gentamicin and silver nanoparticles individually.

Anticancer Activity In Vitro Cytotoxic Activity
The ethanolic extracts of A. glabra (AngTE), A. squamosa (AnsTE), and their nanoparticles, Ang-AgNPs and Ans-AgNPs, respectively, were tested using the SRB assay for their in vitro cytotoxic effects against the PC3, SKOV 3, and HeLa cancer cell lines at doses ranging from 0.01 to 1000 g. All samples tested positive for cytotoxicity against cancer cells in a variety of ways. The cytotoxicity of AnsTE and Ans-AgNPs on HeLa and SKOV 3 cells was the most potent, with IC 50 s ranging from 0.001 ± 0.0001 to 1.6 ± 0.1 g/mL, and the toxicity of Ans-AgNPs on PC3 cells, with IC 50 of 1.7 ± 0.4. Additionally, the AngTE exerted a good potential killing effect against PC3, SKOV 3, and HeLa cancer cells, with IC 50 of 6.8 ± 0.8, 9.1 ± 0.6, and 5.7 ± 0.3 µg/mL, respectively, while that of AnsTE on PC3 cell line was IC 50 3.5 ± 0.1 µg/mL. Moreover, Ang-AgNPs has the strongest cytotoxic effect on PC3 and SKOV 3 with IC 50 2.4 ± 0.3 and 2.8 ± 0.23 µg/mL, respectively ( Figure 8 and Table S3).
The Effect on the Cell Cycle Distribution of Solid Tumor Cells Tracking tumor cell cycle phases was used to investigate anticancer effects. As a result, the effect of AngTE, Ang-AgNPs, AnsTE, and Ans-AgNPs on the distribution of cell cycle phases in PC3, SKOV3, and HeLa cells after 48 h of treatment in cancer cells was studied using flow cytometry as shown in (Figures 9-11), AnsTE and Ans-AgNPs significantly arrested PC3, SKOV3, and HeLa cells in the S phase with a ratio ranging from 44.1 ± 1.5% to 63.5 ± 2.3%, where the effect of Ang-AgNPs on stopping the PC3, SKOV3 cells in the S phase was 49.5 ± 2.2% and 46 ± 2.2%, respectively. In the meantime, after treatment with AngTE and Ang-AgNPs, the percentage of G1 phases in SKOV3 and HeLa cells increased by 42.8 ± 1.6% and 48.5 ± 3.1%, respectively, and treatment with AngTE affected the G1 phase by 38.7 ± 1.1% in SKOV3.
Cell Apoptosis Using Annexin V-FITC In PC3, SKOV3, and HeLa cells, annexin V-FITC/PI staining by a flow cytometer was used to distinguish between cells undergoing apoptosis (programmed cell death) and cells dying via necrosis (non-programmed cell death). There was an increase in the population of apoptotic cells in PC3 cells after treatment with AnsTE of 79.63 ± 3.3%, (Figure 12), as well as in SKOV3 after treatment with AnsTE of 83.1 ± 2.2% and Ans-AgNPs with apoptotic populations of 83.03 ± 2.04, compared to the cell control ( Figure 13). Prolonged exposure of HeLa cells to AnsTE and Ans-AgNPs for 48 h significantly induced the number of apoptotic cells by 80.93 ± 3.4% and 83.49 ± 3.7%, respectively (Figure 14), compared to the control cells. In conclusion, AnsTE used at pre-calculated IC 50 induced a significant apoptotic effect against all tumor cells (PC3, SKOV3, and HeLa).

Docking Study of Bioactive Compounds
At the catalytic site of Bcl-2, the 2D and 3D ligand interaction plots of the identified phytochemical substances are carefully examined. Figure S2 and Table S4 show the Bcl-2 target protein's binding mechanism for the eight identified substances and navitoclax (ABT-263) in PDB: 4LVT. Each simulation's gathered trajectory was statistically analyzed for the protein-ligand contacts, thus identifying the crucial amino acids (ARG143, ASN140, and ASP108) needed for efficient Bcl-2 inhibitors. According to the PDB:4LVT, the identified hit compounds cis-annonacin, motrilin, mosinone A, squamocin L, and squamone primarily form hydrogen-bonding interactions with the binding pocket residues ARG143 and ASN140 through water-bridge interactions. Squamocin G and squamocin, which were shown to be hit compounds, also form a second hydrogen bond with ASP137 and ARG104 in PDB:4LVT, in addition to the first hydrogen bonding contact.
On the other hand, Figure S3 shows the binding mechanism of the Bcl-2 target protein for the eight identified compounds and the navitoclax analogue in PDB: 4MAN. Two interfacing residues that connect hydrophobically are PHE101 and TYR199. The hydrogen bond produced with GLY142, which was also noted in the counterpart of navitoclax, was another significant interaction. In order for cis-Annonacin and squamocin to function as Bcl-2 inhibitors, they had to make essential hydrogen bonds with the GLY142 amino acid residue. The hydrogen bond formed with ARG104, which was seen in cis-annonacin, squamocin, and motrilin, was another important interaction. Both motrilin and squamocin formed a hydrogen connection with the ASP100 amino acid. When compared to other discovered compounds.
Pharmaceuticals 2021, 14, x FOR PEER REVIEW 12 of 24 Figure 9. Effect of the extracts, AngTE, AnsTE, and their nanoparticles, Ang-AgNPs, and Ans-AgNPs, respectively, on cell cycle distributions of PC3 cells. Cell cycle distribution was determined using DNA cytometry analysis after exposure to AngTE, Ang-AgNPs, AnsTE, and Ans-AgNPs for 48 h. Data are presented as the mean ± SD; n = 3; G1, G2 and S are cell cycle phases, One-way ANOVA was used to test for statistical difference (* p < 0.05). Figure 9. Effect of the extracts, AngTE, AnsTE, and their nanoparticles, Ang-AgNPs, and Ans-AgNPs, respectively, on cell cycle distributions of PC3 cells. Cell cycle distribution was determined using DNA cytometry analysis after exposure to AngTE, Ang-AgNPs, AnsTE, and Ans-AgNPs for 48 h. Data are presented as the mean ± SD; n = 3; G1, G2 and S are cell cycle phases, One-way ANOVA was used to test for statistical difference (* p < 0.05).

Cell Apoptosis Using Annexin V-FITC
In PC3, SKOV3, and HeLa cells, annexin V-FITC/PI staining by a flow cytom used to distinguish between cells undergoing apoptosis (programmed cell dea cells dying via necrosis (non-programmed cell death). There was an increase in th lation of apoptotic cells in PC3 cells after treatment with AnsTE of 79.63 ± 3.3%, 12), as well as in SKOV3 after treatment with AnsTE of 83.1 ± 2.2% and Ans-AgN apoptotic populations of 83.03 ± 2.04, compared to the cell control ( Figure 13). Pr exposure of HeLa cells to AnsTE and Ans-AgNPs for 48 h significantly induced t ber of apoptotic cells by 80.93 ± 3.4% and 83.49 ± 3.7%, respectively (Figure 14), co to the control cells. In conclusion, AnsTE used at pre-calculated IC50 induced a sig apoptotic effect against all tumor cells (PC3, SKOV3, and HeLa). Figure 11. Effect of the extracts, AngTE, AnsTE, and their nanoparticles, Ang-AgNPs, and Ans-AgNPs, respectively, on cell cycle distributions of HeLa cells. Cell cycle distribution was determined using DNA cytometry analysis after exposure to AngTE, Ang-AgNPs, AnsTE, and Ans-AgNPs for 48 h. Data are presented as the mean ± SD; n = 3; G1, G2 and S are cell cycle phases, One-way ANOVA was used to test for statistical difference (* p < 0.05).
harmaceuticals 2021, 14, x FOR PEER REVIEW 15 of 24 Figure 12. Apoptosis/necrosis assessment for the extracts, AngTE, AnsTE, and their nanoparticles Ang-AgNPs, and Ans-AgNPs, respectively, against PC3 subjected to previous treatment for 48 h and apoptosis/necrosis quantified using flow cytometry. Data are presented as the mean ± SD; n = 3 One-way ANOVA was used to test for statistical difference (* p < 0.05).

Figure 12.
Apoptosis/necrosis assessment for the extracts, AngTE, AnsTE, and their nanoparticles, Ang-AgNPs, and Ans-AgNPs, respectively, against PC3 subjected to previous treatment for 48 h, and apoptosis/necrosis quantified using flow cytometry. Data are presented as the mean ± SD; n = 3 One-way ANOVA was used to test for statistical difference (* p < 0.05). Pharmaceuticals 2021, 14, x FOR PEER REVIEW 16 of 24 Figure 13. Apoptosis/necrosis assessment of the extracts, AngTE, AnsTE, and their nanoparticles, Ang-AgNPs, and Ans-AgNPs, respectively, against SKOV3 subjected to previous treatment for 48 h, and apoptosis/necrosis quantified using flow cytometry. Data are presented as the mean ± SD; n = 3; One-way ANOVA was used to test for statistical difference (* p < 0.05). Figure 13. Apoptosis/necrosis assessment of the extracts, AngTE, AnsTE, and their nanoparticles, Ang-AgNPs, and Ans-AgNPs, respectively, against SKOV3 subjected to previous treatment for 48 h, and apoptosis/necrosis quantified using flow cytometry. Data are presented as the mean ± SD; n = 3; One-way ANOVA was used to test for statistical difference (* p < 0.05). Apoptosis/necrosis assessment of the extracts, AngTE, AnsTE, and their nanoparticles, Ang-AgNPs, and Ans-AgNPs, respectively, against HeLa subjected to previous treatment for 48 h, and apoptosis/necrosis quantified using flow cytometry Data are presented as the mean ± SD; n = 3; One-way ANOVA was used to test for statistical difference (* p < 0.05).

Docking Study of Bioactive Compounds
At the catalytic site of Bcl-2, the 2D and 3D ligand interaction plots of the identified phytochemical substances are carefully examined. Figure S2 and Table S4 show the Bcl-2 target protein's binding mechanism for the eight identified substances and navitoclax (ABT-263) in PDB: 4LVT. Each simulation's gathered trajectory was statistically analyzed for the protein-ligand contacts, thus identifying the crucial amino acids (ARG143, ASN140, and ASP108) needed for efficient Bcl-2 inhibitors. According to the PDB:4LVT, the identified hit compounds cis-annonacin, motrilin, mosinone A, squamocin L, and squamone primarily form hydrogen-bonding interactions with the binding pocket residues ARG143 and ASN140 through water-bridge interactions. Squamocin G and squamocin, which were shown to be hit compounds, also form a second hydrogen bond with ASP137 and ARG104 in PDB:4LVT, in addition to the first hydrogen bonding contact. Figure 14. Apoptosis/necrosis assessment of the extracts, AngTE, AnsTE, and their nanoparticles, Ang-AgNPs, and Ans-AgNPs, respectively, against HeLa subjected to previous treatment for 48 h, and apoptosis/necrosis quantified using flow cytometry Data are presented as the mean ± SD; n = 3; One-way ANOVA was used to test for statistical difference (* p < 0.05).

Discussion
The chemical profiles of ethanolic extracts of A. squamosa and A. glabra have been analyzed using UPLC-HR-QTOF/MS/MS to highlight the potential active constituents in the extracts. In total, 114 compounds have been identified in the positive ionization mode, including 77 acetogenins, 20 diterpenes, and 17 miscellaneous compounds.
The biosynthesis of silver nanoparticles using A. squamosa and A. glabra fruits extracts displayed high quality, biologically active nano silver particles with substantial cytotoxicities towards PC3 and SKOV3 cell lines. The cytotoxicity of AnsTE and Ans-AgNPs on HeLa and SKOV 3 cells was the most potent, with IC 50 s ranging from 0.001 ± 0.0001 to 1.6 ± 0.1 g/mL, and the toxicity of Ans-AgNPs on PC3 cells, with IC 50 of 1.7 ± 0.4. However, Ang-AgNPs demonstrated the strongest cytotoxic effect on PC3 and SKOV3 with IC 50 2.4 ± 0.3 and 2.8 ± 0.23 µg/mL. In the apoptotic study of the two Annona species and their biogenic nanoparticles. The results obtained using the flow cytometry suggest that there was an increase in the population of apoptotic cells in PC3 cells after treatment with AnsTE of 79.63 ± 3.3% as well as in SKOV3 after treatment with AnsTE of 83.1 ± 2.2% and Ans-AgNPs with apoptotic populations of 83.03 ± 2.04, compared to the cell control.
The microbiological evaluation of the formed nanoparticles indicates that Ans-AgNPs has the highest activity with a lower MIC value (0.093 mg/mL) against C. albicans while Ang-AgNPs has recorded a higher activity against P. aeuroginosa.
These processes were then followed by a molecular docking study to evaluate the role of the identified major active metabolites towards cancer cell lines, where Bcl-2 target protein's binding mechanism for the eight major identified compounds was studied. It was revealed that the mosinone a, cis-annonacin, squamocin g, motrilin, squamocin have higher binding capacity than other compounds with a binding score of (−10.87-10.56.)

Plant Materials and Extract Preparation
At the Mazhar Botanic Garden, fresh fruits of Annona glabra L. and Annona squamosa L. were harvested from trees (Giza, Egypt). The botanical specialist Mrs. Therese Labib of Orman and Qubba Gardens verified the authenticity of the voucher specimens. In a blender, fresh samples were homogenized with 70% ethanol. The extracts were dried in a rotary evaporator (Büchi, Essen, Germany under reduced pressure. For further analysis, the dried ethanolic extracts of A. glabra (AngTE) and A. squamosa (AnsTE) were kept in airtight, dark containers.

Drugs and Chemicals
The highest analytical grades of chemicals and solvents were used, and they were all purchased from Sigma-Aldrich (St. Louis, MO, USA).

UPLC-QTOF-MS/MS Analysis
The LC-MS/MS experiments were carried out at the Faculty of Pharmacy, Fayoum University, using a 6530 Q-TOF LC/MS (Agilent Technologies, Santa Clara, CA, USA) outfitted with an autosampler (G7129A), a quat. pump (G7104C), and a column comp (G7116A). The extracts were separated on an Agilent Technologies Zorbax RP-18 column (150 mm 3 mm, dp = 2.7 m). The flow rate was 0.23 mL/min, and the injection volume was 2 L. The parameters were adjusted as previously described, and mass spectra were acquired using ESI in both ionization modes. A (water with 0.1% formic acid) and B (acetonitrile with 0.1% formic acid) served as the solvents. The following gradient elution times: 0-6 min; 5-50 min; 55-75 min; isocratic elution of 50% A: 50% B; and 75-140 min, were used for the linear gradient from 50% A: 50% B to 100% B.

Green Synthesis of Ang-AgNPs and Ans-AgNPs
To reduce Ag + ions, in a cool, dark setting, 90 mL of a 1 mM AgNO 3 aqueous solution was mixed with 10 mL of each of AngTE and AnsTE in separate flasks. This mixture was then left at room temperature in complete darkness overnight. After 12 h for Ans-AgNPs and 17 h for Ang-AgNPs, an opaque-brown solution developed, indicating AgNP growth. The prepared solutions were immediately subjected to UV testing. Centrifugation at 4000 rpm for 30 min, three washings in double-distilled water, and filtering were all performed to obtain pure AgNPs. FTIR, Zeta potential, XRD, HR-TEM, and SEM were used to further analyze the pure AgNPs [25][26][27].

UV-Vis Spectroscopy
After dilution with distilled water, the UV absorbance of the bio-synthesized AgNPs from AngTE and AnsTE were examined separately using a UV-Vis spectrophotometer (Shimadzu, Kyoto, Japan).

FTIR
The various functional groups of the formed AngTE and AnsTE were measured using an FTIR spectrometer (Jasco, Tokyo, Japan) in the 4000-400 cm −1 range.

HR-TEM
HR-TEM has been used to investigate the particle structure (shape and size) along with SAED. (JEOL-JEM-1011, Kyoto, Japan) and HR-TEM at 200 kV (JEOL-JEM2100, Kyoto, Japan). Three milliliters of the suspended nanoparticles from each sample of Ans-AgNPs and Ang-AgNPs were separately deposited on the copper grid for HR-TEM analysis and were then left to dry at room temperature for 15 min.

Zeta Potential Analysis
A zeta sizer nano ZN was used to investigate zeta potentials and the homogeneity of AngTE and AnsTE (Malvern Panalytical Ltd., Malvern, UK). Prior to the measurements, an aliquot of each Ans-AgNP and Ang-AgNP was combined with ultra-pure water and subjected to sonication for 15 min.

XRD
An XPERT-PRO-PAN-alytical Powder Diffractometer (PAN-alytical B.V., Almelo, The Netherlands) and a monochromatic radiation source Cu-K radiation (θ = 1.5406) at 45 kV and 30 mA at room temperature were used to evaluate metal nanoparticles using XRD analysis as a surface chemical analysis technique. The calculated silver nano-powder intensity data covered a 2 θ range of 4.01 • -79.99 • .

SEM
The morphology of the biosynthesized AgNPs was investigated using SEM (TM 1000, Hitachi, Chiyoda, Japan).

Antibacterial Activity Determination of Antibacterial Activity of the Extracts
The antimicrobial effects were tested against Gram-positive bacteria (Staphylococcus aureus), Gram-negative bacteria (Escherichia coli and Pseudomonas aeruginosa), and the pathogenic fungus (Candida albicans).

Minimal Inhibitory Concentration (MIC) Determination
Using sterile 2 mL 96-well plates, the antimicrobial activities were assessed [28]. Then, 100 µL of double-strength Mueller Hinton agar that had been sterilized was placed in each well. To create a concentration sequence from 33 mg/mL to 0.128 mg/mL in (AngTE and AnsTE) and from 12 till 0.02 mg/mL in (Ang-AgNPs and Ans-AgNPs), an additional 50 µL of a mixture of culture medium and plant extract (dissolved in 10% DMSO) was serially diluted into the wells 1-9. Well 10 served as growth control, while well 11 was the DMSO negative control. The MIC value (the lowest concentration that completely inhibited the bacterial growth) was determined for each sample after 24 h of incubation at 37 • C.

Disc Diffusion Assay
A modified Kirby-Bauer disc diffusion method was used to measure antimicrobial activity. Then, 100 µL of bacteria were removed from 10 mL of fresh culture media and the bacterial inoculums were standardized with respect to McFarland standard no. 0.5 with 1.5 × 108 CFU/ml bacterial density counted until they reached a density of about 108 cells/mL [29]. Onto Mueller Hinton (MH) agar, 100 µL of the microbial suspension was applied. An aliquot (10 µL) of each of the Ans-AgNPs and Ang-AgNPs, which were dissolved in 10% DMSO, was pipetted onto an agar surface using a sterile paper disc (Whatman No. 1, 5.5 mm paper disc). A filter disc coated with 10 µL of DMSO (Sigma, St. Louis, MO, USA) served as a negative control, and an antibiotic disc (gentamicin for Gram negative bacteria and vancomycin for Gram positive bacteria) was coated with Ans-AgNPs and Ang-AgNPs. The plates were inverted and incubated at 37 ºC for 18 h. Inhibition zone was determined by measuring the diameter of the clear zone (inhibition of growth) around each disc and recorded in millimeter according to CLSI [30].

Cytotoxicity Assay
Using the sulphorhodamine B (SRB) assay, the cytotoxicity of the AngTE, AnsTE, Ang-AgNPs, and Ans-AgNPs against (HeLa, PC-3, and SKOV3) human tumor cells was assessed. Before being treated with the AngTE, AnsTE, Ang-AgNPs, and Ans-AgNPs, cells that were growing at 80% confluency, trypsinized and cultured in a 96-well tissue culture plate for 24 h. Cells were subjected to six different doses of each chemical (0.01, 0.1, 1, 10, and 1000 g/mL) with untreated cells added as a control. Before the cells were fixed with TCA (10% w/v) for an hour at 4 • C, they were exposed to the concentrations for 72 h. After multiple washings, cells were stained with a 0.4% (w/v) SRB solution for 10 min in the dark. The surplus stain was eliminated using 1% (v/v) glacial acetic acid. The SRB-stained cells were dissolved in Tris-HCl buffer after drying overnight. A microplate reader was used to gauge the color intensity at 540 nm. Sigma Plot 12.0 software was used to examine the association between each tumor cell line's viability percentage and compound concentrations in order to determine the IC 50 (drug dose that reduces survival to 50%) [31].

Cell Cycle Analysis
The IC 50 values for AngTE, AnsTE, Ang-AgNPs, and Ans-AgNPs were pre-calculated and administered to (HeLa, PC-3, and SKOV3) cells for 48 h. The cells were then fixed in ice-cold 60% ethanol at 40 C and trypsinized before being washed twice in phosphatebuffered saline. After being resuspended, the cells were incubated for 15 min in 500 L of Cell Signaling Technology's (CST) propidium iodide with RNase staining buffer. In order to evaluate the data from 10,000 cells and the distribution of cell cycle phases for each sample,

Protein Preparation
The RCSB Protein Data Bank was used to retrieve the Bcl-2 proteins' structures, which have resolutions of 2.04 and 2.07 (PDB IDs: 4LVT and 4MAN). The system was optimized, the partial charges were computed, and the enzymes were 3D protonated ( Figure S2) [33].

Validation of the Molecular Docking Method
The docking methods and parameters were validated to make sure that the ligand orientations derived from the molecular docking simulations are accurate and that the potential binding modes of ligands are reasonable. Using the (MOE, 2019, Montreal, QC, Canada) software (https://www.chemcomp.com/Products.htm), re-docking for cocrystallized ligands in a target protein's active site after they have been eliminated from Bcl-2 receptor proteins was observed, and a comparison of the results of the final docking (binding mode, hydrogen bond interactions and docking score) was made. The X-ray bioactive conformer of navitoclax and navitoclax analogue was then aligned with the best fitted pose discovered during the docking run [34]. With RMSD = 1.35 and 1.70 for PDB:4LVT and PDB:4MAN (https://www.rcsb.org/structure/4MAN), respectively, the alignment demonstrated good agreement between them, demonstrating the capability of the used docking protocol to retrieve legitimate docking poses. If the returned RMSD value is 3, the method is deemed successful. (Figure S3).

Ligand Preparation
Test ligand structures of eight phytochemical compounds derived from A. squamosa and A. glabra extracts, navitoclax (ABT-263) and navitoclax analog (with Indole) are prepared for docking by applying the following steps: (1) ChemBioDraw Ultra 12.0 (CambridgeSoft Corporation, Cambridge, UK) as used to create the target molecules, which were then copied to MOE. (2) The compounds were protonated in three dimensions. (3) The structures were energy minimized to a gradient of 0.5 by the (MMff94x) Merk Molecular Force Field. (4) The force field partial charges for each molecule were determined. (5) Each molecule underwent a stochastic conformational analysis using default parameters, with the results being saved into a different conformational database. (6) To dock into the active site of the androgen receptor, the most stable conformers of each molecule were recorded in a different database.

Molecular Docking Calculations
The most stable conformers were docked using MOE-dock using flexible ligand-rigid receptors as follows: (1) The receptor atoms are defined as (receptor + solvent). (2) The placement site was defined as (ligand atoms). (3) Ligand atoms were retrieved from the conformational database of the target compounds' least energetic conformers. (4) The placement method was modified to (Alpha triangle). (5) The London dG scoring function was used to estimate the ligand's binding free energy from a given pose. Finally, the best scoring poses of the docked compounds were identified, and the complexes' protein-ligand interactions were investigated.

Conclusions
To the best of our knowledge, this is the first comprehensive study to compare the metabolic profiles of A. squamosa and A. glabra fruits using UPLC-HR-QTOF/MS/MS and correlate them with their cytotoxic and antimicrobial activities using molecular docking simulations. Our results revealed that the silver nanoparticles of both species formed a good synergistic effect in combination with gentamicin against P. aeruginosa and E. coli, respectively. Moreover, the formed nanoparticles potentiated the cytotoxicity of the total extracts against PC3 and SKOV3 cell lines, so the formed nanoparticles are excellent candidates for the treatment of skin infections in patients with ovarian or prostatic cancers.  Figure S3. The alignment between the X-ray bioactive conformer of the navitoclax and its analog (colored in yellow) and the redocked pose of the same compound (colored in cyan) at Bcl-2 receptor binding site proteins. Table S1: Compounds identified in the extracts of A. glabra (AngTE) and A. squamosa (AnsTE) fruits, by UPLC-QTOF-MS/MS (positive mode). Table S2. Results of the minimum inhibitory concentration (MIC) of the extracts. Table S3. The IC 50 (µg) of the ethanolic extracts of A. squamosal and A. glabra, and their nanoparticles against different human solid tumor cells. Table S4. Molecular docking results of major identified compounds from the extracts of A. squamosal and A. glabra.