Methyleugenol Has an Antidepressant Effect in a Neuroendocrine Model: In Silico and In Vivo Evidence

Major depressive disorder is a severe mood disorder characterized by different emotions and feelings. This study investigated the antidepressant activity of the phenylpropanoid methyleugenol (ME) in adult female mice exposed to a stress model induced by dexamethasone. The animals were randomly divided into groups containing eight animals and were pre-administered with dexamethasone (64 μg/kg subcutaneously). After 165 and 180 min, they were treated with ME (25, 50 and 100 mg/kg intraperitoneally) or imipramine (10 mg/kg intraperitoneally) after 45 min and 30 min, respectively; they were then submitted to tests which were filmed. The videos were analyzed blindly. In the tail suspension test, ME (50 mg/kg) increased latency and reduced immobility time. In the splash test, ME (50 mg/kg) decreased grooming latency and increased grooming time. In the open field, there was no statistical difference for the ME groups regarding the number of crosses, and ME (50 mg/kg) increased the number of rearing and time spent in the center. Regarding in silico studies, ME interacted with dopaminergic D1 and α1 adrenergic pathway receptors and with tryptophan hydroxylase inhibitor. In the in vivo evaluation of the pathways of action, the antidepressant potential of ME (50 mg/kg) was reversed by SCH23390 (4 mg/kg intraperitoneally) dopaminergic D1 receptor, Prazosin (1 mg/kg intraperitoneally) α1 adrenergic receptor, and PCPA (4 mg/kg intraperitoneally) tryptophan hydroxylase inhibitor. Our findings indicate that ME did not alter with the locomotor activity of the animals and shows antidepressant activity in female mice with the participation of the D1, α1 and serotonergic systems.


Introduction
Major depressive disorder (MDD) is a mood disturbance influenced by different emotions and feelings.Typically, MDD may feature feelings of sadness and helplessness, reduced or lost interest in daily activities (anhedonia), lack of motivation, depressed mood, irritability, and even suicidal ideation [1,2].
The World Health Organization (WHO) reported that in 2017 there were about 322 million people with depression worldwide (World Health Organization, 2017).The first year of the COVID-19 pandemic intensified the prevalence by 25% and the records show that, The phenylpropanoid ME by i.p. route showed low toxicity, with no evidence of changes in skin, hair, or eyes (a priori in class 5); thus, the LD50 was estimated to be found in values above 2000 mg kg −1 and below 5000 mg kg −1 .In this manner, a range of secure doses could be determined for subsequent in vivo execution.
Pharmaceuticals 2023, 15, x FOR PEER REVIEW 3 of 17 The phenylpropanoid ME by i.p. route showed low toxicity, with no evidence of changes in skin, hair, or eyes (a priori in class 5); thus, the LD50 was estimated to be found in values above 2000 mg kg −1 and below 5000 mg kg −1 .In this manner, a range of secure doses could be determined for subsequent in vivo execution.

In Silico Methodologies
To investigate whether the D1 dopaminergic and α1 adrenergic pathways are involved in the antidepressant impact of methyleugenol, we carried out a study using molecular docking.Through this study, it is possible to assess the probability of interaction between ME and the active site of the selected proteins.The interaction probability can be analyzed through the MolDock score, which calculates the interaction energy, where the higher the stability of the interaction, the lower the energy value.We also analyzed it through intermolecular interactions between ME and amino acid residues present in the binding site of proteins.We validated the docking analysis through redocking and RMSD calculation.

In Silico Methodologies
To investigate whether the D1 dopaminergic and α1 adrenergic pathways are involved in the antidepressant impact of methyleugenol, we carried out a study using molecular docking.Through this study, it is possible to assess the probability of interaction between ME and the active site of the selected proteins.The interaction probability can be analyzed through the MolDock score, which calculates the interaction energy, where the higher the stability of the interaction, the lower the energy value.We also analyzed it through intermolecular interactions between ME and amino acid residues present in the binding site of proteins.We validated the docking analysis through redocking and RMSD calculation.
Table 1 shows the MolDock score energy values for ME of the three proteins, as well as of the molecules SCH 23390 inhibitor of the dopaminergic D1 pathway, prazosin inhibitor of the α1 adrenergic pathway, and PCPA inhibitor of the tryptophan hydroxylase pathway, and also the compounds crystallized in the binding site of each enzyme: fenoldopam (catechol agonist of the D1 dopaminergic pathway), (+)-cyclazosin (an inverse agonist of the α1 adrenergic pathway), and LP533401 (inhibitor of the tryptophan hydroxylase pathway).For a fair comparison between the energy values for every compound within the binding site of each enzyme, the measurements were also standardized according to the value of their molecular masses.Table 1 also provides the RMSD value for the crystallized ligand of each protein.Table 1 shows the MolDock score energy values for ME of the three proteins, as well as of the molecules SCH 23390 inhibitor of the dopaminergic D1 pathway, prazosin inhibitor of the α1 adrenergic pathway, and PCPA inhibitor of the tryptophan hydroxylase pathway, and also the compounds crystallized in the binding site of each enzyme: fenoldopam (catechol agonist of the D1 dopaminergic pathway), (+)-cyclazosin (an inverse agonist of the α1 adrenergic pathway), and LP533401 (inhibitor of the tryptophan hydroxylase pathway).For a fair comparison between the energy values for every compound within the binding site of each enzyme, the measurements were also standardized according to the value of their molecular masses.Table 1 also provides the RMSD value for the crystallized ligand of each protein.

Discussion
In the present study, we investigated and confirmed the antidepressant psychopharmacological profile of the phenylpropanoid methyleugenol using in vivo behavioral and in silico computational methodologies.
The concept that products obtained from natural origin possess important pharmacological properties is widely spread by the scientific community and the population.However, it is important to be aware of the possibility that they present toxic effects that are often overlooked [26].Pharmacological screening qualitatively assesses their actions on the central nervous system (CNS) and autonomic nervous system (ANS) and estimating the acute toxicity of any substance is of utmost relevance to ensure the safety of its use.Methyleugenol at a dose of 2000 mg/kg caused behavioral changes suggestive of drugs with CNS activity but without causing toxic effects and lethality.This estimate of LD50 allowed us to evaluate the doses that were used in the behavioral tests, free of toxic effects.
The antidepressant effect of methyleugenol was investigated in a model of depression induced by dexamethasone applied i.p. [27].Major depressive disorders have been linked to dysregulation of the hypothalamic-pituitary-adrenal (HPA) axis, which involves the secretion of glucocorticoids in excess [28], and persistent activation of its receptors by high levels of glucocorticoids can cause neural damage, with molecular and behavioral alterations, resulting in MDD and other neuropsychiatric disorders [29].Patients treated with glucocorticoids have exhibited psychiatric and cognitive symptoms similar to depression.Already in animal models, chronic and acute glucocorticoid administration caused behavioral modifications similar to those found in humans with the disorder [30,31].Thus, in this animal model of inducing depressive-like behavior, it is possible to mimic symptoms associated with the disease and develop pathophysiological aspects analogous to the human condition where treatment with antidepressant drugs can reverse the changes in the model [32,33].
Besides the choice of model induction that induces depressive-symptomatic behavior, the tests chosen to evaluate this behavior are also important.Indeed, the tail suspension test is one of the most widely used tests for screening antidepressant drugs [34] and is considered sensitive and specific for most classes of antidepressant drug [35].Thus, the animal placed in an inescapable situation (trapped by the tail) develops a posture of immobility after escape attempts, and we perceived this as giving up, which characterizes a depressive state.The antidepressant substances analyzed can reverse such a condition by reducing the time of immobility [36,37].Methyleugenol (at a dose of 50 mg/kg) increased latency in the tail suspension test, showing an effect analogous to that of imipramine and showing the antidepressant-similar activity of the substance.Although little evaluated, the verification of latency contributes to further understanding of the effects of drugs, as certain classes of antidepressants alter this parameter [38].
In 2014, Li and collaborators conducted a study using the Dex protocol, in which Dex was administered during the neonatal period (1 to 3 days) in doses ranging from 0.1 to 0.5 mg/kg body weight.Despite differences in the protocol when compared to the one used in this study, Dex increased the immobility time of the animals in the tail suspension test when compared to the control group, confirming the results obtained in this study.Still validating our data, [39] found the same result as the previous authors when administering a higher dose of Dex (1.5 mg/kg) in rats.Recently, an increase in immobility time was observed in male mice submitted to single-dose injections of Dex (64 µg/kg) in the forced swim test [31,40].
To corroborate our findings in the tail suspension test, we performed the splash test.Animals submitted to this test exhibit a shorter grooming time, which is interpreted as anecdotal behavior.In contrast, administering classical antidepressants increases the time spent on this behavior [41,42].In turn, methyleugenol (intermediate dose 50 mg/kg) decreased the latency to grooming as well as increased the total grooming time, again showing an antidepressant-similar effect in female mice and similar to the standard imipramine.In a study similar to ours, acute treatment of a single dose of dexamethasone increased immobility time in the tail suspension test, inducing depressive behavior in animals, which was counteracted with venlafaxine, a conventional antidepressant [43].
The open field test is widely used for assessing spontaneous locomotor activity and provides clues to the emotional state of rodents [44][45][46].In our study, in the open field test, only the animals that received the positive control (diazepam) and the dexamethasone group reduced the number of crosses.Methyleugenol (50 mg/kg) and animals receiving diazepam increased the number of rearing and dwell time in the center, related to Dex.A methyleugenol dose of 25 mg/kg decreased the latency to rearing.This set of results suggests that methyleugenol did not alter spontaneous locomotor activity and showed an effect indicative of anxiolytic activity, providing clues to the emotional state of the mice.
Once we identified the antidepressant-similar activity of ME in female mice, we needed to understand how this action occurs.Accordingly on the principles of the 3Rs (Replacement, Reduction, and Refinement) in non-clinical research (Russel and Burch, 1992), we conducted the study in silico to evaluate in vivo by which route(s) the ME was exerting the antidepressant effect.
In molecular docking studies, we found interaction probabilities between methyleugenol and receptors of the D1, α1 adrenergic, and tryptophan hydroxylase pathways.The dopaminergic D1 receptor (PDB ID 7CKW), obtained from the PDB Data Bank, is crystallized with the catechol agonist fenoldopam.According to the literature, the active site of D1 receptor agonists and antagonists act on the same active site, and residues Asp103, Trp 321, Phe 289, Phe 288, Trp 285, Asn 292, and Ser 202 are considered the most influential for antagonist activity at the D1 receptor [47,48].
The interactions occurring between the agonist fenoldopam and the active site's amino acid residues of the D1 receptor are stronger than the interactions between the an- Molecular interactions between the compounds (+)-cyclazosin, prazosin, and methyleugenol, and active site's amino acid residues of the α1 adrenergic protein (PDB ID: 7B6W).
These interactions explain the energetic behavior of the connections within the compound and the binding site of the α1 adrenergic protein.Prazosin had the strongest interaction, followed by methyleugenol and (+)-cyclazosin.Imipramine was not evalu-ated for its interaction with any of the proteins in this study as its mode of operation involves inhibiting the reuptake of serotonin and noradrenaline, rather than influencing these proteins.
The dopaminergic D1 receptor (PDB ID 7CKW), obtained from the PDB Data Bank, is crystallized with the catechol agonist fenoldopam.According to the literature, the active site of D1 receptor agonists and antagonists act on the same active site, and residues Asp103, Trp 321, Phe 289, Phe 288, Trp 285, Asn 292, and Ser 202 are considered the most influential for antagonist activity at the D1 receptor [47,48].
The interactions occurring between the agonist fenoldopam and the active site's amino acid residues of the D1 receptor are stronger than the interactions between the antagonists SCH 23390 and methyleugenol.Fenoldopam makes hydrogen-type interactions with the Ser 202; engages in van der Waals contacts with residues Asp 103, Asn 292, Trp 321, Phe 289 and Trp 285; and Pi-Pi interaction with the amino acid residue Phe 288.
Both the control antagonist used in this study (SCH 23390) and the target molecule (methyleugenol) share similarities in interactions with active sites' amino acid residues within the D1 receptor.Both perform van der Waals interactions with residues Asp 103, Asn 292, Phe 289, and Ser 202.The SCH 23390 also performs van der Waals interactions with residues Phe 288, Trp 285, and Trp 321, while methyleugenol performs alkyl or Pi-Pi type interactions with residues Trp 321, Trp 285, and Phe 288.The similarities between the antagonists are also observed in the interaction energy (Table 1).Figure 6 shows the interactions between these molecules and the D2 receptor.The tryptophan hydroxylase protein (PDB ID 7ZIK) obtained from the PDB dataBank is crystallized with the LP53340 inhibitor.Based on the analysis of the active site, the amino acid residues Thr 265, Glu 317, Tyr 235, Ser 336, Ser 337, Gly 333, and Arg 257 were considered the most influential to inhibit the action of the protein [49].
The LP533401 inhibitor crystallized with the protein and the inhibitor used in the in vitro analysis, PCPA, performs hydrogen interactions with amino acid residues Thr 265 and Ser 336, and van der Waals interactions with Glu 317, Ser 337.Methyleugenol interacts The tryptophan hydroxylase protein (PDB ID 7ZIK) obtained from the PDB dataBank is crystallized with the LP53340 inhibitor.Based on the analysis of the active site, the amino acid residues Thr 265, Glu 317, Tyr 235, Ser 336, Ser 337, Gly 333, and Arg 257 were considered the most influential to inhibit the action of the protein [49].
The LP533401 inhibitor crystallized with the protein and the inhibitor used in the in vitro analysis, PCPA, performs hydrogen interactions with amino acid residues Thr 265 and Ser 336, and van der Waals interactions with Glu 317, Ser 337.Methyleugenol interacts with hydrogen only with the amino acid residue Tyr 235, and with the other six amino acid residues considered to be most important, methyleugenol interacts through van der Waals interactions.The LP533401 inhibitor also performs Pi-Pi-type interactions with the Tyr 235 and a salt bridge interaction with the Arg residue 257.In addition, the PCPA inhibitor performs alkyl-type interactions with the Tyr residue 235 and an unfavorable hydrogen interaction with the amino acid residue Arg 257. Figure 7 shows the interactions between these molecules in the active site of the tryptophan hydroxylase protein.We assessed the possibility of interaction through the noradrenergic system in vivo.Because of their involvement in the pathophysiology of depression, some antidepressants increase the availability of noradrenaline at synapses [50].α1-and α2-adrenergic receptors are linked to antidepressant drug responses in behavioral models of depression [51,52].
In the in vivo evaluation, methyleugenol had its antidepressant effect blocked by pretreatment with prazosin (α1-adrenergic antagonist), increasing immobility time in the tail suspension test.Thus, we suggest that ME interacts with the α1-adrenergic receptors of the noradrenergic system to exert its antidepressant-similar effect.Similar to our findings, [53] demonstrated that the structural analog of eugenol (α-asarone) exhibited antidepressant-sympathetic activity reversed by prazosin pretreatment.
Another important system implicated in depression is the serotoninergic system.The increase in the level of 5-hydroxytryptamine (serotonin) would be related to the combatting of depression [54].We used PCPA, a selective inhibitor of serotonin synthesis that inhibits the enzyme tryptophan hydroxylase and depletes the level of serotonin without affecting the levels of noradrenaline and dopamine in the brain.In the present study, PCPA administered alone showed a significant difference regarding immobility time in We assessed the possibility of interaction through the noradrenergic system in vivo.Because of their involvement in the pathophysiology of depression, some antidepressants increase the availability of noradrenaline at synapses [50].α1and α2-adrenergic receptors are linked to antidepressant drug responses in behavioral models of depression [51,52].
In the in vivo evaluation, methyleugenol had its antidepressant effect blocked by pretreatment with prazosin (α1-adrenergic antagonist), increasing immobility time in the tail suspension test.Thus, we suggest that ME interacts with the α1-adrenergic receptors of the noradrenergic system to exert its antidepressant-similar effect.Similar to our findings, [53] demonstrated that the structural analog of eugenol (α-asarone) exhibited antidepressant-sympathetic activity reversed by prazosin pretreatment.
Another important system implicated in depression is the serotoninergic system.The increase in the level of 5-hydroxytryptamine (serotonin) would be related to the combatting of depression [54].We used PCPA, a selective inhibitor of serotonin synthesis that inhibits the enzyme tryptophan hydroxylase and depletes the level of serotonin without affecting the levels of noradrenaline and dopamine in the brain.In the present study, PCPA administered alone showed a significant difference regarding immobility time in the tail suspension test compared to methyleugenol.When the animals were pretreated with PCPA and then with methyleugenol, they also increased the immobility time, thus reversing the anti-immobility effect of methyleugenol.Similar to our results, the study conducted by [55] found that pretreatment with PCPA altered the effect on the immobility time of Carthamus tinctorius L. in the tail suspension test.
To continue the understanding of the mechanism of action of methyleugenol, we evaluated the possibility that its action is still through the dopaminergic system.This is especially involved in mood regulation [56], and it is known that dysregulation of the mesolimbic and mesocortical dopaminergic pathway is related to signs of melancholy and altered cognition in depression [57][58][59].In addition, there are already records of the efficacy of antidepressants that act on the dopaminergic pathway [60].
Thus, involvement of the dopaminergic (D1) system in the antidepressant effect of ME was found, as pretreatment of the mice with SCH23390 (D1 receptor antagonist) blocked the effect induced by methyleugenol, thus increasing the immobility time in the tail suspension test.Only the administration of SCH23390 increased immobility, as was also observed in the findings of [35], who by administering SCH23390 to mice, reversed the antidepressant effect of gallic acid.
Our findings obtained in the present study allow us to conclude that the administration of methyleugenol in animals submitted to the dexamethasone-induced depression model produces a relevant antidepressant effect without altering locomotor activity.Its action on depressive-like behavior seems to be associated with the involvement of the D1 dopaminergic, α1 adrenergic, and serotoninergic systems.
We prepared the substances in their doses minutes before use.The decimal concentrations used allowed the injection of 0.1 mL/10 g of animal weight.Methyleugenol was dissolved in distilled water with the help of tween 80, in doses 25 mg/kg, 50 mg/kg, and 100 mg/kg, for intraperitoneal (i.p.) administration.For the negative control, Tween 80 1% in distilled water (i.p.) was used.

Animals
Female albino Swiss mice (Mus musculus), weighing 25-35 g, approximately 3 months old, were housed in polyethylene cages and kept under controlled temperature (21 ± 1 • C).They had free access to pellet feed (Quimtia ® ) and water available in polyethylene bottles, and were kept under a 12 h light/dark cycle (the light phase from 6 am to 6 pm).We used mice for all experimental protocols.The animals were from the Animal Production Unit (UPA), Institute for Research in Drugs and Medicines (IPeFarM), UFPB.All procedures were conducted according to the guidelines proposed by the National Institutes of Health (NIH) Animal Care and Use Committee and approved by the Animal Ethics Committee (CEUA), UFPB (certificate no.045/2017).

Methyleugenol Acute Toxicity Estimation and Pharmacological Screening
We performed the determination of acute toxicity in the procedures adapted from the Brazilian Health Regulatory Agency [61] and the Organization for Economic Cooperation and Development "Guidelines for testing of chemicals" [62] with modifications.
Groups of Swiss mice (consisting of three mice per group, including the control) were administered a single dose of 2000 mg/kg of ME.The control group received treatment with 0.9% saline + 0.5% Tween (i.p.).Following administration, we conducted observations at intervals of 30 min, 1st hour, 2nd hour, and 4th hour; and subsequently once daily up to the fourteenth day.After the fourteenth day, the animals were weighed and euthanized with an anesthetic overdose injection.From this step, we established the doses from the estimated LD50 value.

Induction of Depressive-like Behavior
To evaluate the antidepressant activity, we initially induced depressive-like behavior in the animals, which was adapted from the studies of Wrobél and colleagues (2014).Dexamethasone (64 µg/kg s.c.) was administered to all groups, except the saline group, 3 h 30 min before we performed the behavioral tests.We ended the experiments by noon because of peak endogenous corticosterone release.

Tail Suspension Test
After the time of administration of dexamethasone (180 min our 165 min), we performed the test according to the protocol described previously [63].Groups of mice were used (n = 8).The groups Vehicle (control group: Tween 80 0.5%), imipramine (30 mg/kg i.p.), and ME (25, 50 and 100 mg/kg i.p.) had their substances administered 30 (imipramine) and 45 (ME) minutes before the beginning of the test.After the treatment period, we suspended the mice 58 cm above the ground with tape placed approximately 1 cm from the tip of the tail.Each animal remained in this position for 6 min.The behavioral parameters observed were latency to immobility and total immobility time.The animal was classified as immobile when it showed no body movement and remained totally still.

Splash Test
To analyze the self-cleaning behavior (grooming) of the animals, a 10% sucrose solution was sprayed on the animal's back.The sucrose solution is thick and induces an unpleasant sensation in the animal, prompting it to engage in body cleaning, a typical self-care action.Animal cleanliness represents a form of self-care behavior, and its decrease is linked to depressive symptoms (10).After administration of dexamethasone, groups of mice were used (n = 8).The groups Vehicle (control group: Tween 80 0.5%), imipramine (30 mg/kg i.p.), and ME (25, 50 and 100 mg/kg i.p.) had their substances administered 30 (imipramine) and 45 (ME) minutes before the beginning of the test.The animal was individually placed inside an acrylic box (9 × 7 × 11 cm) and the latency parameters for the onset of self-cleaning and the time the animal remained in the self-cleaning behavior were observed for 5 min.Between observations, we sanitized the apparatus with 10% alcohol.

Open Field Test
The aim is to observe the locomotion and exploratory behavior of animals in the open field apparatus, similar to that described by an early study [64].The field is manufactured of polyethylene measuring 30 cm × 30 cm, and the floor (50 cm × 50 cm) has 12 equal squares.The groups Vehicle (control group: Tween 80 0.5%), diazepam (DZP) (1 mg/kg i.p.), and ME (25, 50 and 100 mg/kg i.p.) had their substances administered 30 (DZP) and 45 (ME) minutes before the beginning of the test.The mice were positioned in the central region of the arena, and the number of animal crossings in the quadrants, dwell time in the center of the arena, amount of lifting (rearing), and latency to 1st rearing were observed for 5 min.The dwell time in the central area (4 central quadrants) was considered a measure of anxiety.
We cleaned the base of the arena with 10% ethanol between tests to avoid the effects of odors left by previous subjects.

Investigation of the Mechanisms of Action of the Antidepressant Activity of Methyleugenol
The dose of 50 mg/kg of methyleugenol was used because of its good performance in behavioral tests.The experimental groups received dexamethasone (64 µg/kg s.c.) about 3 h 30 min before the beginning of the evaluation; we chose the tail suspension test for this step.
To ascertain the participation of the noradrenergic system in the antidepressant activity of methyleugenol, Prazosin, an α1 receptor antagonist [53,55], was used.Afterward, the different groups received the following treatments: vehicle (Tween 80 1% i.p.), ME (50 mg/kg i.p.), or prazosin (1 mg/kg i.p.).The other group received ME administration 15 min after the prazosin application.After 45 min, we performed the tail suspension test.
To investigate the serotoninergic pathway, p-chlorophenylalanine (PCPA 100 mg/kg i.p.), an inhibitor of tryptophan hydroxylase-serotonin synthesis [55], was used.We also organized two other groups that received the administration of PCPA daily for 5 days; thus, on the day of the experiment, 30 min prior to the test, one group received the methyleugenol (50 mg/kg i.p.) and the other group the saline.After 30 min, we submitted the animals to the tail suspension test.

Video Recording and Analysis
All behavioral tests were engraved, and experienced researchers analyzed the behaviors in a double-blind scheme.

Statistical Analysis
The results obtained during the in vivo experiments passed the Shapiro-Wilk normality test.Parametric tests were used, as the data showed normal distribution.Differences for single-factor group data were ruled out by one-way analysis of variance (ANOVA) followed by Dunnett or Tukey (post hoc).Data were expressed as mean ± SEM (standard error of the mean).A 5% confidence level was adopted.All data analyses were performed with GraphPad Prism 9.0 free version (GraphPad Software, San Diego, CA, USA).

In Silico Methodology Investigation of the Mechanisms of Action of the Antidepressant Activity of Methyleugenol
We performed molecular docking analysis to investigate the mechanisms of the action of methyleugenol via D1 dopaminergic, α1 adrenergic, and tryptophan hydroxylase receptors with antidepressant effects.Both proteins in the complex with the catechol agonist for D1 dopaminergic and inverse agonist for α1 adrenergic and inhibitor LP533401 for tryptophan hydroxylase were obtained from the Protein Data Bank (PDB) at www.rcsb.org(accessed on 1 October 2023): D1 dopaminergic (PDB ID 7CKW, ligand fenoldopam) [65] α1 adrenergic (PDB ID 7B6W, ligand (+)-Cyclazosin) [66] and tryptophan hydroxylase (PDB ID 7ZIK) [49].
We employed redocking as an assessment technique for the docking protocol.In the redocking approach, we utilized the original arrangement of a crystallized ligand bound to the protein and juxtaposed it with the positioning of a ligand affixed to the active site of the identical protein through molecular docking.To gauge the efficacy, we employed the RMSD (root-mean-square deviation) metric, which measures the average spatial disparity between the crystallized ligand and the docked ligand.A docking outcome was considered acceptable if the RMSD measured up to 2.0 Å [67].This evaluation procedure was applied to both proteins under scrutiny in the current investigation.
Methyleugenol underwent molecular docking through Molegro Virtual Docker, version 6.0.1 (MVD).During the docking procedure, elimination of all water molecules was carried out to facilitate the docking process.Preparation of enzyme and compound structures was executed using the default parameter settings within the same software package, adhering to the following configurations: ligand evaluation criteria encompassing internal electrostatics (ES), internal hydrogen bonding, Sp2-Sp2 twists-all of which were verified; a total of 10 runs were conducted; MolDock SE algorithm was employed; interactions capped at 1500; population size limited to 50; a maximum of 300 steps; a neighbor distance factor of 1.00; and retrieval of up to 5 poses.The fitting process employed a GRID with a radius of 10 Å and resolution of 0.30 to encompass the ligand-binding site across both PDB files.
The central point of the target site served as the origin for the grid box, aligning with the ligand, which was then juxtaposed against the remaining trio of proteins, utilizing RMSD as the evaluative yardstick.We generated models incorporating pertinent attributes from both ligands that bore relevance to ligand binding.Employing the MolDock scoring algorithm in conjunction with the MolDock search algorithm [68], we harnessed Molegro Virtual Docker (CLC bio Company, Aarhus, Denmark) to generate five distinct poses for each alkaloid within the active site of every individual protein.Among these, the most stable conformation, characterized by the lowest interaction energy, was cherry-picked and integrated into the Discovery Studio 2021 software for visual scrutiny [69].By employing the MolDock scoring mechanism, we calculated the energy of the crystallized ligand for each protein.Leveraging the MolDock score algorithm, the energy of the ligand from its crystallized configuration was seamlessly converted to the energy scale employed by the algorithm.

Conclusions
The investigation into the antidepressant-like effects of methyleugenol has yielded important results.In this study utilizing female mice, methyleugenol demonstrated a significant improvement in depressive symptoms induced by Dex, without altering locomotor activity.Notably, the observed effects of methyleugenol were associated with modulation of the D1 dopaminergic, α1 adrenergic, and serotoninergic systems.
This suggests a potential therapeutic value of methyleugenol in mitigating depres-sionrelated behaviors in model animals.While these findings hold promise, further research is needed to elucidate the precise molecular mechanisms of methyleugenol's antidepressant action, and explore potential interactions with existing antidepressant treatments.

Figure 3
Figure 3 demonstrates the outcomes derived from the open field assessment.In the evaluation of the number of crossings in the arena (Figure 3A), the group treated with dexamethasone (Dex) (38.8 ± 4.7) and diazepam (DZP) (61.8 ± 5.2) showed a reduction in the number of crossings and differed from the vehicle group (Veh) (86.6 ± 7.4) (F (5, 42) = 9.3; p < 0.0001).

Figure 6 .
Figure 6.Molecular interactions between the compounds fenoldopam, SCH 23390, and methyleugenol and active site's amino acid residues of the dopaminergic protein D1 (PDB ID: 7CKW).

Figure 7 .
Figure 7. Molecular interactions between the compounds LP533401, PCPA, and methyleugenol, and active site's amino acid residues of the tryptophan hydroxylase protein (PDB ID: 7ZIK).

Figure 7 .
Figure 7. Molecular interactions between the compounds LP533401, PCPA, and methyleugenol, and active site's amino acid residues of the tryptophan hydroxylase protein (PDB ID: 7ZIK).

Table 1 .
The MolDock score and the RMSD of the docking molecular.

Table 1 .
The MolDock score and the RMSD of the docking molecular.