Hadieh Monajemi*a,
Mun Hon Cheahb,
Vannajan Sanghiran Leeb,
Sharifuddin Mohd. Zainb and
Wan Ahmad Tajuddin Wan Abdullaha
aDepartment of Physics, Universiti Malaya, 50603 Kuala Lumpur, Malaysia. E-mail: h.monajemi@hotmail.com
bDepartment of Chemistry, Universiti Malaya, 50603 Kuala Lumpur, Malaysia
First published on 3rd January 2014
Boronic acids are well-known for their ability to complex with saccharides and to form cyclic boronic esters or cyclic boronate ions. Increasing the reactivity of boronic acids would lead towards the design of a highly responsive sensor for non-enzymatic glucose monitoring applications. Many studies have been carried out to investigate the reactivity of boronic acids towards diols in different environmental effects. The symmetry around boron in their reactive species, however, is still an open question. In this study, we used computational quantum chemistry calculations to propose a model in which boronic acid is highly reactive towards diols. B3LYP/6-31+G(d,p) model chemistry was used in water to calculate the transition states of several proposed reaction mechanisms and the rates of the reactions were calculated using the transition state theory. The reaction takes place in two steps with the first step known as the rate determining step. The model we proposed in this study includes two different electronegative R-groups attached to boronic acid (i.e. CF3 and CH3) which highly influence the interaction of boron with diols. The kinetic results in our study show that boronate ion with a low electronegative R-group is a better sensor in alkaline medium and boronic acid with a high electronegative R-group is a better sensor in acidic medium. Using the latter in alkaline medium results in a very poor sensor since it is highly reactive and promptly forms a tetrahedral boronate ion. The high electronegative R-group on boronate ion decreases its basic characteristics and results in a low reactive sensor. The same scenario goes for the former i.e. boronate ion attached to a low electronegative R-group in acidic medium. As a result, we can conclude that boronic acid's reactivity is not mainly due to the symmetry around boron, it is due to the characteristics of boron itself which can be affected by both the R-group and the medium.
![]() | (1) |
![]() | (2) |
The method Lorand and Edwards used in their study involved the measurement of the equilibrium constants for the formation of complexes between benzeneboronate ion and hydroxyl containing compounds. Comparing these constants with that of borate, they concluded that the conjugate basic form of boronic acid and the boronate ester is tetrahedral. This resulted in a disagreement on what is the reactive form of boronic acid. It was suggested by Lorand and Edwards that boronate ion is reactive in alkaline medium since its complex formation lowers the pH. Several other studies have also suggested that the tetrahedral boronate ion is predominant since most of the reactions occur in alkaline solutions where sensors are most operative.7–12 On the other hand, there are a few arguments claiming that trigonal boronic acid is the reactive species regardless of the pH of the medium.13–18 The latter group has argued that the high reactivity of boronate ion in alkaline solution in which it exists abundantly is only due to the low concentration of its conjugate acid and does not necessarily mean that the basic boronate ion itself is more reactive than its conjugate acid. Their study shows that “trigonal boronic acid is always a reactive species irrespective of the pH of the solution and the reactivity of its conjugate tetrahedral boronate ion is comparable with or less than that of the boronic acid”.16,18,19 They have also recently validated their conclusion by setting up reaction systems using sensors and ligands with different pKa values, without proton ambiguity and a fully protonated ligand.19 They argue that even in alkaline solution, boronic acid is more reactive than its conjugate basic counterpart.
Determining the reactive species of boronic acid is an important step towards understanding its reaction mechanism in ester formation with diols. Another factor affecting its reactivity is the type of R-group attached to boron as this affects its sensitivity towards saccharides. The majority of studies which have been carried out have used either phenylboronic acid9,11 or boric acid.7,10,13,14 These studies have mainly focused on various combinations of diols and pH values to investigate the selectivity. We on the other hand, have focused on the R-group variety to investigate the reactivity of boronic acid. One way to look at this problem is to study the rate at which boronic acid or boronate ion reacts with diols. There are several limitations in experimental studies to figure out the accurate reaction mechanism and kinetic measurements because of the so called proton ambiguity; the transition structure cannot be clearly identified and the reaction of boronic acid in complexing diols can hardly be measured kinetically.8,14,16–18
In this study, a number of alternative reaction pathways between boronic acid and 1,2-ethanediol have been explored using Density Functional Theory (DFT) calculations. To explore the influence of the R group towards the reactivity of boronic acid, methylboronic acid (CH3B(OH)2) and trifluoromethylboronic acid (CF3B(OH)2) were selected because of the large difference in electronegativity.
![]() | (3) |
While higher levels of theory will produce more accurate energies, the purpose of this study was to discriminate between alternative reaction pathways by comparing the relative energy difference. The relative energy differences between CCSD (Coupled Cluster Single and Double excitation) and B3LYP (Becke three parameter hybrid functional combined with Lee–Yang–Parr correlation functional) methods calculated for the first proposed reaction mechanism in this study are presented in Table 1. A comparison was made of the relative activation energies determined by CCSD and B3LYP methods for the reaction of CH3B(OH)2 and CF3B(OH)2 with diol. It was observed that the DFT method follows the same trend as the coupled cluster method and there is little variation between the two methods. Furthermore, DFT methods have been shown to reproduce fairly high accuracy.24,25 Due to the relatively high computational cost associated with CCSD, DFT was adopted for the purpose of this study.26–29
Structure | Method/basis set | ΔG‡ (kcal mol−1) | k s−1 |
---|---|---|---|
CH3B(OH)2 | CCSD/aug-cc-pVDZ (PCM) | 37.652 | 2.338 × 10−15 |
B3LYP/6-31+G(d,p) (PCM) | 38.110 | 1.078 × 10−15 | |
CF3B(OH)2 | CCSD/aug-cc-pVDZ (PCM) | 28.866 | 6.555 × 10−9 |
B3LYP/6-31+G(d,p) (PCM) | 29.091 | 4.481 × 10−9 |
The basis function we used describes the set of d-functions, added to the heavy atoms as well as a set of p-functions, added to the hydrogen atoms (i.e. d,p). This function is referred to as a polarization function. We also included the diffuse function, representing the portion of atomic orbitals which are distant from the nuclei. The diffuse functions that are denoted by a “+” sign, are particularly important for the anionic structures. In addition for heavy atoms, we also added the diffuse function to light atoms such as hydrogen in the coupled cluster method, indicated by an “aug” sign for the CCSD method.
Five possible reaction paths were considered using protonated and deprotonated forms of both boronic acid and 1,2-ethanediol (mechanisms I–V) leading to two distinct products (products 1 and 2).
![]() | ||
Fig. 1 Five different reaction mechanisms for methylboronic acid in interaction with 1,2-ethanediol in water. All reactions take place in two steps, forming an intermediate in the middle. |
![]() | ||
Fig. 2 Five different reaction mechanisms for trifluoromethylboronic acid in interaction with 1,2-ethanediol in water. All reactions take place in two steps, forming an intermediate in the middle. |
Mechanism I proposes a reaction between the neutral diol and neutral boronic acid. The empty p orbital on the boron atom deshields its nucleus which increases its Lewis acidity. This would result in an interaction between the boron atom in boronic acid and the electron lone pair on the oxygen of the incoming diol.
Although a few studies have shown that the complexation reaction of tetrahedral boronate ion is 3–4 orders of magnitude faster than that of trigonal boronic acid,8–10 it is not wise to conclude that the tetrahedral boronate ion is the reactive species since it would be rather difficult to specify the rate determining step in a two step mechanism. The best way to rationalize this idea is to propose a reaction mechanism in which the trigonal boronic acid in reaction with a diol forms a tetrahedral intermediate. Based on what has been reported above and since the chelation reaction is not the rate determining step, it is more likely for the trigonal boronic acid to be the starting material. This could occur with either a Brønsted basic sensor or a basic incoming ligand. In the former, the diol donates its proton to react with trigonal boronic acid and forms a tetrahedral intermediate (mechanisms IV and V). The latter occurs when the incoming ligand lacks an acidic proton: the negative oxygen tends to react with the Lewis acidic boron (mechanisms II and III). This reaction mechanism could explain the idea pointed out by Pizer and Tihal, in which for a ligand lacking an acidic proton, trigonal boronic acid would be unreactive.8,9 This is only true if the intermediate is trigonal with elimination of a water molecule by removal of hydroxide from the central boron and a proton from the ligand (mechanism II). Otherwise, one proton is enough for the whole reaction to occur with the negative oxygen of the ligand directly attacking the boron atom, forming a tetrahedral intermediate which could easily occur if the starting material is in the trigonal form (mechanism III).
Furthermore, for ligands with no acidic proton, an unchelated tetrahedral intermediate would still form, increasing the coordination number of boron from 3 (sp2) to 4 (sp3). The result of this reaction is a tetrahedral intermediate in which the negatively charged boronate ion has a low first ionization energy. In addition, with electronegative oxygen attached to boron, the B–O bond becomes relatively polar and less stable than that for trigonal boronic acid. Thus, an acidic environment could facilitate the chelate ring closure due to a good leaving hydroxyl group. In this proposed mechanism, a lack of one acidic proton does not have a negative income on the reactivity of boronic acid and makes it rather efficient.
Mechanism/CH3B(OH)2 | Step of reaction | ΔG (kcal mol−1) | ΔS (kcal mol−1 K−1) | ΔG‡ (kcal mol−1) | ΔS‡ (kcal mol−1 K−1) | k s−1 |
---|---|---|---|---|---|---|
I/neutral | Step 1 | 6.621 | −0.031 | 38.110 | −0.047 | 1.078 × 10−15 |
Step 2 | −4.261 | 0.014 | 28.123 | −0.010 | 2.300 × 10−8 | |
Overall | 2.401 | −0.017 | — | — | 1.078 × 10−15 | |
II/acidic | Step 1 | 0.549 | −0.036 | 47.785 | −0.050 | 8.566 × 10−23 |
Step 2 | −8.210 | −0.007 | 17.625 | −0.006 | 1.163 | |
Overall | −7.661 | −0.043 | — | — | 8.566 × 10−23 | |
III/acidic & basic | Step 1 | −4.273 | 0.042 | — | — | — |
Step 2 | −3.396 | −0.001 | 8.035 | 0.003 | 1.269 × 107 | |
Overall | −7.661 | −0.043 | — | — | 1.269 × 107 | |
IV/basic | Step 1 | −12.208 | −0.044 | 20.254 | −0.042 | 1.369 × 10−2 |
Step 2 | −3.396 | −0.001 | 8.035 | 0.003 | 1.269 × 107 | |
Overall | −15.590 | −0.045 | — | — | 1.369 × 10−2 | |
V/basic | Step 1 | −3.081 | −0.032 | 34.819 | −0.047 | 2.805 × 10−13 |
Step 2 | 0.001 | −0.012 | 0.755 | −0.003 | 0.279 × 1013 | |
Overall | −3.071 | −0.044 | — | — | 2.805 × 10−13 |
A few other experimental studies have also reported a decrease in the reactivity of boronic acid as the ligand deprotonates.7,14 Ishihara et al. (1991) reported a lower rate by a factor of 104-fold for 4-isopropyltropolone (Hipt) compared to chromotropic acid (H2cht) as the ligand and Ito et al. (2003) reported a decrease in the reactivity of boronic acid as H2ipt+ deprotonates to Hipt and ipt−. However, both their proposed mechanisms, H2cht in Ishihara et al.'s study and H2ipt+ in Ito et al.'s study result in tetrahedral boronate ion and elimination of a hydronium ion which is different from mechanism I in our study. It is worth mentioning that their studies were carried out in alkaline solution where boronic acid predominantly exists in its basic form and is more reactive with fully protonated diols.
In mechanism III on the other hand, having a deprotonated ligand is an advantage since the first step goes through no activation barrier and the reaction takes place at a rather high rate, forming a tetrahedral intermediate. However, the rate of formation of cyclic boronate ester in mechanism III with a tetrahedral boronate ion as the starting material is rather high compared to the other steps. This could support the idea of a faster reaction rate for tetrahedral boronate ion compared to that for trigonal boronic acid.
Mechanisms IV and V represent the reaction in alkaline medium where boronic acid is present in a trigonal basic form. The second step of mechanism V is similar to the first step of mechanism IV where the proton from the ligand is transferred to the basic oxygen of boronate ion. The rates are however significantly different knowing that the former yields a tetrahedral cyclic boronate ester and the latter yields a Lewis basic tetrahedral intermediate. This indicates the low stability of the trigonal basic form of boronic acid and its tendency to form a tetrahedral Lewis base which further supports the idea of tetrahedral boronate ion as the conjugate basic form of boronic acid. The overall rates of reaction in these mechanisms are relatively high which indicates the importance of a fully protonated ligand for boronate ions.
Mechanism/CF3B(OH)2 | Step of reaction | ΔG (kcal mol−1) | ΔS (kcal mol−1 K−1) | ΔG‡ (kcal mol−1) | ΔS‡ (kcal mol−1 K−1) | k s−1 |
---|---|---|---|---|---|---|
I/neutral | Step 1 | 4.890 | −0.030 | 29.091 | −0.047 | 4.481 × 10−9 |
Step 2 | 0.411 | 0.016 | 26.862 | −0.015 | 1.937 × 10−7 | |
Overall | 5.301 | −0.014 | — | — | 4.481 × 10−9 | |
II/acidic | Step 1 | −4.780 | −0.035 | 32.161 | −0.039 | 2.503 × 10−11 |
Step 2 | −24.160 | 0.001 | 9.601 | −0.003 | 0.900 × 106 | |
Overall | −28.940 | −0.040 | — | — | 2.503 × 10−11 | |
III/acidic & basic | Step 1 | −25.601 | −0.037 | — | — | — |
Step 2 | −3.330 | −0.003 | 26.640 | −0.009 | 2.819 × 10−7 | |
Overall | −28.940 | −0.040 | — | — | 2.819 × 10−7 | |
IV/basic | Step 1 | −5.630 | −0.035 | 24.141 | −0.037 | 1.923 × 10−5 |
Step 2 | −3.330 | −0.003 | 26.640 | −0.009 | 2.819 × 10−7 | |
Overall | −8.960 | −0.038 | — | — | 2.819 × 10−7 | |
V/basic | Step 1 | 4.230 | −0.024 | 36.120 | −0.042 | 3.113 × 10−14 |
Step 2 | −15.010 | −0.007 | 19.460 | −0.011 | 5.238 × 10−2 | |
Overall | −10.780 | −0.032 | — | — | 3.113 × 10−14 |
As for mechanism I, the rate also increased in mechanism II where the diol lacks an acidic proton. The electronegative R-group in this structure slightly increases the acidity of boron and the electrophilicity which increases its tendency to react with the nucleophilic oxygen of diol. This effect is observed more obviously in the second step of mechanism II with a relatively high rate (106 s−1) where the chelate ring closure does not involve deprotonation of the diol since it already lacks a proton. Furthermore, the higher electronegative R-group of CF3–B(OH)2 results in a 105-fold higher rate in this step compared to that for CH3–B(OH)2.
Mechanism III is less likely to occur compared to that for CH3–B(OH)2 due to a relatively high activation barrier in the second step. In this mechanism, the tetrahedral intermediate is in the Lewis base form which makes boron a high electron donor. The high electronegative R-group lowers this characteristic and results in a less electron donating boron which decreases the rate of reaction by a factor of 1014-fold compared to that for CH3–B(OH)2. The same effect was also observed for the Brønsted basic form of boronic acid which lacks a proton on one of the OH groups. The high electronegativity reduces the basic characteristic of boron, resulting in a 103-fold and 10-fold decrease in step 1 of mechanisms IV and V respectively. The second step of mechanism IV is similar to that for mechanism III. However, it is the rate determining step for mechanism IV with CF3 R-group. The higher activation barrier reduces the rate by a factor of 10−2-fold relative to step 1 of mechanism IV.
The effect of the electronegative R-group was further observed in mechanism V where both steps have a relatively higher activation barrier compared to that for the CH3 R-group. Although the TS2V seems rather stable due to a simple proton transfer between the ligand and the basic oxygen of boron and indeed it is for CH3 R-group, it requires a fairly high activation energy with the CF3 R-group to proceed with the chelate ring closure. It is the electronegative R-group which reduces the stability of the transition structure and affects the reaction rate. These results are consistent with those carried out empirically by Watanabe et al.19 Among the R-groups they used were 3-fluorophenylboronic acid (3-FPhB(OH)2) and 3-pyridylboronic acid (3-pyB(OH)2) with the former being more electronegative. The reaction rate for tetrahedral 3-fluorophenylboronate ion is 104-fold lower than that for 3-pyridylboronate ion. Furthermore, the rate of reaction is 103-fold higher for 3-FPhB(OH)2 compared to its tetrahedral Brønsted basic form (i.e. 3-FPhB(OH)3−).
Associated with this is the fact that the more electronegative CF3-group also reduces the basicity of the Brønsted boronate ion (Tables 2 and 3). This results in less nucleophilic basic oxygen on boronic acid which causes a delay in the proton transfer of incoming diol to the boronic acid leaving group. The difference was clearly observed in TS1IV, TS2III&IV, TS1V and TS2V. There was a 103-fold reduction in the reaction rate upon using a CF3 R-group rather than a CH3 R-group in the first step of mechanism IV.
In mechanism V, the neutral diol reacts with boron in a similar fashion as in the first step of mechanism I which also results in a lower rate by a factor of 10-fold for the CF3 R-group. The second step involves a proton transfer from the other end of the diol followed by nucleophilic attack of the resulting negatively charged oxygen towards boron which leads to the formation of a tetrahedral cyclic boronate ester (product 2). Similar to the first step of mechanism IV, the second step of mechanism V involves a proton transfer from the hydroxyl of diol to the Brønsted basic oxygen of boronate ion. The resulting negative oxygen of diol becomes highly reactive with a high tendency to form a tetrahedral Lewis boronate ion. This is clearly understood from the significantly lower activation barrier in mechanism V with a lower electronegative R-group (Fig. 1).
This study made it rather clear that regardless of the properties of the incoming ligand, boronate ions always have Lewis characteristics in forming the cyclic boronate ester. Lewis bases eliminate a water molecule and form a tetrahedral cyclic boronate ester. Brønsted bases initially form a Lewis acid by accepting a proton from the ligand and then the negatively charged oxygen on the ligand forms a tetrahedral intermediate with electrophilic boron on the boronic acid which ultimately leads to the formation of a cyclic boronate ester. For bifunctional ligands lacking an acidic proton, trigonal boronic acid is more reactive due to the empty p orbital on boron; the nucleophilic oxygen of the ligand attacks the electrophilic boron and forms a tetrahedral intermediate without any leaving OH. Chelate ring closure then takes place by proton transfer from the ligand to the leaving OH.
We have evidently observed in this study that boronic acids could be highly reactive in both acidic and basic forms under the effect of different R-groups. An electronegative R-group enhances the reactivity of boronic acids while retarding the reactivity of boronate ions. Therefore, boronic acids attached to highly electronegative R-groups are more reactive in acidic environment, while boronate ions are more reactive in alkaline medium with highly electropositive R-groups. However, using boronic acid in alkaline medium with electronegative R-group results in a weak sensor since it promptly forms a tetrahedral boronate ion upon reacting with the bases in the medium. Now that the structure is basic, being attached to an electronegative R-group further decreases its basicity and as a result, its reactivity. On the other hand, a tetrahedral boronate ion attached to a lower electronegative R-group results in a good sensor due to the fact that boron's basic characteristic and as a result, its reactivity is enhanced. The same scenario is also true for trigonal boronic acids attached to highly electronegative R-groups. Being highly reactive in alkaline medium, they make a better sensor in acidic environment.
This journal is © The Royal Society of Chemistry 2014 |