Guest Editor: C. A. Schalley
Beilstein J. Org. Chem. 2012, 8, 81–89. https://doi.org/10.3762/bjoc.8.8
Received 30 Sep 2011, Accepted 21 Dec 2011, Published 12 Jan 2012
For the condensation of anions such as phosphate and ADP to form ATP and water, nature employs sophisticated supramolecular systems to overcome coulomb repulsion and activation barriers. For an attempt to create a simple, analogous chemical system, the dimerization of vanadate is probably the simplest model. We have investigated Zn-benzylcyclene which favors the dimerization thermodynamically as shown by NMR titration. Moreover, EXSY NMR experiments reveal that the vanadate dimer is also kinetically stabilized with respect to hydrolysis by complexation with Zn-cyclene.
Keywords: EXSY NMR; NMR titration; supramolecular chemistry; vanadate condensation; 51V NMR; Zn-cyclene
Driving endergonic reactions with external energy sources is one of the challenging and so far unsolved problems in supramolecular chemistry. The best known example in nature is probably the condensation of phosphate and ADP to ATP, driven by a proton gradient across a membrane . There are a number of preconditions that have to be met for such an energy driven condensation of inorganic anions. At least three steps (four stages) are required (Figure 1): a) Moving both complexated anions together, b) bond formation (condensation), and c) product release/reactant binding. In the ATP synthase, the energy is put into the system mainly in step c), thus, avoiding product inhibition .
In a first attempt towards the very ambitious goal of creating an artificial ATP synthase model, the a) → b) → c) energy profile has to be carefully designed to operate such a system. Here, we focused on step b). To further simplify the system, we chose the dimerization of vanadate as a model for the phosphate condensation for the following reasons: The condensation of phosphate has a very high activation barrier (disodium phosphate condensates to pyrophosphate at temperatures >250 °C) which is likely difficult to be realized in a prototype artificial system . In contrast to phosphates, mono-, di- and higher vanadates are in equilibrium at room temperature depending on the pH. Moreover, 51V is one of the most sensitive NMR nuclei in the periodic table . The concentration of the different vanadate species can be easily determined even in very diluted solutions by integration of the corresponding NMR peaks . Although vanadates are slightly larger than their corresponding phosphates (V–O bond approximately 0.170 nm, P–O bond approximately 0.152 nm) they exhibit a striking similarity in their chemistry . Vanadates are accepted by a number of phosphatase enzymes . Thus, the plethora of information on both natural and artificial phosphate binding systems can be used to design vanadate coordinating ligands. Vice versa, the vanadate model systems will provide insight into the phosphate condensation, e.g., the formation of ATP .
A large number of thoroughly investigated enzymes hydrolyze phosphates and provide insight into conceivable mechanisms of phosphate or vanadate condensations. In phosphatases, kinases and ATP synthase, the catalyzed transfer of phosphate usually requires the presence of metal cations such as Zn2+, Mg2+ and Mn2+ in the active site . It is generally assumed that the coordination of the metal cations to the phosphate oxygen atoms reduces the Coulomb repulsion in the transition states or in the intermediates, thus, lowering the activation barrier . In tyrosine phosphatase, this role is taken by hydrogen bonding of the phosphate ion to the positively charged arginine side chain . In the key step of the hydrolysis of pyrophosphate by the yeast phosphatase, two Mg2+ ions and arginine H-bonds assist the P–O bond cleavage . A number of artificial enzymes that catalyze phosphate ester hydrolysis have been synthesized and investigated . Kimura et al. demonstrated the acceleration of phosphate ester hydrolysis by Zn2+ complexes of triazacyclododecane and tetraazacyclododecane (cyclene) . The efficiency was increased by tethering two or three Zn-cyclene units in an appropriate distance to allow multiple coordination to the phosphate oxygen atoms . Brown et al. were able to catalyze the cyclization (hydrolysis/intramolecular condensation) of 2-hydroxypropyl aryl phosphate using a tethered dimeric Zn-triazacyclododecane . Lehn et al. catalyzed the condensation of acetyl phosphate and phosphate to pyrophosphate using a protonated macrocyclic polyamine [17-19].
To investigate the fundamental thermodynamic and kinetic effects of metal coordination in inorganic condensation reactions, we focused on the influence of benzyl Zn-tetraazacyclododecane (Zn-benzylcyclene) 1 on the aggregation of vanadate in water.
According to DFT calculations (B3LYP/6-31G*) of a number of different geometries, the monovanadate anion is bound to the Zn-cyclene unit by forming one coordination bond to the Zn2+ ion and two hydrogen bonds to the NH groups of the cyclene ring (Figure 2). The formation of a bidentate complex (two vanadate oxygens coordinating with Zn2+) is considerably less favorable. This is also in agreement with X-ray structures of Zn-cyclene phosphate complexes .
For the analysis of the multicomponent system, 1-D 51V NMR titrations of Na3VO4 with ligand 1 were performed at three different pH values at 25 °C (pH ~9.5 CHES buffer (N-cyclohexyl-2-aminoethanesulfonic acid); pH ~8.5 EPPS buffer (3-[4-(2-hydroxyethyl)-1-piperazinyl]propanesulfonic acid); pH ~7.6 HEPES buffer (2-[4-(2-hydroxyethyl)piperazin-1-yl]ethanesulfonic acid)) (Figure 3, see also Supporting Information File 1). These alkaline pH values were chosen because monovanadate is the main species in this range. Like the monovanadate, the di- and tetravanadate exist primarily in the neutral region.
The 51V NMR titration at pH 9.5 is shown in Figure 3 (for the corresponding spectra recorded at pH 7.6 and 8.5 see Supporting Information File 1). The addition of Zn-benzylcyclene 1 leads to an upfield shift of the monovanadate and the divanadate signal at all pH values that were investigated. This shift indicates that Zn-benzylcyclene 1 is capable of complexing both mono- and divanadate. The monovanadate signal remains the main signal throughout the titrations, however, the signal of the divanadate increases upon addition of Zn-benzylcyclene. By raising the concentration of benzylcyclene 1t from 0 to 7.5 mM in a system with a total concentration of vanadium [V]t of 1.5 mM at a pH value of 9.5, the ratio of vanadium bound in divanadate species rises from 5% to 23%, while the ratio of vanadium bound in monovanadate decreases from 95% to 77%. Hence, Zn-benzylcyclene 1 favors the condensation of monovanadate to divanadate.
To determine the stoichiometry and stability of the most important species in this system, we included all vanadate oxo-anions that are known to exist in aqueous solution [20,21] and the corresponding complexes with ligand 1 in different stoichiometries as basis of a complex formation model. This model was refined with respect to the experimental data.
In the approach mentioned above, the equilibria of the multi component system were described mathematically by mass balance equations. When the equilibrium species are written with H+, VO43−, and Zn-benzylcyclene 1 as components, they are formed according to the general expression (Equation 1) :
The formation constants of the species βpqr are defined using the concentrations of 1, VO43− and of the protons according to Equation 2.
The substitution of the complex species concentration in the mass balances of the components leads to the mass balance equation system (3). The concentrations of the free components are the variables; the stability constants are the parameters to be optimized.
In the 51V NMR titrations for the analysis of the vanadate/Zn-benzylcyclene system the total vanadium concentration [V]t at each pH value was 1.5 mM, and the total concentration of 1 (t) was varied from 0 to 7.5 mM. The NMR data were obtained by careful deconvolution and integration of the signals using routines included in the MestreNova program . The mathematical analysis of the NMR data was performed with the program system LAKE  for chemical equilibrium analysis and with MS-Excel . Data of all titrations (78 titration points in total) were included in the calculation. The formation constants for the vanadium oxo-anions of the VO43−/H+-subsystem , the pKa value of Zn-benzylcyclene 1  as well as the pKa values of the buffers  were adopted from literature. The influence of the ionic medium was taken into account in the calculation by means of the Davies equation . The most important ternary species in the equilibrium model that fit the recorded data are presented in Table 1.
Table 1: Stoichiometry and formation constants of the complex species in the proton/vanadium/Zn benzylcyclene system.
|p, q, ra||formulab||lg βpqrc (3σ)|
|1, 1, 1||HVO41||16.7 (0.3)|
|2, 2, 1||V2O712−||32.2 (0.1)|
|2, 2, 2||V2O712||35.7 (0.8)|
|3, 2, 1||HV2O71−||40.6 (0.3)|
|3, 2, 2||HV2O7122+||43.6 (0.4)|
|8, 4, 2||V4O1212||99.9 (0.7)|
|8, 4, 4||V4O12144+||105.6 (1.8)|
aStoichiometric coefficients according to the formation from the components H+, VO43− and 1 (Equation 1); bformula of the oxo-anion complex. The water formed according to the formation reaction from the components is not given; cionic strength, I = 150 mM, T = 298 K.
According to these data, the association constant for the protonated monovanadate HVO42− with 1 is about 103.8 M−1, while the association constant for V2O74− is Ka = 106.2 M−1. The pKa value for the coordinated divanadate is 8.4 and therewith, 1.9 pKa units lower than for the uncoordinated V2O74− ion. The binding of a second Zn-benzylcyclene unit proceeds with an association constant of 103.5 M−1. Consequently, the affinity of the divanadate for both the proton and the second Zn-benzylcyclene is lowered by the coordination of the first Zn-benzylcyclene unit. A speciation diagram for the titration at pH 9.5 based on these data is given in Figure 4. Fits of calculated and experimental 51V NMR integral data at pH 7.6, 8.5, and 9.5 are given in the Supporting Information File 1.
The formation constants given in Table 1 can be used to calculate the free enthalpy (ΔG10, ΔG20, ΔG30) of the condensation reactions (Equations 4–9):
The analysis reveals that the free reaction enthalpy for the condensation of hydrogen monovanadate is lowered by the complexation of the vanadium oxo-anions with Zn- benzylcyclene. In the absence of the ligand, the dimerization of monovanadate is almost thermoneutral (ΔG10 = 0.9 kJ mol−1). However, the condensation of a free monovanadate with a monovanadate bound to Zn-benzylcyclene is exothermic (ΔG20 = −14.8 kJ mol−1) and the reaction of two complexed monovanadates is less exothermic (ΔG30 = −13.1 kJ mol−1). Hence, the condensation is favored by the complexation with Zn-benzylcyclene and one equiv of the ligand is more effective than two.
Further insight into the stabilization of divanadate by Zn-benzylcyclene is obtained by analyzing the ligand exchange reactions (Equation 10 and Equation 12). While the transfer of a ligand from hydrogenmonovanadate to divanadate is exergonic (ΔG100 = −13.7 kJ mol−1), the transfer of a second ligand from another hydrogen monovanadate is endergonic (ΔG120 = 1.1 kJ mol−1).
Accordingly, a complexed divanadate is destabilized by the coordination of a second ligand.
The influence of the complexation on the rates of condensation and hydrolysis of the involved vanadium oxo-anions was examined by 51V EXSY NMR experiments [27,28]. 51V EXSY NMR data were collected for solutions buffered with EPPS and HEPES and 1.5 mM total vanadium concentration at 25 °C. For each buffer solution, the amount of Zn-benzylcyclene 1 was varied. The mixing time was 1 ms (Figure 5).
The data show a decrease for the pseudo first order rate constant of the condensation (k12 in Equation 14) as well as for the hydrolysis (k21 in Equation 15), while the influence on the rate constant of the hydrolysis is stronger.
In the HEPES buffered solution (pH 7.6), the rate constant for the condensation decreased from 73 s−1 to 68 s−1 upon the addition of Zn-benzylcyclene (3 mM), while the rate constant for the hydrolysis decreased from 206 s−1 to 89 s−1. Compared to that, in the EPPS buffered solution (pH 8.0), the rate constant for condensation decreased from 23 s−1 to 14 s−1 upon the addition of Zn-benzylcyclene 1 (3 mM), while the rate constant for the hydrolysis decreased from 164 s−1 to 70 s−1.
The fact that the condensation of monovanadate is slower in the presence of the ligand is probably due to steric reasons. The nucleophilic attack on the complexed monovanadate is only possible from the non-shielded side. Moreover, the role of the monovanadate acting as the electrophile is defined by complexation which reduces the probability of effective collisions and, thus, decreases the entropy of activation. Both effects, obviously, overcompensate the increased reactivity of the vanadate acting as the electrophile which is induced by complexation (decreased Coulomb repulsion of the two negatively charged reactants). The same arguments are true for the hydrolysis of divanadate (attack of water). However, the thermodynamic stabilization of divanadate further increases the barrier of hydrolysis (provided that the transition state is less effected). Therefore, the rate constant of the divanadate hydrolysis is lowered more strongly by complexation than the rate constant for the condensation of monovanadate.
In conclusion, the condensation of monovanadate to divanadate is considerably favored in the presence of the ligand 1 (Zn-benzylcyclene). In the absence of the ligand at pH 9.5 the condensation is almost thermoneutral (ΔG10 = 0.9 kJ mol−1) whereas upon addition of the ligand the dimerization reaction is exergonic (ΔG20 = −14.8 kJ mol−1). This is mainly due to the fact that the free binding enthalpy of ligand 1 to divanadate is 13.7 kJ mol−1 larger than the complex formation with monovanadate.
Complexation with 1 also favors the formation of divanadate kinetically by decreasing its rate of hydrolysis from 164 s−1 to 70 s−1 upon addition of ligand 1, as determined by 51V EXSY NMR experiments.
Dimerization, however, is less favorable if both reacting monovanadate units are complexed with ligand 1 (ΔG30 = −13.1 kJ mol−1) compared to the reaction of a complexed vanadate with a free monovanadate (ΔG20 = −14.8 kJ mol−1). For the design of an artificial system driving endergonic condensations (Figure 1), we draw the conclusion that only one of the binding sites should provide metal coordination to the vanadate (or phosphate) and the other binding site should associate the nucleophilic anion by neutral hydrogen bonds.
Similar effects obviously hold for natural systems – particularly ATP synthase, phosphatases and kinases catalyzing the condensation, hydrolysis and phosphate transfer in phosphate esters and phosphate oligomers – which are unsymmetric with respect to metal complexation of the substrates. In the P–O bond forming/cleaving step of yeast pyrophosphatase, the phosphate group acting as the electrophile is activated by coordination to Mg2+ as well as by hydrogen bonding to a guanidinium group . The electrophilicity of ADP in ATP synthase is even further increased by electron transfer to Mg2+ forming a radical ion pair . Similar to the ATP synthase, a simplified artificial system should prevent product inhibition by actively releasing the product. This could be realized by a (photo)switchable ligand, see step a) and c) in Figure 1.
1-Benzyl-1,4,7,10-tetraazacyclododecane (benzylcyclene): Benzylchloride (1.00 mL, 8.69 mmol) and cyclene (6.00 g, 34.8 mmol) were dissolved in dry chloroform (240 mL). Under nitrogen atmosphere, triethylamine (1.45 mL, 10.4 mmol) was added and refluxed for 20 h. After cooling, chloroform (200 mL) was added, and the organic phase was washed three times with 1 N sodium hydroxide solution (200 mL) and three times with distilled water (200 mL). The organic phase was dried over potassium carbonate and the solvent was removed under reduced pressure .
For efficient purification, the crude product (2.27 g) was protected with di-tert-butyl dicarbonate (3 equiv) in dry dichloromethane (100 mL). After removal of the solvent the product was purified by column chromatography (cyclohexane:ethylacetate = 1:1, Rf = 0.58). Yield: 4.24 g (87%); 1H NMR (500 MHz, CD2Cl2) δ 7.3–7.2 (m, 5H, CHaryl), 3.65 (s, 2H, Ar-CH2), 3.52 (br s, 4H, CH2), 3.4–3.1 (m, 8H, CH2), 2.60 (br s, 4H, CH2), 1.41, 1.36 (2 s, 27H, CH3) ppm; 13C NMR (125 MHz, CD2Cl2) δ 155.7 (C=O), 137.5 (Caryl,q), 128.5, 127.6, 123.6, (Caryl), 80.0, 79.5, 79.4 (C(CH3)3), 57.5 (Ar-CH2), 50.7, 50.1, 48.2 (CH2), 28.8, 28.6 (CH3) ppm.
To release the free amine, the boc-protected benzylcyclene (3.76 g, 6.69 mmol) was dissolved in dry dichloromethane (30 mL). Trifluoroacetic acid (5 mL) was added under a nitrogen atmosphere and the mixture stirred at room temperature for 12 h. Under ice cooling sodium hydroxide (3 g) and distilled water (30 mL) were added. The organic layer was separated, and the aqueous phase extracted five times with 50 mL dichloromethane. The combined organic phases were dried over potassium carbonate and the solvent was removed under reduced pressure to yield a colorless solid (1.52 g, 86%); 1H NMR (500 MHz, CD2Cl2) δ 7.3–7.2 (m, 5H, CHaryl), 3.64 (s, 2H, Ar-CH2), 2.8–2.4 (m, 16H, CH2) ppm; 13C NMR (125 MHz, CD2Cl2) δ 139.9 (Caryl,q), 129.5, 128.6, 127.4 (Caryl), 60.1 (Ar-CH2), 51.8, 48.1, 46.6, 46.1 (CH2) ppm.
In order to regain the excess of cyclene used in the first step, the aqueous phase was concentrated under reduced pressure and neutralized with 1 N HCl until a white precipitate formed. The precipitate was filtered to yield pure cyclene (2.8 g). 1H NMR (200 MHz, DCl3) δ 2.65 (s, 16H, CH2) ppm.
Zinc-1-benzyl-1,4,7,10-tetraazacyclododecane nitrate (1): 1-Benzyl-1,4,7,10-tetraazacyclododecane (1.00 g, 3.81 mmol) was dissolved in methanol (20 mL). A zinc nitrate solution in water (49.7 mM, 76.6 mL) was added and the solution was stirred for 1 h at room temperature. The solvent was removed under reduced pressure and the colorless solid was dried. Yield: 100 %; 1H NMR (500 MHz, D2O) δ 7.4–7.3 (m, 5H, CHaryl), 3.94 (s, 2H, CH2), 3.2–3.0 (m, 4H, CH2), 2.9–3.0 (m, 12H, CH2) ppm; 13C NMR (125 MHz, D2O) δ 138.5 (Caryl,q), 131.3, 129.1, 128.7 (Caryl), 56.7 (Ar-CH2), 49.1, 44.5, 43.6, 42.2 (CH2) ppm; ESI m/z (%): 363.1 (100), 364.1, 365.1, 366.1, 367.1, 368.1 [C15H26N4Zn(OH2)2], 326.1 (78), 327.1, 328.1, 329.1, 330.1, 331.1 [C15H25N4Zn].
For the NMR titrations three different buffer stock solutions were used: pH 7.6: HEPES (4-(2-hydroxyethyl)piperazine-1-ethanesulfonic acid), pH 8.5: EPPS (4-(2-hydroxyethyl)-1-piperazinepropanesulfonic acid) and pH 9.5: CHES (2-(cyclohexylamino)ethanesulfonic acid). The stock solutions had a total buffer concentration of 100 mM and 50 mM of sodium hydroxide. They were prepared by dissolving the appropriate amount of buffer and sodium hydroxide in 90% deionized water and 10% deuterium oxide. To obtain the vanadate stock solution (50 mM) sodium orthovanadate (92.0 mg) was dissolved in the respective buffer solution (10 mL). For the receptor stock solution (10 mM) the receptor 1 (45.2 mg) was dissolved in 10 mL of the respective buffer solution.
All measurements were carried out at a 1.5 mM total vanadate concentration and a 100 mM buffer concentration. The samples of the 51V NMR titrations were prepared by mixing the appropriate amount of the vanadate stock solution and the receptor 1 stock solution in an NMR tube. Finally, the volume of each sample was brought to a volume of 600 µL by addition of the appropriate amount of the respective buffer stock solution. All NMR spectra were recorded on a Bruker FT-NMR-spectrometer DRX 500. Vanadium chemical shifts are reported relative to the external reference VOCl3 (0 ppm).
The EXSY measurements were performed at two different pH values (pH 7.6 and pH 8.0) and at a total vanadate concentration of 1.5 mM. The concentration of 1 was 3.0 mM. For all stock solutions deionized water was used as solvent. The samples were prepared in the following way: The buffer stock solutions had a total buffer concentration of 500 mM (HEPES for pH 7.6, EPPS for pH 8.0) and 250 mM of sodium hydroxide. The vanadate stock solution contained of 50 mM sodium orthovanadate, and the concentration of the Zn-benzylcyclene 1 stock solution was 10 mM. For each experiment 60 μL of the vanadate stock solution, 400 μL of the respective buffer stock solution, and 600 μL of the Zn-benzylcyclene 1 stock solution were mixed with 100 μL deuterium oxide and 840 μL deionized water. The NMR tubes were then filled with 600 μL of the mixture. In all experiments a mixing time of 1 ms was applied.
|Supporting Information File 1: 51V NMR spectra and data fitting results.|
|Format: PDF||Size: 256.0 KB||Download|