Benzimidazole is an important heterocyclic fragment, present in many biologically active compounds with a great variety of therapeutic purposes. Most of the benzimidazole activities are explained through the existence of 1,3-tautomeric equilibrium. As the binding affinity of each tautomer to a protein target depends on an established bioactive conformation, the effect of tautomers on the ligand protein binding mechanism is determinant.
1. Introduction
The tautomerism phenomenon is considered as a dynamic equilibrium between interconvertible structural isomers, named tautomers, with the migration of one atom or group. When a hydrogen atom migrates, the phenomenon is known as prototropy; however, other groups such as alkyl, aryl, acyl, cyano, halogens, amines, and nitro, as well as methals can migrate; Scheme 1
[1].
Scheme 1. Tautomeric equilibria in organic compounds.
Nuclear magnetic resonance (NMR) spectroscopy is one of the most useful techniques to study tautomerism. The signals observed in the spectra depend on the activation energy between the tautomers, which determines their lifetimes and chemical shift difference (Δδ). When lifetimes are long, compared with 1/Δδ, a slow exchange regime gives rise to separate narrow signals for each of the tautomers. In this case, integration of the 1H-NMR signal intensities is the method of choice to study the tautomerism. For shorter lifetimes, the exchange is raised, leading to line broadening. Then, high magnetic fields and lower temperatures can be used to achieve slow exchange conditions. Interpolation is the method of choice when lifetimes are shorter than 1/Δδ, and fast exchange and signal coalescence predominate to give averaged narrow signals, but the NMR chemical shifts of the individual tautomers remain unknown. Four techniques have been used to solve this problem: (1) the use of blocked derivatives of individual tautomer, replacing the tautomeric proton by a methyl group, and performing a correction for the substituent effect; (2) the use of model compounds that exclusively exist in one tautomeric form; (3) the use of the properties measured in the solid state where only one tautomer exists, but phase effects should be considered; and (4) the use of theoretically calculated properties as GIAO absolute shielding, however, solvent effects are difficult to estimate.
Benzazoles (BZs) are constituted by a benzene ring fused to an oxazole (BO), thiazole (BT), or imidazole (BI) ring; left in Figure 1. A typical kind of tautomerism in BI is the relocation of a proton, known as annular tautomerism (a). BZs containing an exocyclic heteroatom, where a proton can migrate from cyclic nitrogen to an exocyclic heteroatom, retaining the aromaticity, represent an exocylic tautomerism (b); right in Figure 1.
Figure 1. (a) Intracyclic and (b) exocyclic tautomerism in benzazoles (BZs).
In commercial databases, the presence of tautomeric duplicates has been found, and both tautomers are offered as different products
[2]. On the other hand, different names for tautomers based on the IUPAC rules have been used
[3]. In general, compounds with well-known tautomerism are named as the predominant tautomer at equilibrium. For example, 2-mercaptobenzimidazole (MBI) is named as benzimidazolidinethione (BIT) or 1,3-dihydro-2
H-benzo[d]imidazole-2-thione (X = NH, Y = S);
Figure 1.
The tautomeric interconversion with a low-energy barrier, in general, displays averaged signals in the NMR spectrum. To slow down the interconversion rate in solution, the following strategies are used: lowering temperature and/or the use of hexamethylphosphoramide-d
18, which prevents the formation of hydrogen bonds. Thermodynamic and kinetic aspects of the tautomeric equilibrium have been explored using dynamic NMR studies
[4][5][6] Because of the wide range of chemical shifts and high sensitivity,
15N-NMR spectroscopy has been used to study the tautomerism of nitrogen in heteroaromatic compounds
[7]. Solid-state NMR experiments
[4] 4, X-ray crystallography
[8][9][10], and
15N-NMR studies in DMSOd
6 for both BIT and 1-MeBIT revealed that the thione form is the predominant tautomer
[11]. Recently, Pandey et al. theoretically studied 5-MeOBIT using B3LYP methods with a 6-311++G (d, p) basis set. A comparison between the experimental and calculated structure as well as the calculated and experimental
1H and
13C chemical shifts showed a good correlation with the thione isomer
[12].
Among the BZs, benzimidazole (BI) is the most known and the one whose tautomeric equilibria have been widely studied
[4][5][6][7]. It has been stablished that, when a large substituent is bonded at the N1 position in 2-heterosubstituted BIs, two isomeric compounds exist because of the annular tautomerism. BI is an important heterocycle, present in several biologically active compounds. It has been found in progesterone receptor antagonists
[13], luteinizing hormone-releasing hormone antagonists (leuprolide, goserelin, triptorelin)
[14][15], antiviral (enviradine)
[16][17], antiprotozoal
[18][19], antimicrobial
[20][21][22], analgesic and anti-inflammatory
[23], anticonvulsant, antidiabetics
[24], anthelmintics (albendazole, mebendazole, and thiabendazole), proton pump inhibitors (omeprazole, lansoprazole, and pantoprazole), antihistaminic (astemizole), and antihypertensives (candesartan, cilexitil, and telmisartan), among others.
The binding affinity of each tautomer to a protein target depends on an established bioactive conformation. Therefore, the method proposed herein to calculate in solution the prototropic ratio in BZs, on the basis of the electronic effect of pyrrole like atom (NH = Npr) and pyridine like atom (Npd) on C4 and C7 chemical shift resonances, is very valuable. On the other hand, the method is also helpful to correctly assign the carbon atoms on the benzene ring.
2. Estimation of the Anular Tautomerism on Benzimidazol
The chemical shift reference values for the estimation of the tautomeric equilibrium were set using the 1-methyl-benzimidazole (1-MeBI, 4) as a model compound. In case 4, the tautomeric equilibrium is non-existent; therefore, a set of seven narrow signals are observed in CDCl3. The δC4 appears at 120.4 ppm, characteristic of an Npd effect, whereas the δC7 at 109.5 ppm is characteristic of an Npr effect. In DMSOd6, the δC4 and δC7 values are 119.2 and 110.1, respectively. Then, a good approximation for 100% of the pyridine like character of δC4 and 100% of the pyrrole like character of δC7 is 120.0 and 110.0 ppm, respectively. Furthermore, the tautomeric proportion in BIs can be calculated on a base equal to 120.0 ppm for δC4ref (Npd) and 110.0 ppm for δC7ref (Npr) as reference values (Figure 2).
The tautomerism in BIs can be calculated considering that the observed chemical shift δ
obsC4/C7 is the result of the averaged C4 and C7 chemical shifts weighted by the respective molar fraction contributions x
pd (pyridine character) and x
pr (pyrrolic character). Equations (1)–(3) can be used to calculate the molar fractions of each tautomer;
Figure 2.
Figure 3. 13C chemical shifts as references to calculate the tautomeric proportion of pyrrole–pyridine nitrogen atoms.
From the reported
13C data for BI
3, as well as their protonated and deprotonated derivatives, in DMSOd
6 as solvent,
Table 1 [25], the tautomeric proportion can be calculated. For BI, the δ
obs value is 116.3, then the estimation is 63% of the pyridine character (+N3H). In the case of δ
obs for protonated BI (118.3 ppm), a small effect of positive charge on nitrogen atoms shifts this signal to high frequencies in approximately 2.0 ppm. On the contrary, any effect of the negative charge on nitrogen in deprotonated BI was observed.
Table 1. Protonation and deprotonation of BMZ 13C in DMSOd6.
Comp. |
C2 |
C3a |
C4 |
C5 |
C6 |
C7 |
C7a |
3a/3b |
142.5 |
139.0 |
116.3 |
122.3 |
122.3 |
116.3 |
139.0 |
[3c/3d]+ |
143.6 |
133.8 |
118.3 |
130.0 |
130.0 |
118.3 |
133.8 |
[3d/3f]− Na+ |
153.1 |
146.9 |
116.4 |
116.7 |
116.7 |
116.4 |
146.9 |
[3d/3f]− Li+ |
153.5 |
147.0 |
116.8 |
117.0 |
117.0 |
116.8 |
147.0 |
This approximation can be applied to 2-aminomethylbenzimidazole (2-AMBI);
Scheme 2. In the work of Sierra-Centeno et al.
[26], the acid–base equilibrium constants of 2AMBI were determined in aqueous solutions at 25 °C through
13C NMR spectroscopy, potentiometric and spectrophotometric techniques, as well as through theoretical methods.
13C NMR spectra were recorded for solutions of 0.1 M of 2-AMBI at pH values from 1.1 to 13.2, containing 10 %
v/v of D
2O. During the titration process, only five signals were observed. The sets of C5 and C6 signals were shifted approximately 7.0 ppm to lower frequencies and C2 and C3a/C7a in 10.0 and 15.5 ppm, respectively, to higher frequencies, whereas the set of C4/C7 signals remained almost constant at 114.5 ppm. Therefore, the BI heterocycle retains its 45% pyridine character independently from the pH.
Scheme 2. Titration process of 2-(aminomethyl)benzimidazole dihydrochloride.