Synthesis, Physicochemical Characteristics and Plausible Mechanism of Action of an Immunosuppressive Isoxazolo[5,4-e]-1,2,4-Triazepine Derivative (RM33)

Previous studies demonstrated strong anti-inflammatory properties of isoxazolo[5,4-e]-1,2,4-triazepine (RM33) in vivo. The aim of this investigation was to describe synthesis, determine physicochemical characteristics, evaluate biological activities in murine and human in vitro models, as well as to propose mechanism of action of the compound. The compound was devoid of cell toxicity up to 100 μg/mL against a reference A549 cell line. Likewise, RM33 did not induce apoptosis in these cells. The compound stimulated concanavalin A (ConA)-induced splenocyte proliferation but did not change the secondary humoral immune response in vitro to sheep erythrocytes. Nevertheless, a low suppressive effect was registered on lipopolysaccharide (LPS)-induced splenocyte proliferation and a stronger one on tumor necrosis factor alpha (TNFα) production by rat peritoneal cells. The analysis of signaling pathways elicited by RM33 in nonstimulated resident cells and cell lines revealed changes associated with cell activation. Most importantly, we demonstrated that RM33 enhanced production of cyclooxygenase 2 in LPS-stimulated splenocytes. Based on the previous and herein presented results, we conclude that RM33 is an efficient, nontoxic immune suppressor with prevailing anti-inflammatory action. Additionally, structural studies were carried out with the use of appropriate spectral techniques in order to unequivocally confirm the structure of the RM33 molecule. Unambiguous assignment of NMR chemical shifts of carbon atoms of RM33 was conducted thanks to full detailed analysis of 1H, 13C NMR spectra and their two-dimensional (2D) variants. Comparison between theoretically predicted chemical shifts and experimental ones was also carried out. Additionally, N-deuterated isotopologue of RM33 was synthesized to eliminate potentially disturbing frequencies (such as NH, NH2 deformation vibrations) in the carbonyl region of the IR (infrared) spectrum to confirm the presence of the carbonyl group.


Introduction
Compounds containing an isoxazole ring [1] constitute an important source of valuable drugs, designed to treat inflammatory states and diseases of different etiology. Syntheses, classification, mechanisms of action and therapeutic application of registered isoxazole derivatives, as well as these under preclinical investigations, were recently described in detail in several reviews [2][3][4].
Among immunosuppressive isoxazole derivatives synthesized by our research team, an isoxazolo [5,4-e]-1,2,4-triazepine derivative (denoted as RM33) has attracted special interest due to its decisive immunosuppressive properties in in vivo models with a potency compared with that of cyclosporine in inhibition of the immune response. We showed that RM33 administered intraperitoneally (i.p.) before immunization significantly decreased number of antibody forming cells to sheep red blood cells (SRBCs) in mice [27] and carrageenaninduced foot pad inflammation in rats [28]. The compound also inhibited an inductive and effectual phase of delayed type hypersensitivity to ovalbumin (OVA) [27]. RM33 was also effective when administered per os. A suppressive effect of RM33 in generation of the delayed type hypersensitivity to OVA could be also achieved when the compound was admixed with Freund's complete adjuvant for a subcutaneous immunization of mice with OVA. RM33 also reduced serum level of TNFα induced by intravenous administration of LPS [28]. In vitro experiments, performed in parallel, unexpectedly revealed that RM33 did not inhibit proliferative response of mitogen-induced peripheral blood mononuclear cells and LPS-induced TNFα production.
This apparent discrepancy between in vivo and in vitro actions of RM33 prompted us to investigate the mechanism of action of the compound in a mouse in vitro model by determination of mitogen-induced splenocyte proliferation, secondary humoral immune response to SRBC and cytokine production by peritoneal cells. In addition, we investigated cytotoxic actions of the compound towards cell lines and cells from the lymphoid organs. In order to find a plausible mechanism of action of RM33, we analyzed effects of the compound on expression of signaling molecules in mouse bone marrow cells, thymuses and splenocytes, as well as cell lines representing major cell types such as T, B cells, macrophages and myelocytic precursors. In addition, we evaluated effects of RM33 on viability and apoptosis of resident cells from the lymphoid organs and J744 macrophage cell line. The function of macrophages was evaluated by RM33 effects on LPS-induced TNFα production in peritoneal exudates cell culture and pinocytic activity of J774 cell line. Finally, we evaluated effect of RM33 on the levels of cyclooxygenases in LPS-stimulated splenocytes. The biological part of the investigation is accompanied by description of RM33 synthesis and analysis of its identity and physicochemical properties by application of various spectroscopic techniques. It is important to note that the isoxazolo [5,4-e]-1,2,4-triazepine system has been so far little studied, both physicochemically and biologically [29].
Hydrazide 4 was obtained in a reaction of 5-amino-3-methylisoxazole-4-carboxylic acid azide with methylhydrazine, according to the method described in [25]. The final RM33 was synthesized in two steps through intermediate product, i.e., N-acylhydrazone 4a, which was formed in the nucleophilic addition reaction of terminal NH 2 group of hydrazide 4 to carbonyl carbon of acetaldehyde. In second step intermediate 4a was cyclized as a result of reaction of intramolecular nuclephilic addition of isoxazole amine group (aza-Michael type reaction) to hydrazone carbon catalyzed by indium (III) trifluoromethanesulfonate.

Proof and Analysis of the Structure of RM33 Compound
The oxazolo [5,4-e]-1,2,4-triazepine system is relatively little known and only a few its derivatives have been described. As, up to date, no proof of a structure of the considered compound, RM33, has been conducted, it is necessary to present reliable research results confirming its structure. Firstly, the elemental analysis and MS spectrum of the investigated compound point out that its molecular formula is C 8 H 12 N 4 O 2 .
In order to prove the presence of a carbonyl group in the considered molecule of RM33, its ATR FT-IR spectrum was recorded. This spectrum ( Figure S1 in Supplementary Materials) shows that the strong band of the carbonyl region occurs at an extremely low wavenumber (1607 cm −1 ), as is the case for a carbonyl of tertiary cyclic amides (lactams). Assuming the presence of a carbonyl group in the tested compound, it should be taken into account that stretching vibration bands (mostly strong or of very strong intensity) of carbonyl groups in the IR spectrum can be often covered or perturbed with intensive bands of deformation vibrations of N-H (mainly scissoring vibration of NH 2 ) groups. For example, 5-amino-3-methyl-4-isoxazolecarbohydrazide (without an N-methyl group), structurally similar to compound 4, exhibits the intense bending (scissoring) vibration of the NH 2 amino group with a maximum at about 1616 cm −1 and a weak scissoring vibration of the terminal NH 2 of the hydrazide group at about 1684 cm −1 [35]. Therefore, for this reason, the N-deuteratation of RM33 was carried out to remove potentially disturbing and covering bands to unambiguously assign the appropriate band to the carbonyl group (see Figure 1). Generally, the results show that the structure of bands and the wavenumber values in the vibration region of the carbonyl groups do not change much after N-deuteration (the visualization of the spectra, Figures S2-S4). After deuteration, this band is still strong and visible with maximum at 1603 cm −1 , and thus the red-shift of this band is slight, which proves that was not covered and disturbed by bending vibrations of bonds NH (NH 2 ) and really coming from stretching vibration of carbonyl group. Additionally, stretching vibrations of the bonds between ring atoms of an aromatic isoxazole occur below 1600 cm −1 , i.e., 1587 and 1564 cm −1 , so they also cannot cover, disturb or be mistakenly taken for the considered carbonyl band. For comparison, for the 5-amino-3-methyl-4-isoxazolecarbohydrazide (without an N-methyl group) analogous to compound 4, they were observed at 1562 and 1523 cm −1 as very intense bands of stretching vibrations, C=N and C=C, respectively [35]. Therefore, the conclusion is that such a low frequency of stretching vibration of carbonyl group of considered compound RM33 is the result of strong coupling of this group with the aromatic system of the isoxazole ring through resonance effect in conjugated system. It additionally proves a low ring tension, which confirms the presence of seven membered ring of triazepine forming by indium(III) catalyzed cyclization step. However, N-deuteration of the RM33 made the disappearance of intense bands associated with NH bond stretching vibrations visible in the 3100-3400 cm −1 range, which are shifted to a lower frequency and appear as intense stretching vibrations of ND bonds in the 2300-2500 range due to the mass of deuterium being double that of a proton ( Figures S2 and S3). The unchanged but better visible bands in the 2900-3000 cm −1 range (maxima at 2932 and 2988 cm −1) after deuteration come from the C-H stretching vibrations of the methyl and methanetriyl (methin) groups.
Since standard one-dimensional (1D) 1 H and 13 C NMR spectroscopy is insufficient to unambiguously assign chemical shifts to the corresponding hydrogen and carbon atoms ( Figure 2), two-dimensional (2D) NMR methods were used. Admittedly, on the basis of the 1 HNMR spectrum in DMSO-d 6 (see Figures S5A,B in Supplementary Materials), it can be assumed that there is a coupling between protons of the methyl group with a chemical shift of 1.26 ppm (broadened doublet, J ≈ 5 Hz) and a single proton with a chemical shift of 4.41 ppm (multiplet). However, due to the fact that the 4.41 ppm signal (multiplet) is heavily distorted and broadened, it is impossible to determine unquestionably its constant coupling, confirming the scalar coupling of this proton to the methyl group at 1.26 ppm (broadened doublet, J ≈ 5 Hz). Moreover, it is also impossible to clearly determine whether the protons with chemical shifts of 5.72 and 8.54 ppm are coupled to this proton at 4.41 ppm. For this reason, the spectrum of 1 H-1 H 2D-COSY45 was recorded, which dispels this doubt by indicating that they (i.e., N-H) are indeed scalar coupled to the considered proton (4.41 ppm, multiplet) (see Figure S9A in Supplementary Materials). It also confirms the scalar coupling between the protons of the methyl group at 1.26 ppm (broadened doublet) and the proton of the CH group (multiplet at 4.41 ppm). 1 H-13 C 2D HSQC spectroscopy, which reveals the coupling between protons and carbons through one bond, was used for correct assignment of chemical shifts of the carbon atoms directly bound to protons. The HSQC spectrum (see Figure S9B in Supplementary Materials) shows that carbon with a chemical shift 37.4 ppm is directly bound to three protons with a chemical shift of 2.94 ppm and it means that this methyl group is connected to nitrogen atom (position 5 of the system). In addition, the HSQC spectrum shows that the methyl group protons with a chemical shift of 1.26 ppm (position 7 of the system) are coupled to a carbon with a chemical shift of 20.9 ppm, which means that three protons (isoxazole methyl) with a chemical shift of 2.21 ppm are bound to a carbon with a shift of 11.7 ppm. Due to the fact that both protons connected to nitrogen atoms (N-H at positions 6 and 8 of the system) are scalarly coupled through three bonds with the proton from position 7, it is not possible to assign them unambiguously to the appropriate chemical shifts by means of the COSY45 spectroscopy. For this reason, the NOESY spectrum was recorded, which shows the spatial couplings through the space between distant protons (up to 4 Å). The NOESY spectrum showed that a proton with a chemical shift of 5.94 ppm is coupled with the protons (2.94 ppm) of the N-methyl group at position 5 of the system (see Figure S9C). On the other hand, there is no visible spatial coupling between the proton with chemical shift of 8.64 ppm and the protons of the N-methyl group, which means there is too large a spatial distance. Therefore, the conclusion is that the proton with chemical shift 5.94 ppm is a proton at position 6 of the system, which necessarily places the proton with a chemical shift of 8.64 ppm at position 8. In general, the proton at position 8 should be more acidic than the proton at position 6 due to the proximity of the oxygen atom of the isoxazole ring and an electron withdrawing effect through this ring. Additionally, the theoretically predicted chemical shifts for these protons indicate that the proton at position 6 will be more shielded by electrons than the proton at position 8 (see Section 2.1.2, Theoretical Prediction of Chemical shifts of 1 H NMR and 13 C NMR Spectra of RM33). To unambiguously assign quaternary carbon atoms to the appropriate values of their chemical shifts, the 2D 1 H-13 C Heteronuclear Multiple Bond Correlation (HMBC) method of NMR spectroscopy is particularly useful. In the HMBC spectrum (in DMSO-d 6 ), long range couplings (through 2 to 4 bonds) between protons and carbons are visible (see Figure S9D): a.
Strong long range coupling through two bonds between protons of methyl group (2.21 ppm) at position 3 of system and the ring carbon 3 (161.1 ppm); b.
Slightly weaker through three bonds between protons of methyl group (2.21 ppm) at position 3 of system with the bridged ring carbon 3a (88.0 ppm); c.
Very weak, but noticeable, through four bonds between protons of methyl group (2.21 ppm) at position 3 of system with the bridged carbon 8a (166.2 ppm); d. Long range coupling through three bonds between protons of methyl group (2.94 ppm) at position 5 of the system and the carbonyl carbon atom (165.7 ppm) are also clearly visible. It should also be noted that the chemical shifts of the carbons of the carbonyl group and bridged carbon 8a of the isoxazolotriazepine system differ only by about 0.5 ppm and strongly depend on the solvent used. In the strongly polar protic methanol, the carbon nucleus of the carbonyl group is the most deshielded of all the carbon atoms of the isoxazolotriazepine system, and in the slightly less polar aprotic DMSO-d 6 , it is carbon 8a (Figure 2 and Figures S7, S8 and S9D,E in Supplementary Materials).
Finally, the HR-ESI-MS technique was used, the major advantage of which is that it is very selective since it measures the exact mass of a compound, allowing even minor changes in the structure to be distinguished. For RM33, a mass accuracy error was below 5 ppm (i.e., 0.736 ppm), which also confirms the assumed structure (see Section 4.1.1 in Materials and Methods).

Theoretical Prediction of Chemical Shifts of 1 H NMR and 13 C NMR Spectra of RM33
For comparison, experimental chemical shifts with theoretically predicted appropriate calculations were carried out using methods described in detail in Sections 4.1.3 and 4.2. The obtained results (Table 1) for proton chemical shift simulation indicate that the good agreement between the calculated chemical shifts and the experimental values is only for the protons bound to carbon atoms. The biggest discrepancy between theoretically predicted and experimental chemical shifts in the 1 H NMR spectrum is for the N-H acidic protons. Their calculated electron shielding values are much bigger (i.e., lower chemical shifts) than the experimental values. Additionally, the calculated chemical shift for the nitrogen bonded proton in position 6 is smaller than for the carbon bonded proton in position 7, which is inconsistent with the experimental data, where it is the reverse. In turn, the experimental 13 C NMR spectrum and theoretically predicted results consistently show that the carbon C7 atom is the most shielded by elelctrons (of course, not counting the carbons of the methyl groups) of all ring carbon atoms of isoxazolo[5,4-e]-1,2,4-triazepine system. In the aromatic isoxazole ring of the isoxazolotriazepine system, the quaternary bridged carbon C3a is the most shielded by electrons. It is caused by substantial contribution of resonance structures where high electron density appears, which is consistent with previous observations regarding 5-amino-3-methyl-4-isoxazolecarbohydrazide [35]. Addi-tionally, the calculated chemical shifts for the 13 C NMR spectrum show that the quaternary carbon of the carbonyl group is most deshielded, followed by the C3 quaternary and in the order of the bridged quaternary C8a carbon, which is inconsistent with the experimental (in DMSO-d 6 ) results, which show that the most deshielded is C8a quaternary carbon, followed by the carbonyl, and the least is C3. However, in CD 3 OD, the most deshielded is carbonyl carbon, then C8a, and the least is C3. A very strong solvent effect can be seen here. This shows that, although predicted theoretical chemical shifts are useful in interpreting the NMR spectra, unambiguously certain assignments of chemical shifts to the respective atoms can only be achieved with appropriate NMR experimental techniques.  * Acquisition and process parameters used during recording 1 H and 13 C NMR spectra are included in Supplementary Materials; experimental chemical shifts were corrected relative to chemical shift of the solvent; ** SD-the standard deviations of the chemical shifts were estimated by measuring the NMR spectra of three independent samples of the compound under the same conditions.

Activity of RM33 in Immunological Models In Vitro
Previous studies on RM33 activity in vivo (Table 2) revealed significant immunosuppressive and anti-inflammatory properties of the compound in several mouse and rat experimental models. These actions regarded both antigen-specific and nonspecific immune responses, of both humoral and cellular types. In addition, the compound affected both inductive and effectual phase of the immune response. The anti-inflammatory actions of RM33 were reflected by respective changes in cytokine release and amelioration of histopathological alterations. Table 2. The summary of effects of RM33 on innate and adaptive immune responses in vivo.

Effects of RM33 Administration Reference
Humoral immune response Inhibition of the inductive phase [27] Delayed type hypersensitivity Inhibition of the inductive and effectual phase [27] Edema induced by complete Freund's adjuvant Decrease in footpad edema [28] LPS-induced endotoxemia Significant inhibition of serum TNFα, small inhibitory effect on IL-6 and no effect on IL-10 serum level [28] Carrageenan reaction Decreased serum TNFα, reduced: mast cell and macrophage infiltration and edema of the connective tissue [28] In a preliminary experiment, RM33 was tested for its potential toxicity against a reference A549 cell line. Appropriate dilutions of the solvent (DMSO) served as controls.
The results show no sign of toxicity up to 100 µg/mL concentration of the compound as evaluated by MTT colorimetric method and cell morphology ( Figure S13). A small toxic effect of RM33 was only observed at 200 µg/mL concentration ( Figure 3). To determine effects of RM33 on metabolism and number of viable cells upon culture with RM33, the cells from mouse lymphatic organs were incubated overnight with 20 µg/mL of RM33 and cell viability was determined by MTT colorimetric method in comparison to appropriate DMSO control culture. The results (Table S1)   In addition, the results presented in Table 3 revealed no induction of apoptosis by the compound, at 2-50 µg/mL concentration range, with regard to these cell populations. Table 3. Effect of RM33 on apoptosis of cells from mouse lymphoid organs. Cells isolated from the lymphoid organs were cultured 106/well overnight with RM33 (2, 10 or 50 µg/mL) and appropriate dilutions of DMSO. For DNA fragmentation evaluation cells were stained with propidium iodide (50 µg/mL) and analyzed by FACS. The results are presented as increase in % of apoptotic cells vs. DMSO control.

Organ
Increase in % of Apoptotic Cells RM33 (µg/mL) 2 10 50 Bone marrow 0 0 0 Thymus 1 3 6 Spleen 1 1 1 Figure 5 presents effects of RM33 on ConA-induced proliferation of splenocytes at a concentration range of 1.25-5 µg/mL. Cyclosporin A (CsA) served as a reference antiproliferative drug. The effects of RM33 were compared with appropriate concentrations of DMSO. The results showed stimulatory effects of RM33 in the investigated concentration range. A very strong suppressive effect was registered using CsA. Nevertheless, in some experiments, this stimulation was not registered (Table S2) and was absent with proliferation of human peripheral blood lymphocytes stimulated by PHA (Table S3).

Figure 5.
Effect of RM33 on ConA-induced mouse splenocyte proliferation. Splenocytes were distributed into tissue culture plates at a density of 2 × 10 5 /100 µL/well. An amount of 2.5 µg/mL ConA was added to induce cell proliferation. RM33 and CsA were added to the cultures at 1.25-5 µg/mL concentration range. After a 3-day incubation, the cell proliferation was determined using the colorimetric MTT assay [36]. The results are presented as the mean OD ± SE from quadruplicate determinations, * p < 0.05 (ANOVA; vs. DMSO dilution equivalent).
The compound did not affect the magnitude of the secondary humoral immune response to SRBC in vitro (Table S4). Nevertheless, RM33 moderately inhibited, in a dose dependent manner, LPS-induced splenocyte proliferation (Figure 6), where B cells are the targets [37]. In addition (Figure 7), the compound inhibited LPS-induced TNFα production by rat peritoneal cells at concentration of 5 µg/mL (3581 vs. 13,162 pg/mL). Of note, RM33 elevated the spontaneous production of TNFα by 2-fold. IL-6 production was not significantly affected. Figure 6. Effect of RM33 on LPS-induced mouse splenocyte proliferation. Splenocytes were distributed into 96-well flat-bottom tissue culture plates at a density of 2 × 10 5 /100 µL/well. LPS (50 µg/mL) was added to induce cell proliferation. RM33 and CsA were added to the cultures at 1.25-5 µg/mL concentration range. After a 3-day incubation, the cell proliferation was determined using the colorimetric MTT assay. The results are presented as the mean OD ± SE from quadruplicate determinations, * p < 0.05 (ANOVA; vs. DMSO dilution equivalent).

Figure 7.
Effect of RM33 on LPS-induced cytokine production by rat peritoneal cells. Peritoneal cells were incubated overnight in 24-well culture plates containing 5 × 10 6 cells/mL. LPS was applied at concentration of 1 µg/mL and RM33 at 1 and 5 µg/mL. Control cultures contained appropriate dilutions of DMSO, corresponding to DMSO dilutions at respective RM33 doses. The determination of TNFα and IL-6 was performed by bioassays [38,39]. The results are presented in ng/mL.
The effects of the compound on expression of signaling molecules in mouse bone marrow cells, thymocytes and splenocytes are presented in Tables 4 and 5. RM33 caused differential changes in the expression of MAP kinases in cells from the respective organs (Table 4). In bone marrow cells, significant changes included block of ERK1 and decreases in p38γ and p38δ subunit expression as well as significant increases in p38β and JNK expression. In the splenocytes, significant augmentation of signaling protein expression included all p38 subunits and JNK with an exceptional increase (256×) in p38δ. In the thymocytes, expressions of p38β and p38δ, as well JNK, were elevated. Table 4. Effect of RM33 on expression of MAP kinases in the mouse lymphoid organs. Cells from the organs were incubated overnight with 10 µg/mL of RM33. Total RNA isolation, reverse transcription and quantitative analysis of gene expression by Real Time PCR were performed as described in the Materials and Methods. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as a housekeeping gene to normalize values obtained for the target genes. The values presented in the table reflect relative fold changes in response of the cells to RM33 compound, when the GAPDH gene expression was a control and assumed as equal to one. In the bone marrow cells, the compound elevated expression of caspase 9 and Fas and moderately NFκB1 (Table 5). A stronger stimulation of caspases 9 and 3 was observed in thymocytes, accompanied by an increase in Bcl-2 and NFκB1 expression. In the splenocytes, RM33 induced comparable levels of caspases 3 and 9 as in thymocytes, and significant increases in expression of Fas and NFκB1. Table 6 presents effects of the compound on expression of MAP kinases in Jurkat cells. The most notable changes were observed among expression of p38 subunits, i.e., a very deep inhibition of p38γ and to a lesser degree of p38α, and 3.7-fold elevation of the p38δ subunit. On the other hand, no changes (not shown) in expression of caspases 3, 7, 8, 9, Bcl2, Fas, NFκB1 and p53 were registered. The effects of RM33 on the expression of signaling molecules in macrophage J774 cell line were strikingly strong (Table 7). Among signaling molecules of the MAP kinase family, strong elevations in expression were observed for ERK-2, p38β, p38γ, p38δ and JNK proteins. Relatively lower increases were observed with regard to caspase 3, Bcl-2 and Fas molecules, except NFκB1, whose expression was strongly upregulated. Although effects of RM33 on the pinocytic activity of J774 cells were not significant ( Figure S14), this effect seemed to be dose dependent-i.e., 100 µg/mL concentration of the compound inhibited and a concentration of 25 µg/mL tended to stimulate pinocytosis. The changes in expression of signaling molecules in immature B cell WEHI 231 line and monocyte/neutrophil precursors HL-60 cell line were negligible and are shown in the Supplementary Materials (Tables S5-S7).
In the last stage of our investigations on the mechanism of action of RM33, the compound was added to the culture of splenocytes stimulated with 20 µg/mL of LPS and on day 3 of the incubation the cellular concentration of cyclooxygenase 2 was measured by ELISA assay. The results shown in Table 8 revealed a small increase in COX-2 level in nonstimulated control culture and a 38.9% increase in the culture of LPS-stimulated cells.

Discussion
In this report, we attempted to elucidate the intriguing phenomenon of divergent activities of RM33 in vivo and in vitro. RM33 demonstrated significant suppressive effects on the humoral and cellular response in vivo, as well as on inflammatory responses [27,28], but negligible actions on mitogen-induced proliferative response of splenocytes and secondary humoral immune response in vitro, as found in this investigation (Table S4). The discrepancies between strong immunosuppressive actions of RM33 in vivo and lack of significant effects on humoral immune response and lymphocyte proliferation in vitro suggested that the compound could act as a prodrug and should be metabolized in the liver to become effective, as in the case of leflunomide [5]. However, since RM33 inhibited in vitro LPS-inducible TNFα production and splenocyte proliferation and regulated cyclooxygenase 2 production, such a possibility could be excluded.
In the majority of experiments, we studied actions of RM33 at relatively low concentrations in nonstimulated cells and cell lines, so we did not expect increases in expression of signaling molecules associated with apoptosis. In fact, RM33 did not show cytotoxic effects on the reference cell line and resident cells from lymphatic organs. Eventually, the experiment with iodine propidine excluded a proapoptotic property of the compound (Table 3).
We consistently observed changes in expression of signaling molecules from MAP kinase family responsible for cell proliferation/differentiation or apoptosis, indicating antiapoptotic action of RM33. This was particularly evident in cases of low or lack of p38α expression, initiating a proapoptotic signal [40] and increases in p38β and p38δ, the molecules mediating cell activation or differentiation [41][42][43]. Therefore, the increases in p38δ expression in Jurkat cells, accompanied by the lack of changes in expression of caspases, Fas, NFκB1 and p53, involved in cell apoptosis (Tables 5 and 6), could explain stimulation of ConA-induced splenocyte proliferation ( Figure 5). In addition, the tremendous increases in NFκB1 expression in the macrophage J774 cell line (Table 7) seem to be correlated with stimulation of the spontaneous TNFα production in rat peritoneal cells (Figure 7) since TNF α production is dependent on NFκB1 expression [44]. At the same time, both the metabolism and pinocytic activity of J774 cells were not significantly affected ( Figure S14, Table S8) demonstrating no negative effect of RM33 on cell function. On the other hand, relevant to the postulated anti-inflammatory RM33 action, rat peritoneal cells, stimulated with LPS, produced less TNFα upon treatment with RM33 ( Figure 7). In addition, LPS-induced proliferation of B cells from mouse splenocytes was also inhibited ( Figure 6). Thus, both actions of RM33, observed in cell cultures, support anti-inflammatory properties of the compound in vivo against proinflammatory stimuli, such as LPS, Freund's complete adjuvant and carrageenan [27,28]. Of importance, all the studied proinflammatory stimuli [45][46][47] use both TLR2 and TLR4 receptors, respectively.
The suppressive property of RM33 on the humoral immune response to SRBC was limited only to the inductive phase in the in vivo system. The lack of its effects in vitro (Table S4) could be explained by differences in the experimental models, since in the model of the secondary immune response in vitro antigen-specific T and B cells are transferred from antigen-primed mice to the culture additionally immunized with the antigen. This model already expresses a low level SRBC-specific background response and antigenspecific B cells may act as both antigen presenting and accessory cells to develop a strong secondary immune response [48,49].
Although the inhibition by RM33 of LPS-induced IL-8 in human whole blood cultures was small (Table S9), the inhibition of IL-8 production could have a significance in the in vivo suppression of carrageenan-induced inflammation [28] by lowering neutrophil migration [50] to the site of inflammation. In the rat model of carrageenan-induced inflammation, IL-8 plays a major role [51] and the signaling pathway induced by carrageenan involves TLR4 [45], as in the case of LPS [47].
The effect of RM33 on cell signaling in bone marrow cells (Tables 4 and 5) may be difficult to interpret since the composition of the bone marrow cell population is heterologous and contains pluripotential cells and a big reservoir of granulocytes and myelocytic precursors in rodents [52], where the HL-60 cell line could serve as a representative cell type. Therefore, the changes in expression of signaling molecules may represent a resultant effect of RM33 on various cell populations contained in the bone marrow.
In the case of HL-60 cells, the increases in expression of MAP kinases were modest, except ERK-2 (5.2× increase), and thus are not conclusive (Table S7). Nevertheless, the overall effect of RM33 on signaling pathways in bone marrow cells suggests cell activation, since all members of MAP kinases were strongly stimulated (ERK-2 16×, p38β 138× and JNK 31×) ( Table 4). The changes in cell signaling were positively correlated with the increase in cell viability by 17.8% (Figure 4).
The discovery that RM33 affects production of COX-2 (Table 8) significantly contributed to understanding its immunosuppressive actions towards humoral and cellular immune responses in vivo [27] since PGE2 is involved in suppression of these responses [53]. Stimulation of COX-2 production by RM33 could also have a role in suppression of skin inflammatory reactions such as adjuvant inflammation or carrageenan reaction [28] where CCL27 chemokine is a mediator [54]. Nevertheless, at this stage of investigation, it is difficult to propose definite and complete mechanism of action of the compound.
In conclusion, given the anti-inflammatory properties of RM33, accompanied with the lack of cytotoxicity at high doses and its ability to regulate cyclooxygenase activities, the compound may find broad potential application in therapy of inflammatory disorders. Our unpublished data revealed a good potency of RM33, administered topically in ointment, in decreasing manifestations of contact sensitivity to oxazolone, indicating a possibility of treating skin inflammatory disorders.

Chemistry
Melting points were determined on the Büchi M560 melting point apparatus (BÜCHI Labortechnik AG, Meierseggstrasse 40, CH-9230 Flawil, Sankt Gallen, Switzerland) and were uncorrected. Thin-layer chromatography method (TLC) was applied to monitor the reaction progress as well as to confirm the purity of the obtained compounds. Polygram SIL G/UV254 plates (Mocherey-Nagel, Düren, Germany) for TLC were used, eluting medium was chloroform/methanol (9:1), detection of the compounds on the chromatograms was carried out with UV light. The Attenuated Total Reflectance IR (ATR-FT-IR) spectra (4000-450 cm −1 ) were recorded on a Nicolet iS50 FT-IR spectrophotometer (Thermo Fisher Scientific Inc., Waltham, MA, USA) using clean solid forms of the compounds (on the diamond crystal surface, 32 scans, resolution: 1 cm −1 , measurement temperature: 20-25 • C). IR Spectra were recorded with ATR intensity correction. It is known that the relative intensity of bands in an ATR spectrum increases with wavelength, causing the distortion of relative peak intensities in comparison with the classical transmission experiment [55]. Omnic Specta software (Thermo Fisher Scientific Inc., Waltham, MA, USA) was used for IR spectra analysis. Baselines of all IR spectra were corrected with autocorrection (fit order 2) and then all the spectra were normalized. Frequencies are reported in cm −1 . The samples were applied as solids. 1  150 mg of RM33 was refluxed with 5 mL of deuterated methanol (CD 3 OD) in a moisture-proof flask for 48h. Upon cooling the solution to room temperature, the product crystallized, which was filtered off to obtain 45.1 mg of dry compound after thorough drying in a desiccator. M.p. = 197-200 • C (dec.) (measured in sealed capillary with the Büchi M560 melting point apparatus). In order to confirm the exchange of two NH protons to N-D deuterium, the 1 H NMR spectrum (in CD 3 OD) was performed, which confirmed the complete conversion of N-protons to N-deuteriums. Then, the ATR-FTIR spectrum (with ATR correction) of the isotopologue RM33D sample was measured. ATR-FTIR (T = 298 K): ν max (cm −1 ) 1603 (CO stretching vibration) (Figures S2-S4). 1 Figure S10). 13 C NMR (CD 3 OD, 75.5 MHz, T = 299.6 K): δ (ppm) 12.1 (carbon of methyl group at position 3 of the isoxazolotriazepine system), 21.2 (carbon of methyl group at position 7 of the system, broadened signal due to coupling with neighboring deuterium atom), 38.2 (carbon of N-methyl group at position 5 of the system), 69.1 (C7 carbon at position 7 of the system broadened signal due to C-D coupling), 89.5 (quaternary bridgehead carbon C3a of the system), 163.4 (quaternary C3 carbon at position 3 of the system), 168.3 (quaternary C5 carbon at position 5 of the system), 168.8 (quaternary carbonyl carbon at position 4 of the system) ( Figure S11).

Computational Details about Theoretical Prediction of Chemical Shifts of 1 H NMR and 13 C NMR Spectra of RM33
In the beginning, the geometry optimization of RM33 and tetramethylsilane (TMS) molecules was performed on the basis of ab initio quantum mechanical DFT (density functional theory) method using the B3LYP [56,57] hybrid density functional with the following basis functions: Pople type basis sets augmented with diffuse functions [58][59][60], and Dunning's correlation consistent basis set-augmented with diffuse functions aug-cc-pVTZ [61]. The Conductor-like Polarizable Continuum Model (CPCM) [62] was used for simulation of DMSO environment with Gaussian 16 default set of parameters (scaling factor for van der Waals radius for all atoms-i.e., alpha = 1.1). All the calculations were computed using the Gaussian 2016 revision C.01 software [63]. To possibly get most accurate approximation the following keywords, i.e., Fopt = (Tight,CalcAll), SCF = (Direct,VeryTight) and Integral(Grid = SuperFineGrid) were used in the calculations. The absence of imaginary wavenumbers in the result of the frequency calculations (CalcAll or Freq keywords) proved a stationary point found (a minimum on the potential energy surface). The calculations of isotropic chemical shielding for 1 H and 13 C nuclei and then chemical shifts were performed similarly to a way described in [64,65]. Isotropic chemical shielding of NMR were calculated using Gauge-Independent Atomic Orbital (GIAO) method [66] with following combinations functional/basis set: B3LYP/6-311+G(d,p), B3LYP/6-311++G(df,pd) and B3LYP/aug-cc-pVTZ with application of CPCM [61] with default prarameters for simulation of solvent environment. Following parameters were used nmr (All, PrintEigenvectors, mixed, Susceptibility). Nondeuterated DMSO, as the solvent was set and its appropriate parameters (such as electric permittivity) used instead of its deuterated form in CPCM calculations. TMS was set as a reference. Hence, chemical shifts were calculated for 1 H and 13 C NMR. The calculated shifts and comparison with experimental values are collected in Table 1. The results show that the least deviations between calculated and experimental values are obtained for combination of B3LYP functional and medium 6-31+G(d,p) basis set. Larger basis sets, such as 6-311++G(df,pd) and aug-cc-pVTZ, gave worse results in this case.  Table S10. RM33 was synthesized as described in Section 4.1.1. RM33 was initially dissolved in DMSO (1 mg/100 µL) and subsequently in the culture medium.

Propagation of Cell Lines
J774E cells were maintained in a culture medium consisting of RPMI-1640 containing 10% FCS and antibiotics, at a density of 2 × 10 5 cells/mL. HL-60 cells were maintained in Iscove's Medium Dulbecco's Medium containing 10% FCS and antibiotics, at a density of 1 × 10 5 cells/mL. Jurkat cells were cultured in a culture medium consisting of RPMI-1640 and 10% FCS and antibiotics, at a density of 1 × 10 5 cells/mL. WEHI 231 cells were cultured in Dulbecco's Modified Eagle's Medium (DMEM) and 10% FCS and antibiotics, at a density of 1 × 10 5 cells/mL. The cell cultures were maintained at 37 • C in 5% CO 2 and 95% humidified atmosphere.

Cell Toxicity Test
RM33 was initially dissolved in DMSO, further dilutions of the compound were performed in RPMI-1640 medium supplemented with 2% FCS. RM33 was tested at 0.78-200 µg/mL concentration range. The evaluation of a potential cytotoxic action of RM33 was performed in a monolayer culture of epithelial lung cancer cell line A549. The cells at density of 5 × 10 4 /well were incubated for 24 h in a cell culture incubator. After the incubation, the culture supernatants were removed and to the monolayer cultures of cells appropriate dilutions of RM33 in the culture medium (200 µL/well) were added and incubated for additional 48 h. Control cultures contained corresponding dilutions of DMSO. The cell viability was determined by MTT colorimetric assay [36].

Isolation of Cells from the Lymphoid Organs
Mice were sacrificed by cervical dislocation. Thymuses, spleens and femurs were isolated and placed in disposable Petri dishes containing sterile, ice-cold PBS. The cells were released from the lymphatic organs by passing them through a nylon mesh and separating by centrifugation (200× g). The cells were then collected, washed with icecold PBS supplemented with 1% BSA and re-suspended in culture medium at density of 1 × 10 6 cells/mL. The viability of the cell suspension, determined by trypan blue dye exclusion assay, was 90-98%.

Evaluation of Viability of Cells from the Lymphoid Organs
The cells, isolated as described above, were diluted 20× in 0.1% trypan blue and counted in a hemocytometer. The number of viable (unstained) cells was determined.

Isolation of Peritoneal Exudates Cells and Cytokines Determination
Rats were sacrificed in a CO 2 atmosphere. The peritoneal cavities were washed with 20 mL of cold Hanks' medium containing 5 units/mL of heparin. The cells were centrifuged at 200× g for 10 min, washed once with Hanks' medium and counted in a hemocytometer. The cells were re-suspended in the culture medium in 1 mL aliquots in 24-well culture plates (Nunc) containing 5 × 10 6 cells/mL. The cells were incubated overnight. LPS was applied at concentration of 1 µg/mL and RM33 at 1 and 5 µg/mL. Control cultures contained appropriate dilutions of DMSO, corresponding to DMSO dilutions at respective RM33 doses. Determination of TNFα and IL-6 were performed by bioassays [38,39].

Proliferation Tests
The spleens were pressed against a plastic screen into 0.83% NH 4 Cl solution to lyze erythrocytes (5 min incubation at room temperature). The cells were then washed twice with Hanks' medium, passed through a glass wool column to remove debris, and re-suspended in the culture medium, referred to below as the culture medium, consisting of RPMI-1640, supplemented with 10% FCS, L-glutamine, sodium pyruvate, 2-mercaptoethanol and antibiotics. The cells were then distributed into 96-well flat-bottom tissue culture plates (Nunc) at a density of 2 × 10 5 /well. 2.5 µg/mL ConA was added to induce cell proliferation. To induce B-cell proliferation LPS at concentration of 50 µg/mL was applied. RM33 was added to the cultures at doses of 1.25, 2.5 and 5 µg/mL. After a 3-day incubation, the cell proliferation was determined using the colorimetric MTT assay (see Section 4.2.10). The results are presented as the mean optical density (OD) at 550/630 nm ± standard error (SE) from quadruplicate determinations.

Colorimetric MTT Assay
The cell proliferation was determined using the colorimetric MTT assay according to Hansen et al. [36]. Briefly, 25 µL of 3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazoliumbromide (MTT, Sigma-Aldrich, St. Louis, MO, USA). From stock solution (5 mg/mL), 25 µL was added per well at the end of cell incubation period and the plates were incubated for additional 3 h in a cell culture incubator. Then, 100 µL of the extraction buffer (20% SDS with 50% DMF, pH 4.7) was added. After an overnight incubation, the OD was measured at 550 nm with the reference wavelength of 630 nm.

Determination of Apoptosis
Cells from the lymphoid organs were isolated as described above and placed in 24-well plates at a concentration of 10 6 /well in the culture medium. The cells were cultured overnight with RM33 (2, 10 or 50 µg/mL) and appropriate dilutions of the solvent (DMSO). Next, the nonadherent cells were transferred to FACS tubes and the adherent cells were treated for 30 min with 100 µL of 0.1% trypsin solution containing EDTA, pH 8.0, to regain remaining adherent cells. Adherent cells were then washed with PBS containing 2.5% FCS, centrifuged and added to nonadherent cells, centrifuged, and the cell pellets were eventually re-suspended in 1 mL of PBS containing 2.5% FCS and placed in an ice bath. For DNA fragmentation evaluation, cells were fixed with 70% ethanol (POCh, Gliwice, Poland) and subsequently stained with propidium iodide (50 µg/mL) and RNA-se (0.02 mg/mL) (Sigma-Aldrich) as described [67].

Total RNA Isolation
Total RNA isolation was performed with TRIzol Reagent (Ambion) accordingly to manufacturer's recommendations. The cell pellet (2 × 10 6 cells) was re-suspended in 1 mL of TRIzol reagent, shaken, incubated for 10 min at room temperature (RT), supplemented with 0.2 mL of chloroform, shaken vigorously for 15 s, incubated for 3 min at RT and centrifuged at 12,000× g for 15 min at 4 • C. The water phase was collected, transferred to a new tube, supplemented with 0.5 mL of isopropanol, incubated at RT for 10 min and centrifuged at 12,000× g for 10 min at 4 • C. The RNA pellet was washed with 1 mL of 75% ethanol, dried in air and dissolved in 20-30 µL of sterile diethylpyrocarbonate-treated Mili-Q water. RNA samples were stored at −20 • C.

Reverse Transcription
Single stranded complementary DNA (cDNA) was synthesized with oligo (dT)12-18 primers from 5 µg of total RNA using Novazym VerteKit, according to the manufacturer's instruction. The list of primers used in the study is presented in Table S10.

Quantitative Analysis of Gene Expression by Real Time PCR
The expression of studied genes, i.e., caspases 3, 7, 8 and 9, Bcl-2, Fas, NFκB1 and p53, was determined with AmpliQ 5× HOT EvaGreen ® qPCR Mix Plus (noROX) (Novazym). The change in expression of each particular gene in samples treated with RM33 was calculated by comparison to such expression in control culture after earlier normalization of total mRNA content by GADPH expression.

Determination of Cyclooxygenases in LPS-Stimulated Mouse Splenocytes
Splenocytes from DBA/1 mice were incubated in the culture medium for 72 h in 24-well culture plates and a density of 5 × 10 6 /mL. The cells were stimulated with LPS (20 µg/mL) and RM33 was used at 10 µg/mL. The control cultures contained concentration of DMSO in the culture medium corresponding to the concentration of DMSO in cultures containing RM33.
After the incubation, the cells were detached by treatment with 0.1% trypsin solution for 10 min followed by 3× wash in PBS. The dry cell pellets were suspended in 0.5 mL of lysis buffer consisting of 1 mM EDTA, 0.5% Triton X-100 and containing 10 mg/mL leupeptin, 10 mg/mL pepstatin, and 3 mg/mL aprotinin in PBS, pH 7.4. Lysis was carried out on the ice for 15 min. Lysates were preserved by being deep-frozen (−80 • C). Before use, the lysate samples were centrifuged at 2000× g for 5 min, and the supernatant was transferred to a new tube. Sample total protein concentration was determined by the spectrophotometric method for equilibrating lysates dilutions. The content of cyclooxygenases in cells was determined by ELISA kits (Human/Mouse Total COX-2 DuoSet IC ELISA, R&D Systems, Minneapolis, MN, USA) according to the manufacturer's instructions.

Statistics
Each experimental group consisted of four wells (determinations) for in vitro tests. The results were subjected to statistical analysis using analysis of variance (one-way ANOVA) in STATISTICA 7 for Windows. Brown-Forsyth's test was used to determine the homogeneity of variance between groups. Due to nonconstant variance, the data were analyzed using the nonparametric the Kruskal-Wallis' analysis of variance, followed by Dunn's test to estimate the significance of the difference between groups. Significance was determined at p < 0.05. The results are presented as mean values ± SE.

Data Availability Statement:
The data presented in the current study are available from the corresponding author upon reasonable request.