Alkaloid-Based Isoxazolylureas: Synthesis and Effect in Combination with Anticancer Drugs on C6 Rat Glioma Model Cells

Alkaloid-based urea derivatives were produced with high yield through the reaction of anabasine and cytisine with isoxazolylphenylcarbamates in boiling benzene. Their antitumor activity, in combination with the commonly used five anticancer drugs, namely cyclophosphane, fluorouracil, etoposide, cisplatin, ribomustine with different mechanisms of action, was investigated. Based on the quantum chemical calculations data and molecular docking, hypotheses have been put forward to explain their mutual influence when affecting C6 rat glioma model cells.


Introduction
Cancer is one of the leading causes of death worldwide.Chemotherapy has been widely used as the primary therapy, particularly against inoperable cancer or as an adjunct therapy before or after another treatment.However, it is usually accompanied by negative side effects.The combination of classical anticancer drugs with synergistic adjuvants makes it possible to reduce the therapeutic dose and thereby reduce the toxicity of chemotherapy.Therefore, the approach of using synergistic adjuvants is being seriously considered.This strategy is based on the susceptibility of different molecular pathways involved in the genesis of a particular disease to the different mechanism of action of each individual drug to improve treatment efficacy and reduce the development of drug resistance [1,2].
The main goal when developing a drug combination is to achieve an effect significantly greater than the additive effect of the individual drugs, thereby reducing dosage and toxicity.On the other hand, many types of complementary medicines may be ineffective.Worse, they can produce unintended adverse side effects or alter the actions of conventional medical treatments by inhibiting or amplifying their effects to dangerous levels [3,4].Therefore, the search for new effective chemotherapeutic drug combinations, and the study of the mechanisms underlying their mutually influencing processes, seem to be an urgent task of medicinal chemistry.
The isoxazole ring is often used as a central core of compounds with anticancer effects acting as inhibitors of numerous targets including aromatase, tyrosine kinase, thymidylate, ERα, PLA2, HDAC, HER2, HSP90 inhibition and apoptosis induction [5][6][7].The structural features of isoxazole make it possible for multiple non-covalent interactions, especially hydrogen bonds (hydrogen bond acceptor N and O), π-π stacking (unsaturated five-membered ring) and hydrophilic interactions [8].The authors of this article have been conducting research on molecular design of structures, synthesis and biotesting of new synergistic adjuvants in compositions with first-line anticancer drugs for a number of years.As a result, it was found that some isoxazole derivatives are promising candidates for a synergist's role [9][10][11].Previously conducted quantum chemical calculations of structural and electronic changes occurring in the cisplatin-adjuvant system based on morpholinium and 4-methylpiperazinium 4,5-dichloroisothiazol-3-carboxylates allowed us to assume that the synergistic effect appeared due to the conjugation of cisplatin with adjuvant through the relocation of frontier molecular orbitals, as well as an increase in the conjugate's dipole moment, which led to a change in the interaction character with the cellular target and an increase in the bioactivity of the system [12].
Continuing our research in this area, we expanded the scope of such studies to include other chemotherapy drugs that are used in modern oncology, as well as new promising 1,2-azole derivatives possessing synergistic effects.

Chemistry
We previously discovered that some azolylcarbamides and carbamates can exhibit a synergistic effect in compositions with antitumor substances [13].In this work, to obtain 5-arylisoxazolyl ureas with alkaloid fragments, an approach we used in [11], was chosen.(Scheme 1).The starting 5-phenylisoxazolylcarbamate was synthesized through the reaction of 5-phenylisoxazole-3-carbonyl azide with phenol, as described by us before [11].Isoxazolylphenylcarbamates were further introduced in the reaction with cytisine and anabasine to obtain 1-3 in 87-94% yields.NMR, MS, IR spectra of synthesized ureas are available in the Supplementary Materials.
Molecules 2024, 29, x FOR PEER REVIEW 2 of 17 can produce unintended adverse side effects or alter the actions of conventional medical treatments by inhibiting or amplifying their effects to dangerous levels [3,4].Therefore, the search for new effective chemotherapeutic drug combinations, and the study of the mechanisms underlying their mutually influencing processes, seem to be an urgent task of medicinal chemistry.The isoxazole ring is often used as a central core of compounds with anticancer effects acting as inhibitors of numerous targets including aromatase, tyrosine kinase, thymidylate, ERα, PLA2, HDAC, HER2, HSP90 inhibition and apoptosis induction [5][6][7].The structural features of isoxazole make it possible for multiple non-covalent interactions, especially hydrogen bonds (hydrogen bond acceptor N and O), π-π stacking (unsaturated fivemembered ring) and hydrophilic interactions [8].The authors of this article have been conducting research on molecular design of structures, synthesis and biotesting of new synergistic adjuvants in compositions with first-line anticancer drugs for a number of years.As a result, it was found that some isoxazole derivatives are promising candidates for a synergist's role [9][10][11].Previously conducted quantum chemical calculations of structural and electronic changes occurring in the cisplatin-adjuvant system based on morpholinium and 4methylpiperazinium 4,5-dichloroisothiazol-3-carboxylates allowed us to assume that the synergistic effect appeared due to the conjugation of cisplatin with adjuvant through the relocation of frontier molecular orbitals, as well as an increase in the conjugate's dipole moment, which led to a change in the interaction character with the cellular target and an increase in the bioactivity of the system [12].
Continuing our research in this area, we expanded the scope of such studies to include other chemotherapy drugs that are used in modern oncology, as well as new promising 1,2azole derivatives possessing synergistic effects.

Chemistry
We previously discovered that some azolylcarbamides and carbamates can exhibit a synergistic effect in compositions with antitumor substances [13].In this work, to obtain 5arylisoxazolyl ureas with alkaloid fragments, an approach we used in [11], was chosen.(Scheme 1).The starting 5-phenylisoxazolylcarbamate was synthesized through the reaction of 5-phenylisoxazole-3-carbonyl azide with phenol, as described by us before [11].Isoxazolylphenylcarbamates were further introduced in the reaction with cytisine and anabasine to obtain 1-3 in 87-94% yields.NMR, MS, IR spectra of synthesized ureas are available in the supplementary Materials.

Evaluation of Antitumor Activity
Determination of anticancer activity of isoxazolylureas 1-3 drugs and their binary mixtures on rat glioma cell line was carried out using flow cytometry analysis.C6 cell line is considered to be a safe and popular glioma model in the literature, providing a good simulation of glioblastoma multiforme [14].
At the beginning, the antitumor effect of synthesized isoxazolylureas was studied.The compounds 1-3 are found to have a little effect on the viability of C6 glioma cells (Figure 1a).Suppression of cell proliferation was observed: the drugs reduced the number of cells in the samples by 19, 23 and 6% compared to the control (Figure 1b).The total number of live cells in the samples decreased by 19, 25 and 8%, respectively (Figure 1c).Thus, the anticancer effect of isoxazolylureas 1-3 is determined mainly by the antiproliferative action of the compounds.Alkaloids anabasine and cytisine in doses of 100-500 µM had a weak inhibitory effect on cells or have no effect.

Evaluation of Antitumor Activity
Determination of anticancer activity of isoxazolylureas 1-3 drugs and their binary mixtures on rat glioma cell line was carried out using flow cytometry analysis.C6 cell line is considered to be a safe and popular glioma model in the literature, providing a good simulation of glioblastoma multiforme [14].
At the beginning, the antitumor effect of synthesized isoxazolylureas was studied.The compounds 1-3 are found to have a little effect on the viability of C6 glioma cells (Figure 1a).Suppression of cell proliferation was observed: the drugs reduced the number of cells in the samples by 19, 23 and 6% compared to the control (Figure 1b).The total number of live cells in the samples decreased by 19, 25 and 8%, respectively (Figure 1c).Thus, the anticancer effect of isoxazolylureas 1-3 is determined mainly by the antiproliferative action of the compounds.Alkaloids anabasine and cytisine in doses of 100-500 µM had a weak inhibitory effect on cells or have no effect.
Hence, even at such a fairly high dose, compounds 1-3 exhibited weak cytotoxicity along with an antiproliferative effect on tumor cells.This makes them suitable candidates for the role of low-toxic synergists of anticancer drugs.A study of the synergistic effects of the five commonly used anticancer drugs, namely cyclophosphamide (CP), fluorouracil (FU), etoposide (ET), cisplatin (CPt), ribomustine (RM), with different mechanisms of action, was undertaken (Figure 2).
The antitumor activity of combinations at fixed drug and isoxazolylurea concentrations was assessed to identify the most promising binary mixture.We used drugs and isoxazolylureas at concentrations of 5 and 200 µM, respectively, and compared their combined effect to a tenfold increased drug dose (50 µM).Hence, even at such a fairly high dose, compounds 1-3 exhibited weak cytotoxicity along with an antiproliferative effect on tumor cells.This makes them suitable candidates for the role of low-toxic synergists of anticancer drugs.

Evaluation of Antitumor Activity
Determination of anticancer activity of isoxazolylureas 1-3 drugs and their binary mixtures on rat glioma cell line was carried out using flow cytometry analysis.C6 cell line is considered to be a safe and popular glioma model in the literature, providing a good simulation of glioblastoma multiforme [14].
At the beginning, the antitumor effect of synthesized isoxazolylureas was studied.The compounds 1-3 are found to have a little effect on the viability of C6 glioma cells (Figure 1a).Suppression of cell proliferation was observed: the drugs reduced the number of cells in the samples by 19, 23 and 6% compared to the control (Figure 1b).The total number of live cells in the samples decreased by 19, 25 and 8%, respectively (Figure 1c).Thus, the anticancer effect of isoxazolylureas 1-3 is determined mainly by the antiproliferative action of the compounds.Alkaloids anabasine and cytisine in doses of 100-500 µM had a weak inhibitory effect on cells or have no effect.
Hence, even at such a fairly high dose, compounds 1-3 exhibited weak cytotoxicity along with an antiproliferative effect on tumor cells.This makes them suitable candidates for the role of low-toxic synergists of anticancer drugs.A study of the synergistic effects of the five commonly used anticancer drugs, namely cyclophosphamide (CP), fluorouracil (FU), etoposide (ET), cisplatin (CPt), ribomustine (RM), with different mechanisms of action, was undertaken (Figure 2).
The antitumor activity of combinations at fixed drug and isoxazolylurea concentrations was assessed to identify the most promising binary mixture.We used drugs and isoxazolylureas at concentrations of 5 and 200 µM, respectively, and compared their combined effect to a tenfold increased drug dose (50 µM).The antitumor activity of combinations at fixed drug and isoxazolylurea concentrations was assessed to identify the most promising binary mixture.We used drugs and isoxazolylureas at concentrations of 5 and 200 µM, respectively, and compared their combined effect to a tenfold increased drug dose (50 µM).
The results of the study are summarized in the Table 1 and Figure 3.As can be seen from the experiment, compound 1 in combination with all drugs except ribomustine showed the effect of additive synergy, demonstrating an activity level approximately equal to that of the sum of each compound separately.Compositions with isoxazolylurea 2 showed a clear example of antagonism in combination with the drugs cyclophosmamide, fluorouracil and etoposide, inhibiting the pharmacological activity of each other or at least one of the compounds, and resulting in 2-3 fold decrease of antiprofilerative action compared to the sum of the effects of the two compounds.At the same time, the combination of compound 2 with cisplatin and ribomustine led to the opposite result: the synergistic effect was realized in a significant increase of the drug's action.The cytostatic activity of the binary mixture 1+RM (5 µM) exceeded the effect of the drug at a 10-fold increased dose (50 µM) by 22%.Generally, the most pronounced synergistic effect was observed in the case of ribomustine mixtures with all synthesized compounds.Σ is the sum of the effects of a single compound and drug.
The results of the study are summarized in the Table 1 and Figure 3.As can be seen from the experiment, compound 1 in combination with all drugs except ribomustine showed the effect of additive synergy, demonstrating an activity level approximately equal to that of the sum of each compound separately.Compositions with isoxazolylurea 2 showed a clear example of antagonism in combination with the drugs cyclophosmamide, fluorouracil and etoposide, inhibiting the pharmacological activity of each other or at least one of the compounds, and resulting in 2-3 fold decrease of antiprofilerative action compared to the sum of the effects of the two compounds.At the same time, the combination of compound 2 with cisplatin and ribomustine led to the opposite result: the synergistic effect was realized in a significant increase of the drug's action.The cytostatic activity of the binary mixture 1+RM (5 µM) exceeded the effect of the drug at a 10-fold increased dose (50 µM) by 22%.Generally, the most pronounced synergistic effect was observed in the case of ribomustine mixtures with all synthesized compounds.The anabasine derivative 3 showed a synergistic effect with fluorouracil, cisplatin and ribomustin increasing the effectiveness of the drug by 8-12% and distinct antagonism with cyclophosphamide reducing the drug's activity by 15%.In the combination with etoposide, additive action is noted.The anabasine derivative 3 showed a synergistic effect with fluorouracil, cisplatin and ribomustin increasing the effectiveness of the drug by 8-12% and distinct antagonism with cyclophosphamide reducing the drug's activity by 15%.In the combination with etoposide, additive action is noted.
Generally, a synergistic effect is observed to varying degrees for all ureas 1-3 when using binary mixtures with alkylating chemotherapy drugs cisplatin and ribomustine.Antagonism and the tendency to antagonism are noted for the remaining drugs, except for the combination of isoxazolylurea 3 with fluorouracil.
Antagonism is known to occur through various means.The most common type is receptor antagonism, where one drug acts as an antagonist to block the effects of another drug by binding to the same receptor.This prevents the agonist drug from activating the receptor and producing its intended effect.Some drugs follow chemical antagonism, which involves the coupling of two drugs to form an inactive product, etc. [15].
To interpret the observed experimental picture, quantum chemical modeling of binary drug-urea mixtures and molecular docking of isoxazolylureas was carried out, taking into account targets that determine the mechanism of drugs action.

Quantum Chemical Modeling
Based on the assumption that the effect of the combined action of drugs can be mediated by the formation of non-covalent drug-adjuvant conjugates, as was shown before [12], we determined structural and electronic changes that occur in the system of two conjugated molecules.DFT/B3LYP-D3/cc-pvdz/LanL2DZ(Pt) level of theory was used for the calculation of molecules' optimal geometry, dipole moment, localization and energy characteristics of frontier molecular orbitals (FMO).Calculations were carried out both for individual compounds and their conjugates, with consideration given to aqueous medium, which simulates situation in living cells.
The optimized geometry of the drug-urea conjugates is presented in Figure 4, and energy characteristics in Tables 1 and 2. The energy characteristics of conjugates formed from the initial components were calculated using formula (Table 2): It follows from calculation data that drug molecules (FU and CPt) and urea derivatives 1-3 form conjugates with shortened interatomic distances due to non-covalent interactions caused by hydrogen and van der Waals bonds (Figure 4a-f).As can be seen from Figure 4, the structure of isoxazolylurea 2 differs from the other two in its most linear spatial arrangement, which results in the varying degrees of conjugates' stability due to the participation of different atoms in intermolecular bond formation.Thus, in the case of fluorouracil, the most stable complex is observed with compound 2 with ∆E f = −16.88kcal/mol (Table 2).Hydrogen bonds are expressed between the oxygen and nitrogen atoms of the urea fragment, and between the nitrogen atom of isoxazole heterocycle and the N-H and C=O fragments of fluorouracil C=O urea  2).In contrast to conjugate 2+CPt (e), the formation of hydrogen bonds involved the nitrogen atom of the pyridine heterocycle of anabasine and the oxygen atom of the cytisine C=O group, in addition to the N, O atoms of the isoxazole heterocycle and the urea fragment.The above allows us to conclude that the formation of stable complexes leads to the manifestation of antagonistic effects in the drug-urea system.The energies of frontier molecular orbitals (FMO) are often considered as important descriptors for determining the biological activity of the molecule.One of the key characteristics in FMO theory is the difference between the energies of HOMO and LUMO (energy gap ∆E).The obtained values calculated using DFT/B3LYP-D3/cc-pvdz/LanL2DZ(Pt) level of theory are given in Table 2. Generally, the formation of all conjugates is accompanied by a significant decrease in ∆E: when going from fluorouracil FU to conjugates, the average value of decrease is 0.88 eV, from cisplatine to its conjugates-0.37eV.This indicates an increase in reactivity of the system, which could mean more active binding with biological sites.It is worth mentioning that the decrease of ∆E for cisplatin conjugates occurs due to an increase in the HOMO energy by 0.47 eV, which contributes to the realization of their nucleophilic potential.
Global reactivity descriptors such as electronegativity (χ), chemical potential (µ), chemical softness-hardness (η and S) and electrophilic index (ω) are highly successful in predicting stability properties and reactivity trends of the molecular systems.The formal definitions of all these descriptors and working equations for their computation have been described in [16].From theoretical calculations, it was found that chemical hardness (softness) value decreases (increases) when moving from drugs to conjugates, which indicates an increase in their chemical reactivity.The values of ω for conjugates shows that the electrophilicity of binary systems with fluorouracil increases, while with cisplatin, on the contrary, it decreases (Table 3).The work in [12] provides data on the correlation between an increase in the potentiating activity of cisplatin complexes with heterocyclic compounds (having a rigid conformational structure) and an increase in the values of their calculated dipole moments, in comparison to pure cisplatin.In this work, acylurea derivatives have the possibility of conformational rotation of C-C, C-N and N-N bonds, and no such correlation is observed.
Determining the localization of FMO is important for establishing preferred directions and regions of the molecule for attack by nucleophiles and electrophiles.Results of calculating the localization of HOMO and LUMO are shown in the form of 3D isosurfaces (Figure 5).
The formation of conjugates with fluorouracil 1+FU, 2+FU, 3+FU leads to the displacement of LUMO from the drug molecule and its localization on isoxazole heterocycles (Figure 5).Since LUMO determines interaction with target nucleophilic sites, this type of binding will be realized through the isoxazole heterocycle.Binding to electrophilic protein sites determines the localization of HOMO in conjugate.Calculations for conjugate molecules showed that HOMO in all fluorouracil conjugates is completely localized on isoxazole heterocycle or cytisine fragments, but not on the drug molecule.This means that the first act of interaction with the electrophilic and nucleophilic regions of protein sites is realized mainly due to the adjuvant molecule.In addition, based on previous findings, the 2+FU complex turned out to be the most stable, which likely interferes with the drug molecule's own action and makes the resulting conjugate the least effective against biological targets.
interferes with the drug molecule's own action and makes the resulting conjugate the least effective against biological targets.
As for cisplatin conjugates 1+СPt and 2+CPt, LUMO are uniformly distributed between the isoxazole heterocycle and the cisplatin molecule, and HOMO are located on the cytisine fragment.Meanwhile the cisplatin molecule is responsible for interaction with nucleophilic sites in the 3-CPt conjugate, and LUMO are concentrated here on the drug molecule.In contrast to 1+CPt and 2+CPt conjugates, HOMO in 3-CPt do not affect the alkaloid fragment, being dispersed between cisplatin and the isoxazole heterocycle.All the given facts may indicate a coordinated mechanism of action of cisplatin and the adjuvant in binary mixtures towards biological targets.As for cisplatin conjugates 1+CPt and 2+CPt, LUMO are uniformly distributed between the isoxazole heterocycle and the cisplatin molecule, and HOMO are located on the cytisine fragment.Meanwhile the cisplatin molecule is responsible for interaction with nucleophilic sites in the 3-CPt conjugate, and LUMO are concentrated here on the drug molecule.In contrast to 1+CPt and 2+CPt conjugates, HOMO in 3-CPt do not affect the alkaloid fragment, being dispersed between cisplatin and the isoxazole heterocycle.All the given facts may indicate a coordinated mechanism of action of cisplatin and the adjuvant in binary mixtures towards biological targets.

Docking Studies
Thymidylate synthase is one of the key enzymes in carcinogenesis and, for this reason, is a critical target for cancer chemotherapy.As previously noted, the main mechanism of action of fluorouracil is its conversion into 5-fluoro-2-deoxy-5-monophosphate (FdUMP), resulting in the inhibition of TS and DNA synthesis.In order to study the possibility of thymidylate synthase inhibition by isoxazolyl ureas 1-3 as a process competing with the binding of fluorouracil and its metabolite to the same target, molecular docking was carried out using software packages AutoDock/Vina 4.2.6 and CHIMERA 1.16, Biovia Discovery Studio 2024 for the visualization of the docked pose.The crystal structure of target protein TS (mouse thymidylate synthase, PDB ID: 5FCT) was downloaded from the Protein Data Bank.
The efficiency of interactions was evaluated based on the changes observed in predicted binding energies and bonds formed with the active site residues of TS.
As can be seen from Figure 6 and Table 4, all isoxazolyl ureas demonstrate a good binding degree to the thymidylate synthase protein site, forming an extensive system of residue interactions.The binding energy of 1-3 with tymidylate synthase is −9.3, −9.9, −10.0 kcal/mol, respectively, which even exceeds that of fluorouracil (−5.1 kcal/mol).Generally, ureas 1-3 form stable complexes with thymidylate synthase protein, mainly due to the hydrophobic interactions of electron-rich π systems of isoxazole heterocycle and aromatic rings with aminoacid residues, rather than hydrogen bonds, as occurs with fluorouracil and its metabolite.It was observed that ureas 2, 3, as well as 5-FU and FdUMP, occupied similar protein site-forming bonds with ASN212, CYS189, ARG169 and ARG44 amino acid residues (Table 4), which confirms the possibility of competing binding processes with the thymidylate synthase site.Meanwhile, for urea 1, the weakest binding was observed and it also interacts with other amino acid residues of the target protein.
Considering the high binding energy with TS of all ureas, molecular docking data cannot explain such a different manifestation of fluorouracil in the composition with adjuvants (both antagonism and synergism are observed).Ureas can probably be multitarget ligands that attack other proteins as well.In addition, summing up the data on energy characteristics of the fluorouracil-urea conjugates formation and docking data on binding energy with TS, we can assume that adjuvants will instead bind the FU into a stable complex, preventing the drug from exerting its effect (∆E f ranges from −6.19 to −16.88 kcal/mol); this especially applies to isoxazolylurea 2.
binding degree to the thymidylate synthase protein site, forming an extensive system of residue interactions.The binding energy оf 1-3 with tymidylate synthase is −9.3, −9.9, −10.0 kcаl/mоl, respectively, which even exceeds that of fluorouracil (−5.1 kcal/mol).Generally, ureas 1-3 form stable complexes with thymidylate synthase protein, mainly due to the hydrophobic interactions of electron-rich π systems of isoxazole heterocycle and aromatic rings with aminoacid residues, rather than hydrogen bonds, as occurs with fluorouracil and its metabolite.

General Chemistry Section
UV spectra were recorded on a Varian Cary 300 spectrophotometer using quartz cuvettes with l = 1 cm.The concentration of the studied compounds in methanol was 4 × 10 −5 -1 × 10 −4 M. IR spectra were registered on a Thermo Nicolet Protege 460 FT-IR spectrometer (Nicolet, Waltham, MA, USA) in KBr pellets. 1 H and 13 C NMR spectra were acquired on a Bruker Avance 500 spectrometer (Bruker, Bremen, Germany) in CDCl 3 .The residual solvent signals (CDCl 3 , δH 7.26, δC 77.2 ppm) were used as an internal standard.The assignment of signals in the 13 C NMR spectra was performed using the DEPT technique.
Liquid chromatography-mass spectrometry spectra were recorded on an Agilent 1200 LC-MS system with an Agilent 6410 Triple Quad Mass Selective Detector (Agilent Technologies Inc., Santa Clara, CA, USA) with electrospray ionization in the positive ion registration mode (MS2 scanning mode).An Agilent ZORBAX Eclipse XDB-C18 (4.6 × 50 mm, 1.8 µm) column (Agilent Technologies Inc., USA) was used.The mobile phase was MeCN-H 2 O + 0.05% HCO 2 H, with gradient elution from 40 to 90% MeCN in 10 min.A flow rate of 0.5 mL/min was used.
The optical activity of the compounds was measured on a polarimeter MCP 100 (Anton Paar, Graz, Austria).
Reagents and solvents used were of analytical grade, with the content of the main component being more than 99.5%.Phenyl (5-arylisoxazol-3-yl)carbamate was synthesized according to the previously described procedures [11].
(-)-Anabasine (colorless viscous liquid, turning yellow in air and in light; bp 276 General Procedure for the Synthesis of Anabasine and Cytisine Ureas 1-3 A mixture containing 1.6 g of anabasine (or 1.9 g of cytisine, 10 mmol) and 10 mmol of phenyl (5-arylisoxazol-3-yl)carbamate in 30 mL of benzene was boiled and stirred for 6 h.The solvent was removed in vacuum.Diethyl ether (50 mL) was added to the oily residue, and the mixture was boiled until complete crystallization had occurred.The crystalline product was separated by filtration, washed with ether and dried in air at room temperature.

In Vitro Biological Assays
Cell lines.The work was carried out on C6 glioma cell lines (rat) from the collection of the Republican Research and Practical Center for Epidemiology and Microbiology (RRPCEM, Minsk, Belarus).
Sample preparation.Test samples were dissolved in dimethyl sulfoxide (DMSO) to a concentration of 0.1 M. Then the resulting solution was diluted in an isotonic solution (0.9% sodium chloride solution) to a concentration of compounds of 2000 µM.
Solutions of the test compounds were added to the wells with cells in a volume of 10% of the total volume of the culture medium (in a ratio of 1 (test compounds): 9 (medium with cells)).The final concentration of the test compounds was 200 µM, DMSO was 0.2%.
Conducting experiments.Cells were seeded in wells of 96-well plates (Costar, Corning Inc., Corning, NY, USA) in DMEM (Merck KGaA, Darmstadt, Germany), supplemented with 10% fetal bovine serum (Sigma, St. Louis, MO, USA) and antibiotics (penicillin, streptomycin and neomycin, Biological Industries Beit Haemek, Israel).The test compounds were added to the wells at a final concentration of 200 µM, and/or antitumor drugs at a final concentration of 5 and 50 µM.In the control, instead of the tested compounds, a solvent dimethyl sulfoxide (DMSO), diluted in 0.9% sodium chloride, at a final concentration of 0.2%, was added.The number of samples in each series was n = 16 when studying the effects of 1, 2, 3, and n = 8 when studying the effects of combined use of 1, 2, 3 with antitumor drugs.Cells were cultured for 48 h at 37 • C and 5% CO 2 .Then the cells were removed from the plate and cell samples were prepared.To assess cell viability, cells were stained with 7-aminoactinomycin (7-AAD) or propidium iodide.To determine the total number of cells and the number of viable cells in the samples (as an end point characterizing the antitumor effect), FLOW-COUNT™ reference fluorospheres (Beckman Coulter, Inc., Brea, CA, USA) were used according to the manufacturer's instructions.Cellular samples were analyzed using a BD FACSCanto II flow cytometer (Becton Dickinson, Franklin Lakes, NJ, USA).To evaluate the effect of the combination against a tenfold increase in the dosage of the drug, a concentration of 5 µM was taken where a weak effect and low toxicity would appear, and this was increased to 50 µM for evaluating the effect of the drug-adjuvant combination against a tenfold increase in the drug dosage.
Statistical processing.The resulting digital material was processed by methods of variation statistics using the Excel 2016 and Statistica 7 software packages.The data are presented as the mean and its standard error.Differences between series were considered significant at a significance level of p < 0.05, according to Student's t test.

Quantum Chemical Methods
The Gaussian 16 package software was used for DFT calculation, and for the visualization of the results of the DFT calculation, Gauss View 6.0.16 software was used.The calculations were performed within the framework of the DFT-D3 dispersion correction method [17], which, as shown in a number of works [18][19][20], is successfully used to take into account the inter-and intramolecular long-range dispersion interactions.Becke's three-parameter exchange functional [21], in combination with the Lee et al. correlation functional [22] (B3LYP), were used.This was because this functional gives a high accuracy compared to experimental values for structures with intermolecular interactions, as demonstrated in [23,24].
Dunning's correlation consistent polarized valence double zeta basis set with adding diffuse functions aug-cc-pVDZ was used.Diffuse functions were added for describing long-range interactions, such as Van der Waals forces, and noncovalent interactions, such as hydrogen bonding [25][26][27].The LANL2DZ (Los Alamos National Laboratory 2 double zeta) basis set with effective core potential (ECP) was used for Pt [28].
The polarizable continuum model (PCM solvent is considered as a continuous dielectric medium) was used for solvent phase calculations [29].The PCM model implements a self-consistent reaction field (SCRF) approach and defines solvent polarization in terms of electrostatic potential.The following discussions are based on this method, if not noted otherwise.No symmetric constraints were imposed during geometrical optimizations.The energy minima were identified by subsequent frequency calculations.

Molecular Docking
Molecular docking software packages AutoDock/Vina 4.2.6 and CHIMERA 1.16 were employed to calculate the protein-ligand binding interactions.The 3D crystal structure of the mouse thymidylate synthase (TS) was downloaded from the RCSB Protein Data Bank (PDB ID: 5FCT, 1.55 Å resolution).The protein was parameterized using AMBER f14SB force field, and the Gasteiger charges and hydrogen atoms were added to the examined receptor with native ligands.All water molecules and cofactors were removed before docking.The binding pockets of protein were identified using the CASTp package (http://cast.engr.uic.edu,accessed on 15 March 2024), and the grid box had the following parameters: ∆X 24 Å, ∆Y 32 Å, ∆Z 29 Å centered at (5.0; −5.5; −5.0).The docking calculations were performed with exhaustiveness of 8 and an energy range equal to 3 kcal/mol.
The validation of the procedure involved docking the X-ray ligand conformation rigidly to the active site following the protein preparation step.The 3D structure of UFP ligand was isolated from the crystal UFP-C2F-TS complex.The results of a comparative analysis of complexes of UFP inhibitor with TS, constructed using the methods of X-ray crystallography and molecular docking (Figure 7), indicate that the calculation protocol we used provided sufficient accuracy (RMSD < 2 Å, Figure 7) in predicting the orientation of the ligand in the catalytic site of the enzyme.This suggests that the data we obtained for the isoxazolylurea adequately describe the main geometric and energetic characteristics of their complexes with thymidylate synthase.Intermolecular interactions between ligand and receptor protein were analyzed by BIOVIA Discovery Studio Visualizer 2024. kcal/mol.
The validation of the procedure involved docking the X-ray ligand conformation rigidly to the active site following the protein preparation step.The 3D structure of UFP ligand was isolated from the crystal UFP-C2F-TS complex.The results of a comparative analysis of complexes of UFP inhibitor with TS, constructed using the methods of X-ray crystallography and molecular docking (Figure 7), indicate that the calculation protocol we used provided sufficient accuracy (RMSD < 2 Å, Figure 7) in predicting the orientation of the ligand in the catalytic site of the enzyme.This suggests that the data we obtained for the isoxazolylurea adequately describe the main geometric and energetic characteristics of their complexes with thymidylate synthase.Intermolecular interactions between ligand and receptor protein were analyzed by BIOVIA Discovery Studio Visualizer 2024.

Conclusions
Taking into account the combination of theoretical calculations, molecular docking and test results, we can conclude: In vitro tests of alkaloid-based isoxazolylureas on the C6 rat glioma model showed that a synergistic effect was observed to varying degrees for all studied compounds when using binary mixtures with alkylating chemotherapy drugs cisplatin and ribomustine, whereas antagonism and the tendency to antagonism are noted for the remaining drugs, except for the combination with fluorouracil, where conflicting results were observed.

Conclusions
Taking into account the combination of theoretical calculations, molecular docking and test results, we can conclude: In vitro tests of alkaloid-based isoxazolylureas on the C6 rat glioma model showed that a synergistic effect was observed to varying degrees for all studied compounds when using binary mixtures with alkylating chemotherapy drugs cisplatin and ribomustine, whereas antagonism and the tendency to antagonism are noted for the remaining drugs, except for the combination with fluorouracil, where conflicting results were observed.
The clear antagonistic effect of the combination of urea with cyclophosphamide, fluorouracil and etoposide, resulting in a 2-3 fold decrease of antiproliferative action, can be explained based on the results of quantum chemical modeling of the spatial structure of the adjuvants.The structure of isoxazolylurea 2 differed from the other two in its most linear spatial arrangement, which could result in the varying degrees of conjugate stability.In addition, analysis of FMO localization in fluorouracil conjugates showed that HOMO and LUMO, in contrast to cisplatin conjugates, were completely localized on isoxazole heterocycle or alkaloid fragments, but not on drug molecules.This means that the interaction with the electrophilic and nucleophilic regions of protein sites is realized mainly due to the adjuvant molecule.
Based on the molecular docking data, all isoxazolylureas demonstrated a good degree of binding to the thymidylate synthase protein site, forming an extensive system of residue interactions mainly due to the hydrophobic interactions of electron-rich π systems of isoxazole heterocycle and aromatic rings with amino acid residues.This probably confirms the possibility of competing binding processes with the thymidylate synthase site, but does not explain the synergy observed in the case of the composition of fluorouracil with anabasine isoxazolylurea 3. It cannot be ruled out that ureas could be multitarget ligands that attack other proteins as well, and observed discrepant behavior of fluorouracil compositions should instead be explained by varying degrees of stability of drug-adjuvant complexes with the formation of an inactive product.
Isothiazolylurea adjuvants shaped conjugates with cisplatin with non-covalent interactions between molecules due to hydrogen and van der Waals bonds generally forming more stable complexes than with fluorouracil.FMOs here were uniformly distributed between the adjuvant and the cisplatin molecule.This may indicate a coordinated mechanism of action cisplatin and the adjuvant in binary mixtures towards biological targets, which was confirmed by the experimentally observed synergistic effect.
The most pronounced synergistic effect was observed in the case of ribomustine mixtures with all synthesized compounds, which encourages further research of the dosedependent characteristics for this composition.

Figure 2 .
Figure 2. Anticancer drugs used in composition with isoxazolylurea adjuvants.Figure 2. Anticancer drugs used in composition with isoxazolylurea adjuvants.

Figure 2 .
Figure 2. Anticancer drugs used in composition with isoxazolylurea adjuvants.Figure 2. Anticancer drugs used in composition with isoxazolylurea adjuvants.

Figure 3 .
Figure 3. Decrease in viable C6 rat glioma cell concentration for binary mixtures and sum of the effects of an individual compound and drug (C1-3 = 200 µM, Cdrug = 5 µM; number of repetitions for drug-urea composition in each experiment n = 8).Additivity-the effect of drugs combination equals the sum of their effects when given alone; antagonism-drug combinations show adverse

Figure 3 .
Figure 3. Decrease in viable C6 rat glioma cell concentration for binary mixtures and sum of the effects of an individual compound and drug (C 1-3 = 200 µM, C drug = 5 µM; number of repetitions for drug-urea composition in each experiment n = 8).Additivity-the effect of drugs combination equals the sum of their effects when given alone; antagonism-drug combinations show adverse results of a decrease in antitumor activity compared to the sum of individual effects; synergy-agents given together produce higher antitumor activity than the sum of individual effects.

Figure 4 .
Figure 4. Calculated optimized structures with consideration given to aqueous medium for conjugates 1-3 with fluorouracil and cysplatine.

1+FUFigure 5 .
Figure 5. Localizations of HOMO and LUMO in molecules of conjugates in the form of 3D isosurfaces, with consideration given to water.

Table 1 .
Decrease in viable C6 rat glioma cell concentration (%) for antitumor drugs (c = 5 and 50 µM), compounds 1-3 (c = 200 µM) and their composition.Σ is the sum of the effects of a single compound and drug.

Table 2 .
Energy characteristics of the conjugates formation process from the initial components 1

Table 3 .
Theoretical electronic parameters (energies of HOMO and LUMO, energy gap ∆E,) and global reactivity descriptors for compounds 1-3, flourouracil, cisplatin and their conjugates.