Propranolol-Loaded Trehalosome as Antiproliferative Agent for Treating Skin Cancer: Optimization, Cytotoxicity, and In Silico Studies

This study aims at preparing propranolol-loaded trehalosomes (a trehalose-coated liposome) to be used as an antiproliferative agent for treating skin cancer. A factorial design was used to select the optimum formula, where trehalose, lecithin, and Tween 80 levels were studied. A total of 24 runs were prepared and characterized according to size, charge, entrapment efficiency, and release after 3 h to select the optimum formula. The optimized formula was investigated using TEM, DSC, and FTIR. Cell studies were carried out against the human melanoma cell line to measure cytotoxicity, apoptosis/necrosis, and cell cycle arrest. In silico studies were conducted to understand the interaction between propranolol and the influential receptors in melanoma. The results showed the selected formula consisted of trehalose (175 mg), lecithin (164 mg), and Tween 80 (200 mg) with a size of 245 nm, a charge of −9 mV, an EE% of 68%, and a Q3 of 62%. Moreover, the selected formula has good cytotoxicity compared to the free drug due to the synergistic effect of the drug and the designed carrier. IC50 of free propranolol and the encapsulation of propranolol were 17.48 μg/mL and 7.26 μg/mL, respectively. Also, propranolol and the encapsulation of propranolol were found to significantly increase early and late apoptosis, in addition to inducing G1 phase cell cycle arrest. An in silico virtual study demonstrated that the highest influential receptors in melanoma were the vitamin D receptor, CRH-R1, VEGFR 1, and c-Kit, which matches the results of experimental apoptotic and cell cycle analysis. In conclusion, the selected formula has good cytotoxicity compared to the free drug due to the synergistic effect of the drug and the designed carrier, which make it a good candidate as an antiproliferative agent for treating skin cancer.


Introduction
One of the biggest challenges in cancer research is understanding how cancer immunetolerance mechanisms work. Melanoma is a highly aggressive tumor with low five-year survival rates. It has been discovered that β2-adrenoceptors play a role in cancer and are a new target in treating melanoma. Recent studies have shown that β3-adrenoceptors have a pleiotropic effect on the melanoma microenvironment, leading to cancer progression, but the mechanisms are still not well-understood. Beta-blockers are widely used therapeutic agents and have shown potential anticarcinogenic effects in patients with prostate cancer, breast cancer, and melanoma. However, it should also be noted that β-blockers can cause or aggravate various skin conditions such as psoriasis, lichen planus-like drug eruptions (LDE), acrocyanosis, and alopecia [1].
Propranolol hydrochloride was purchased from Sigma-Aldrich (Saint Louis, MO, USA). Trehalose was purchased from Oxford Lab fine chem LLP (Mumbai, India). Lecithin from soyabean and Tween 80 were purchased from Alpha Chemika Company (Mumbai, India). Potassium dihydrogen phosphate and disodium hydrogen phosphate were supplied from El Nasr Pharmaceutical Chemicals (Cairo, Egypt). DMEM media, streptomycin, penicillin, heat-inactivated fetal bovine serum (FBS), and DNAse/RNAse-free water were purchased from Thermo Fisher Scientific (Cleveland, OH, USA).

Preparation of PPL-Loaded Trehalosomes (PTHs)
PPL (10 mg), lecithin, and Tween 80 were dissolved in ethanol using a magnetic stirrer at 60 • C. Trehalose was dissolved in distilled water. A dropwise addition of the organic phase to the aqueous phase was conducted with the aid of a stirrer for 30 min at 60 • C, then left at room temperature to cool. PTHs dispersions were sonicated for 2 min (5 s on/off; Ultrasonic LC 60 H, Elma, Singem, Germany). The prepared PTHs were kept in a refrigerator for further investigation [13].

Experimental Design
An experimental design was conducted with the aid of Design-Expert ® Software version 13 (Stat-Ease, Inc., Minneapolis, MN, USA) using 2 3 factorial design. In the study, the independent variables were trehalose (X1), lecithin (X2), and Tween 80 (X3). PTHs were prepared formulations as demonstrated in Table 1.
Each formulation was prepared 3 times according to the levels demonstrated in Table 2 to calculate the mean and standard deviation. The prepared formulations were characterized, where the dependent variables were particle size, charge, entrapment efficiency percent, and percentage of drug released after 3 h, as shown in Table 3. The optimized formula was selected according to the lowest size and the highest charge (absolute value), EE %, and 3 h dissolution, as shown in Table 3.  The particle size and surface charge of the prepared PTHs were determined using Malvern Nano-ZS (Malvern Instruments Ltd., Worcestershire, UK) after appropriate dilution with distilled water at a temperature of 25 • C [14].

Entrapment Efficiency (EE)
Entrapment efficiency (EE) was determined by indirect technique, where PTHs were separated from the dispersion by an ultracentrifugation method at 15,000 rpm for 1 h at 4 • C. The collected pellets were washed 3 times and added to the supernatant. Finally, the PPL concentration was measured spectrophotometrically (Shimadzu UV-1600, Kyoto, Japan) in the supernatant at a wavelength of 290 nm. The entrapment efficiency was calculated according to the following equation: EE% = (Total PPL − Unentraped PPL)/(Total PPL) × 100 (1)

In Vitro Release Studies
The dissolution profile of PPL was determined using the dissolution apparatus (Dr. Schleuniger Pharmaton, type Diss 6000, Solothurn, Switzerland). An amount of 1 mL of PTH was added to a dialysis tube (12 kDa, Sigma-Aldrich, USA), attached to a dissolution paddle, and then placed into a dissolution vessel containing 100 mL phosphate buffer. The solubility was tested prior to the test to confirm the maintenance of the sink condition. The release condition was stirred at 100 rpm and the temperature was maintained at 37 • C. Samples were collected each hour and replaced by fresh buffer. The concentration of the released PPL was measured by spectrophotometer at 290 nm and the cumulative release profile was determined.

Transmission Electron Microscopy (TEM)
The morphological examination of the optimum PTHs was studied by transmission electron microscopy (Joel JEM 1230, Tokyo, Japan). The PTH sample was dried on a carbon-coated grid and stained with phosphotungstic acid (2%). After the drying of the phosphotungstic acid, the images were acquired using TEM at 200 kV.

Differential Scanning Calorimetry (DSC)
DSC thermograms of PPL, lecithin, trehalose, Tween 80, and the selected formula were generated by differential scanning calorimetry equipment (DSC50, Shimadzu, Japan) connected to TA-501 thermal analyzer. A total of 5 mg of each sample was heated in a sealed aluminum pan at a temperature range of 25 to 300 • C at a heating rate of 10 • C/min, under a nitrogen flow of 20 mL/min.

Fourier-Transform Infrared Spectroscopy (FTIR)
FTIR spectra of PPL, lecithin, trehalose, Tween 80, and the selected formula were obtained by FTIR spectrophotometer (FTIR-8400, Shimadzu, Japan). Then, 3 mg of each sample was mixed with KBr and compressed into a disk, followed by scanning at a range from 400 to 4000 cm −1 .

Cytotoxicity Cell Study
The cell viability was measured using Sulforhodamine B (SRB) assay. Human melanoma cells A375 (ATCC ® CRL-1619™) were used to test unloaded and loaded formula, PPL, and cisplatin as control. Aliquots of 100 µL cell suspension (5 × 103 cells) were incubated in complete media for 24 h. The cells were treated with 100 µL media containing different samples at various concentrations (five concentrations prepad in cells media were used). After 72 h of drug exposure, the cells were fixed by replacing media with 150 µL of 10% Trichloroacetic acid (TCA) and incubated at 4 • C for 1 h. The TCA was removed, and the cells were washed 5 times with distilled water. An amount of 70 µL of 0.4% w/v SRB solution was mixed and incubated in a dark place at room temperature for 10 min, followed by washing 3 times with 1% acetic acid and drying overnight. Then, 150 µL of 10 mM of Tris base solution (TRIS) was added to dissolve the protein-bound SRB stain. The absorbance was measured at 540 nm using a BMG LABTECH ® -FLUOstar Omega microplate reader (Ortenberg, Germany) [15]. The cytotoxicity was expressed as a Half-maximal inhibitory concentration (IC50). The IC50 was determined by plotting dose-response curves, where IC50 values were obtained from dose-response curves and regression analysis by GraphPad PRISM version 6 program (GraphPad Software, Inc., San Diego, CA, USA).

Flow Cytometric Analysis
Apoptosis/necrosis and cell cycle distribution were investigated by flow cytometric analysis (ACEA Biosciences Inc., San Diego, CA, USA) as reported [15]. Human melanoma cells (A375) were maintained in DMEM media, then 100 mg/mL of streptomycin, 100 units/mL of penicillin, and 10% of heat-inactivated fetal bovine serum in a humidified, 5% (v/v) carbon dioxide atmosphere were added at 37 • C.

Annexin-V-FITC Apoptosis Assay
Apoptotic cells were determined using an Annexin V-FITC apoptosis detection kit (Abcam Inc., Cambridge Science Park, Cambridge, UK). Cells were treated with different samples for 48 h. According to the manufacturer protocol, 105 cells were collected by trypsinization, washed 2 times with ice-cold PBS (pH 7.4), and incubated in the dark with 0.5 mL of Annexin V-FITC/PI solution for 30 min at room temperature. After staining, the cells were injected via a flow cytometer and analyzed for fluorescent signals FITC (FL1, λex/em 488/530 nm) and PI (FL2, λex/em 535/617 nm). For each sample, 12,000 events were acquired and positive FITC and/or PI cells were quantified by quadrant analysis and calculated using ACEA NovoExpress™ software version 1.3 (ACEA Biosciences Inc., San Diego, CA, USA) [15].

Cell Cycle Analysis
The cell distribution was analyzed by cell treatment with test compounds for 48 h using cisplatin as a positive control. The cells (105 cells) were collected by trypsinization and washed twice with ice-cold PBS (pH 7.4). The cells were resuspended in 2 mL of 60% ice-cold ethanol and incubated at 4 • C for 1 h for fixation. The fixed cells were washed twice with PBS and re-suspended in 1 mL of PBS containing 50 µg/mL RNAase A and 10 µg/mL propidium iodide (PI). After 20 min of incubation in the dark at 37 • C, the cells were analyzed for DNA contents using flow cytometry analysis using an FL2 (λex/em 535/617 nm) signal detector (ACEA Novocyte™ flowcytometer, ACEA Biosciences Inc., San Diego, CA, USA). For each sample, 12,000 events were acquired. The cell cycle distribution was calculated using ACEA NovoExpress™ software (ACEA Biosciences Inc., San Diego, CA, USA) [16].
2.6. In Silico Studies 2.6.1. Ligand File Preparation Propranolol, a beta-blocker drug, was downloaded as sdf to be in its minimized energy form from the Pubchem database (https://pubchem.ncbi.nlm.nih.gov/ (accessed on 20 January 2023)). This step was carried out to prepare the drug file in the required form to be used in the docking studies with different enzymes.

Preparation of Enzyme Files
Different melanocyte receptors involved in melanoma were investigated and then downloaded from either the Swiss model (https://swissmodel.expasy.org/ (accessed on 21 January 2023)) [17,18] or the protein database bank (https://www.rcsb.org/ (accessed on 20 January 2023)) in the form of pdb files. Then, the files were refined by removing the corresponding co-crystal ligand from each enzyme file to obtain the enzyme-only files in pdb format using PyMol software version 2.2.3.

Interaction of Propranolol and Different Melanocyte Receptors
The interaction between each receptor individually and propranolol was visualized using Discovery Studio-19 software version 19.1.0.18287. The studied interactions included conventional H-bonding, aromatic ring center, hydrophobic interaction, covalent bond, perpendicular and parallel pi-stacking, charge center, etc. The receptor-propranolol binding interactions were assessed through calculations of binding free energy (DG).

Statistics
All experiments were conducted in triplicates and all data were reported as the mean ± SD. The GraphPad Prism 6 software (San Diego, CA, USA) was used to analyze the statistically significant difference using one-way analysis of variance (ANOVA) and Tukey post-test (* p < 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001).  Table 3, the formula was significantly affected by the levels of different factors and the interaction between them. The design used shows a promising fit to the design space, as all adequate precision values exceed four. Additionally, the adjusted R 2 and predicted R 2 are similar, which displays the consistency of the model. The impact of the factors on the responses was analyzed, and the findings are presented in 3D diagrams as shown in Figure 1. conventional H-bonding, aromatic ring center, hydrophobic interaction, covalent bond, perpendicular and parallel pi-stacking, charge center, etc. The receptor-propranolol binding interactions were assessed through calculations of binding free energy (DG).

Statistics
All experiments were conducted in triplicates and all data were reported as the mean ± SD. The GraphPad Prism 6 software (San Diego, CA, USA) was used to analyze the statistically significant difference using one-way analysis of variance (ANOVA) and Tukey posttest (* p < 0.05, ** p < 0.01, *** p < 0.001 and **** p < 0.0001).

Analysis of Responses
Different formulations were prepared according to the composition demonstrated in Table 1. In total, 23 factorial designs were created by Design Expert software. The experimental runs were 24 in order to consider the effect of three factors at different levels. As presented in Table 3, the formula was significantly affected by the levels of different factors and the interaction between them. The design used shows a promising fit to the design space, as all adequate precision values exceed four. Additionally, the adjusted R 2 and predicted R 2 are similar, which displays the consistency of the model. The impact of the factors on the responses was analyzed, and the findings are presented in 3D diagrams as shown in Figure 1. The color scale started from blue as a low value and increased gradually through green till reached red color for high value.

Particle Size
The prepared PTH formulations were evaluated via different techniques in which particle size was considered as one of the critical factors for the optimization of PTHs. Smaller particle size usually provides a good surface area to encapsulate drugs such as propranolol [21]. The data show that the obtained PTH sizes ranged from 65 up to 386.3 nm (Figure 1ac). Three factors have a significant effect on particle size. Increasing trehalose and lecithin led to an increase in particle size, while increasing Tween 80 led to a decrease in the particle size, as shown in Table 4, where the p values were <0.0001. Moreover, the interaction of the second and third factors caused a significant increase in particle size (p values < 0.0001), as shown in Figure 1a-c. It was reported by Ochi et al. that lecithin hydrophobicity enlarges the vesicle dimension [22]. The incorporation of Tween 80 in PTHs causes a decrease in particle size due to steric repulsion rendered by surfactant molecules, which prevents or minimizes the aggregation of the vesicles [23].

Particle Size
The prepared PTH formulations were evaluated via different techniques in which particle size was considered as one of the critical factors for the optimization of PTHs. Smaller particle size usually provides a good surface area to encapsulate drugs such as propranolol [21]. The data show that the obtained PTH sizes ranged from 65 up to 386.3 nm (Figure 1a-c). Three factors have a significant effect on particle size. Increasing trehalose and lecithin led to an increase in particle size, while increasing Tween 80 led to a decrease in the particle size, as shown in Table 4, where the p values were <0.0001. Moreover, the interaction of the second and third factors caused a significant increase in particle size (p values < 0.0001), as shown in Figure 1a-c. It was reported by Ochi et al. that lecithin hydrophobicity enlarges the vesicle dimension [22]. The incorporation of Tween 80 in PTHs causes a decrease in particle size due to steric repulsion rendered by surfactant molecules, which prevents or minimizes the aggregation of the vesicles [23].

Surface Charge
The zeta potential is another parameter that affects the stability of nanocarriers, where high negativity or positivity of the surface charge can form a stable colloidal system [21,24]. It was found that the surface charge was not influenced by the three factors alone, but the interaction between trehalose and lecithin beside trehalose and Tween 80 were significant where p-values were 0.0121 and 0.0242, respectively, as shown in Table 4. The surface charges of the PTH formulation values ranged from −5.9 mV up to −12.8 mV, as shown in Figure 1d-f. Their surface's negative charge indicates that most of the samples are physically stable where the negative charge will prevent aggregation. Furthermore, the presence of Tween 80 could prevent the agglomeration of nanoparticles through the steric hindrance effect [25]. The combination of small particle size and low negative surface charge enabled high dermal penetration by trehalosomes, which is appropriate for treating skin cancer [26].

Entrapment Efficiency
Entrapment efficiency is usually measured to determine the encapsulation ability of the nanocarrier. The EE percentage ranged from 57% up to 92.5%, as shown in Figure 1g-i. It has been found that increasing trehalose and lecithin have a significant influence on the EE% results, as shown in Table 4, with p-values of 0.0013 and <0.0001, respectively. Furthermore, the interaction between trehalose and the two other factors negatively affects the EE%. According to Table 4, increasing the Tween 80 level showed an insignificant effect on EE%. Therefore, the formulations with a low level of Tween 80 showed high EE%.

Dissolution after 3 h
The release profile of PPL was also conducted to study the ability of the nanocarrier system to control drug release. The three factors significantly affected the Q3, as shown in Figure 1j-l and Table 4. The release after 3 h ranged from 12% to 100%, depending on the level of the three factors. As displayed in Table 4, trehalose has a significant (p < 0.0001) positive effect on the drug released. This could be attributed to the high solubility of trehalose, which could provide channels in the bilayer structure of trehalosomes, allowing the drug release. On the other hand, lecithin and Tween 80 had a significant negative effect (p < 0.0001) on the drug released.

Optimization
Numerical and graphical analyses were conducted to determine the best formula for achieving the minimum size and maximum surface charge, EE%, and Q3 ( Table 3). The results of both optimizations were similar, indicating the reliability of the chosen model ( Figure 2). The best optimum formula has a desirability of 0.468. This optimum formula consisted of trehalose (175 mg), lecithin (164 mg), and Tween 80 (200 mg), as shown in Table 5. According to program analysis, this formula should have a size of 245 nm, a charge of −9 mV, an EE% of 68%, and a Q3 of 62%. This formula was then prepared and characterized to confirm the accuracy of the predictions through experimental methods as shown in Table 5.

Characterization of Optimum Formula
A TEM micrograph of the optimum formula showed spherical non-aggregated vesicles with particle sizes ranging from 280 to 310 nm, as shown in Figure 3a. In addition, the results indicate that the particle size of the developed vesicles is consistent with the size predicted by the design. To confirm the drug encapsulation, DSC thermograms were measured, as shown in Figure 3b. The propranolol thermogram showed an endothermic peak at 163.8 °C, corresponding to its melting point. This finding was in agreement with the previously reported melting point of propranolol [27]. The melting point of trehalose was around 100, which was in agreement with the previously reported melting point [28]. The melting point of the prepared formula showed the disappearance of the drug melting point, which confirms the drug diffusion in the prepared vesicles [29,30]. Furthermore, Figure 3c demonstrated the FTIR spectrum of the prepared formula. The drug showed a characteristic absorption peak at 2965 cm −1 due to the presence of a secondary amine group and 3283 cm −1 due to the hydroxyl group (secondary); the aryl alkyl ether displayed a stretching band at 1267.27 cm −1 and a peak at 798 cm −1 due to a-substituted naphthalene [31]. The spectrum of the prepared formula showed the characteristic absorption peaks of the drug and vesicle components, which emphasize the success of the trehalosome assembly and drug encapsulation.

Characterization of Optimum Formula
A TEM micrograph of the optimum formula showed spherical non-aggregated vesicles with particle sizes ranging from 280 to 310 nm, as shown in Figure 3a. In addition, the results indicate that the particle size of the developed vesicles is consistent with the size predicted by the design. To confirm the drug encapsulation, DSC thermograms were measured, as shown in Figure 3b. The propranolol thermogram showed an endothermic peak at 163.8 • C, corresponding to its melting point. This finding was in agreement with the previously reported melting point of propranolol [27]. The melting point of trehalose was around 100, which was in agreement with the previously reported melting point [28]. The melting point of the prepared formula showed the disappearance of the drug melting point, which confirms the drug diffusion in the prepared vesicles [29,30]. Furthermore, Figure 3c demonstrated the FTIR spectrum of the prepared formula. The drug showed a characteristic absorption peak at 2965 cm −1 due to the presence of a secondary amine group and 3283 cm −1 due to the hydroxyl group (secondary); the aryl alkyl ether displayed a stretching band at 1267.27 cm −1 and a peak at 798 cm −1 due to a-substituted naphthalene [31]. The spectrum of the prepared formula showed the characteristic absorption peaks of the drug and vesicle components, which emphasize the success of the trehalosome assembly and drug encapsulation.

Characterization of Optimum Formula
A TEM micrograph of the optimum formula showed spherical non-aggregated vesicles with particle sizes ranging from 280 to 310 nm, as shown in Figure 3a. In addition, the results indicate that the particle size of the developed vesicles is consistent with the size predicted by the design. To confirm the drug encapsulation, DSC thermograms were measured, as shown in Figure 3b. The propranolol thermogram showed an endothermic peak at 163.8 °C, corresponding to its melting point. This finding was in agreement with the previously reported melting point of propranolol [27]. The melting point of trehalose was around 100, which was in agreement with the previously reported melting point [28]. The melting point of the prepared formula showed the disappearance of the drug melting point, which confirms the drug diffusion in the prepared vesicles [29,30]. Furthermore, Figure 3c demonstrated the FTIR spectrum of the prepared formula. The drug showed a characteristic absorption peak at 2965 cm −1 due to the presence of a secondary amine group and 3283 cm −1 due to the hydroxyl group (secondary); the aryl alkyl ether displayed a stretching band at 1267.27 cm −1 and a peak at 798 cm −1 due to a-substituted naphthalene [31]. The spectrum of the prepared formula showed the characteristic absorption peaks of the drug and vesicle components, which emphasize the success of the trehalosome assembly and drug encapsulation.

Cytotoxicity Cell Study
The study evaluated the cytotoxicity of free propranolol, blank and loaded formula, and cisplatin against an A 375 human melanoma cell line using the SRB assay. After 72 h of treatment with different concentrations (0.1-1000 µg/mL) of the four test groups, the cell viability was measured as shown in Figure 4a. The results showed that the cell viability of the A 375 human melanoma cell line decreased in a dose-dependent manner in response to treatment by either free propranolol, blank and loaded formula, or cisplatin. Treatment with cisplatin resulted in a greater decrease in the cell viability percentage, with IC50 0.59 µg/mL compared to treatment with other groups, which could be attributed to the potency of cisplatin as a standard drug. However, cisplatin still has serious side effects, which motivates researchers to find alternative effective drugs with lower side effects [32]. On the other hand, propranolol as a beta-blocker drug was previously reported as a promising anticancer drug. Beta-blockers, particularly propranolol, have been found to significantly improve the outlook of melanoma patients. Propranolol works by reducing the growth of blood vessels in the tumor, as well as the proliferation, invasiveness, and immune suppression of tumor cells. Studies have shown that only β3-adrenoceptors, rather than β2-adrenoceptors, were increased in peripheral blood mononuclear cells under hypoxia. These receptors were selectively found in immune cell sub-populations, including Treg, MDSC, and NK cells, which increased the number and cytotoxicity of NK and CD8 cells. Catecholamines may slow down the progression of melanoma, and β-blockers may have untapped potential as a treatment for melanoma, both as a preventive measure and as an adjuvant to other targeted and immune therapies [1,33]. In the current study, the IC50 of free propranolol was 17.48 µg/mL, as shown in Figure 4b. The encapsulation of propranolol showed a significant decrease (p < 0.0001) in IC50 to reach 7.26 µg/mL. This confirmed that the required dose of propranolol could be decreased almost to half when being incorporated within trehalosomes instead of being administered as a free drug. The blank formula achieved IC50 136.7 µg/mL, which confirmed that the synergistic effects of both propranolol and trehalosmoes led to a decrease in the IC50. To understand the cytotoxicity mechanism of free and encapsulated propranolol, further investigations were conducted, as shown in the next section.

Cytotoxicity Cell Study
The study evaluated the cytotoxicity of free propranolol, blank and loaded formula, and cisplatin against an A 375 human melanoma cell line using the SRB assay. After 72 h of treatment with different concentrations (0.1-1000 µg/mL) of the four test groups, the cell viability was measured as shown in Figure 4a. The results showed that the cell viability of the A 375 human melanoma cell line decreased in a dose-dependent manner in response to treatment by either free propranolol, blank and loaded formula, or cisplatin. Treatment with cisplatin resulted in a greater decrease in the cell viability percentage, with IC50 0.59 µg/mL compared to treatment with other groups, which could be attributed to the potency of cisplatin as a standard drug. However, cisplatin still has serious side effects, which motivates researchers to find alternative effective drugs with lower side effects [32]. On the other hand, propranolol as a beta-blocker drug was previously reported as a promising anticancer drug. Beta-blockers, particularly propranolol, have been found to significantly improve the outlook of melanoma patients. Propranolol works by reducing the growth of blood vessels in the tumor, as well as the proliferation, invasiveness, and immune suppression of tumor cells. Studies have shown that only β3-adrenoceptors, rather than β2-adrenoceptors, were increased in peripheral blood mononuclear cells under hypoxia. These receptors were selectively found in immune cell sub-populations, including Treg, MDSC, and NK cells, which increased the number and cytotoxicity of NK and CD8 cells. Catecholamines may slow down the progression of melanoma, and β-blockers may have untapped potential as a treatment for melanoma, both as a preventive measure and as an adjuvant to other targeted and immune therapies [1,33]. In the current study, the IC50 of free propranolol was 17.48 µg/mL, as shown in Figure 4b. The encapsulation of propranolol showed a significant decrease (p < 0.0001) in IC50 to reach 7.26 µg/mL. This confirmed that the required dose of propranolol could be decreased almost to half when being incorporated within trehalosomes instead of being administered as a free drug. The blank formula achieved IC50 136.7 µg/mL, which confirmed that the synergistic effects of both propranolol and trehalosmoes led to a decrease in the IC50. To understand the cytotoxicity mechanism of free and encapsulated propranolol, further investigations were conducted, as shown in the next section.

Annexin-V-FITC Apoptosis Assay
To determine how different treatments inhibit the proliferation of human melanoma cell lines, an Annexin V/propidium iodide assay kit was utilized. The assay measured early apoptosis (Q4), late apoptosis (Q2), necrosis (Q1), and cell viability (Q3) after a 48 h incubation period, as shown in Figure 5a-c. The experiment included untreated human melanoma cell lines as a control, as well as cells treated with free and encapsulated propranolol. In

Annexin-V-FITC Apoptosis Assay
To determine how different treatments inhibit the proliferation of human melanoma cell lines, an Annexin V/propidium iodide assay kit was utilized. The assay measured early apoptosis (Q4), late apoptosis (Q2), necrosis (Q1), and cell viability (Q3) after a 48 h incubation period, as shown in Figure 5a-c. The experiment included untreated human melanoma cell lines as a control, as well as cells treated with free and encapsulated propranolol. In Figure 5d, the two treatments were found to significantly increase early and late apoptosis compared to the control, which is consistent with previous reports of the antiproliferative activity of propranolol [34]. The cytotoxicity of free propranolol was mainly due to the additive effect of early and late apoptosis, while necrosis played a minor role. Similarly, propranolol-loaded trehalosomes showed cytotoxicity due to the additive effect of early and late apoptosis, with a lower contribution from necrosis. Interestingly, apoptosis was similar in the case of free and encapsulated propranolol, indicating the trehalosomes' ability to induce programmed death in cancer cells [15]. To learn more about the nature of the cytotoxicity mechanism of the prepared formula, a cell cycle analysis was conducted in the next section.
Pharmaceutics 2023, 15, x FOR PEER REVIEW 12 of 24 Figure 5d, the two treatments were found to significantly increase early and late apoptosis compared to the control, which is consistent with previous reports of the antiproliferative activity of propranolol [34]. The cytotoxicity of free propranolol was mainly due to the additive effect of early and late apoptosis, while necrosis played a minor role. Similarly, propranolol-loaded trehalosomes showed cytotoxicity due to the additive effect of early and late apoptosis, with a lower contribution from necrosis. Interestingly, apoptosis was similar in the case of free and encapsulated propranolol, indicating the trehalosomes' ability to induce programmed death in cancer cells [15]. To learn more about the nature of the cytotoxicity mechanism of the prepared formula, a cell cycle analysis was conducted in the next section.

Cell Cycle Analysis
This experiment investigated different phases of the cell cycle, including G0-G1 (apoptosis or non-dividing cells), S (DNA synthesis phase), and G2-M (proliferation or mitosis). An untreated human melanoma cell line was used as a control, and cell cycle analysis of free and encapsulated propranolol was conducted using nuclear propidium iodide staining and flow cytometry as shown in Figure 6a-c. The results show that free propranolol mainly suppressed the human melanoma cell line in the G1 phase, with 61.8% of cells accumulating in this phase. Propranolol also demonstrated the ability to suppress cancer cells in the S and G2 phases, with 25.7% and 15.2% accumulation, respectively. This is consistent with previous reports of propranolol's ability to de-regulate cell cycle progression and induce caspase threeand nine-mediated apoptosis [35]. Similarly, propranolol-loaded trehalosomes suppressed the human melanoma cell line in the G1 phase, followed by the S phase and then the G2

Cell Cycle Analysis
This experiment investigated different phases of the cell cycle, including G0-G1 (apoptosis or non-dividing cells), S (DNA synthesis phase), and G2-M (proliferation or mitosis). An untreated human melanoma cell line was used as a control, and cell cycle analysis of free and encapsulated propranolol was conducted using nuclear propidium iodide staining and flow cytometry as shown in Figure 6a-c. The results show that free propranolol mainly suppressed the human melanoma cell line in the G1 phase, with 61.8% of cells accumulating in this phase. Propranolol also demonstrated the ability to suppress cancer cells in the S and G2 phases, with 25.7% and 15.2% accumulation, respectively. This is consistent with previous reports of propranolol's ability to de-regulate cell cycle progression and induce caspase three-and nine-mediated apoptosis [35]. Similarly, propranolol-loaded trehalosomes suppressed the human melanoma cell line in the G1 phase, followed by the S phase and then the G2 phase, with 67.8%, 21.1%, and 15.7% accumulation, respectively. This can be attributed to the presence of trehalose within the nanocarrier structure [10]. There was a marked buildup in the cells in the G1 phase by nearly 1.2 times for free propranolol and propranolol-loaded trehalosome-treated cells when compared to the control, as shown in Figure 6d. Moreover, a low concentration of cells transitioned to the S and G2 phases. Moreover, there is a noticeable increase in the sub G1 phase in cells treated with propranolol and propranolol-loaded trehalosomes by two folds due to apoptosis induction, as displayed in Figure 6d. These results demonstrated the G1 phase cell cycle arrest for the propranolol and propranolol-loaded trehalosome. Furthermore, the current results confirmed the role of trehalosomes in the prevention of malignant melanoma cell proliferation, which is attributed to both propranolol and propranolol-loaded trehalosomes [10].  [10]. There was a marked buildup in the cells in the G1 phase by nearly 1.2 times for free propranolol and propranolol-loaded trehalosome-treated cells when compared to the control, as shown in Figure 6d. Moreover, a low concentration of cells transitioned to the S and G2 phases. Moreover, there is a noticeable increase in the sub G1 phase in cells treated with propranolol and propranolol-loaded trehalosomes by two folds due to apoptosis induction, as displayed in Figure 6d. These results demonstrated the G1 phase cell cycle arrest for the propranolol and propranolol-loaded trehalosome. Furthermore, the current results confirmed the role of trehalosomes in the prevention of malignant melanoma cell proliferation, which is a ributed to both propranolol and propranolol-loaded trehalosomes [10].

Molecular Docking between Propranolol and Different Melanocytes Receptors
A CB-Dock docking server (http://clab.labshare.cn/cb-dock/php/index.php (accessed on 2 March 2023)) [19,20] was used to calculate Gibbs free energy. This was carried out to study the individual strengths of the interactions between propranolol on the one hand and different studied melanocyte receptors on the other hand. Different melanocyte receptors that have been reported in the literature to be involved in melanoma were selected [36]. It was observed that propranolol has a high binding affinity to all of the studied enzymes, as the Gibbs free energy ranged from-6.8 in the case of endothelin and the frizzled receptors

Molecular Docking between Propranolol and Different Melanocytes Receptors
A CB-Dock docking server (http://clab.labshare.cn/cb-dock/php/index.php (accessed on 2 March 2023)) [19,20] was used to calculate Gibbs free energy. This was carried out to study the individual strengths of the interactions between propranolol on the one hand and different studied melanocyte receptors on the other hand. Different melanocyte receptors that have been reported in the literature to be involved in melanoma were selected [36]. It was observed that propranolol has a high binding affinity to all of the studied enzymes, as the Gibbs free energy ranged from-6.8 in the case of endothelin and the frizzled receptors up to-8.6 in the case of the vitamin D3 receptors, as shown in Table 6. This was found to be in accordance with what has been reported previously in the literature that states that beta-blockers possess anti-cancer activity in melanoma patients [33,37,38]. Table 6.
Gibbs free energy shows the interaction between propranolol and different melanocyte receptors.

Propranolol Enzyme
Gibbs Free Energy

Visualization of Interaction between Propranolol and Different Melanocyte Receptors
The individual interactions between propranolol and different melanocyte receptors were visualized using Discovery Studio-19 software to determine their implications in melanomas.

G-Protein Coupled Receptors
Interaction visualizations are demonstrated in the following sections.

Corticotrophin-Releasing Hormone Receptor (CRH-R1)
Both melanocytes and melanoma cells have been reported to express CRH-R1 responding to CRH peptides via activating cAMP, IP3, and Ca-mediated pathways in order to modify the phenotype of melanocytes. CRH-R1 regulates cell proliferation and differentiation. As shown in Figure 7 and Table 7, propranolol could interact with CRH-R1 via eight interactions: one pi-pi stacking with PHF (B: 203), two hydrophobic (alkyl-sulfur) interactions with MET (B:206), two hydrophobic (alkyl-alkyl) interactions with LEU (B: 280) and LEU (B: 320), and three hydrogen bonds, two with ASN (B: 283) and one with PHE (B: 284). These interactions suggest that propranolol can be a good candidate as a CRH-R1 antagonist, and consequently as a cell proliferation inhibitor in malignant tumors, including melanoma. This is found to be in good agreement with what has been reported in the literature concerning CRH-R1 antagonists as anti-cancer agents, including melanoma [36,39].

Frizzled receptor
Frizzled receptors initiate signals resulting in β-catenin stabilization upon Wnt binding. This leads to regulating the genes responsible for tumorigenesis, e.g., cyclin D1, c-Myc, and matrix metalloproteinase. Hence, blocking frizzled receptors would enhance malignancy inhibition. As shown in Figure 8 and Table 8, propranolol could interact with frizzled receptors through ten interactions: two hydrophobic (alkyl-anion) interactions with GLU (B: 112) and ASP (C: 60), two hydrophobic (alkyl-alkyl) interactions with LYS (B: 115) and ALA (B: 116), one hydrophobic (alkyl-sulfur) with CYS (B: 135), three hydrogen bonds with ASN (C: 58), GLU (C: 64), and ASN (D: 143), and two hydrophobic (alkyl-cation) interactions with ARG (D: 108). These interactions suggest that propranolol can be a good candidate as a frizzled antagonist and consequently as a melanoma suppressor. This is found to be in concurrence with what has been reported in the literature concerning frizzled antagonists as anti-melanoma agents [36,40,41].   •

Frizzled receptor
Frizzled receptors initiate signals resulting in β-catenin stabilization upon Wnt binding. This leads to regulating the genes responsible for tumorigenesis, e.g., cyclin D1, c-Myc, and matrix metalloproteinase. Hence, blocking frizzled receptors would enhance malignancy inhibition. As shown in Figure 8 and Table 8, propranolol could interact with frizzled receptors through ten interactions: two hydrophobic (alkyl-anion) interactions with GLU (B: 112) and ASP (C: 60), two hydrophobic (alkyl-alkyl) interactions with LYS (B: 115) and ALA (B: 116), one hydrophobic (alkyl-sulfur) with CYS (B: 135), three hydrogen bonds with ASN (C: 58), GLU (C: 64), and ASN (D: 143), and two hydrophobic (alkyl-cation) interactions with ARG (D: 108). These interactions suggest that propranolol can be a good candidate as a frizzled antagonist and consequently as a melanoma suppressor. This is found to be in concurrence with what has been reported in the literature concerning frizzled antagonists as anti-melanoma agents [36,40,41].

Frizzled receptor
Frizzled receptors initiate signals resulting in β-catenin stabilization upon Wnt binding. This leads to regulating the genes responsible for tumorigenesis, e.g., cyclin D1, c-Myc, and matrix metalloproteinase. Hence, blocking frizzled receptors would enhance malignancy inhibition. As shown in Figure 8 and Table 8, propranolol could interact with frizzled receptors through ten interactions: two hydrophobic (alkyl-anion) interactions with GLU (B: 112) and ASP (C: 60), two hydrophobic (alkyl-alkyl) interactions with LYS (B: 115) and ALA (B: 116), one hydrophobic (alkyl-sulfur) with CYS (B: 135), three hydrogen bonds with ASN (C: 58), GLU (C: 64), and ASN (D: 143), and two hydrophobic (alkyl-cation) interactions with ARG (D: 108). These interactions suggest that propranolol can be a good candidate as a frizzled antagonist and consequently as a melanoma suppressor. This is found to be in concurrence with what has been reported in the literature concerning frizzled antagonists as anti-melanoma agents [36,40,41].   •

Melanocortin receptor
Melanocortin receptors are involved in gene polymorphism, which is responsible for pigmentation and sensitization toward UV radiation. Overexpression of these receptors would consequently lead to the propagation of melanoma. Hence, the blockage of these receptors would help in minimizing melanoma propagation. As shown in Figure 9 and Table 9, propranolol could interact with melanocortin receptors through eight interactions: two hydrogen bonds with ARG (A: 218) and THR (N: 104); three hydrophobic (alkyl-alkyl) interactions, one with ARG (A: 266) and two with ALA (N: 101); and three hydrophobic (alkyl-anion) with ASP (B: 228), ASP (B: 290), and ASP (N: 106). These interactions suggest that propranolol can be a good candidate as a melanocortin antagonist and consequently as a melanoma suppressor. This is found to be in alignment with what has been reported in the literature concerning melanocortin antagonists as anti-melanoma agents [36,42].

Melanocortin receptor
Melanocortin receptors are involved in gene polymorphism, which is responsible for pigmentation and sensitization toward UV radiation. Overexpression of these receptors would consequently lead to the propagation of melanoma. Hence, the blockage of these receptors would help in minimizing melanoma propagation. As shown in Figure 9 and Table 9, propranolol could interact with melanocortin receptors through eight interactions: two hydrogen bonds with ARG (A: 218) and THR (N: 104); three hydrophobic (alkyl-alkyl) interactions, one with ARG (A: 266) and two with ALA (N: 101); and three hydrophobic (alkyl-anion) with ASP (B: 228), ASP (B: 290), and ASP (N: 106). These interactions suggest that propranolol can be a good candidate as a melanocortin antagonist and consequently as a melanoma suppressor. This is found to be in alignment with what has been reported in the literature concerning melanocortin antagonists as anti-melanoma agents [36,42].

Melatonin receptor
Melatonin agonists have been reported for their efficiency in producing melatonin as a protectant against external stress factors that may lead to melanoma. Hence, the interaction between propranolol and melatonin receptors could lead to participation in suppressing melanoma. As shown in Figure 10 and Table 10, propranolol could interact with melanocortin receptors through six interactions: one hydrophobic (alkyl-sulfur) interaction with MET (A: 107), a pi-sigma bonding and hydrophobic (alkyl-alkyl) interaction with VAL (A: 111), a donor-donor interaction with GLN (A: 181), and a hydrogen bonding and a donor-donor interaction ASN (A: 255). This suggests that propranolol can be a good activator for melatonin receptors. This is found to be similar to what has been reported in the literature concerning the role of melatonin receptors in melanoma treatment [43,44].

• Melatonin receptor
Melatonin agonists have been reported for their efficiency in producing melatonin as a protectant against external stress factors that may lead to melanoma. Hence, the interaction between propranolol and melatonin receptors could lead to participation in suppressing melanoma. As shown in Figure 10 and Table 10, propranolol could interact with melanocortin receptors through six interactions: one hydrophobic (alkyl-sulfur) interaction with MET (A: 107), a pi-sigma bonding and hydrophobic (alkyl-alkyl) interaction with VAL (A: 111), a donor-donor interaction with GLN (A: 181), and a hydrogen bonding and a donor-donor interaction ASN (A: 255). This suggests that propranolol can be a good activator for melatonin receptors. This is found to be similar to what has been reported in the literature concerning the role of melatonin receptors in melanoma treatment [43,44]. Table 10. A summary of types of interactions detected between propranolol and melatonin receptor.

Metabotropic glutamate receptor
Metabotropic glutamate receptor 1 (GRM1) is overexpressed in the incidence of melanoma. Hence, its inhibition suppresses melanoma progression. As shown in Figure 11 and Table 11, propranolol could interact with metabotropic glutamate receptors through seven interactions: two hydrogen bonds with THR (B: 190) and MET (B: 463), two hydrophobic (alkyl-alkyl) interactions with ALA (B: 326) and ARG (B: 465), two pi-pi stacking interactions with HIS (B: 328), and one hydrophobic (alkyl-sulfur) interaction with MET (B: 391). This is found to be in conformity with what has been reported in the literature concerning metabotropic glutamate antagonists as anti-melanoma agents [36].

• Metabotropic glutamate receptor
Metabotropic glutamate receptor 1 (GRM1) is overexpressed in the incidence of melanoma. Hence, its inhibition suppresses melanoma progression. As shown in Figure 11 and Table 11, propranolol could interact with metabotropic glutamate receptors through seven interactions: two hydrogen bonds with THR (B: 190) and MET (B: 463), two hydrophobic (alkyl-alkyl) interactions with ALA (B: 326) and ARG (B: 465), two pi-pi stacking interactions with HIS (B: 328), and one hydrophobic (alkyl-sulfur) interaction with MET (B: 391). This is found to be in conformity with what has been reported in the literature concerning metabotropic glutamate antagonists as anti-melanoma agents [36].

Metabotropic glutamate receptor
Metabotropic glutamate receptor 1 (GRM1) is overexpressed in the incidence of melanoma. Hence, its inhibition suppresses melanoma progression. As shown in Figure 11 and Table 11, propranolol could interact with metabotropic glutamate receptors through seven interactions: two hydrogen bonds with THR (B: 190) and MET (B: 463), two hydrophobic (alkyl-alkyl) interactions with ALA (B: 326) and ARG (B: 465), two pi-pi stacking interactions with HIS (B: 328), and one hydrophobic (alkyl-sulfur) interaction with MET (B: 391). This is found to be in conformity with what has been reported in the literature concerning metabotropic glutamate antagonists as anti-melanoma agents [36].   Growth Factors C-Kit (CD117) is a growth factor for melanocytes that affects their proliferation, melanogenesis, and survival. Thus, the overexpression of c-kit in melanoma should be targeted with antagonists to suppress that propagation. As shown in Figure 12 and Table 12, propranolol could interact with C-Kit via eight interactions: four hydrophobic (alkyl-alkyl) interactions, two with LEU (A: 644), one with ILE (A: 653), and one with LEU (A: 783); and four hydrogen bonds, one with CYS (A: 809) and three with ASP (A: 810). This is found to be in good agreement with what has been reported in the literature concerning C-Kit antagonists as anti-melanoma agents [36,45]. Growth Factors C-Kit (CD117) is a growth factor for melanocytes that affects their proliferation, melanogenesis, and survival. Thus, the overexpression of c-kit in melanoma should be targeted with antagonists to suppress that propagation. As shown in Figure 12 and Table 12, propranolol could interact with C-Kit via eight interactions: four hydrophobic (alkyl-alkyl) interactions, two with LEU (A: 644), one with ILE (A: 653), and one with LEU (A: 783); and four hydrogen bonds, one with CYS (A: 809) and three with ASP (A: 810). This is found to be in good agreement with what has been reported in the literature concerning C-Kit antagonists as anti-melanoma agents [36,45].

Endothelin (EDN) Receptors
In melanoma, EDN has been found to stimulate the down-regulation of E-cadherin expression and the up-regulation of N-cadherin, leading to an increase in tumor proteolytic activity, which in turn provokes melanoma cell adhesion, proliferation, and migration. Hence, blocking EDN results in decreasing melanoma cell viability as well as DNA synthesis while increasing melanoma apoptosis. As shown in Figure 13 and Table 13

Endothelin (EDN) Receptors
In melanoma, EDN has been found to stimulate the down-regulation of E-cadherin expression and the up-regulation of N-cadherin, leading to an increase in tumor proteolytic activity, which in turn provokes melanoma cell adhesion, proliferation, and migration. Hence, blocking EDN results in decreasing melanoma cell viability as well as DNA synthesis while increasing melanoma apoptosis. As shown in Figure 13 and Table 13,

Death Receptors
Death receptors such as Fas (CD95) can be stimulated to trigger the cellular death of melanocytes, including melanoma. Hence, agents that are capable of triggering these death receptors have proven their efficiency in the treatment of metastatic melanoma. As shown in Figure 14 and Table 14, propranolol could interact with Fas receptors via eight interactions: three hydrophobic (alkyl-alkyl) interactions, one with LYS (E: 299) and two with LEU (E: 303); one hydrophobic (alkyl-cation) interaction with GLU (E: 337); two hydrophobic (alkyl-anion) interactions with GLU (F: 179); one hydrogen bond with SER (G: 325); and one pi-pi stacking with PHF (G: 327). This is found to be in alignment with what has been reported in the literature concerning death receptor activators as anti-melanoma agents [36,47].

Death Receptors
Death receptors such as Fas (CD95) can be stimulated to trigger the cellular death of melanocytes, including melanoma. Hence, agents that are capable of triggering these death receptors have proven their efficiency in the treatment of metastatic melanoma. As shown in Figure 14 and Table 14, propranolol could interact with Fas receptors via eight interactions: three hydrophobic (alkyl-alkyl) interactions, one with LYS (E: 299) and two with LEU (E: 303); one hydrophobic (alkyl-cation) interaction with GLU (E: 337); two hydrophobic (alkyl-anion) interactions with GLU (F: 179); one hydrogen bond with SER (G: 325); and one pi-pi stacking with PHF (G: 327). This is found to be in alignment with what has been reported in the literature concerning death receptor activators as anti-melanoma agents [36,47].

Vascular Endothelial Growth Factors
Vascular endothelial growth factors such as the VEGFR-1 receptor possess a tyrosine kinase activity capable of activating signaling pathways that could alter patterns of gene expression in order to enhance melanocyte proliferation. Hence, the blockage of these fac-

Vascular Endothelial Growth Factors
Vascular endothelial growth factors such as the VEGFR-1 receptor possess a tyrosine kinase activity capable of activating signaling pathways that could alter patterns of gene expression in order to enhance melanocyte proliferation. Hence, the blockage of these factors is considered a promising target in suppressing melanoma. As shown in Figure 15 and Table 15

Vascular Endothelial Growth Factors
Vascular endothelial growth factors such as the VEGFR-1 receptor possess a tyrosine kinase activity capable of activating signaling pathways that could alter patterns of gene expression in order to enhance melanocyte proliferation. Hence, the blockage of these factors is considered a promising target in suppressing melanoma. As shown in Figure 15 and Table 15, propranolol could interact with VEGFR-1 via eleven interactions: five hydrophobic (alkyl-alkyl) interactions, two with VAL (A: 841), two with ALA (A: 859), and one with VAL (A: 909); one hydrophobic (alkyl-cation) with LYS (A: 861); four hydrogen bonds, one with GLU (A: 878), one with ILE (A: 1038), and two with ASP (A: 1040); and one hydrophobic (alkyl-sulfur) interaction with CYS (A: 1039). This is found to be in conformity with what has been reported in the literature concerning VEGFR receptor inhibitors as anti-melanoma agents owing to their antiproliferative and anti-angiogenic activities [35,47,48].    Besides their involvement in calcium metabolism, vitamin D3 receptors play a role in melanocyte regulation and apoptosis. Hence, it has been reported that vitamin D derivatives or vitamin D3 receptor triggers could aid in melanoma treatment. As shown in Figure 16 and Table 16, propranolol could interact with vitamin D3 receptors via eight interactions: one hydrophobic (alkyl-cation) interaction and one pi-sigma interaction with TYR (A: 143), two hydrophobic (alkyl-alkyl) interactions with LEU (A: 233), one positivepositive interaction with ARG (A: 274), one hydrogen bond with SER (A: 278), and two pi-pi stacking interactions with TRP (A: 286). This is found to be in alignment with what has been reported in the literature concerning vitamin D3 receptor-triggering agents as anti-melanoma agents [35,49,50].
Besides their involvement in calcium metabolism, vitamin D3 receptors play a role in melanocyte regulation and apoptosis. Hence, it has been reported that vitamin D derivatives or vitamin D3 receptor triggers could aid in melanoma treatment. As shown in Figure 16 and Table 16, propranolol could interact with vitamin D3 receptors via eight interactions: one hydrophobic (alkyl-cation) interaction and one pi-sigma interaction with TYR (A: 143), two hydrophobic (alkyl-alkyl) interactions with LEU (A: 233), one positive-positive interaction with ARG (A: 274), one hydrogen bond with SER (A: 278), and two pi-pi stacking interactions with TRP (A: 286). This is found to be in alignment with what has been reported in the literature concerning vitamin D3 receptor-triggering agents as anti-melanoma agents [35,49,50].  It is inferred from the in silico studies that the highest influential receptors in melanoma are vitamin D receptors, CRH-R1, VEGFR 1, and c-Kit, in descending order. This has been detected by assessing the binding energy that showed −8.6, −8.5, −8.4, and −8.3, respectively. This was also confirmed through detecting high bond interactions through discovery studio visualization. Interestingly, the results of in silico studies were found to match and support the results of experimental apoptotic and cell cycle analysis reported in Sections 3.2.2 and 3.2.3, respectively.
For instance, the vitamin D receptor, which is responsible for melanocyte regulation and apoptosis, could be triggered using propranolol, leading to melanoma apoptosis. This is found to be in good agreement with apoptotic results through caspase enzyme inhibition. In addition, CRH-R1, which was reported for its proliferative activity, has also been found to be strongly inhibited by propranolol. This is found to be in good accordance with what has been reported in cell cycle analysis. Furthermore, in silico studies showed that both c-Kit and VEGFR 1, which are responsible for mutagenesis, can be blocked using propranolol, resulting in anti-proliferation activity. This also supports the experimental results above. Finally, experimental and in silico studies collectively suggest that propranolol can be used as a promising antiproliferative and apoptotic agent in melanoma treatment.  It is inferred from the in silico studies that the highest influential receptors in melanoma are vitamin D receptors, CRH-R1, VEGFR 1, and c-Kit, in descending order. This has been detected by assessing the binding energy that showed −8.6, −8.5, −8.4, and −8.3, respectively. This was also confirmed through detecting high bond interactions through discovery studio visualization. Interestingly, the results of in silico studies were found to match and support the results of experimental apoptotic and cell cycle analysis reported in Sections 3.2.2 and 3.2.3, respectively.
For instance, the vitamin D receptor, which is responsible for melanocyte regulation and apoptosis, could be triggered using propranolol, leading to melanoma apoptosis. This is found to be in good agreement with apoptotic results through caspase enzyme inhibition. In addition, CRH-R1, which was reported for its proliferative activity, has also been found to be strongly inhibited by propranolol. This is found to be in good accordance with what has been reported in cell cycle analysis. Furthermore, in silico studies showed that both c-Kit and VEGFR 1, which are responsible for mutagenesis, can be blocked using propranolol, resulting in anti-proliferation activity. This also supports the experimental results above. Finally, experimental and in silico studies collectively suggest that propranolol can be used as a promising antiproliferative and apoptotic agent in melanoma treatment.

Conclusions
Propranolol-loaded trehalosomes were developed using the ethanol injection method. A factorial design was used to select the optimum formula, where trehalose, lecithin, and Tween 80 levels were studied to achieve minimum size and maximum charge, EE%, and Q3. The optimized formula was investigated using TEM, DSC, and FTIR. Cell studies were carried out against the human melanoma cell line to measure cytotoxicity, apoptosis/necrosis, and cell cycle arrest. The results show that the selected formula has good cytotoxicity compared to the free drug due to the synergistic effect of the drug and the designed carrier, which make it a good candidate as an antiproliferative agent for treating skin cancer. Therefore, the future plan is to study the efficacy of the combination of the selected formula with the traditional chemotherapeutic agent using an in vivo animal model to decrease the chemotherapeutic agent dose, which in turn will decrease the side effects.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.