- Research article
- Open Access
Synthesis, structure and toxicity evaluation of ethanolamine nitro/chloronitrobenzoates: a combined experimental and theoretical study
© The Author(s) 2017
- Received: 25 April 2017
- Accepted: 10 November 2017
- Published: 6 December 2017
Nitroaromatic and chloronitroaromatic compounds have been a subject of great interest in industry and recently in medical-pharmaceutic field. 2-Chloro-4-nitro/2-chloro-5-nitrobenzoic acids and 4-nitrobenzoic acid are promising new agents for the treatment of main infectious killing diseases in the world: immunodeficiency diseases and tuberculosis.
New ethanolamine nitro/chloronitrobenzoates were synthesized and characterized by X-ray crystallography, UV–vis, FT-IR and elementary analysis techniques. The toxicity of the compounds prepared and correspondent components was evaluated using Hydractinia echinata as test system. A significant lower toxicity was observed for nitro-derivative compared with chloronitro-derivatives and individual components. Crystallographic studies, together with the chemical reactivity and stability profiles resulted from density functional theory and ab initio molecular orbital calculations, explain the particular behavior of ethanolamine 4-nitrobenzoate in biological test.
The experimental and theoretical data reveal the potential of these compounds to contribute to the design of new active pharmaceutical ingredients with lower toxicity.
- Nitrobenzoic and chloronitrobenzoic acids and derivatives
- Single crystal X-ray diffraction
- Chemical reactivity
Nitroaromatic and chloronitroaromatic compounds are versatile precursors, the vast majority synthetic and frequently employed as important intermediates for the synthesis of industrial chemicals, dyes, pigments and pharmaceutical drugs [1–3]. The functional groups nitro and chloride provide chemical and structural diversity and a significant impact on properties and reactivities of chemicals, making these compounds attractive in different research fields over the past decades. Many pharmaceuticals have their chemical origins in nitro- and chloronitroaromatic compounds. They are used to treat a wide variety of diseases categories: Parkinson, angina, insomnia and parasitic infection (e.g. Giardiasis, Amebiasis, Trichomoniasis) [4–6]. Recently, 4-nitrobenzoic acid (4-NO2BA), 2-chloro-4-nitrobenzoic acid (2-Cl-4-NO2BA) and 2-chloro-5-nitrobenzoic acid (2-Cl-5-NO2BA) have been used as active ingredients in the main infectious killing diseases in the world. Therefore, 4-NO2BA is used as inhibitor agent for identification of Mycobacterium tuberculosis complex [7–9] and 2-Cl-4-NO2BA/2-Cl-5-NO2BA used in a novel therapy for immunodeficiency diseases, including the human immunodeficiency virus (HIV) infection [10, 11].
Taking into account the prospects of nitro- and chloronitrobenzoic acids and the fact that about a half of all active pharmaceutical ingredients are used today as salts due to their improved drug’s physicochemical properties, we have focused to prepare new nitro/chloronitroderivatives with lower toxicity. Limited information on the experimental toxicity profile of benzoic agents have been founded in literature [12–15]. Experimental and/or theoretical toxicity studies on different organisms are required for new derivatives of an active pharmacological substance. The safety profile of a salt depends significantly on the cation nature and the alkyl side chain. The majority of hydrogen-bonded complexes of nitro-, respectively chloronitrobenzoic acids use heterocyclic amines (e.g. pyridine, piperazine, morpholine), considered strong mutagens [16–18].
This study proposes the development of new compounds with dual biological activity, based on biologically active components: ethanolammonium as cation and 2-chloro-4-nitrobenzoate, 2-chloro-5-nitrobenzoate and respectively 4-nitrobenzoate as anions. Ethanolamine (EA) is a naturally occurring component and a suitable model of alkanolamines, which is safer, economical and commercially available, used as base chemical in the production of pharmaceuticals. It is an essential component of cell membranes, present in the synthesis of membrane lipids, such as phosphatidylethanolamine and phosphatidylcholine . Recent studies introduce EA as a new therapeutic agent in the treatment of age-associated human diseases, being involved in autophagy regulation .
Besides the plethora of properties noted above, both EA and substituted benzoic acids have the ability to establish strong and directional hydrogen bonds, forming supramolecular systems with different topologies, important in crystal engineering [21–23]. Continuing our interest in multicomponent organic crystal with dual biological activity, some of the major objectives of this paper are to obtain crystalline ethanolamine salts of 2-Cl-4-NO2BA (1), 2-Cl-5-NO2BA (2) and 4-NO2BA (3), to characterize them physico-chemically and structurally, and to investigate their supramolecular synthons and molecular packing. A comparative analysis on toxicity of the compounds studied is also presented, using Hydractinia echinata, previously demonstrated to be an excellent test system [12, 24]. In order to explain the particular behavior of ethanolamine salts of nitro- and chloronitrobenzoic acids in biological test, a theoretical study regarding the chemical reactivity and stability profiles was described in correlation with their biological activity (toxicity) and X-ray structures.
All the chemicals used for the synthesis were of analytical grade and purchased from Fluka AG (Buchs SG). EA was freshly distilled before any use. Melting points were determined on Boetius melting point apparatus and are uncorrected. FT-IR spectra of new synthesized compounds were recorded as KBr pellet on a JASCO—FT/IR-4200 spectrometer, in the range 4000–400 cm−1, with a resolution of 4.0 cm−1 and a scanning speed of 16 mm s−1. The optical properties were examined by using a UV–vis spectrophotometer at room temperature. UV–vis spectra in the 190–800 nm range were measured on PERKIN-ELMER LAMBDA 12 UV–vis spectrometer.
Synthesis and characterization of compounds 1–3
C9H12ClN2O5 (262.58), m.p. 112–115 °C; λmax = 277.93 nm; FT-IR spectra (KBr pelet, cm−1): 3194 (νNH3 +), 1613 (νasCOO−), 1520 (δNH3 +), 1394 (νsCOO−), 1065 (νC–O), calcd. (%): C 41.13, H 4.57, Cl 13.52, N 10.66; found (%): C 41.04, H 4.41, Cl 13.39, N 10.51;
C9H12ClN2O5·H2O (280.58), m.p. 116–118 °C; λmax = 279.81 nm; FT-IR spectra (KBr pelet, cm−1): 3166 (νNH3 +), 1621 (νasCOO−), 1519 (δNH3 +), 1404 (νsCOO−), 1087 (νC–O), calcd. (%): C 38.49, H 4.99, Cl 12.65, N 9.98; found (%): C 38.39, H 4.82, Cl 12.48, N 9.84;
C9H12N2O5 (228.20), m.p. 173–175 °C; λmax = 273.55 nm; FT-IR spectra (KBr pelet, cm−1): 2955 (νNH3 +), 1648 (νasCOO−), 1515 (δNH3 +), 1392 (νsCOO−), 1018 (νC–O), calcd. (%): C 47.33, H 5.25, N 12.27; found (%): C 47.18, H 5.13, N 12.09;
Single crystal X-ray structure analysis
Crystal data and structure refinement parameters for compounds 1–3
Unit cell dimensions
D calc (g cm−3)
Crystal size (mm)
0.22 × 0.14 × 0.07
0.15 × 0.12 × 0.08
0.35 × 0.4 × 0.6
θ range for data collection (°)
− 10 ≤ h ≤ 1
− 8 ≤ h ≤ 8
0 ≤ h ≤ 7
− 7 ≤ k ≤ 0
− 4 ≤ k ≤ 11
0 ≤ k ≤ 11
− 26 ≤ l ≤ 26
− 11 ≤ l ≤ 11
− 25 ≤ l ≤ 25
Reflections with [I > 2σ(I)]
Goodness-of-fit on F2
R 1, wR 2 [I > 2σ(I)]
R 1, wR 2 (all data)
Largest difference in peak and hole (e Å−3)
Crystallographic data were deposited in the Cambridge Crystallographic Data Centre (CCDC Numbers 853484, 853485, 257807).
The toxicity of compounds 1–3 and the correspondent components was evaluated using H. echinata as test system. This simple and common invertebrate living in European and North American coastal areas is rapidly reproductible and considered sustainable for in vivo experiments. The used method was identical with the one described in previous articles [12, 24]. Dishes with 30 H. echinata larvae were exposed to 3 mL seawater (980 mosmol, pH 8.2, 18 °C) containing 20 mM CsCl and the test compounds were added for 3 h. The concentration (mol L−1) at which the frequency of metamorphosis induction larvae to polyp was reduced by 50% with respect to control was determined after 24 h and noted as MRC50 (metamorphosis reducing concentration) and correspond to EC50 (50-effective concentration) in literature. The measured value of toxicity (M) was showed as logarithm of reciprocal value of MRC50 (M = log 1/MRC50). Triplicate experiments were performed for each concentration assessment of title compounds and each experiment was repeated twice.
Synthesis and characterization
Tabulated ΔpKa and hydrogen-bonding geometry (Å, °) for compounds 1–3
Symmetry transformation for H-acceptor
x, y, z
x + 1, y + 1, z + 1
− x, − y + 1, − z + 1
− x, − y + 1, − z + 1
x, y, z
− x, − y, − z
− x + 1, − y, − z + 1
x − 1, y, z
− x + 2, − y, − z + 1
x, y, z
x, y, z
x, y, z
− x − 1, + y − 1/2, − z + 1/2
x + 1, y, z
The IR spectra provide the evidence of salt formation by the presence of absorption bands in the regions 1650–1540 and 1450–1360 cm−1, arising from the asymmetric and symmetric vibrations of the COO− group and by the absence of bands at 1710–1680 cm−1, corresponding to carbonyl stretch (νC=O) and 1320–1210 cm−1 characteristic for νC–OH vibrations in a COOH group (Additional file 1: Figure S2) [35, 36]. The appearance of C–O vibrations at 1100–1000 cm−1, belonging to –CH2OH from EA was a supplementary proof of salt formation. The asymmetric –NH3 + stretching vibrations are observed in the region 3200–2800 cm−1 and weak bands of symmetric stretching –NH3 + near 2600 and 2100 cm−1. A strong –NH3 + deformation band is observed at 1550–1485 cm−1, which almost overlaps with asymmetric vibration of –NO2 group at 1570–1485 cm−1.
The anions in all compounds studied form a non-planar systems, more evident in compounds 1 and 2, where the dihedral angles between the least-square planes of the phenyl rings C(3)C(4)C(5)C(6)C(7)C(8) and the least-square planes of the –COO− and –NO2 groups are equal to 82.5, 11.0° and 47.6, 10.4°, while for compound 3 these values are 5.8 and 4.2° respectively. The cation adopts the—Syn-Clinal conformation, the N(1)C(1)C(2)O(1) torsion angles in compounds 1 and 2 being equal to − 63.2° and 51.5°, respectively, towards compound 3 where is 77.4°.
In compound 1, the system of hydrogen bridges is formed only by the –NH3 + and –OH groups from ethanolamine and carboxylate oxygen atoms. In this crystal, two cations and two anions are held together by two N–H…O hydrogen bonds (N(1)…O(2) 2.747(2) Å) and two O–H…O bonds (O(1)…O(3) 2.682(2) Å) (Fig. 2a). The R 4 4(18) synthon is stabilized by two intermolecular hydrogen bonds N–H…O (N(1)…O(1) 2.868(2) Å). These synthons form infinite chains by four hydrogen bonds N–H…O (N(1)…O(2) 2.815(2) Å) (Fig. 2b), which are further consolidated into 2-D layers through C–H…O hydrogen bonds (where O is one oxygen atom from –NO2 group of benzoate) developed along y direction (Fig. 2c). It is notable that a chlorine atom and the other oxygen atom from the NO2-group of benzoate are not involved in the crystalline structure.
Compound 2 differs from compound 1 by the position of the substituent atoms, and this compound crystallizes as a hydrate, which changes the system of hydrogen bonds (Table 2). Thus, in addition to the system of hydrogen bonds formed by the –NH3 + and –OH groups from ethanolamine and carboxylate oxygen atoms, a water molecule is involved (Fig. 1). Water molecule acts as donor in the hydrogen bond with anion and as acceptor in hydrogen bond with cation. So, the nets contain R 3 3(11) synthons formed by three hydrogen bonds (N(1)…O(1W) 2.858(2) Å), O(1)…O(3) 2.716(2) Å and O(1W)…O(2) 2.754(2) Å) linked between them by the R4 2(8) synthons formed by two O–H…O bonds (O(1W)…O(2) 2.754(2) Å and O(1W)…O(2) 2.857(2) Å). The cations, anions, and water molecules form hydrogen bonded chains, which are further hydrogen bonded to one another by pairs of water molecules, to form layers (Fig. 3a, b). The layers formed in this way stack along the y axis (Fig. 3c).
In the crystal structure of compound 3, the anion and cation are held together by two charge-assisted hydrogen bonds and form the infinite helix-like chains along b direction . These chains are joined by glide plane in double chains through the C(1)–H…O(5) H-bonds (check labeling). Thus in the crystal packing of compounds 1 and 2, the anions and cations are self-assembled via N–H…O, O–H…O hydrogen bonds to form the chains. These ones are consolidated into 2-D layers by C–H…O H-bonds and water molecules respectively, while the crystal packing of three adopts the chain-like structure. The anions in compounds 1 and 2 form essentially non-planar systems in comparison with that in compound 3.
Experimental determination and theoretical investigation of toxicity behaviour
The influence of the compounds studied against the transformation from larva to polyp of marine organism H. echinata was evaluated. The measured values of toxicity were found to be 3.22 logarithm unites (log u.), respectively 3.10 log u. in the case of chloronitro-compounds (1, 2) and 1.78 log u. for nitro-compound (3). This situation is of particular interest, because compound 3 shows a much lower toxicity (1.78 log u.) even if compared to individual components (4-NO2BA: 3.11 log u. and EA: 2.67 log u.).
With this purpose we proceed to investigate the molecular structures of the three compounds using ab initio methods given that the quantum chemical descriptors have shown good correlations with biological activity since the later is dependent on the nature of the compound. Using the frontier orbitals values, the energy gap and several global reactivity descriptors, such as chemical hardness (η), chemical potential (µ), electronegativity (χ) and electrophilicity index (ω), were calculated to evaluate the chemical reactivity and the stability of these compounds in order to elucidate a possible mechanism of action for this type of compounds.
The energy gaps between HOMO and LUMO orbitals (Additional file 1: Figure S3) offer information about chemical reactivity of a molecule, a small gap suggesting an easy electronic transition and therefore a higher chemical reactivity. Chloronitro-compounds 1 and 2 have a larger energy gap than nitro-compound 3, which means that the electron transfer between the HOMO and LUMO orbitals is more possible to occur in the case of compound 3 than in the case of compounds 1 and 2.
Chemical reactivity descriptors
Chemical stability descriptors
Heat of formation (kcal mol−1)
On the contrary, compound 3, which is less toxic, has a lower stability but a higher reactivity, although toxicity is usually associated with large values of reactivity . In this light, probably the mechanism of action involves an accumulation of compounds 1 and 2 in the system since they do not react easily and are very stable.
Anion stability and reactivity
Ionization potentials and acidic constants for anions A1–A3
Ionization potential (eV)
The electron density of the anions is determined by the electronic effects of substitutents, so indirectly, if an acid is stronger, its conjugate base (anion) is weaker and more stable. Anion A3 is less stable, since the pKa of its conjugate acid is higher . All three anions contain a nitro group with a resonance effect − R, that decreases the electron density at the aromatic ring, and thus it increases its inductive effect against carboxyl group, which can free easily the proton. 4-NO2BA and 2-Cl-4-NO2BA are stronger acids compared to 2-Cl-5-NO2BA due to the presence of nitro group in para position, which implies a greater closeness of the carboxyl group to the positive charge which occurs on the aromatic ring through conjugation with nitro group. If the positive charge is closer to the carboxyl, the electrons are shifted closer to the oxygen atom and the anion is more stable. Chlorine atom containing pairs of non-bonding electrons has a resonance effect + E against the aromatic ring, functioning as a strong electron donor. Implicitly, the − R effect of nitro group is mitigated by the presence of chlorine per core, so the anion of compound 3 is a stronger acid and therefore has greater stability than those with chlorine.
Biological active salts with ethanolammonium as cation and 2-chloro-4-nitro/2-chloro-5-nitrobenzoate or 4-nitrobenzoate as anion were structurally, chemically and toxicologically investigated. Single crystal X-ray diffraction confirmed the proton transfer and ionic hydrogen bonds formation between the components. Therefore, in compounds 1 and 2, the components are linked only by one HB, while in compound 3 two charge-assisted HB are formed. Chloronitro compounds (1 and 2) adopt layered structure and nitrocompound (3) chain-like structure. The presence of chlorine atom in benzene rings has led to rotation of carboxyl groups with respect to this ring which forms the essentially non-planar anions of 1 and 2 in comparison with compound 3. The position of nitro-substituent leads to change in the system of hydrogen bonds between anion and cation in compound 2 which crystallizes as hydrate. The measured value of toxicity indicates a very low order of toxicity for compound 3 (1.78 log u.), almost half of the value of compound 1, respectively 2 (3.22 log u., 3.10 log u.) and lower compared to the individual components (4-NO2BA: 3.11 log u. and EA: 2.67 log u.). The theoretical study regarding chemical reactivity and stability profiles explains the experimental values of toxicity. Compound 3 showed a higher reactivity and a lower stability compared to its compound 1 and 2, which is in agreement with the lowest toxicity value measured in biological assay. In conclusion, toxicity test on H. echinata in relation with density functional theory, ab initio molecular orbital calculations and crystallographic study leads to a better understanding of nitro/chloronitro substituent effect on toxicity, contributing to the design of new compounds with low toxicity and practical applicability.
MC performed the synthesis and characterization of the compounds. LH realized the theoretical calculations. YC and PB did the crystallographic studies. SAC performed the toxicity test. All authors have contribution in write-up. All authors read and approved the final manuscript.
This work was developed through a bilateral project between Romania and Moldova, CCCDI-UEFISCDI, PN3-P3-217/24BM/19.09.2016 (Romania), 16.80013.5007.04/RO (Moldova).
The authors declare that they have no competing interests.
Ethics approval and consent to participate
No human subjects are involved in this research.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ono N (2002) The nitro group in organic synthesis. Wiley, New YorkGoogle Scholar
- Tripathi KD (2009) Essentials of medical pharmacology. Jaypee Brothers Medical Publishers (P) Ltd, New DelhiGoogle Scholar
- Zaragoza DF (2012) Nitro compounds. Wiley-VCH Verlag GmbH & Co, WeinheimGoogle Scholar
- Truong DD (2009) Tolcapone: review of its pharmacology and use as adjunctive therapy in patients with Parkinson’s disease. Clin Interv Aging 4:109–113View ArticleGoogle Scholar
- Raether W, Hanel H (2003) Nitroheterocyclic drugs with broad spectrum activity. Parasitol Res 90(Supp 1):19–39Google Scholar
- Wilkinson SR, Christopher B, John MK, Belinda SH (2011) Trypanocidal activity of nitroaromatic prodrugs: current treatments and future perspectives. Curr Top Med Chem 11:2072–2084View ArticleGoogle Scholar
- Singla D, Tewari R, Kumar A, Raghava G et al (2013) Designing of inhibitors against drug tolerant Mycobacterium tuberculosis (H37Rv). Chem Cent J 7:49View ArticleGoogle Scholar
- Shakoor S, Ahsan T, Jabeen K, Raza M, Hasan R (2010) Use of p-nitrobenzoic acid in 7H10 agar for identification of Mycobacterium tuberculosis complex: a field study. Int J Tuberc Lung Dis 14:1644–1646Google Scholar
- Giampaglia CMS, Martins MC, Chimara E, Oliveira RS, de Oliveira Vieira GB et al (2007) Differentiation of Mycobacterium tuberculosis from other mycobacteria with p-nitrobenzoic acid using MGIT960. Int J Tuberc Lung Dis 11:803–807Google Scholar
- Kinchinqton D, Ng T, Mathews N, Tisdale M, Devine D, Ayuko WO (1997) Tcell costimulation by derivatives of benzoic acid. Antivir Chem Chemother 8:121–130View ArticleGoogle Scholar
- Lemmerer A (2012) Covalent assistance to supramolecular synthesis: modifying the drug functionality of the antituberculosis API isoniazid in situ during co-crystallization with GRAS and API compounds. CrystEngComm 14:2465–2478View ArticleGoogle Scholar
- Chicu SA, Grozav M, Kurunczi L, Crisan M (2008) SAR for amine salts of carboxylic acids to Hydractinia echinata. Rev Chim 59:582–587Google Scholar
- Isayev O, Rasulev B, Gorb L, Leszczynski J (2006) Structure-toxicity relationships of nitroaromatic compounds. Mol Divers 10:233–245View ArticleGoogle Scholar
- Keshavarz MH, Pouretedal HR (2013) Simple and reliable prediction of toxicological activities of benzoic acid derivatives without using any experimental data or computer codes. Med Chem Res 22:1238–1257View ArticleGoogle Scholar
- Sun Y, Li Z, Yan X, Wang L, Meng F (2009) Study on the quantitative structure–toxicity relationships of benzoic acid derivatives in rats via oral LD50. Med Chem Res 18:712–724View ArticleGoogle Scholar
- Hemamalini M, Loh WS, Quah CK, Fun HK (2014) Investigation of supramolecular synthons and structural characterisation of aminopyridine-carboxylic acid derivatives. Chem Cent J 8:31View ArticleGoogle Scholar
- Chen Z, Peng M (2011) Supramolecular architectures constructed from piperazine and substituted benzoic acids. J Chem Crystallogr 41:137–142View ArticleGoogle Scholar
- Ishida H, Rahman B, Kashino S (2001) Morpholinium 2-chloro-4-nitrobenzoate, 2-chloro-5-nitrobenzoate and 4-chloro-3-nitrobenzoate. Acta Crystallogr Sect C Cryst Struct Commun 57:1450–1453View ArticleGoogle Scholar
- Gibellini F, Smith TK (2010) The Kennedy pathway—de novo synthesis of phosphatidylethanolamine and phosphatidylcholine. IUBMB Life 62:414–428View ArticleGoogle Scholar
- Rockenfeller P, Koska M, Pietrocola F, Minois N, Knittelfelder O et al (2015) Phosphatidylethanolamine positively regulates autophagy and longevity. Cell Death Differ 22:499–508View ArticleGoogle Scholar
- Chumakov Y, Simonov Y, Grozav M, Crisan M, Bocelli G et al (2006) Hydrogen-bonding network in the organic salts of 4-nitrobenzoic acid. Cent Eur J Chem v4:458–475Google Scholar
- Crisan M, Bourosh P, Chumakov Y, Petric M, Ilia G (2013) Supramolecular assembly and Ab initio quantum chemical calculations of 2-hydroxyethylammonium salts of para-substituted benzoic acids. Cryst Growth Des 13:143–154View ArticleGoogle Scholar
- Crisan ME, Bourosh P, Maffei ME, Forni A, Pieraccini S, Sironi M, Chumakov YM (2014) Synthesis, crystal structure and biological activity of 2-hydroxyethylammonium salt of p-aminobenzoic acid. PLoS ONE 9:e101892View ArticleGoogle Scholar
- Chicu SA, Herrmann K, Berking S (2000) An approach to calculate the toxicity of simple organic molecules on the basis of QSAR analysis in Hydractinia echinata (Hydrozoa, Cnidaria). Quant Struct Act Relat 19:227–236View ArticleGoogle Scholar
- Dorosencu M, Rad R, Grozav M, Neamtiu I (2001) The quantitative determination of some quaternary ammonium salts of benzoic acids. Ann W U T-Ser Chem 10:337–342Google Scholar
- Sheldrick GM (1997) SHELX 97: programs for crystal structure analysis. University of Göttingen, Germany Google Scholar
- Burnett MN, Johnson CK. ORTEP-III: Oak Ridge thermal ellipsoid plot program for crystal structure illustrations, Oak Ridge National Laboratory Report ORNL-6895, 1996Google Scholar
- Bochevarov AD, Harder E, Hughes TF, Greenwood JR, Braden DA et al (2013) Jaguar: a high-performance quantum chemistry software program with strengths in life and materials sciences. Int J Quantum Chem 113:2110–2142View ArticleGoogle Scholar
- Schrödinger L (2015) New York, NY, Schrödinger Release 2015-3: Jaguar, 8.9 edGoogle Scholar
- Schrödinger L (2015) New York, NY, Schrödinger Release 2015-3: MacroModel, version 10.9 edGoogle Scholar
- Schrödinger L (2015) New York, NY, Schrödinger Release 2015-3: Maestro, version 10.3 edGoogle Scholar
- Koopmans T (1934) Über die Zuordnung von Wellenfunktionen und Eigenwerten zu den Einzelnen Elektronen Eines Atoms. Physica 1:104–113View ArticleGoogle Scholar
- Cruz-Cabeza AJ (2012) Acid-base crystalline complexes and the pKa rule. CrystEngComm 14:6362–6365View ArticleGoogle Scholar
- Delori A, Galek PTA, Pidcock E, Patni M, Jones W (2013) Knowledge-based hydrogen bond prediction and the synthesis of salts and cocrystals of the anti-malarial drug pyrimethamine with various drug and GRAS molecules. CrystEngComm 15:2916–2928View ArticleGoogle Scholar
- Pavia DL, Lampman GM, Kriz GS (2009) Introduction to spectroscopy: a guide for students of organic chemistry. W. B. Saunders Co., PhiladelphiaGoogle Scholar
- Silverstein RM, Webster FX (1998) Spectrometric identification of organic compounds. Wiley, New YorkGoogle Scholar
- Cronin MTD, Bajot F, Enoch SJ, Madden JC, Roberts DW, Schwöbel J (2009) The in chemico-in silico interface: challenges for integrating experimental and computational chemistry to identify toxicity. ATLA 37:513–521Google Scholar
- Jover J, Bosque R, Sales J (2008) QSPR prediction of pKa for benzoic acids in different solvents. QSAR Comb Sci 27:563–581View ArticleGoogle Scholar