Introduction
Although substituted N-arylbenzylimines have been investigated as antibacterial and antifungal molecules [1-3], little is known about their possible mechanism of action for inhibiting bacterial growth. According to the World Health Organization (WHO), it is important to design and develop compounds directed at specific targets in bacteria and determine the antibacterial mechanism of action involved [4]. Thus, the priority is to find molecules that can inhibit the bacteria most harmful to the world population and that do so with mechanisms of action distinct from those of current clinical drugs [5-8].
The bacteria responsible for the most lethal communicable diseases are Acinetobacter, Pseudomonas, Klebsiella, Escherichia coli, Serratia, Staphylococcus aureus and Proteus [9, 10]. They trigger local and systemic diseases such as soft-tissue infections, pneumonia, meningitis, sepsis, septic shock, and multiple organ failure. Multidrug-resistant Pseudomonas aeruginosa, responsible for 70% of invasive infections, induces hospital infections and sepsis [11].
An imine, or Schiff base, is formed by the reaction of ammonia or a primary amine with an aldehyde or ketone [12]. Similar to alkenes, Schiff bases can express syn/anti ( Z/E) isomerism in the C=N double bond [13, 14], which is a useful property for catalytic processes [15, 16], organic synthesis, and the formation of metal complexes [17].
The aim of the current study was to synthesize seven N-arylbenzylimines, evaluate their antibacterial activity in vitro on P. aeruginosa, and carry out a quantitative structure- activity relationship (QSAR) analysis, considering lipophilicity as well as the electronic and steric effects proposed by Hansch-Fujita. The N-arylbenzylimines synthesized and tested against P. aeruginosa are structurally very similar. The common skeleton of N-aryl-arylmethyleneimines [18-23] was strategically substituted in the para position with different electron donating and withdrawing groups, and then characterized spectroscopically. NMR spectra and X-ray diffraction showed that all seven synthesized compounds had an E-configuration. Biological activity was examined in vitro by the agar diffusion method (Kirby-Bauer), and the MIC was calculated with the Hill method.
Experimental
In silico study
Evaluation of the ADMET properties
By utilizing the admetSAR (version 2.0) [24, 25], SwissADME [26-29], SMARTCyp [30-32] and MarvinSketch 6.0.0 [33] bioinformatic server, the physicochemical ADMET properties (absorption, distribution, metabolism, excretion and toxicity) of the seven N-arylbenzylimines herein synthesized were theoretically calculated and related to their pharmacodynamic behavior. An in-depth comparison of the structural characteristics of the test compounds was made by running experimental values obtained from similar known compounds in the corresponding software, and acute and chronic toxicity were estimated [24, 25].
Materials and methods
All reagents and solvents were used as received from Sigma- Aldrich®. Thin-layer chromatography (TLC) was performed with Merck-DC-F254® silica gel on aluminum plates and revealed with a short-wave UV lamp at λ= 254 nm. Melting points were determined on a Mel- Temp apparatus (Laboratory Devices Inc., USA). IR spectra were recorded on a Pike Technologies GladiATR spectrophotometer and UV- Vis spectra on a Beckman Coulter DU 650 spectrophotometer. 1H, 13C, NOE and HSCQ NMR spectra were acquired on a Varian VNMRS-500 spectrometer at 500 MHz and a Bruker Avance spectrometer at 600 MHz (125.8 and 150.9 MHz for 13C, respectively) in CDCl3, using TMS as internal reference. Chemical shifts (δ) are reported in parts per million (ppm). Coupling constants (J) are given in Hz. High-resolution electrospray ionization mass spectrometry (ESI-MS) was performed with an Agilent 6545 Q -TOF LC/MS instrument. Chemical purity was examined by high-performance liquid chromatography (HPLC) with an Agilent 1200 Infinity Series apparatus.
Synthesis and characterization
Equimolar amounts of the corresponding para -substituted anilines (1a-g) and benzaldehyde (2) diluted in 60 mL of ethanol were put in a round-bottom flask. The reaction was carried out in an acidic medium (with concentrated HCl as the catalyst) and left under reflux for 12 h. Subsequently, recrystallization was achieved by evaporating the solvent and placing the reaction crude in an ice bath [34-35]. The product was filtered under vacuum and purified by sublimation (Scheme 1).
(E)-1-(N-phenyl)-1-phenylmethanimine (3a)
A yellow solid was obtained in 80% yield; mp 48 -50 °C; Rf 0.86 (n-hexane/AcOEt 7: 3); UV-Vis (MeOH) λmax = 267.8 ± 3.5 nm, ε = 135500 mol-1·dm3·cm-1; IR = 1574.97 (-C=N-), 3059.43, 1480.92, 1448.20 (aromatic) cm-1; 1H NMR (CDCl3, 600 MHz,): δ= 8.53 (- HC=N- , 1H, s), 7.99 (H2', H6', 2H, d, 3 JH-H = 8.4 Hz), 7.56-7.55 (H3', H4', H5', 3H, m), 7.48 (H2, H6, 2H, d, 3 JH-H = 7.8 Hz), 7.32-7.29 (H3, H4, H5, 3H, m); 13C NMR (CDCl3, 150 MHz,): δ = 160.4 (-HC=N-), 152.1 (C-1), 136.2 (C-1’), 131.4 (C- 4’), 129.1, (C-2, C-6), 128.8 (C-2’, C-6’), 128.8 (C-3’, C-5’), 125.9 (C-4), 120.9 (C-3, C-5); EIMS m/z [M+H]+ 182.09 (calcd. 181.09).
(E)-1-[N-(4-iodophenyl)]-1-phenylmethanimine (3b)
bright grey solid was provided in 40% yield; mp 74-76 °C; Rf 0.50 (n-hexane/AcOEt 7:3); UV-Vis (MeOH): λmax = 269.25 ± 2.85 nm, ε= 2002900 mol-1·dm3·cm-1; IR = 1569.90 (-C=N-), 3053.87, 1620.94, 1475.68, 1452.13 (aromatic), 529.55 (-C-I), 1958.53, 1903.60, 823.99 (para-substituted) cm-1; 1H NMR (CDCl3, 600 MHz,): δ= 8.47 (-HC=N-, 1H, s), 7.95 (H2', H6', 2H, d, 3 JH-H = 7.80 Hz), 7.56-7.52 (H3', H4', H5', 3H, m), 7.76 (H3, H5, 2H, d, 3 JH-H = 7.80 Hz), 7.02 (H2, H6, 2H, d, 3 JH-H = 8.40 Hz); 13C NMR (CDCl3, 150 MHz,): 160.8 (-HC=N-), 151.7 (C-1), 138.1 (C-3, C-5), 137.9 (C-1’), 135.9 (C-4), 131.7 (C-4’), 128.9 (C-2’, C-6’), 128.8 (C-3’, C-5’), 123.0 (C-2, C-6); EIMS m/z [M+H]+ 307.99 (calcd. 306.99).
(E)-4-(N-methylenephenylimino)benzoic acid (3c)
white solid was afforded in 60% yield; mp 194-196 °C; Rf 0.63 (n-hexane/AcOEt 7:3); UV-Vis (MeOH): λmax = 274.3 ± 0.5 nm, ε= 161740 mol-1·dm3·cm-1; IR = 1677.21 (-C-O- acid), 1574.97 (-C=N-), 3063.92, 1595.42, 1427.65 (aromatic), 851.15 (para-substituted) cm-1; 1H NMR (DMSO, 500 MHz,): δ= 12.68 (-COOH, 1H, br), 8.63 (-HC=N-, 1H, s), 7.98 (H3, H5, 2H, d, 3 JH-H = 9.0 Hz), 7.95 (H2', H6', 2H, d, 3 JH-H = 8.0 Hz), 7.57-7.51 (H3', H4', H5', 3H, m), 7.316 (H2, H6, 2H, d, 3 JH-H = 8.50 Hz); 13C NMR (DMSO, 125 MHz,): δ=167.9 (-COOH), 167.4 (-HC=N-), 155.8 (C-1), 136.1 (C-1’), 135.0 (C-3, C-5), 131.0 (C-4’), 129.3 (C-2’, C-6’), 129.3 (C-3’, C-5’), 128.4 (C-2, C-6), 121.4 (C-4); EIMS m/z [M+H]+ 226.08 (calcd. 225.08).
(E)-1-[N-(4-bromophenyl)]-1-phenylmethanimine (3d)
bright grey solid was produced in 50% yield; 66-68 °C; Rf 0.80 (n-hexane/AcOEt 7:3); UV-Vis (MeOH): λmax = 268.8 ± 1.8 nm, ε= 259520 mol-1·dm3·cm-1; IR = 1574.19 (-C=N-), 1619.96. 1480.92, 1399.13 (aromatic), 523.99 (-C-Br), 1888.62, 818.43 (para-substituted) cm-1; 1H NMR (CDCl3, 600 MHz,): δ=8.48 (-HC=N-, 1H, s), 7.95 (H2', H6', 2H, d, 3 JH-H = 7.8 Hz), 7.56-7.54 (H3', H4', H5', 3H, m), 7.56 (H3, H5, 2H, d, 3 JH-H=8.4 Hz), 7.15 (H2, H6, 2H, d, 3 JH-H = 9.0 Hz); 13C NMR (CDCl3, 150 MHz,) δ= 160.7 (-HC=N-), 151.0 (C-1), 135.9 (C-1’), 132.2 (C-3, C-5), 131.6 (C-4’), 128.9 (C-2’, C-6’), 128.8 (C-3’, C-5’), 122.6 (C-2, C-6), 119.3 (C-4); EIMS m/z [M+H]+ 260.00 (calcd. 259.00).
(E)-1-[N-(4-chlorophenyl)]-1-phenylmethanimine (3e)
bright yellow solid was furnished in 50% yield; 64-66 °C; Rf 0.83 (n-hexane/AcOEt 7:3); UV-Vis (MeOH): λmax = 268.5 ± 3.9 nm, ε= 661250 mol-1·dm3·cm-1; IR = 1570.88 (-C=N-), 1619.96, 1468.65, 1452.29, 1399.13 (aromatic), 683.48 (-C-Cl), 1889.86, 818.43 (para-substituted) cm-1; 1H NMR (CDCl3, 500 MHz,) δ= 8.34 (-HC=N-, 1H, s), 7.806 (H2', H6', 2H, d, 3 JH-H = 8.0 Hz), 7.40-7.38 (H3', H4', H5', 3H, m), 7.26 (H3, H5, 2H, d, 3 JH-H = 9.0 Hz), 7.06 (H2, H6, 2H, d, 3 JH-H = 9.0 Hz); 13C NMR (CDCl3, 125 MHz,) δ= 160.9 (-HC=N-), 150.7 (C-1), 136.1 (C-1’), 131.8 (C-4’), 131.6 (C-4), 129.4 (C-3, C-5), 129.1 (C-2’, C-6’), 129.0 (C-3’, C-5‘), 122.1 (C-2, C-6); EIMS m/z [M+H]+ 216.05 (calcd. 215.05).
(E)-1-[N-(4-fluorophenyl)]-1-phenylmethanimine (3f)
bright white solid was obtained in 50% yield; 54-56°C; Rf 0.90 (n-hexane/AcOEt 7:3); UV-Vis (MeOH): λmax = 266.2 ± 1.6 nm, ε= 150470 mol-1·dm3·cm-1; IR= 1570.88 (-C=N-), 1619.96, 1485.01, 1448.20 (aromatic), 753.00 (-C-F), 830.70 (para-substituted) cm-1; 1H NMR (CDCl3, 600 MHz,): δ= 8.53 (-HC=N-, 1H, s), 8.00-7.98 (H2', H6', 2H, m), 7.56-7.55 (H3', H4', H5', 3H, m), 7.30 (H2, H6, 2H, d, 3 JH-H = 7.0 Hz), 7.47 (H3, H5, 2H, d, 3 JH-H = 7.0 Hz); 13C NMR (CDCl3, 150 MHz,): δ= 160.4 (-HC=N-), 152.1 (C-1), 136.2 (C-1’), 131.4 (C-4’),129.1 (C-2’, C-6’), 128.80 (C-3’, C-5’), 128.85 (C-2, C-6), 125.9 (C-4), 120.9 (C-3, C-5); EIMS m/z [M+H]+ 200.08 (calcd. 199.22).
(E)-1-[N-(4-methoxyphenyl)]-1-phenylmethanimine (3g)
bright grey solid was provided in 75% yield; 70-72 °C; Rf 0.76 (n-hexane/AcOEt 7:3); UV-Vis (MeOH): λmax = 285.3 ± 49.75 nm, ε= 132680 mol-1·dm3·cm-1; IR = 2953.11 (CH3, -CH2-), 1574.97 (-C=N-), 1497.27 (aromatic), 1104.69 (-C-O-CH3), 834.79 (para-substituted) cm-1; 1H NMR (CDCl3, 600 MHz,): δ= 8.55 (-HC=N-, 1H, s), 7.97-7.96 (H2', H6', 2H, m), 7.53-7.53 (H3', H4', H5', 3H, m), 7.32 (H2, H6, 2H, d, 3 JH-H = 9.0 Hz), 7.01 (H3, H5, 2H, d, 3 JH-H = 8.4 Hz), 3.89 (-CH3, 3H, s); 13C NMR (CDCl3, 150 MHz,): δ= 158.3 (-HC=N-), 144.9 (C-1), 136.4 (C-1’), 131.1 (C-4’),128.7 (C-2’, C-6’), 128.6 (C-3’, C-5’), 122.3 (C-2, C-6), 116.4 (C-4), 114.4 (C-3, C-5), 55.5 (-OCH3); EIMS m/z [M+H]+ 212.10 (calcd. 211.10).
Single-crystal X-ray crystallography
Recrystallization of N-arylbenzylimine 3f from CH2Cl2/n-hexane led to the formation of colorless crystals, one of which was mounted on a glass fiber. Crystallographic measurements were carried out at room temperature on an Oxford Diffraction Xcalibur S diffractometer with Mo Ka radiation (graphite monochromator, λmax= 0.71073 Å) and a CCD detector. Empirical absorption corrections were applied. The structure was solved with SHELXT and refined with SHELXL [36], both running within the WinGX environment [37]. Anisotropic temperature factors were introduced for all non -hydrogen atoms. Hydrogen atoms were placed in idealized positions, and their atomic coordinates refined by using unit weights. The plot in Fig. 2 was made with PLATON [38, 39]. Data for N-arylbenzylimine 3f are summarized in Table 1. Supplementary crystallographic data have been deposited at the Cambridge Crystallographic Data Centre (CCDC 2046921) [40].
Empirical formula | C13H10FN |
Molecular weight | 199.22 |
Crystal size (mm) | 0.49×0.26×0.16 |
Crystal system | Monoclinic |
Space group | P 21/n |
Unit cell parameters (Å, °) | a = 5.7528(3), α = 90 b = 25.1417(11), β = 90.198(5) c = 7.2892(4), γ = 90 |
Volume (Å3) | 1054.27(9) |
Z | 4 |
Density (calcd, Mg/m3) | 1.255 |
Absorption coefficient (µ, mm−1) | 0.086 |
Ɵ range for data collection (°) | 3.231−32.461 |
Reflections collected | 6489 |
Independent reflections | 3203 |
Observed reflections | 1932 |
Goodness-of-fit on F2 | 1.028 |
Final R indices [I > 2σ(I)] | R1 = 0.0610, wR 2 = 0.1236 |
General method for the preparation of the solutions: positive controls and (E) N-arylbenzylimines
The solutions of the positive control were prepared by taking successive 1:10 aliquots from the stock solution (14.6743 mg aztreonam, 337 µM) in order to construct the standard curve and calculate the MIC [41-43]. Successive solutions of the N-arylbenzylimines were prepared to obtain the gradual concentration-response curves and calculate the MIC values. From the stock solution (4.0 mM), the dilutions were elaborated at 400, 200, 175, 150, 125, 100, 75, and 50 µM. The compounds were added according to the general method to furnish 4.0 mM, requiring 7.2496 mg of 3a, 15.3565 mg of 3b, 9.0100 mg of 3c, 13.0065 mg of 3d, 8.6272 mg of 3e, 7.9692 mg of 3f and 8.4504 mg of 3g. Methanol served as the negative control [41-43].
The N-arylbenzylimines derivatives showed stability under the conditions in which the experiments were conducted. The stock solutions did not show any change during preparation, whether in appearance, coloration, or the formation of precipitates. They remained translucent, which allowed for the preparation of successive dilutions with the same characteristics.
In vitro antibacterial activity
Bacterial strains
All the materials used for culturing the bacterium and conducting the in vitro microbiological assays were autoclaved at 121 °C and 15 psi for 20 min. Bacterial lyophilizate, taken from P. aeruginosa ATCC 27853, was activated in nutritive broth containing 1.0 L pancreatic digest comprised of 5.0 g gelatin and 3.0 g of an extract of bovine meat (at pH 7.3 ± 0.2). It was seeded by cross striation in the culture medium of trypticase soy agar (TSA) containing 1.0 L pancreatic digest constituted by 15.0 g casein and papaic digest of 5.0 g soy flour, 5.0 g NaCl and 15.0 g agar (at pH 7.3 ± 0.2). Primary biochemical evaluations (the gram stain, oxidase and catalase tests) were employed to assure the purity of the strains and confirm their identity [44].
Quantitative assessment of biological activity
Antibacterial activity was examined in vitro by the agar diffusion method (Kirby -Bauer). Live cells of P. aeruginosa were harvested from the TSA medium at 18-24 h of growth. They were suspended in 0.85% (m/V) saline with a density of 1.5x106 colony forming units (CFU)/mL (0.5 McFarland units), measured at an absorbance of 625 nm in a range of 0.08- 0.10 AU on a Model CS-200 PC Globe spectrophotometer and a Model Vortex-Mixer Fisher Scientific shaker. Subsequently, this suspension (1.5x106 CFU/mL) was seeded with sterile swabs in the Müeller-Hinton agar culture medium, which was composed of 2.0 g beef extract, 17.5 g hydrolyzed casein acid, 1.5 g starch and 17.0 g agar (at pH 7.3 ± 0.2).
After placing 6 Sensi discs in each Petri dish, 10 µL of one of the compounds was added (positive control, negative control, and 50, 75, 100, 125, 150, 175, 200, and 400 µM of the corresponding N-arylbenzylimine) . The plates were inverted and incubated for 24 h at 37 °C, and then the inhibition halos were measured with a Vernier [41-43]. The MICs were determined for the reference and test compounds [45-47]. The aztreonam-induced inhibition was taken as 100%, and the inhibitory activity of the test compounds was expressed as a relative percentage. Each N-arylbenzylimine was examined by constructing a gradual concentration-response curve. MIC values, calculated by the Hill method, were considered equivalent to the CE95 value [48], according to the following equation:
The QSAR was performed under QSAR- 2D, proposed by Corwin Hansch and Toshio Fujita [49, 50]. Concerning the physicochemical properties of the N-arylbenzylimines, size was measured by molar refractivity (Advanced Chemistry Development, Inc.) and lipophilicity by the partition coefficient [51, 52]. These parameters were determined on ACD/ChemSketch 2015 and CS ChemDraw Pro v.6 software, respectively [53].
The lipid solubility descriptor was established by the following equation:
where, π is the lipid solubility descriptor, and Px and PH are the partition coefficients of a given substituted compound and the reference drug, respectively.
The Hammett para substituent constant (σp) was utilized as the criterion of electronic effects [54, 55]. The steric descriptor was calculated by means of the following equation:
Where Es is the steric descriptor proposed by Taft [56, 57], and MRx and MRH are the molar refractivity of a given substituted compound and the reference drug, respectively.
Statistical analysis
Data were examined by using simple linear regression analysis with the least-squares technique and one -way ANOVA. The slopes were ordered to the origin and linear correlation coefficients were compared to the controls by the Student’s t-test. Significance was considered at p<0.05 in all cases. The QSAR model was analyzed by utilizing multiple regression analysis with the determinant technique and one-way ANOVA. Statistical evaluations were run on Sigma Stat 4.0 (Systat Software Inc., San Jose, CA, USA). The external prediction of the QSAR models was assessed by the predictive squared correlation coefficient (Q2) [58], calculated by the following equation:
Where PRESS is the index related to the predictive residual sum of squares, and SS is the activity sum of squares.
Results and discussion
In silico study
Determination of the ADMET properties
The physicochemical characteristics of the seven N-arylbenzylimines (3a-g) are described by their ADMET properties (Table 2).
Physicochemical properties | Compound | ||||||
3a (-H) | 3b (-I) | (-COOH) | 3d (-Br) | 3e (-Cl) | 3f (-F) | 3g (-OCH3) | |
Molecular weight g/mol | 181.09 | 306.99 | 225.08 | 259.00 | 215.05 | 199.22 | 211.10 |
Num. heavy atoms | 14 | 15 | 17 | 15 | 15 | 15 | 16 |
Num. aromatic heavy atoms | 12 | 12 | 12 | 12 | 12 | 12 | 12 |
Fraction Csp3 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.07 |
Num. rotatable bonds | 2 | 2 | 3 | 2 | 2 | 2 | 3 |
Num. H-bond acceptors | 1 | 1 | 3 | 1 | 1 | 2 | 2 |
Num. H-bond donors | 0 | 0 | 1 | 0 | 0 | 0 | 0 |
Molar Refractivity | 60.14 | 72.86 | 67.10 | 67.84 | 65.15 | 60.10 | 66.63 |
TPSA Ų | 12.36 | 12.36 | 49.66 | 12.36 | 12.36 | 12.36 | 21.59 |
Lipophilicity | |||||||
Log Po/w (iLOGP) | 2.50 | 2.86 | 2.06 | 2.90 | 2.77 | 2.63 | 2.81 |
Log Po/w (XLOGP3) | 2.82 | 3.47 | 2.35 | 3.51 | 3.45 | 2.92 | 2.79 |
Log Po/w (WLOGP) | 3.44 | 4.04 | 3.14 | 4.20 | 4.09 | 4.00 | 3.45 |
Log Po/w (MLOGP) | 3.26 | 4.05 | 2.72 | 3.93 | 3.80 | 3.67 | 2.87 |
Log Po/w (SILICOS- IT) | 3.90 | 4.83 | 3.26 | 4.54 | 4.50 | 4.28 | 3.88 |
Consensus Log Po/w | 3.18 | 3.85 | 2.70 | 3.81 | 3.72 | 3.50 | 3.16 |
ADMET predicted profile A (Absorption) | |||||||
Blood-brain barrier (BBB) penetration | 0.9962 | 0.9926 | 0.9434 | 0.9926 | 0.9911 | 0.9918 | 0.9775 |
Human intestinal absorption (HIA) | 0.9730 | 0.9319 | 0.9474 | 0.9649 | 0.9749 | 0.9734 | 0.9825 |
Human oral bioavailability (HOB) | 0.5857 | 0.5000 | 0.6571 | 0.5429 | 0.5857 | 0.5714 | 0.5000 |
Caco-2-permeability | 0.9070 | 0.8085 | 0.8062 | 0.8494 | 0.8565 | 0.9207 | 0.9367 |
P-glycoprotein inhibitor | 0.9815 | 0.9774 | 0.9646 | 0.9670 | 0.9682 | 0.9692 | 0.9365 |
P-glycoprotein substrate | 0.9879 | 0.9761 | 0.9822 | 0.9840 | 0.9899 | 0.9560 | 0.9695 |
D (Distribution) | |||||||
Subcellular location: Mitochondria | 0.4784 | 0.5888 | 0.6476 | 0.5111 | 0.5950 | 0.6152 | 0.6539 |
M (Metabolism) Cytochrome P450 (CYP450) substrate, inhibitor | |||||||
CYP3A4 substrate | 0.8087 | 0.7411 | 0.7660 | 0.7322 | 0.6478 | 0.7202 | 0.6432 |
CYP2C9 substrate | 0.5856 | 0.5878 | 0.8032 | 0.5878 | 0.6261 | 0.6017 | 0.5878 |
CYP2D6 substrate | 0.7217 | 0.7395 | 0.8785 | 0.7383 | 0.7585 | 0.7344 | 0.6843 |
CYP3A4 inhibitor | 0.9568 | 0.9372 | 0.9578 | 0.9518 | 0.9392 | 0.9443 | 0.7476 |
CYP2C9 inhibitor | 0.9269 | 0.9077 | 0.9252 | 0.8505 | 0.9106 | 0.9100 | 0.8828 |
CYP2C19 inhibitor | 0.6873 | 0.7239 | 0.8843 | 0.7872 | 0.8198 | 0.6449 | 0.5612 |
CYP2D6 inhibitor | 0.6535 | 0.5092 | 0.8502 | 0.6208 | 0.6064 | 0.6950 | 0.8435 |
CYP1A2 inhibitor | 0.7705 | 0.9484 | 0.6455 | 0.9158 | 0.9414 | 0.9047 | 0.8552 |
CYP inhibitory promiscuity | 0.6829 | 0.6997 | 0.9130 | 0.7371 | 0.7180 | 0.6651 | 0.7467 |
Drug transporters in pharmacokinetics | |||||||
OATP2B1 inhibitor | 1.0000 | 1.0000 | 1.0000 | 1.0000 | 1.0000 | 1.0000 | 1.0000 |
OATP1B1 inhibitor | 0.9430 | 0.9016 | 0.9297 | 0.9038 | 0.9082 | 0.9274 | 0.9040 |
OATP1B3 inhibitor | 0.9528 | 0.9509 | 0.9464 | 0.9508 | 0.9516 | 0.9519 | 0.9631 |
MATE1 inhibitor | 1.0000 | 1.0000 | 0.9800 | 1.0000 | 0.9000 | 1.0000 | 0.9600 |
OCT2 inhibitor | 0.7500 | 0.7750 | 0.9750 | 0.7750 | 0.8000 | 0.7750 | 0.7250 |
BSEP inhibitor | 0.5759 | 0.6161 | 0.6979 | 0.6664 | 0.5721 | 0.5120 | 0.6124 |
T (Toxicity) Organ toxicity | |||||||
Hepatotoxicity | 0.6250 | 0.5250 | 0.7500 | 0.5750 | 0.6000 | 0.5750 | 0.6500 |
Human ether-a-go- go-related gene (hERG) inhibition | 0.5723 | 0.6390 | 0.7617 | 0.5099 | 0.4685 | 0.4220 | 0.6873 |
Eye corrosion | 0.9123 | 0.8324 | 0.9298 | 0.9387 | 0.9418 | 0.9294 | 0.9119 |
Eye irritation | 0.9966 | 0.8651 | 0.9599 | 0.9503 | 0.9048 | 0.9069 | 0.9702 |
Estrogen receptor binding | 0.8236 | 0.9015 | 0.9394 | 0.9052 | 0.9160 | 0.8763 | 0.8815 |
Androgen receptor binding | 0.8062 | 0.6558 | 0.5248 | 0.4929 | 0.5591 | 0.6429 | 0.7434 |
Thyroid receptor binding | 0.5858 | 0.5349 | 0.6916 | 0.5143 | 0.5160 | 0.5627 | 0.5395 |
Glucocorticoid receptor binding | 0.6120 | 0.4737 | 0.6554 | 0.6016 | 0.5485 | 0.5523 | 0.5807 |
Aromatase binding | 0.8636 | 0.8324 | 0.8750 | 0.8454 | 0.8233 | 0.7086 | 0.7297 |
PPARγ | 0.8093 | 0.7090 | 0.7858 | 0.9066 | 0.8889 | 0.7891 | 0.5638 |
Toxicogenomics | |||||||
Ames mutagenicity | 0.7100 | 0.6600 | 0.5400 | 0.6200 | 0.5300 | 0.7300 | 0.6500 |
Carcinogenicity (binary) | 0.5336 | 0.5336 | 0.6626 | 0.5194 | 0.5479 | 0.5336 | 0.6075 |
Carcinogenicity (ternary) | 0.5455 | 0.6711 | 0.6198 | 0.5808 | 0.6727 | 0.6390 | 0.5216 |
Micronuclear | 0.7299 | 0.7032 | 0.6900 | 0.7032 | 0.7632 | 0.6332 | 0.5600 |
Eco-toxicity | |||||||
Biodegradation | 0.9250 | 0.9500 | 0.8000 | 0.9250 | 0.9500 | 0.9500 | 0.9500 |
Honey bee toxicity | 0.6907 | 0.6535 | 0.7249 | 0.6372 | 0.6509 | 0.6431 | 0.5000 |
Aquatic toxicity (crustaceans) | 0.8600 | 0.8800 | 0.7100 | 0.8700 | 0.9200 | 0.8600 | 0.9100 |
Aquatic toxicity (fish) | 0.9489 | 0.9877 | 0.9416 | 0.9890 | 0.9876 | 0.9878 | 0.8363 |
ADMET predicted profile-Regressions | |||||||
Water solubility (logS) | -3.803 | -4.633 | -3.153 | -4.515 | -4.471 | -3.657 | -3.889 |
Plasma protein binding | 0.439 | 0.731 | 0.680 | 0.650 | 0.732 | 0.583 | 0.554 |
Acute oral toxicity (mol/kg) | 1.914 | 1.866 | 1.866 | 2.187 | 2.088 | 2.123 | 1.729 |
Tetrahymena pyriformis pIGC50 (µg/L) | 0.559 | 1.826 | 0.496 | 1.875 | 1.684 | 1.322 | 1.144 |
Absorption: The degree of fat solubility of the N-arylbenzylimines allows them to cross intestinal cell membranes and therefore reach sufficient concentrations in the central compartment. On the other hand, these compounds do not undergo any change in their structure at pH values greater than 4.5, and therefore could be highly absorbable in the duodenum. The exception is compound 3c (-COOH), which would have a lower absorption rate because of being a conjugate base.
Distribution: Once the compounds reach the general circulation, they have a probability of binding to plasma proteins in the range of 43.0-74.0%. The lowest binding percentage is for 3a (-H ) and the highest for 3e (-Cl) . Thus, the free fraction of the N-arylbenzylimines should be adequately distributed in the body, since they have high enough log P values to be able to cross the cell membranes of some organs and tissues, such as the blood-brain barrier (BBB).
Metabolism: N-arylbenzylimines are susceptible to metabolic biotransformation by CYP450. Considering the oral route (the most common route of administration), these compounds undergo adequate absorption. However, they are susceptible to the first-pass effect, since their oral bioavailability (Foral) ranges from 0.50-0.66.
Overall, the N-arylbenzylimines show a greater tendency to inhibit the CYP isoforms 3A4, 2C9 and 2D6 than act as substrates of the same. The inhibition of CYP450 enzymes is important from the clinical point of view, because they represent a potential risk of drug-drug interactions, which occurs for example with some groups of antibiotics such as macrolides, fluoroquinolones and antiamoebic agents from the group of metronidazole [59].
In case N-arylbenzylimines are capable of acting as a substrate for CYP450, there is the potential risk of drug-drug interactions, as well as the formation of inactive, active and/or toxic metabolites. The greatest probability of biotransformation reactions mediated by CYP3A4 and CYP2D6 are aromatic hydroxylation and O-dealkylation. The former could possibly involve all derivatives (3a-g), while the latter would apply only to 3g (Table 3) . The greatest activity of the CYP3A4 isoform is on the disubstituted ring and the imine nitrogen, while that of the CYP2D6 isoform is the oxidation of the monosubstituted aromatic ring (Table 3).
Elimination: The oxidation of N-arylbenzylimines by CYP450 generates more polar molecules, which would mainly lead to their renal excretion due to their high water solubility and low molecular weight (M<500 g/mol).
According to the results of the in silico study, the N-arylbenzylimines do not violate any of the acceptability standards of Lipinski’s rule of five. Their fat solubility, size and polarity values indicate a capacity to easily cross cell membranes, such as the BBB and enterocytes, allowing them to reach plasma concentrations adequate for producing the desired pharmacological effect. The antibacterial activity found for the present N-arylbenzylimines is in agreement with their adsorption properties. The biological membranes they can cross include those of Gram-negative bacteria.
Toxicity: In relation to the acute toxicity of N-arylbenzylimines 3a-g, the LD50 are greater than 5000 mg/kg, indicating safe substances. Regarding toxicogenomics, these compounds have moderate activity. However, the harmful effects would be caused by chronic exposure.
Synthesis and characterization
Seven para-substituted N-arylbenzylimines were obtained in solid form with medium polarity and moderate to high yields (≥50%). As for the NMR determination of their relative configuration, the nuclear Overhauser effect (NOE) experiment conclusively demonstrated the E-configuration for all seven molecules. When the hydrogen corresponding to the formation of the N-arylbenzylimines was irradiated, 3b was displaced at approximately 8.4 ppm (Fig. 1 ). The hydrogens that undergo this displacement are H2 and H6, giving signals at 6.9 and 7.0 ppm, respectively.
After the selective irradiation of the imine hydrogen, there was a disturbance in the intensity of the signals of the H2 and H6 hydrogens, which indicates nuclear spin polarization transfer in the space between these hydrogens. Clearly, this could only occur for the E-configuration of the N-arylbenzylimines, which should be favored on the basis of the steric effects of the aromatic rings.
Single-crystal X-ray crystallography
The E-configuration of N-arylbenzylimine 3f was further supported by single-crystal X-ray diffraction analysis (Fig. 2).
In vitro antibacterial activity
Regarding macroscopic morphology, the cultured bacteria was seen in small colonies with irregular and translucent edges [60]. The microscopic morphology of the bacteria corresponds to isolated Gram-negative bacilli. The primary biochemical evaluations found were the following: Gram-, catalase+ and oxidase+ [61, 62].
On the other hand, most of the compounds (3a- e and g) generated good inhibitory activity on P. aeruginosa, in the range of 90-95% (relative to aztreonam). Since all the relevant N-arylbenzylimines displayed similar slopes of the concentration-response curve, they likely have the same mechanism of action. Contrarily, the slope of aztreonam was different (Fig. 3).
The activity of the test compounds on P. aeruginosa is shown by the MIC values (Table 4), which are lower for 3b-d (198.47, 249.60 and 290.20 µM, respectively) than aztreonam (337.39 µM).
ND: Not determined
N-arylbenzylimines with substituents generally produced a better effect than the lead compound with hydrogen (MIC = 790.10 µM) . The derivative bearing the -I substituent (3b) (MIC = 198.47 µ M) was the most active of the test compounds, giving another indication that the antibacterial mechanism of action of the N-arylbenzylimines is probably different from that of the monobactam control drug, known to inhibit penicillin-binding proteins [63-67]. On the other hand, the MIC value for N-arylbenzylimine with substituent -F (3f) could not be determined because activity was only observed at the highest levels (200 and 400 µM) of the concentration range selected (50-400 µM).
QSAR/SAR analysis
The QSAR analysis was conducted on the six relevant N-arylbenzylimines (3a-e and g) evaluated in vitro against P. aeruginosa. Their inhibitory activity was examined with the Hansch-Fujita model, which considers lipid solubility as well as electronic and steric effects. The substituents in molecules with the greatest antibacterial activity are bulky, electron-withdrawing and highly polar. Size was the most significant descriptor, as evidenced by the coefficient of 7.326.
Where π corresponds to lipid solubility, σ to the electronic effects in the para position, and Es to the steric effects. The model provided adequate predictability and a good correlation value on P. aeruginosa. The calculated slope is very close to 1.0, and the linear correlation coefficient is r =0.998 (Fig. 4).
The in vitro experiment with 3f afforded inhibition halos at concentrations of 200 µM and 400 µM. The size of the halos evidenced concentration dependence and suggested that the MIC value is probably greater than 400 µM. Interestingly, the general QSAR model predicts a MIC of 724.56 µM for 3f, a value in accordance with the present results.
For the SAR analysis, the compounds were subclassified into G1 with resonance effects and G2 with inductive effects, taking 3a as the point of reference. The values of the lipid solubility descriptor (π) for the compounds in G1 (3c and 3g) exhibited an inverse correlation with biological activity (the more polar substituents gave better antibacterial activity). For the compounds in G2 (3b, 3d and 3e), contrarily, a lower polarity of the substituent is associated with greater activity, resulting in a positive slope of the concentration-response curve (Fig. 5). Since 3f bears a fluorine atom with inductive effects (similar to the other halogens), it would definitely fit into the G2 group.
Regarding the descriptor of electronic effects, both the G1 and G2 sets showed a linear correlation with positive slope. In other words, the electron-withdrawing substituents of N-arylbenzylimines conferred greater antibacterial activity, with ρ values of 0.345 and 1.900, respectively. The difference in the magnitude of ρ evidences the sensitivity of the biological response to electronic effects. Better antibacterial activity was found for the inductive effects conferred by the halogens in G2 (3b, 3d and 3e) than the resonance effects of the other compounds in G1 (3c and 3g) (Fig. 6). The inductive effects of 3f would make it a halogenated compound belonging to the G2 group, though with the least antibacterial activity (MIC> 400 µM) of the molecules in this group.
Concerning the steric descriptor (Es), antibacterial activity was correlated with the size of the molecule for both G1 and G2. The data plotted for each group forms a straight line with a positive slope. Hence, the biological effect increased with the size of the molecule, and molecular size was dependent on the substituent. The slope was very similar for the G1 and G2 sets (8.1423 and 7.2107, respectively), indicating the great susceptibility of biological activity to steric effects (Fig. 7). The molecules with the smallest substituents were 3a (the lead compound) with a hydrogen atom and 3f with a fluorine atom, which coincides with their relatively low biological activity.
The following interpretation is provided for how the N-arylbenzylimines crossed the bacterial barriers of P. aeruginosa. For G1, better activity was found for the highly polar compounds, 3c and 3g (especially the former). The -COOH substituent in compound 3c generates a partial positive charge on the ring that delocalizes the positive charge (ρ=0.345), causing it to be repelled due to the positive charge on the outside of the outer membrane. Therefore, the molecule requires a transporter or a porin of the bacterium to pass through the membrane into the periplasmic space, where it strongly interacts with acidic peptidoglycan. For the G2 set (3b, 3d, 3e and 3f), antibacterial activity was favored by high lipid solubility, suggesting the penetration of the external, cytoplasmic and cell wall membranes by simple diffusion. The polarizability and the soft base properties of the halogens give better interaction and affinity with the phospholipids of the membranes and with the peptidoglycan layer, thus explaining why they can cross by simple diffusion (Fig. 8). Rather than delocalizing any type of charge, these compounds have inductive electron-withdrawing effects.
Conclusion
Seven N-arylbenzylimines were derived from para -substituted anilines and obtained in moderate to high yields with a simple synthetic methodology. They all showed an E-configuration, evidenced by NMR experiments and definitively confirmed by X-ray diffraction. These compounds displayed relevant inhibitory activity against P. aeruginosa compared to the corresponding reference drug, aztreonam. The microbiological data and the appraisal of the curves suggest that the mechanism of action is very likely the same for the seven N-arylbenzylimines, but different from that of the positive control. According to the QSAR analysis, biological activity is significantly dependent on steric effects, while lipid solubility and electronic effects play a vital role in the capacity of the molecules to cross the barriers of the bacterium. A greater capacity to penetrate these barriers was found for the inductive effect of halogenated compounds than the resonance effect of the other compounds. The current physicochemical results are very important for the design and development of new molecules directed against multi-drug resistant P. aeruginosa, a highly infective and lethal pathogen.