Possible Mechanisms Underlying the Antispasmodic, Bronchodilator, and Antidiarrheal Activities of Polarity–Based Extracts of Cucumis sativus L. Seeds in In Silico, In Vitro, and In Vivo Studies

Apart from the nutritional value, Cucumis sativus L. has also been used in the traditional medicine of Iran, Pakistan, and India. Its seeds are used by herbalists to treat gastrointestinal, respiratory, and urinary problems. However, more investigations are required to explain its mechanisms for treating GI, respiratory, and urinary diseases. Accordingly, the aim of the present work was to investigate the antispasmodic, bronchodilator, and antidiarrheal activities of C. sativus seeds extracts and the underlying mechanisms of action. For this purpose, sequential extracts of C. sativus seeds were prepared in n-hexane, dichloromethane, ethanol, and water. Bioactive compounds in C. sativus seed extracts were identified and quantified by utilizing LC ESI–MS/MS and HPLC. Moreover, network pharmacology and molecular docking were employed to examine the antispasmodic and bronchodilator effects of the bioactive substances in the extracts. In vitro and in vivo experiments were also conducted to validate the mechanistic insights gained from the in silico analysis. Results indicated the presence of kaempferol with a concentration of 813.74 µg/g (highest concentration) in the seed extract of C. sativus, followed by quercetin (713.83 µg/g), narcissin (681.87 µg/g), and orientin (676.19 µg/g). In silico investigations demonstrated that the bioactive chemicals in C. sativus seeds inhibited the expression of the target genes involved in smooth muscle contraction and calcium-mediated signaling. Sequential seed extracts of C. sativus caused a dose-dependent relaxant response for spasmolytic reaction and resulted in a relaxation of K+ (80 mM) spastic contraction. In animal models, C. sativus seed extracts exhibited partial or complete antiperistalsis, antidiarrheal, and antisecretory actions. By modulating the contractile response through calcium-mediated signaling target proteins, C. sativus seeds generated bronchodilator, antispasmodic, and antidiarrheal therapeutic effects.


Introduction
Cucumis sativus L., often called cucumber, is a member of the Cucurbitaceae family extensively cultivated as a culinary vegetable in tropical Asia, South America, and Africa. [1,2]. Because of its high vitamin, protein, carbohydrate, and fatty acid content, it is frequently served as a green appetizer [1]. On the other hand, herbalists in Pakistan and Pharmaceuticals 2022, 15  The target genes of the ethanol extract's bioactive compounds were significantly portrayed in numerous GO pathways, with a particular emphasis on neurotransmitter regulation, blood circulation regulation, chemical synaptic transmission modulation, transsynaptic signaling control, homeostasis of cellular calcium ions, homeostasis of cellular divalent inorganic cations, peptidyl threonine phosphorylation, and concentration regulation of cytosolic calcium ions. In addition, the ethanol bioactive compounds modulate the following KEGG signaling pathways: cholinergic synapse; calcium signaling pathway; GnRH signaling system; dopaminergic synapse; lipid and atherosclerosis, and ErbB signaling; contraction of vascular smooth muscle; the production, secretion and action of growth hormones; cardiomyocyte adrenergic signaling; gastric acid secretion; oxytocin signaling; diabetic cardiomyopathy; cancer proteoglycans; and cAMP signaling.
Similarly, the numerous GO terms exhibited a predominance of target genes for aqueous extracts' bioactive compounds clearly emphasizing chemical synaptic transmission modification, trans-synaptic signaling, muscular system regulation, metal ion transport regulation, neurotransmitter levels control, amyloid-beta response, cellular response to chemical stress, and peptidyl-threonine. The route of calcium signaling, cholinergic synapse, proteoglycans in cancer, ErbB cascade pathways, GnRH cascade systems, HIF-1 cascade systems, oxytocin cascade systems, lipids and atherosclerosis, cAMP signaling routes, cGMP-PKG signaling pathways, dopaminergic synapses, vascular smooth muscle contractions, cardiomyocyte adrenergic signaling, and leukocyte transendothelial migration are the key KEGG signaling pathways affected by active substances.

Network Construction
All retained pathogenic targets were subjected to the STRING plugin of Cytoscape for protein-protein interaction (PPI). The target genes for the gastrointestinal and respiratory tracts contained 28 nodes and 79 edges. Between the bioactive chemicals from C. sativus and disease target genes, a C-T-D network with 103 nodes and 265 edges was built in Cytoscape. According to network analysis, kaempferol has the most favorable impact on the target genes with a degree of 45, after quercetin (degree = 43), luteolin (degree = 38), apigenin (degree = 38), and linoleic acid (degree = 23) (Table S9). Based on its interactions with gene ontology biological process (GO-BP) terms (Table S10) and KEGG signaling pathways, a compound target pathway (C-T-P) network was constructed for each C. sativus seed extract ( Figure 2).
In the C-T-P of GO biological process terms, n-hexane extract had 41 nodes and 154 edges. The network study of n-hexane revealed that linoleic acid and palmitic acid, with degrees of 33.6 and 16.8, had stronger, more positive impacts on the disease target genes. The following GO biological process terms had a strong positive interaction with n-hexane bioactive constituents: protein-coupled receptor signaling pathway coupled to cyclic nucleotide second messenger (degree = 76.8); trans-synaptic signaling (degree = 18.9); modulation of chemical synaptic transmission (degree = 18.9); regulation of neurotransmitter levels (degree = 18.9); adenylate cyclase-modulating G protein-coupled receptor signaling pathway (degree = 16.8); and neurotransmitter transport (degree = 16.8). The n-Hexane had 45 nodes and 133 edges in the KEGG pathways C-T-P. The network analysis revealed that linoleic and palmitic acids, with degrees of 44.8 and 11.2, had more powerful and positive impacts on the disease target genes. The following KEGG pathways had a strong positive interaction with n-hexane bioactive constituents: cholinergic synapse (degree = 47.6); neuroactive ligand-receptor interaction (degree = 44.4); signaling route in cAMP (degree = 36); calcium-dependent signaling (degree = 11.2); cancer proteoglycans (degree = 9.6); estrogen signaling pathway (degree = 9.6); and signaling route for relaxin (degree = 9.6).   DCM extract had 39 nodes and 134 edges in the C-T-P of GO biological process terms. β-Sitosterol (degree = 21), umbelliferone (degree = 18.9), and stigmasterol (degree = 18.9) exhibited more powerful therapeutic impacts on the disease targeted genes, according to network analysis. Peptidyl-serine modification (degree = 18.9), peptidyl-serine phosphorylation (degree = 18.9), regulation of neurotransmitter levels (degree = 18.9), regulation of DNA-binding transcription factor activity (degree = 16.8), and modulation of chemical synaptic transmission (degree = 16.8) were the biological processes with the most target genes and compounds. There were 36 nodes and 121 edges in the KEGG pathways C-T-P. β-sitosterol (degree = 12.8), umbelliferone (degree = 12.8), and stigmasterol (degree = 11.2) exhibited more powerful positive impacts on disease target genes, according to the network analysis. Cholinergic synapse (degree = 56), adrenergic signaling in cardiomyocytes (degree = 29.6), the chemokine signaling route (degree = 21.2), and the calcium signaling system (degree = 11.2) were the KEGG pathways that had the highest target genes and chemicals.
Aqueous extract had 56 nodes and 199 edges in the C-T-P of GO biological process terms. Kaempferol (degree = 94.5), kaempferol 3-O-glucoside (degree = 23.1), and ferulic acid (degree = 21) exhibited more powerful therapeutic impacts on disease-targeted genes, according to the network analysis. The muscular system process (degree = 25.2), control of trans-synaptic signaling (degree = 25.2), modulation of chemical synaptic transmission (degree = 25.2), and regulation of metal ion transport (degree = 23.1) were the biological processes with the most target genes and chemicals. The KEGG pathways C-T-P has 52 nodes and 199 edges. According to the results of the network analysis, kaempferol (degree = 65.6), kaempferol 3-O-glucoside (degree = 16), and ferulic acid (degree = 16) have more potent therapeutic effects on disease target genes. The calcium signaling route (degree = 54.4), cholinergic synapse (degree = 36), proteoglycans in cancer (degree = 19.2), and cAMP signaling pathway (degree = 17.6) were the KEGG pathways with the highest target genes and chemicals.

Physicochemical Characteristics
ProtParam was employed to determine the physicochemical characteristics of voltagegated calcium channel sequences, muscarinic acetylcholine receptor, myosin light chain kinase, and phospholipase C-gamma-1. The amino acid profile of the voltage-gated calcium channel revealed a molecular weight of 243,482.46 and 2169 total number of amino acids with a maximum of 10.3% of Leu, 7.7% of Ser, and 7.5% of Ala residues. The negatively charged (Asp + Glu) residues were 229, whereas the positively charged residues (Arg + Lys) were 95. The PI value of 6.60 indicates that the protein is almost neutral, while the instability index of 50.36 shows that it is more unstable than 40. This instability was anticipated owing to the presence of specific dipeptides that are absent from stable proteins. A higher aliphatic index (91.34) indicates that the protein's thermostability is significant, whereas a lower GRAVY score (−0.059) indicates that the protein might interact better with water. At 280 nm, the extinction coefficient is 230,910 M −1 cm −1 for the Cys, Trp, and Tyr concentrations. This coefficient helps quantify the protein-protein and protein-ligand interactions in solution.
The amino acid profile of the muscarinic acetylcholine receptor revealed a molecular weight of 66,065.76 and 589 total number of amino acids with a maximum of 12.8% Ser, 7.1% Ala, and 5.3% Asn residues. The negatively charged (Asp + Glu) residues were 45 and positively charged (Arg + Lys) 62. The PI value of 9.22 indicates that the protein is alkaline, while the instability index is 52.69. A high aliphatic index (86.45) shows that the protein is thermostable, whereas a lower GRAVY score (−0.129) indicates that the protein might interact better with water. At 280 nm, the extinction coefficient is 105,310 M −1 cm −1 for the Cys, Trp, and Tyr concentrations.
The amino acid profile of myosin light chain kinase revealed a molecular weight of 65,816.22 and 610 total number of amino acids with a maximum of 9.2 percent Ala residues, followed by 8.7 percent Leu residues, and 8.5 percent Glu residues. The negatively charged (Asp + Glu) residues were 89, and the positively charged (Arg + Lys) residues were 67. The PI value of 5.05 indicates that the protein is acidic, while the instability index is 48.27. A higher aliphatic index (73.00) showed that the protein is thermostable, whereas a lower GRAVY score (−0.497) indicated that the protein might interact better with water. At 280nm, the extinction coefficient is 33,920 M −1 cm −1 for the Cys, Trp, and Tyr concentrations.
The amino acid profile of phospholipase C-gamma-1 revealed a molecular weight of 148,548.09 and 1290 total number of amino acids with a maximum of 9.2% Leu, followed by 8.8% Glu, and 7.4% Ser residues. The negatively charged (Asp + Glu) residues were 187, and the positively charged (Arg + Lys) residues were 158. The PI value of 5.67 indicates that the protein is acidic, and the instability index is 45.12. A high aliphatic index (74. 47) shows that the protein is thermostable, whereas a low GRAVY score (−0.549) indicates that the protein might interact better with water. At 280 nm, the extinction coefficient is 179,000 M −1 cm −1 for the Cys, Trp, and Tyr concentrations.
The Ramachandran plot was used to verify the model when it was checked using PROCHECK, and this plot reflected the amino acid distribution in phi and psi angles. ( Figure 3B). The calculated values helped in determining the structure. Results showed that the amino acid residues in VGCC, muscarinic M3, MLCK, and PLCγ-1 had 86.9%, 82.9%, 82.2%, and 82.2%, respectively, in the favored regions and 0.9%, 1.3%, 1.5%, and 0.9% in the disallowed regions. When the alignment is completed correctly, the low sequence similarity may give a decent model that uses molecular docking research when completed.

Molecular Docking
Molecular docking calculations may be used to determine the position of the ligand inside the binding pockets of the target proteins. The physical energy terms such as solvation energy are combined with a suitable force field for docking calculation accuracy [33,34]         In spontaneously contracting jejunum preparations, the antispasmodic action of n-hexane, DCM, ethanol, and aqueous seed extracts of C. sativus were studied. All extracts exhibited dose-dependent spasmogenic and spasmolytic effects, and a dose-dependent partial or complete relaxation of K + (80 mM) and K + (25 mM) elicited contractions. When both n-hexane and DCM were exposed to rhythmic or periodic jejunal contractions, they displayed dose-dependent contractile responses ( Figure 6). The spasmogenic responses of n-hexane and DCM extracts were attenuated and a dose-dependent spasmolytic effect was established in an atropine (1 µM) pretreatment jejunum preparation, with EC 50 values of 0.3572 mg/mL (0.2340-0.5689; 95% CI) and 0.8343 mg/mL (0.6196-1.141; 95% CI). On the other hand, ethanol and aqueous seed extracts diminished the spontaneous rhythmic contraction in a dose-dependent pattern, with EC 50 values of 0.3313 mg/mL (0.2237-0.5074; 95% CI) and 0.5625 mg/mL (0.4072-0.7868; 95% CI), respectively.
In addition, the n-hexane, DCM, ethanol, and aqueous seed extracts of C. sativus were evaluated in jejunum preparations against contractions of K + (80 mM) and K + (25 mM). The n-hexane extract showed a dose-dependent partial inhibition towards the sustained contractions of K + (80 mM) with an EC 50 Figure 6).
The CRCs of calcium were established in the presence and absence of extracts in the jejunal preparations that were entirely devoid of cytosolic calcium to demonstrate the VGCC antagonistic activity of C. sativus seed extracts. Ethanol and aqueous seed extracts produced non-parallel rightward shifting and calcium CRCs suppression at doses of 0.3 and 1 mg/mL. Nevertheless, both extracts significantly depressed the calcium concentrationresponse curves at 1 mg/mL. (Figure 6). In contrast, n-hexane and DCM extracts did not reduce calcium CRCs but produced a competitive rightward shift at 3 and 5 mg/mL. Moreover, n-hexane and DCM extracts suppressed calcium CRCs in atropine (1 µM) pretreated tissues. They caused a non-parallel rightward shift of calcium CRCs and suppressed them in atropinized tissue, suggesting that the cholinomimetic components mediated the contractile response in the absence of atropine ( Figure 6) [35,36]. Atropine had no impact on the calcium CRCs [36]. These results suggest that the calcium ion channel antagonists are abundant in the polar region of the C. sativus seed. However, in n-hexane and DCM extracts, the effect was less pronounced. Results were validated and compared to a calcium channel blocker called verapamil to see whether they had any calcium antagonistic effects. Verapamil relaxed the spontaneous concentrations, and K + (80 mM) and K + (25 mM) provoked contractions in jejunal preparations with EC 50 values of 0.05872 µM (0.04739-0.07325; 95% CI), 0.1163 µM (0.1003-0.1354; 95% CI), and 0.09815 µM (0.06016-0.1678; 95% CI), respectively ( Figure 6). In addition, the CRCs for calcium were created and produced a full blockage of calcium CRCs at a concentration of 0.3 µM, comparable to Cucumis sativus extracts (Figure 6), showing that Cucumis sativus extracts exhibit antagonistic action towards calcium ions channels.

Effect of C. sativus Seed Extracts on Isolated Rat Ileum Preparations
A bolus dose of n-hexane and DCM extracts was employed to induce spasmogenic activity in isolated rat ileum preparations. Doses of 0.3, 1, and 3 mg/mL of n-hexane and DCM extracts elicited contractile responses comparable to the cholinergic agonist acetylcholine (0.1, 0.3, and 1 µM). In the pretreatment of atropine (1 µM), the contractile response of n-hexane and DCM was significantly (p < 0.001 vs. Ach 1 µM) reduced, equivalent to cholinergic contraction (Figure 7).  A bolus dose of n-hexane and DCM extracts was employed to induce spasmogenic activity in isolated rat ileum preparations. Doses of 0.3, 1, and 3 mg/mL of n-hexane and DCM extracts elicited contractile responses comparable to the cholinergic agonist acetylcholine (0.1, 0.3, and 1 µM). In the pretreatment of atropine (1 µM), the contractile response of n-hexane and DCM was significantly (p < 0.001 vs. Ach 1 µM) reduced, equivalent to cholinergic contraction (Figure 7). Values are expressed as the mean ± SD (n = 4). The data were evaluated using one-way ANOVA followed by Dunnett's test in comparison to Ach. Results were considered significant at p < 0.05. (* p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001).

Effect of C. sativus Seed Extracts on Isolated Rabbit Tracheal Preparations
Our findings revealed that the n-hexane, DCM, ethanol, and aqueous C. sativus seed extracts exhibit dose-dependent relaxant effects when tested against persistent contrac-  Values are expressed as the mean ± SD (n = 4), and data were evaluated using dose-response curves. for CCh were also produced for verapamil, which resulted in the CRCs for CCh being suppressed at a dose equivalent to that of C. sativus extracts (Figure 8).

Maximum Tolerated Doses for C. sativus Seed Extracts
C. sativus seed extracts were investigated in rats for a maximum tolerated dosage. During the 28-day test period, no changes in body weight, behavioral alterations, clinical symptoms of pain and distress, or death of the rats were observed.

Effect of C. sativus Seed Extracts on GI Charcoal Meal Intestinal Transit
Pretreatment with ethanol and aqueous C. sativus seed extracts induced a significant (p < 0.001 vs. the negative control) reduction in charcoal meal propulsive movement. The peristalsis index of the ethanol extract was 22.6% at 150 mg/kg and 16.30% at 300 mg/kg, whereas the peristalsis index of the aqueous extract was 32.1% at 150 mg/kg and 18.4% at 300 mg/kg. Loperamide and verapamil were used as the positive controls and showed similar results. On the other hand, the n-hexane and DCM extracts failed to reduce peristaltic motions, which differs somewhat from the negative control (71.4%). At 150 mg/kg, the peristalsis index for n-hexane and DCM extracts was 49.6% and 38.6%, respectively, and at 300 mg/kg, it was 64.6% and 52.1%, respectively ( Figure 10).

Effect of C. sativus Seed Extracts on Castor Oil-Induced Diarrhea
This investigation showed that pretreatment of animals with ethanol and aqueous C. sativus seed extracts causes a significant (p < 0.001 vs. the negative control) inhibitory effect on castor oil-induced diarrhea. At 150 mg/kg, the ethanol and aqueous extracts showed 66.6% and 60.4% protection against castor oil-induced diarrhea, respectively, and caused 81.3% and 77.6% protection at 300 mg/kg. That was equivalent to the positive controls of loperamide and verapamil in terms of diarrhea prevention (p < 0.001 vs. the negative control). However, unlike the negative control, the n-hexane and DCM extracts failed to suppress diarrhea (4.3%). At 150 mg/kg, the n-hexane and DCM extracts showed 25.5% and 52% protection, respectively, while at 300 mg/kg, they exhibited 11% and 32.3% protection, which is non-significant against diarrhea ( Figure 10).

Discussion
Herbal therapy has been globally recognized as a safe and effective approach for treating a wide range of chronic disorders. Consumption of non-nutrient secondary metabolites [37] from a high vegetable, fruit, and grain diet can prevent chronic illnesses [38][39][40]. Within this context, local healers and practitioners have long used Cucumis sativus L. seeds to treat various diseases [11,16,25,41]. The present study was designed to explore the possible therapeutic potential of C. sativus seed extracts in digestive and respiratory ailments. Our findings related to phytochemical screening revealed the presence of sterols, anthraquinones (anthracene), alkaloids, flavones, and flavonoids in the seed extracts of C. sativus [11,15,38,42]. In addition, published research indicated that C. sativus seeds are rich in polyphenols such as quercitrin, rutin, and orientin [27], along with phytosterols (β-sitosterol, stigmasterol) [25]. In this respect, β sitosterol and polyphenols such as rutin and orientin exert antispasmodic activity. These phytochemicals in plants, fruits, or seeds may be antispasmodic [43,44].
In silico analyses of bioactive chemicals isolated from C. sativus seed extracts anticipated that these compounds would have antispasmodic action via interfering with calcium-mediated signaling and cholinergic synapses (Figure 1). The GO and KEGG pathways investigation revealed that the target genes of bioactive substances have a direct or indirect impact on smooth muscle contraction via calcium-mediated and cholinergic synapse signaling. The bioactive chemicals in the C. sativus seeds extracts interact with both pathway members In molecular docking studies, these compounds were shown to have a The data are presented as the mean ± SD (n = 5). Data were analyzed using one-way ANOVA followed by Dunnett's test when compared to the control group (normal saline or castor oil); differences were considered significant at p ≤ 0.05. (* p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001). N.S.: Normal saline (10 mL/kg); Lop.: Loperamide (10 mg/kg); Ver. Verapamil (10 mg/kg).

Effect of C. sativus Seed Extracts on Intestinal Fluid Accumulation
When compared to the control group, oral administration of castor oil led to a significant increase (p < 0.001 vs. the negative control) in fluid accumulation (136 ± 1.9 g) in the intestinal fluid accumulation assay (112 ± 1.8 g). When mice were pretreated with ethanol and aqueous C. sativus seed extracts, the intestinal secretions were decreased considerably (p < 0.001 vs. castor oil group). At 150 mg/kg, the weights of fluid for ethanol and aqueous extracts were 96.6 ± 1.15 and 101.3 ± 2.06 g, respectively, whereas, at 300 mg/kg, the weights were 83.7 ± 1.2 and 88.1 ± 0.93 g, respectively, which were similar to the positive control groups treated with loperamide (75.4 ± 1.2 g) and verapamil (78.0 ± 0.13 g). In contrast, the n-hexane and DCM extracts failed to cause an antisecretory effect. At a dosage of 150 mg/kg, the n-hexane and DCM extracts raised fluid accumulation to 109 ± 2.6 g and 105 ± 2.1 g, respectively, whereas at 300 mg/kg, these extracts increased fluid accumulation to 122 ± 1.7 g and 118 ± 2.3 g, respectively ( Figure 10).

Discussion
Herbal therapy has been globally recognized as a safe and effective approach for treating a wide range of chronic disorders. Consumption of non-nutrient secondary metabolites [37] from a high vegetable, fruit, and grain diet can prevent chronic illnesses [38][39][40]. Within this context, local healers and practitioners have long used Cucumis sativus L. seeds to treat various diseases [11,16,25,41]. The present study was designed to explore the possible therapeutic potential of C. sativus seed extracts in digestive and respiratory ail-ments. Our findings related to phytochemical screening revealed the presence of sterols, anthraquinones (anthracene), alkaloids, flavones, and flavonoids in the seed extracts of C. sativus [11,15,38,42]. In addition, published research indicated that C. sativus seeds are rich in polyphenols such as quercitrin, rutin, and orientin [27], along with phytosterols (β-sitosterol, stigmasterol) [25]. In this respect, β sitosterol and polyphenols such as rutin and orientin exert antispasmodic activity. These phytochemicals in plants, fruits, or seeds may be antispasmodic [43,44].
In silico analyses of bioactive chemicals isolated from C. sativus seed extracts anticipated that these compounds would have antispasmodic action via interfering with calcium-mediated signaling and cholinergic synapses (Figure 1). The GO and KEGG pathways investigation revealed that the target genes of bioactive substances have a direct or indirect impact on smooth muscle contraction via calcium-mediated and cholinergic synapse signaling. The bioactive chemicals in the C. sativus seeds extracts interact with both pathway members In molecular docking studies, these compounds were shown to have a strong affinity for the L type voltage-gated calcium ion channel and phosphoinositide phospholipase C, myosin light chain kinase (MLCK), and muscarinic M3 receptor. Thus, it can be assumed that the strong antispasmodic activity of C. sativus seed extracts was due to the high affinity of compounds for target proteins, which may inhibit signal transduction responsible for contraction.
The gastrointestinal motor tone is regulated by several physiological mediators [45]. These physiological mediators have a role in the gut's stimulatory actions. In this respect, an increase in the cytosolic calcium concentration, caused by either extracellular calcium influx or cytosolic calcium release from endoplasmic reticulum storage, invokes depolarization and disrupts the resting membrane's action potential and spontaneous contractions of jejunum preparations. Any medication that inhibits this pathway is a recognized antispasmodic [18,46] and an active therapeutic agent in hyperactive gastrointestinal illnesses [32,45,47]. Similarly, the antidiarrheal and antispasmodic effects of C. sativus seed extracts were tested against spontaneous contractions of jejunal preparations, and the findings of the polar and nonpolar extracts were opposite. The ethanol and aqueous C. sativus seed extracts reduced involuntary contractions by preventing calcium entry into the cells. However, n-hexane and DCM extracts revealed a contractile response at low dosages and a spasmolytic reaction at higher concentrations. As previously stated, physiological mediators elicit contractile response binding with their receptors to speed up cytosolic calcium influx and trigger depolarization. The spasmogenic response of n-hexane and DCM extracts was reduced by atropine (1 µM), demonstrating the existence of cholinergic components in both extracts [48].
C. sativus seed extracts were evaluated for their spasmolytic activity in jejunum preparations to see whether it was mediated by calcium channel blockage or potassium channel opening [49]. Results indicated that K + (80 mM) promotes depolarization in the cell through calcium influx, thus increasing the calcium ions. In jejunum preparations, this intense depolarization of the membrane action potential induced a long-lasting sustained contractile response [50]; the smooth muscle relaxes after repolarization. In this investigation, the ethanol and aqueous C. sativus seed extracts relaxed K + (80 mM), elicited contractions, inhibited the provoked depolarization, and repolarized the membrane potential by inhibiting the calcium current inflow. In addition, the following series of events were inhibited [45,51]: (1) cytosolic calcium ions concentration would be decreased; (2) due to a lack of or the inadequate number of cytosolic calcium ions, calcium ions did not interact with regulatory protein phosphokinase C, and the calcium-calmodulin complex would not be formed; (3) as a result, MLCK activation was blocked, and MLC phosphorylation was attenuated; and (4) actin and the MLCs did not interact because the MLCs were not phosphorylated, hence no contractile response was produced in tissue [32].
At a concentration of 10 mg/mL, the DCM extract relaxed, while the n-hexane exhibited partial relaxation in response to contractions caused by K + (80 mM). Similar to verapamil, pretreatment with ethanol and aqueous C. sativus seed extracts resulted in the suppression and rightward shift of calcium ion channel CRCs on jejunum preparations [52], which might support previous findings that ethanol and aqueous C. sativus seed extracts acted as calcium ion channel blockers and repolarized the membrane action potential [50]. On the other hand, the n-hexane and DCM extracts failed to diminish the calcium CRCs, although both extracts showed a shift to the right of the CRCs with the reduction in atropine pretreated tissues. Findings showed that the n-hexane and DCM extracts contain cholinomimetic and CCBs compounds [36]. Cholinomimetic constituents stimulate muscarinic receptors in the absence of atropine, counteracting calcium CRCs suppression; however, in the pretreatment of atropine, muscarinic receptors were occupied and rendering the cholinergic compounds of both n-hexane and DCM extracts to be impotent [35,36]. Although drugs that relax the K + (25 mM) provoke contractions and are referred to as K + channel openers, CCBs have an equal affinity to inhibit both K + (80 mM) and K + (25 mM) provoked contractions [46]. Thus, C. sativus seed extracts and verapamil inhibited K+ (80 mM) entirely, while K+ (25 mM) induced contractions [53]. The spasmogenic effect of n-hexane and DCM extracts was investigated using rat ileum, which is advantageous for assessing spasmogenic response due to its quiescent nature [48]. When n-hexane or DCM extract was administered in bolus dosages, it induced a dose-dependent contractile response similar to acetylcholine. Atropine (1 µM) inhibited the contractile response of both n-hexane and DCM extracts, suggesting the presence of cholinomimetic drugs capable of initiating the activation of M3 muscarinic receptors [36].
C. sativus seed extracts were tested for antispasmodic activity by depolarizing tracheal and urinary bladder preparations with persistent contractions of CCh (1 µM) and K + (80 mM) [54]. By establishing a repolarization state, C. sativus seed extracts relaxed the CCh (1 µM) induced contractions. However, C. sativus seed extracts demonstrated the same findings for K + (80 mM)-induced contractions on the tracheal and urinary bladder, as previously found in jejunal preparations, that ethanol and aqueous extracts relaxed the contractions. However, the DCM extract relaxed the contractions at 10 mg/mL, whereas the n-hexane extract could not attain complete relaxation. Based on the EC 50 value of CCh (1 µM), C. sativus seed extracts act on the muscarinic receptors in tracheal and urinary bladder preparations. A rightward shift in carbachol CRCs on tracheal preparations similar to verapamil revealed the anticholinergic effect of C. sativus seed extracts [55]. The activation of the M3 muscarinic receptor stimulates the phosphoinositide phospholipase C, which hydrolyzes phosphatidylinositol 4,5-bisphosphate into two secondary messengers; inositol 1,4,5-trisphosphate (IP 3 ) and diacylglycerol (DAG). IP 3 stimulates inositol 1,4,5trisphosphate receptors (IP3R) on the sarcoplasmic reticulum to release calcium ions, increasing the cytosolic calcium ion levels. Diacylglycerol and calcium ions activate a regulatory protein phosphokinase C (PKC), through which phosphorylation of calmodulin occurs to form a calcium/calmodulin complex. This calcium/calmodulin complex activates MLCK, which causes the phosphorylation of MLCs. The phosphorylated MLCs and actin form an interaction network to produce a contractile response [50,56]. The anticholinergic activity of C. sativus seed extracts was confirmed by a rightward shift in the carbachol CRCs on tracheal preparations such as verapamil [55]. Consequently, C. sativus extracts behaved as bronchodilators to reduce cytosolic calcium release from the sarcoplasmic reticulum, by inhibiting the muscarinic receptor path of contractile response.
On the other hand, the ethanol and aqueous extracts significantly inhibited peristalsis, diarrhea, and electrolyte imbalance (p < 0.05 vs. control), whereas the n-hexane and DCM analogs failed to suppress them. Moreover, the ethanol and aqueous seed extracts produced an antispasmodic response, thus preventing the charcoal meal from traveling the distance and blocking the gut's rhythmic contractions. Castor oil hydrolysis in the gut releases an active ingredient called ricinoleic acid [57,58], which may raise intestinal fluid content and induce electrolyte and water transport imbalance, resulting in massive contractions in the transvers distal colon. Furthermore, the ethanol and aqueous extracts inhibited diarrhea and fluid secretion produced by castor oil in animals [57]. As with loperamide and verapamil, the antiperistalsis, antisecretory, and antidiarrheal effects of the ethanol and aqueous extracts were mediated by the calcium channel blocking activity [59,60]. In contrast, n-hexane and DCM seed extracts showed no substantial suppression of the peristaltic movements, castor oil-induced fluid secretion, and castor oil-induced diarrhea. These findings suggest the existence of certain constituents that prevent the n-hexane and DCM extracts from inhibiting spasms.

Preparation of Extracts
The seeds of C. sativus were collected from fresh fruits grown locally in Multan, Punjab, Pakistan, from May to July 2018; these seeds were then depulped to remove the fruit pulp. Dr. Zafarullah Zafar, a taxonomist at Bahauddin Zakariya University Multan, Pakistan, identified and authenticated the seeds and fruit of C. sativus. A voucher specimen (PKF 1916/30) was deposited at the herbarium located at the same institute. After removing the husks, the inner part of the seeds was ground with an herbal grinder into a granular powder. Approximately 120 g of the pulverized material was subjected to successive extraction (fractions) in the Soxhlet apparatus with solvents of varying polarity. Solvents were evaporated at a reduced pressure using rotary evaporation to yield oils (n-hexane and DCM) and thick brownish yellowish tinted residues (ethanol and aqueous) of C. sativus seed extracts. The overall yield of C. sativus seeds extracts was 53.9%, whereas the n-hexane, DCM, ethanol, and aqueous yields were 28 g, 16.1 g, 12 g, and 8.6 g, respectively.
The extracts of C. sativus seed were placed in an amber-colored container and kept at a temperature of −20 ± 3 • C. The n-hexane and DCM extracts were freshly prepared in 1.5 percent v/v Tween 80 on the experiment day, whereas the ethanol or aqueous extract was dissolved in distilled water or 0.9% normal saline solution. For Tween 80, we first dissolved both n-hexane and DCM extracts in respective solvents n-hexane and dichloromethane, then added the 1.5 percent v/v Tween 80 and thoroughly mixed it on the ultra-sonicator. Then we evaporated the respective solvent at 30 ± 2 • C and then dissolved it in distilled water or 0.9 percent normal saline solution. All preparations were diluted to the desired concentration.

Sample Preparation for HPLC and LC-ESI-MS/MS
The sequential polarity-based C. sativus seed extracts were diluted in 1 mL of 100% methanol for HPLC and LC-ESI-MS/MS analysis, and the supernatant was collected after centrifugation for 12 min at 14,000 rpm. The supernatant was filtered using a 0.22 µm syringe filter.

LC ESI-MS/MS Analysis
LC ESI-MS/MS analyses were performed on C. sativus seed extracts to identify bioactive constituents [18,61]. We employed direct sample injections into an ESI probe configured for positive (M+) and negative (M-) ionization modes.
Binary Eluent system: Different solvent systems were employed for C. sativus seed extracts due to the different nature of extracts. The mobile phase for separating the constituents of the ethanol and aqueous extracts was solvent A, made up of 0.1 percent v/v formic acid (FA) in methanol, and solvent B, which was made up of 0.1 percent v/v FA in aqueous and ACN (80:20). Similarly, the mobile phase was used to separate the constituents of Cu.Hexane and Cu.DCM was solvent A, which comprised 0.1 percent v/v FA in methanol, and solvent B, which comprised 0.1 percent v/v FA in ACN. For solvent A (v/v), the elution gradient profile was modified as follows: 98%; 98-80%; 80-70%; 70-10%; and 10-98% from 0-4.0; 4.0-8.0; 8.0-15.0; 15.0-17.0; and 17.0-20.0 min, respectively [49].
Optimal chromatographic conditions: The optimal chromatographic parameters for separation of the constituents of C. sativus seed extracts were: the flow-rate of the mobile phase was 0.4 mL/min; the temperature of the column was 25 • C; the volume of sample-injection was 8 µL; the range for mass-spectrum (m/z) was 50 to 1000; and the size of the column (C-18) was 2.1 mm × 100 mm, 3.0 µm. The optimal mass-spectrometry conditions were: the positive and negative mode ion spray voltage was 5500 V and −5500 V, respectively; the pressure of curtain gas (CUR), nebulizer gas (GS1), auxiliary gas (GS2) were 25, 55, 55 psi, respectively; the ion source temperature for mass-spectrometry was 400 • C. Nitrogen was used as a CUR and GS2. A collision-induced dissociation (CID) energy range of 10 to 45 was used for ion fragmentation based on parent ion stability. The ESI-MS/MS data were collected and processed using XCalibur 3.0.
Identification of Compounds: Bioactive constituents of C. sativus seed extracts were identified by comparing their mass spectral and chromatographic data to those previously published and reference libraries [49]. ChemDraw 18.0 was used to predict structure fragmentation.

RP-HPLC Method Optimization and Validation
We performed RP-HPLC experiments to establish the validity and quantification of bioactive compounds [18,49].
Binary Optimal chromatographic conditions: The optimal chromatographic parameters for separation of constituents of C. sativus seed extracts were: the flow-rate of the mobile phase was 0.8 mL/min; the temperature of the column was 25 • C; the volume of sample-injection was 8 µL; the size of the column (C-18) was 4.6 mm × 150 mm, 5.0 µm; and the appropriate wavelengths for peak detection were 250; 270; 280; and 320 nm. The results of the C. sativus seed extracts were compared with the peaks and retention time (Rt) of external standards.

Validation of an Analytical Method
A validation method was developed under the International Conference on Harmonization (ICH) principles, as outlined in our prior work [18].
Preparation of external standards: A total of 150 µg of external standards (umbelliferone, stigmasterol, campesterol, β-sitosterol, epicatechin, scopoletin, ellagic acid, luteolin, apigenin, orientin, 1,4-dicaffeoylquinic acid, narcissin, quercetin, naringenin, luteolin 7-O-glucoside, kaempferol, ferulic acid, p-coumaric acid, kaempferol-3-O-glucoside, and hesperidin) were dissolved in 1 mL of 100 percent methanol to prepare the stock solution of external standards, then centrifuged at 14,000 rpm for 10 min and the supernatant was collected and stored at −20 • C. On the experiment day, the stock solution of external Linearity, limits of detection, and quantification: The dilutions (7.81-500 µg/mL) of each external standard were prepared for linearity validation of HPLC, and a linear regression standard was developed for each external standard, based on the concentration and peak area of each standard. The detection (LOD) and quantification (LOQ) limits were calculated by using the following formula: where S is the slope of the linear regression equation, and σ is the intercept of the standard deviation. Specificity: An injection volume of 8 µL of C. sativus seed extracts and external standards samples were injected at a flow rate of 0.8 mL/min into RP-HPLC using the above ideal chromatographic settings. To measure the peak purity of the C. sativus seed extracts, the peaks and retention time (Rt) of the C. sativus seed extracts were compared with external standards. The peaks were quantified using a linear regression standard for each external standard.
Precision and repeatability: Under the above-mentioned appropriate chromatographic parameters, the precision and repeatability of this method was studied for each external standard on inter-day (n = 4) and intra-day (n= 4 for 3 days), and the relative standard deviation (RSD) was determined.
Accuracy and recovery: The external standards were recovered from the C. sativus seed extracts to assess the accuracy and recovery of this method. Three different concentrations (50, 100, and 150 µg/mL) of external standards were used for this approach. Furthermore, 100 mg of C. sativus seed extract and 50, 100, and 150 µg of the external standards were mixed thoroughly in 0.5 mL of methanol and then centrifuged at 14,000 rpm for 10 min, and the supernatant was collected, and the volume was made up to 1 mL with methanol, to get a final concentration of external standards of 50, 100, and 150 µg/mL. An injection volume of 0.8 µL of this mixture was injected into the RP-HPLC utilizing the optimal chromatographic parameters, as described above. The sample recovery percentage of an external standard was calculated using the following equation % Recovery = (observed concentration ÷ actual concentration) × 100

Animals and Housing Conditions
Sprague Dawley rats (170 ± 20 g), mice (28 ± 3 g), and locally bred rabbits (1.2-1.4 kg) of either sex were used throughout this investigation. The temperature of these animals was maintained at 23-25 • C in dark/light cycled animal houses, with free access to a standard diet and water. They were maintained according to the Animals By-Laws, which were authorized by the Ethical Committee of the Department of Pharmacology (vide No. EC /04-PhDL/S2018), and the research was carried out in compliance with the Commission of Laboratory Animal Resources [62]. Rabbits were sacrificed for in vitro experiments, whereas rats and mice were sacrificed through cervical dislocation for in vivo experiments.

Network Pharmacology Analysis
We adopted the approach provided by Xiao et al. [63] for network pharmacology. Screening of potential protein targets: The potential protein targets were retrieved from DrugBank (http://drugbank.com, accessed on 10 July 2021) and Swiss target prediction (http://www/swisstargetpredication.ch, accessed on 10 July 2021) for the bioactive constituents of C. sativus seed extracts. The potential protein targets for diseases were obtained from GeneCards (https://www.genecards.org, accessed on 16 July 2021), Dis-GeNET (http://www.disgenet.org/web/DisGeNET, accessed on 16 July 2021), Pubmed (https://pubmed.ncbi.nlm.nih.gov, accessed on 16 July 2021), and Online Mendelian Inher-itance in Man (OMIM) (https://www.omim.org, accessed on 16 July 2021) with keywords "diarrhea"; "chronic diarrhea"; "chronic constipation"; "constipation"; "irritable bowel syndrome (IBS)"; "asthma"; "obstructive pulmonary disorder "; and "coughing." The VarElect tool (https://ve.genecards.org/, accessed on 16 July 2021) was used to select the top 150 disease target genes for gastrointestinal and tracheal disorders under clinical phenotype and genetic association. Using the package "ggvenn" in RStudio, we conducted an intersection evaluation between the bioactive chemicals of C. sativus seed extracts and disease target genes, and a Venn Diagram was constructed to illustrate the relationship between the potential protein targets of bioactive chemicals and disease. Further, we evaluated the KEGG pathway and gene ontology (GO) enrichment for intersected potential protein targets.
Network construction and pathway analysis: We used RStudio with packages "BiocManager", "Clusterprofiler", and "DOSE" and utilized the database "org.Hs.eg.db" for KEGG pathway and GO enrichment analysis of potential target genes. The compound target disease proteins (C-T-D) network, the Protein-Protein Interaction (PPI) network, and the Compounds-Targets-Pathways (C-T-P) network were all constructed and analyzed in Cytoscape 3.8.0.

Protein Homology Modeling
Homology modeling identifies the three-dimensional structure of a protein, based on UniProt sequence similarities to one or more PDB templates. This strategy is consistent with the notion that a target structural conformation is more conserved than its amino acid sequence. The ExPASy ProtParam tool (http://www.web.expasy.org/protparam, accessed on 23 July 2021) was used to identify different physical and chemical qualities [63], such as molecular weight, amino acid composition, theoretical isoelectric point, extinction coefficient, the total number of negative and positive charged residues, instability index, aliphatic index, and the grand average of hydropathicity (GRAVY).
The following stages are involved in homology modeling (Maestro v11.8, Schrodinger suite 2018-4): template identification from UniProt sequence; BLAST homology search (NCBI PDB Database); sequence alignment (ClustW); model construction; and then validation of the model. The protein model failed to pass the Protein Reliability report and the Ramachandran plot. We proceeded further for the refinement of the build model, which included loops refinement and energy reduction using the Prime Refine Loops module and the Prime minimize module, choosing OPLS3e force field and VSGB solvation model. The minimized protein model was then evaluated in PROCHECK for the Ramachandran plot on the SAVES version 6 server (https://saves.mbi.ucla.edu/, accessed on 23 July 2021) and the PDBSum server (http://www.ebi.ac.uk/thornton--srv/databases/cgi--bin/pdbsum/, accessed on 23 July 2021) [64,65].

Molecular Docking
We employed the method described by Sirous and colleagues [33] for the molecular docking studies and used different modules of Maestro v11.8, Schrodinger suite 2018-4, and Discovery Studio. The binding forces and stability of the target proteins were investigated using prominent ligands or active compounds in network pharmacology. The LigPrep module optimized, minimized, and ionized the ligands. PubChem (https://pubchem.ncbi. nlm.nih.gov, accessed on 19 July 2021) was used to obtain three-dimensional structures of ligands in mol 2 format. In addition, the force field of OPLS3e was utilized to optimize and minimize the energy of ligands to increase the precision and effectiveness of the proteinligand model with the lowest energy conformer. The Epik tool was utilized for putative promoters, and ionization states were created at the desired pH (7.0 ± 2.0). The tautomeric states were created for each ionized ligand, and the ligands with the most atoms were kept. The stereoizer tool was used to generate stereoisomers of ligands, while retaining data on stereogenic chirality from the input file.
In addition, the protein preparation module was used to assign hydrogen bonds, minimize het states, ionize, and optimize the three-dimensional crystal structures of the built protein model. The built protein model structure was modified to eliminate duplicate binding sites and determine if the protein-ligand interaction is a dimer or a multimer. The chemical component dictionary (CCD) database was utilized to reorder the het groups and non-templated amino acid residues. The formal charges assigned to metal ions and surrounding atoms had been altered, and covalent bonds between metal ions and neighboring atoms had been reduced to zero-order. If two sulfur atoms are within 3.2 Å of each other in the protein structure, disulfide bonds are formed between them. Additionally, the prime tool had been used to validate and fill in the structural gaps in the protein, such as those caused by missing side-chain atoms. Finally, the Epik tool was used to produce the protonation and metal charge states of the het groups for cofactors and metals with a pH of (7.0 ± 2.0) within a range of 5.0 Å. The hydrogen atoms were substituted for water molecules and het groups in the protein structure within a range of 5.0 Å. The PROPKA tool had been used to optimize the assignment of hydrogen bonds in protein structures at a pH of 7.0 and the OPLS3e force field for constraint energy minimization and protein structural optimization.
The binding pockets were identified for molecular docking and receptor grid generation modules. In this scenario, a literature search and sitemap module were used to identify the coordinates of the cubic grid box to facilitate the ligand binding in the protein pockets. The extra precision Glide was utilized to perform molecular docking with the previously prepared compounds, proteins, and receptor grid files. For the Van der Waals radii, the partial charge was adjusted to 0.15 with a scaling factor of 0.80. The ligand sampling was optimized for flexible docking and includes nitrogen (non-ring) inversions and ring conformations to generate conformers. Additionally, the torsions around the bond were sampled for bias. With the Epik tool, penalties were applied to the docking score. The ligand-protein complex binding energies were computed using the Prime MM-GBSA module with VSGB solvation and the force field OPLS3e.
Inhibition Constant (Ki): The inhibition constant was derived using the ligand's binding free energy, using the following equation: where ∆G is the ligand's binding free energy, R is the gas constant (cal mol −1 K −1 ), and T is the room temperature (298 Kelvin).

Isolated Tissue Experimentation
We have used the methods of Saqib and Janbaz [50], and Wahid et al., [18] for the in vitro experiments.

Isolated Rabbit Jejunum Preparations
Animals were sacrificed and dissected to open the abdomen. Tissues from the adjacent mesenteries were removed, and sections of jejunum tissue about 2 to 3 cm in length were produced. To keep the jejunum preparation alive and ready for use, they were removed and stored in a freshly prepared Tyrode's solution, aerated with carbogen gas (95% O 2 and 5% CO 2 ). A tissue organ bath of 15 mL was used to fix each jejunum preparation filled with Tyrode's solution, aerated with carbogen gas, and maintained at 37 ± 0.5 • C with a circulating thermoregulator. Before assessing the drug, a force of 1 ± 0.1 g was preloaded onto the jejunal preparation. Spontaneous rhythmic contractions were observed after allowing the jejunal preparation to be equilibrated for 25 ± 5 min with Tyrode's solution changes every 8 min. The tissue physiological responses were recorded using an isotonic transducer (MLT0015), which were amplified by a data-collecting device, Power Lab ® (4/25), and shown in Lab Chart Pro (Version 7). The physiological response was determined using the contraction percentage obtained immediately before the addition of the test substance, to test whether n-hexane and DCM, ethanol, and aqueous C. sativus seed extracts had spasmolytic and spasmogenic effects. A tissue organ bath materialized the jejunal preparation [55].
The antispasmodic activity of C. sativus seed extracts was assessed by either inducing smooth muscle contractions through a calcium ion channel blockade with K + (80 mM) or potassium ion channel activation with K + (25 mM). The extracts were applied cumulatively to get a dose-dependent inhibitory response of C. sativus seed extracts. K + (80 mM) may affect the smooth muscle contraction by depolarizing the cell and causing an influx of calcium ions. Due to the extract's ability to inhibit contractions induced by K + (80 mM), it was designated a calcium channel blocker [32].
Further, the C. sativus seed extracts were tested to see whether they have a calcium antagonistic effect. To diminish intracellular calcium concentrations, the jejunal preparation was incubated with K + (80 mM) three times and then equilibrated in an EDTA containing calcium-free Tyrode solution for 25 ± 5 min, followed by 50 ± 5 min incubation in calciumfree and potassium-rich Tyrode's solution. After incubation, calcium was added to the jejunum preparation in a cumulative fashion to produce calcium concentration-response curves (CRCs). When calcium CRCs were made as a control, jejunal preparations were washed and incubated with C. sativus seed extracts for 50 ± 5 min. These calcium CRCs were made with various doses of C. sativus seed extracts and compared to control calcium CRCs, to examine whether any calcium antagonistic activity was present [32].

Isolated Rat Ileum Preparation
The mesenteries encircling the ileum tissue of rats were removed, and 2 to 3 cm long ileum tissue segments were obtained. Each ileum preparation was processed in the same manner as previously described in the jejunum preparation section. The ileum preparation was incubated in a tissue organ bath to determine the cholinergic response of n-hexane and DCM seed extracts [43]. Bolus administrations of n-hexane and DCM extracts were used to assess the contractile response.

Isolated Rabbit Tracheal Preparations
A rabbit tracheal tissue was dissected, and 3-4 mm length rings were formed by excising adjacent fatty and adhesion materials. Smooth muscles were positioned between cartilaginous segments after laterally incising the tracheal rings to produce a tracheal strip. A tissue organ bath of 15 mL was used to fix each tracheal preparation, filled with Krebs solution, aerated with carbogen gas, and maintained at 37 ± 0.5 • C with a circulating thermoregulator. Before assessing the drug, a force of 1 ± 0.1 g was preloaded onto the tracheal preparation after allowing the tracheal preparation to be equilibrated for 50 ± 5 min, with Krebs solution changes every 8 min. The physiological responses of the tissue were recorded using an isometric transducer (FORT100), which were amplified by a data-collecting device, Power Lab ® (4/25), and shown in Lab Chart Pro (Version 7). The physiological response was determined using the contraction percentage obtained immediately before the addition of the test substance. To test whether n-hexane and DCM, ethanol, and aqueous C. sativus seed extracts had bronchodilation activity [55], they were exposed to persistent contractions of carbachol (CCh 1 µM), K + (80mM), and K + (25mM).
Further, C. sativus seed extracts were tested to see whether they have a muscarinic receptor antagonism effect. The tracheal preparation was incubated with CCh (1µM) three times and followed by 50 ± 5 min of incubation in Krebs solution. After incubation, carbachol was added to the tracheal preparation in a cumulative fashion to produce carbachol (CCh) CRCs. When the CCh CRCs were made as a control, tracheal preparations were washed and incubated with C. sativus seed extracts for 50 ± 5 min. These CCh CRCs were made with various doses of C. sativus seed extracts and compared to the control CCh CRCs to examine whether any muscarinic receptor antagonistic activity was present [32].

Isolated Urinary Bladder Preparations
The urinary bladder tissue of a rabbit was dissected, and the rings of the urinary bladder were incised laterally to form a urinary bladder strip. A tissue organ bath of 15 mL was used to fix each tracheal preparation, filled with Krebs solution, aerated with carbogen gas, and maintained at 37 ± 0.5 • C with a circulating thermoregulator. Before assessing the drug, a force of 1 ± 0.1 g was preloaded onto the urinary bladder preparation after allowing the urinary bladder preparation to equilibrate for 50 ± 5 min with Krebs solution changes every 8 min. The physiological responses of tissue were recorded using an isometric transducer (FORT100), which were amplified by a data-collecting device, Power Lab ® (4/25), and shown in Lab Chart Pro (Version 7). The physiological response was determined using the contraction percentage obtained immediately before the addition of the test substance. To test whether n-hexane and DCM, ethanol, and aqueous C. sativus seed extracts had antispasmodic activity [55], they were exposed to persistent contractions of carbachol (CCh 1 µM), K + (80 mM), and K + (25 mM).
Further, C. sativus seed extracts were tested to see whether they have a calcium antagonistic effect. To diminish intracellular calcium concentrations, the urinary bladder preparation was incubated with K + (80 mM) three times and then equilibrated in an EDTA containing calcium-free Krebs solution for 25 ± 5 min, followed by 50 ± 5 min incubation in calcium-free and potassium-rich Krebs solution. After incubation, calcium was added to the jejunum preparation in a cumulative fashion to produce calcium CRCs. When calcium CRCs were made as a control, urinary bladder preparations were washed and incubated with C. sativus seed extracts for 50 ± 5 min. These calcium CRCs were made with various doses of C. sativus seed extracts and compared to control calcium CRCs to examine whether any calcium antagonistic activity was present [32].

Evaluation of Maximum Tolerated Dose
The maximum tolerated dose of C. sativus seed extracts was determined by randomly dividing the rats into six groups (five rats of either gender per group). For 28 days, one group received saline, while the other groups received 50, 100, 150, 200, and 300 mg/kg of C. sativus seed extracts every day. The extract was injected into the stomach via the esophagus using a stainless-steel needle linked to a plastic syringe. Bodyweight, behavioral and mood changes, clinical symptoms of pain or discomfort, and animal mortality were all tracked for 28 days [66].

Charcoal Meal GI Transit Test
This test was conducted using the method described by Saqib and Janbaz [50]. The antiperistalsis activity of C. sativus seed extracts was determined by randomly dividing the mice into 12 groups (five mice of either gender per group). Before the trial, all groups were allowed unrestricted access to water, but the food was removed for an overnight period.
Group I was marked as a negative control, orally given 10 mL/kg of 0.9% normal saline; Group II was marked as a positive control, orally given normal saline (10 mL/kg, p.o.), then after 15 min, a charcoal meal (10 mL/kg) was given orally; Group III, as standard control, was orally given loperamide (10 mg/kg); Group IV, as standard control, was orally given verapamil (10 mg/kg); Groups V and VI were orally given 150 and 300 mg/kg of n-hexane C. sativus seed extract, respectively; Groups VII and VIII were orally given 150 and 300 mg/kg of DCM C. sativus seed extract, respectively; Groups IX and X were orally given 150 and 300 mg/kg of ethanol C. sativus seed extract, respectively; Groups XI and XII were orally given 150 and 300 mg/kg of aqueous C. sativus seed extract, respectively.
After 30 min of treatment, all groups except Group I received an oral dose of castor oil (10 mL/kg). After 30 min, mice were sacrificed by cervical dislocation, their abdomens were dissected, and both ends of the small intestine pylorus and caecum were ligated, then the whole small intestine was removed so that the weight of the small intestine with and without intestinal fluid could be calculated. The findings were determined using the formula (Pi/Pm) × 1000, where Pm represents the mouse's body weight, and Pi represents the intestine's weight. The weight of fluid collected in the intestine was used to calculate the amount of fluid accumulated (g).

Statistical Analysis
The results are expressed as the mean ± standard deviation (SD), and 95% confidence intervals (CI) for the half-maximal effective concentration (EC 50 ) values were calculated. Using logarithmic sigmoidal dose-response graphs, concentration-response curves were plotted. Data were subjected to one-way analysis of variance (ANOVA). Statistical analysis was performed with the aid of the Dunnett test for significance with the aid of GraphPad Prism-8 (GraphPad Software, San Diego, CA, USA); differences were considered significant at p ≤ 0.05.

Conclusions
In conclusion, Cucumis sativus L. seed extracts, mainly ethanol and aqueous, exhibited medicinal effects that could be used to treat disorders such as asthma and diarrhea by regulating the contractile response via the L voltage-gated calcium ion channel, M3 muscarinic receptor, or other proteins, such as PLC, to the repolarize membrane action potential. Our results demonstrated that the dichloromethane and n-hexane extracts elicited a spasmogenic reaction, with higher dosages eliciting a spasmolytic response. This might be due to the cholinomimetic phytoconstituents in jejunum preparation, which cause the spasmogenic reaction. In vitro and in silico studies of C. sativus seed extracts validate these suppression phenomena of phytoconstituents in the same extract. Finally, C. sativus seed extracts include several bioactive chemicals that might be utilized to cure diseases and illnesses. Furthermore, our results may explain why C. sativus is used as an antidiarrheal medicinal plant, which explains its usage in traditional medicine to treat diarrhea, dysuria, and urge incontinence.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/ph15050641/s1, Table S1: Precision validation of analytical method of C. sativus sequential extracts; Table S2: Accuracy validation of analytical method of C. sativus seed sequential extracts through % recovery method; Table S3: Gastrointestinal and respiratory pathogenic genes of sequential seed extracts of C. sativus; Table S4: Gastrointestinal and respiratory pathogenic genes retrieved for compounds; Table S5: Top 15 GO Biological process of bioactive compounds of sequential extracts of C. sativus for gastrointestinal and respiratory target genes; Table S6: Top 15 KEGG pathways of bioactive compounds of sequential extracts of C. sativus for gastrointestinal and respiratory target genes; Table S7: A network analysis of pathogenic target genes interaction with phytoconstituents; Table S8: A network analysis of pathogenic target genes interaction with phytoconstituents and GO biological process (BP); Table S9: A network analysis of pathogenic target genes interaction with phytoconstituents and KEGG pathway. Table S10. Network analysis of pathogenic target genes interaction with phytoconstituents and GO biological process (BP). Table S11. Network analysis of pathogenic target genes interaction with phytoconstituents and KEGG pathway. Figure S1. Negative and positive ESI-MS/MS spectra for putative chemicals of sequential seed extracts of C. sativus. Umbelliferone, stigmasterol, campesterol, β-sitosterol, epicatechin, scopoletin, ellagic acid, luteolin, orientin, 1,4-dicaffeoylquinic acid, narcissin, quercetin, kaempferol, apigenin, naringenin, p-coumaric acid, luteolin 7-O-glucoside, kaempferol 3-O-glucoside, ferulic acid, hesperidin, and linoleic acid. Figure S2. RP-HPLC chromatograms of sequential C. sativus seed extracts at different wavelengths and chromatograms of bioactive components confirmed by standard addition. The HPLC chromatograms of ethanol extract (A) revealed epicatechin, scopoletin, ellagic acid, luteolin, and orientin at 280 nm. (B) showed 1,4-dicaffeoylquinic acid, narcissin, quercetin, apigenin, and naringenin at 320 nm. The HPLC chromatograms of DCM extract (C) revealed umbelliferone, stigmasterol, campesterol, and β-sitosterol at 270 nm. The HPLC chromatograms of aqueous extract (D)showed kaempferol-3-O-glucoside and hesperidin at 320 nm.