Daniele
Del Giudice‡
ab,
Emanuele
Spatola‡
ab,
Matteo
Valentini
a,
Cecilia
Bombelli
b,
Gianfranco
Ercolani
*c and
Stefano
Di Stefano
*ab
aDipartimento di Chimica, Università degli Studi di Roma “La Sapienza”, P.le A. Moro 5, I-00185 Rome, Italy. E-mail: stefano.distefano@uniroma1.it
bISB-CNR Sede Secondaria di Roma – Meccanismi di Reazione c/o Dipartimento di Chimica, Università degli Studi di Roma “La Sapienza”, P.le A. Moro 5, I-00185 Rome, Italy
cDipartimento di Scienze e Tecnologie Chimiche, Università di Roma Tor Vergata, Via della Ricerca Scientifica, 00133 Rome, Italy. E-mail: ercolani@uniroma2.it
First published on 21st April 2021
In this report it is shown that nitroacetic acid 1 (O2NCH2CO2H) can be conveniently used to control the pH of a water solution over time. Time-programmable sequences of the kind pH1(high)–pH2(low)–pH3(high) can be achieved, where both the extent of the initial pH jump (pH1(high)–pH2(low)) and the time required for the subsequent pH rising (pH2(low)–pH3(high)) can be predictably controlled by a judicious choice of the absolute and relative concentrations of the reagents (acid 1 and NaOH). Successive pH1(high)–pH2(low)–pH3(high) sequences can be obtained by subsequent additions of acid 1. As a proof of concept, the method is applied to control over time the pH-dependent host–guest interaction between alpha-cyclodextrin and p-aminobenzoic acid.
Since biochemistry has evolved in water, one of the ancestral tools used by nature to control the chemical properties of a compound, or a group of compounds, is pH. For this very reason, the operation of many artificial systems and materials inspired by nature, and aspiring to have life-like properties, is based on pH changes affecting protonation or deprotonation of the involved molecular structures.3 DNA-based molecular machines whose motions are guided by pH changes are appropriate examples.4 In this case, a time-predictable change of pH would guarantee a full control on the molecular motions of the machine over time. The same holds for chemical architectures or aggregates, which assemble and disassemble under the influence of pH variations.5 In one case, the time controlled assembly/disassembly of a micellar system was elegantly achieved through programmed variations of pH (up–down sequences) obtained by a judicious choice of the hydrolysis conditions of a series of ester compounds of abiotic nature.6 In other cases, the time controlled assembly/disassembly was made possible by the presence of enzymes, which allow the required pH variations.7 Alternatively, pH variations effecting aggregation/disaggregation phenomena have also been achieved by bubbling gaseous CO2 in solution.8 More in general, a number of systems including dissipative systems9,10 could take advantage from time-programmable control of pH.
Lately, some activated carboxylic acids such as 2-cyano-2-phenylpropanoic acid,10,11 and its derivatives,12 or trichloroacetic acid,13 have been successfully used to promote time controlled cycles of motions of molecular machines, both switches and motors, in organic media. In all cases, the movements within the molecular machine were due to sequential protonation and deprotonation of a nitrogen base present in the machine skeleton. Initially, the carboxylic acid donates the proton to the nitrogen base of the molecular machine that passes from state A to state B with a consequent large amplitude motion. The subsequent decarboxylation of the resulting carboxylate generates a strong base, which is able to take back the proton from the protonated molecular machine. As a consequence, the state A of the latter is restored with a second large amplitude movement that closes the cycle of motion.9
Taking inspiration from these studies, which were all carried out in organic solvents, we wanted to devise a chemical system able to modulate the pH in water over time. After some preliminary search,14 we resorted to nitroacetic acid15 (O2NCH2CO2H) as putative carboxylic acid to achieve such pH control. Hereafter we report the results of our investigation.
In an ideal, prototypal experiment, a 0.010 M water solution of NaOH is prepared. To this, solution one or more mol equiv. of nitroacetic acid are added so that the pH rapidly changes from basic to acidic. At this point, the mono-anion 2 begins to irreversibly decarboxylate, and the resulting nitromethide ion, 4, (pKa = 10.28 at 20 °C and μ = 0.50 M)16 takes back the proton from the solvent to eventually raise the pH of the solution back to basic values. As illustrated in Scheme 1, however, the pH is affected not only by the aforementioned processes, but also by the formation of carbon dioxide. If the experiment is carried out by bubbling Ar (as well as N2 or air) into the solution by means of a capillary connected to a gas cylinder, one can avoid carbon dioxide oversaturation of the solution. Under this condition, the concentration of carbon dioxide remains at the constant value of its solubility in water, i.e. 1.63 × 10−5 M at 20 °C, during the entire course of the experiment.§17 It is known, that more than 99% of CO2 exists in solution as the dissolved gas and less than 1% as carbonic acid H2CO3, which partially dissociates to give HCO3− and CO32− ions. Owing to the difficulty of differentiating between CO2 and H2CO3 in solution, the apparent first ionization constant, defined as [HCO3−][H3O+]/([CO2] + [H2CO3]), is generally considered. Accordingly, carbon dioxide can be treated as a diprotic acid with an apparent pKa1 = 6.35, and a pKa2 = 10.33 at 25 °C.17
By assuming that all the equilibria in Scheme 1 are much faster than the process of decarboxylation, all the kinetic and equilibrium constants necessary to investigate the dynamics of the system are known. However, owing to its formidable complexity, the system is not amenable to a simple analytic solution. For this very reason, in order to forecast the behavior of the system, we recurred to numerical integration using the program COPASI (COmplex PAthway SImulator), which is an open-source software application for simulation and analysis of chemical networks and their dynamics.18 We simulated two experiments: in the first experiment (a), 1 equiv. of nitroacetic acid 1 is added to an initial 0.010 M NaOH solution (pH = 12), while, in the second experiment (b), 2 equiv. of nitroacetic acid 1 are added to the same initial NaOH. The kinetics of pH changes for the two experiments are reported in Fig. 1, whereas the kinetics of the concentration changes of all the species appearing in Scheme 1 are reported in the ESI (pages S3–S7).†
The kinetics in Fig. 1 are readily understood. In experiment (a), nitroacetic acid, 1, is completely neutralized by NaOH, and the resulting nitroacetate, 2, solution shows an initial pH of 5.43. However, this solution is not buffered, and thus the progress of the decarboxylation reaction causes a steep pH jump. After the initial jump, the pH changes more slowly because of the buffer pairs present in solution (2/3, 5/4, CO2 aq/HCO3−, and HCO3−/CO32−) up to the final pH value of 9.05 which is determined by the latter three buffer pairs. The final dominant species in solution are nitromethane, 5, and the hydrogen carbonate anion. In experiment (b), nitroacetic acid, 1, is partially neutralized by NaOH, however since 1 is a relatively strong acid, initially the ratio [2]/[1] is larger than 1 (actually the ratio is about 5) and the resulting pH = 2.18. Until the undissociated acid 1 is present, the solution is buffered and the progress of the decarboxylation reaction causes a limited pH increase. When the acid 1 is completely neutralized, the solution is no more buffered and a steep pH jump ensues. After the jump, the pH increases more slowly as in experiment (a) and for the same very reasons. It is noteworthy that the final pH value is very similar in the two experiments. Comparing the two kinetic profiles, it is immediately apparent that the experiment with excess acid benefits from a larger pH jump, and therefore the presence of excess acid is more suitable for making molecular machines or other dissipative systems work. Moreover, Fig. 1 suggests that by increasing the excess acid, the time necessary to reach the pH jump can be modulated at will. Thus, a system based on nitroacetic acid could be exploited for chemically programming subsequent pH changes of the kind pH1(high)–pH2(low)–pH3(high). The choice of the absolute and relative initial concentrations of NaOH and nitroacetic acid dictates values of pH1 and pH2.
Fig. 2 shows the real version of one of the prototypal experiments previously designed. Here, the pH of an initial 0.010 M NaOH solution is followed over time (T = 20 °C, μ = 0.50 M with NaCl¶) by means of a glass microelectrode. At t = 0 pH is 12. When 2 mol equiv. of 1 are added, the pH immediately drops to about 2.05 and, subsequently, starts to rise up again drawing a sigmoidal curve. After 30 min the curve tends to sit down and at t = 60 min a pH of 9.0 is read (the kinetic profile described by the experimental points is highly reproducible). The initial and the final pH values are in good agreement with those predicted by the simulation experiment, however the experimental profile is less steep than the theoretical one (compare profiles in Fig. 2 and 1b). This behavior very likely depends on the difficulty for the glass microelectrode to follow the pH changes in solution in real time, probably because of the formation of gas microbubbles over its surface.||
Interestingly, the initial and final pH values of two runs carried out under the same conditions, but bubbling either argon or air in solution, were the same (see Fig. S8†), supporting the assumption that in both cases the solution is saturated by CO2 throughout the experiment.
On the other hand, different time-programmable evolutions after the initial pH jump were achieved by fixing the initial concentration of the base and varying the amount of added acid 1. In Fig. 3b the case in which the [NaOH]o is fixed to 0.010 M and the amount of added acid 1 is varied (1 mol equiv. black trace, 2 mol equiv. blue trace, 3 mol equiv. red trace and 4 mol equiv. green trace) is reported. It is clear that the time required to reach the plateau value, round about 9 in the different cases, increases on increasing the initial concentration of acid 1, which can be varied at will. It is noteworthy that once the excess acid is over, all the kinetic profiles are practically superimposable. Interestingly, although the time needed to reach the plateau value is more or less the same, when 1 or 2 mol equiv. of acid 1 are added to 0.010 M NaOH (black and blue traces of Fig. 3b, respectively), the pH time-profiles of the two solutions predictably and strongly diverge for the reasons discussed in the previous section.**
Eventually, there is a chance to adjust the concentrations of NaOH and acid 1 in order to obtain exactly the same pH (pH2(low)) after addition of 1, and regulate the time necessary to reach the final pH. This kind of control is illustrated in Fig. 3c that reports two experiments with specific concentrations of NaOH and 1 (0.010 M NaOH + 0.026 M 1 and 0.036 M NaOH + 0.060 M 1, red and blue trace, respectively) aimed at obtaining the same pH2(low) = 2.
The effectiveness of nitroacetic acid 1 in promoting a time-programmed pH1(high)–pH2(low)–pH3 = 9 sequence of a water solution can also be visually illustrated using classical pH indicators. Fig. 4 shows selected pictures taken at different times of a 0.010 NaOH solution before and after addition of 0.020 M acid 1 in the presence of bromocresol green (pKa = 4.8). The corresponding movie can be found in the ESI page S10,† together with the movie related to the same experiment carried out in the presence of methyl red (pKa = 5.1), see ESI page S11.†
Fig. 4 Selected pictures of a 0.010 M NaOH solution before and after addition of 0.020 M acid 1 in the presence of bromocresol green (pKa = 4.8). Second pictures from the left (t = 0) were taken immediately after addition of acid 1 (see ESI page S10 for related movie†). |
Fig. 5 shows that it is possible to repeat more times the sequences pH1(high)–pH2(low)–pH3 = 9 by successive additions of nitroacetic acid. After every two shots, it is convenient to reset the initial pH (pH1) by adding the required amount of NaOH in order to have a full efficiency of the system.††
Scheme 2 State of charge of 7 at different pHs and corresponding affinity for alpha-cyclodextrin 6. Relative emission intensities are also schematically shown. |
Fig. 6a shows the fluorescence and the corresponding pH variations as a function of time when acid 1 (2.00 mM) is added to a solution containing NaOH (1.00 mM), 6 (2.00 mM) and 7 (0.0050 mM). Before addition of acid 1, the solution pH is 11.0 (t = 0, red trace), 7 is in its negative form and is largely outside the cavity of cyclodextrin, and consequently the fluorescence emission is low (t = 0, blue trace). Immediately after the addition of acid 1 (t = 4 min), the pH drops to 3.0 and 7, which is now mostly in its positive form, is more efficiently bound by 6, resulting in an immediate strong enhancement of the fluorescence.
Fig. 6 (a) Time dependence of pH and fluorescence emission intensity (λexc = 282 nm, λem = 338 nm) of 0.0050 mM aminobenzoic acid 7 in the presence of 2.0 mM alpha-cyclodextrin 6 (1.0 mM NaOH, μ = 0.50 M; T = 20 °C) before and after the addition (at t = 4 min) of 2.0 mM nitroacetic acid 1. At t = 84 min the amount of NaOH required to restore pH 11 was added. (b) Corresponding I/Iovs. pH from t = 4 min to t = 84 min obtained from graphs in Fig. 6a. Red and green curves are a guide to the eye. |
Then, decarboxylation takes place and pH begins to increase. After 16.5 min from start, at pH 4.9, the concentration of the zwitterion is at its maximum value, concentration of the inclusion complex 6·7 is at its maximum too, and, correspondingly, the fluorescence reaches its highest value (see also Fig. 6b). From 16.5 min (pH 4.9) onwards, the pH still rises, and the zwitterionic form of 7 is partially transformed into the negatively charged one. The above complex partially dissociates and the fluorescence decreases again. Interestingly, the fluorescence variation shows an opposite behavior in the absence of cyclodextrin 6 (see Fig. S12†). For a technical difficulty due to the instrumental set up, the experiment was carried out without argon bubbling. This is the reason why the final pH does not rise to the usual value of 9 but reaches 7.4 at 84 min from start. Furthermore, the pKa of the second dissociation of p-aminobenzoic acid 7 is strongly enhanced by inclusion into cyclodextrin,20 and thus at the end of the monitoring (pH = 7.4), 7 is still partially in its zwitterionic form.§§ However it is evident that an association–disassociation cycle of the kind 6 + 7 → 6·7 → 6 + 7 has been achieved via a nitroacetic acid driven time-programmed pH1(high)–pH2(low)–pH3(high) sequence (namely, 11–3–7.4).
Footnotes |
† Electronic supplementary information (ESI) available: Details on kinetic simulations, experimental procedures, control experiments and multimedia files. See DOI: 10.1039/d1sc01196k |
‡ These authors contributed equally. |
§ The concentration of CO2 in the atmosphere is 415 ppm (value taken from the website https://www.co2.earth on January 18th, 2021), meaning that its partial pressure is 4.15 × 10−4 atm. Its solubility in water at 20 °C is 0.0393 mol per kg water, at a partial pressure of 1 atm CO2 (see ref. 17). At such low pressure Henry's law is followed (see ref. 17), and thus the molar concentration of CO2 in water at 20 °C can be calculated as the product of its solubility times its partial pressure in the atmosphere. |
¶ Reaction time and profile do not change if ionic strength is decreased by a factor of 5 (NaCl 0.100 M). |
|| A delayed answer of the glass electrode was effectively proved by following the pH evolution of a solution in which 0.010 M NaOH and 0.020 M acid 1 were reacted in the presence of bromocresol green. The end point was visually detected by the color change of the indicator a non-negligible time before it was recorded by the electrode pH reading. |
** In an additional series of runs, 0.015 M 1 was added to NaOH of varying concentration, namely, 0.0075, 0.015, 0.030 and 0.050 M. The results obtained, which are in line with expectation, are reported at page 8 of the ESI (Fig. S7).† |
†† If NaOH is not added after the second pulse of acid 1, the system loses efficiency and doesn't reach the final pH = 9 anymore, probably due to the accumulation of reaction products (see Scheme 1) that someway affect the decarboxylation process (see Fig. S10 in the ESI†). |
‡‡ pKa1 and pKa2 of p-aminobenzoic acid are 2.39 and 4.91, respectively (M. S. K. Niazi and J. Mollin, Bull. Chem. Soc. Jpn., 1987, 60, 2605–2610), and are expected to considerably change in the presence of excess cyclodextrin. For example, complexation into beta-cyclodextrin, which is a weaker binder than alpha-cyclodextrin, sensibly enhances the pKa2 (ΔpKa = 1.5) of p-aminobenzoic acid. |
§§ At t = 84 min, addition of the required amount of NaOH to reset the pH solution to 11 restored approximately the initial value of the fluorescence. |
This journal is © The Royal Society of Chemistry 2021 |