I’ve started so I’ll finish. Kinetic isotope effect models for a general acid as a catalyst in the protiodecarboxylation of indoles.

Earlier I explored models for the heteroaromatic electrophilic protiodecarboxylation of an 3-substituted indole, focusing on the role of water as the proton transfer and delivery agent. Next, came models for both water and the general base catalysed ionization of indolinones. Here I explore general acid catalysis by evaluating the properties of two possible models for decarboxylation of 3-indole carboxylic acid, one involving proton transfer (PT) from neutral water in the presence of covalent un-ionized HCl (1) and one with PT from a protonated water resulting from ionised HCl (2).

The original study[1] noted that the rate of decarboxylation fitted well to the kinetic expression: rate = {a + b[L3O+]/(1 + c[L3O+])}[indole], where L can be H or D. Experimentally, [L3O+] is controlled by adding a strong general acid such as HCl, which when the appropriate number of water molecules are added[2] fully ionizes to H3O+.OH. Now for B3LYP+D3/Def2-TZVPD/SCRF=water calculations:

  • Model takes the pure water model and adds HCl (blue above) via hydrogen bonding to the H2O that is transferring the proton to the indole ring. Three water molecules are hydrogen bonding to the carboxylate oxygens to create a bicyclic network in which a ring of either 8 or 10 atoms can act as the proton relay structure. The question now arises whether the proton relay takes the longer (red) route or the slightly shorter green route.
  • Isomeric model 2 uses H3O+ for proton transfer, with an adjacent Cl to complete the ion-pair.
Model ΔG298 (0.044M) DataDOIs kH/kD[3]
1 27.4 [4],[5],[6],[7] 5.69
2 16.8 (14.8) [5],[8],[9] 2.45

Reactant as a non-ionised covalent HCl. reactant as an isomeric ionized H3O+.Cl–  beng 2.0 kcal/mol higher in energ within this solvation model. Note added in proof. A significantly lower form of the reactant has subsequently been located which increases the free energy barrier to 22.1 kcal (vis 22.0 actually measured!). Diacussion of this can be seen in the associated post here.

  1. An IRC for Model 1 shows that the proton relay takes the red path, whereas without the HCl the green path is followed.

     

     

     

    The transition state free energy however is ..

  2. 10.6‡ or 12.6 kcal/mol higher than model (click on the image below to load a 3D model). The general acid catalysed model is therefore preferred. The difference in free energy between the two models corresponds to a rate acceleration of >106, which is indeed similar to that observed[1].

The clincher comes with calculation[3] of the kinetic isotope effects (KIE). For general acid catalysis, they were measured as kH/kD ~2.5.[1]

  • For model 1, using an un-ionised reactant and un-ionised transition state, the calculated KIE is 5.69 (very similar to that calculated for the water catalysed reaction, 5.66) but not a good fit to experiment.
  • For model 2, using the same un-ionised reactant but an ionised transition state, KIE = 2.04, a much better fit.
  • For model 2, using ionised reactant AND transition state, KIE = 2.45, an even better fit to experiment.

So we now have a model for the general acid catalysed decarboxylation of a 3-indole carboxylate which agrees with both the kinetic behaviours and the isotope effects measured for this reaction. Since the barrier is a relatively large one, proton tunnelling may play a lesser role in this interpretation, and the stage is set to use this model to e.g. explore how isotope effects are indeed influenced by tuning the reactivity using ring substitutents, the original purpose of my researches all those years ago. Perhaps the catch phrase I’ve started so I’ll start is now more apposite.


Acknowledgments

This post has been cross-posted in PDF format at Authorea.

References

  1. B.C. Challis, and H.S. Rzepa, "Heteroaromatic hydrogen exchange reactions. Part 9. Acid catalysed decarboxylation of indole-3-carboxylic acids", Journal of the Chemical Society, Perkin Transactions 2, pp. 281, 1977. http://dx.doi.org/10.1039/P29770000281
  2. A. Vargas‐Caamal, J.L. Cabellos, F. Ortiz‐Chi, H.S. Rzepa, A. Restrepo, and G. Merino, "How Many Water Molecules Does it Take to Dissociate HCl?", Chemistry – A European Journal, vol. 22, pp. 2812-2818, 2016. http://dx.doi.org/10.1002/chem.201504016
  3. Henry Rzepa., "Ionic model for general acid catalysed decarboxylation", 2016. http://dx.doi.org/10.14469/hpc/204
  4. Henry S Rzepa., "C 9 H 16 Cl 1 N 1 O 6", 2016. http://dx.doi.org/10.14469/ch/191792
  5. Henry S Rzepa., "C 9 H 16 Cl 1 N 1 O 6", 2016. http://dx.doi.org/10.14469/ch/191795
  6. Henry S Rzepa., "C 9 H 16 Cl 1 N 1 O 6", 2016. http://dx.doi.org/10.14469/ch/191794
  7. Henry S Rzepa., "C9H16ClNO6", 2016. http://dx.doi.org/10.14469/ch/191767
  8. Henry S Rzepa., "C 9 H 16 Cl 1 N 1 O 6", 2016. http://dx.doi.org/10.14469/ch/191790
Henry Rzepa

Henry Rzepa is Emeritus Professor of Computational Chemistry at Imperial College London.

View Comments

Recent Posts

Internet Archeology: reviving a 2001 article published in the Internet Journal of Chemistry.

In the mid to late 1990s as the Web developed, it was becoming more obvious…

1 month ago

Detecting anomeric effects in tetrahedral carbon bearing four oxygen substituents.

I have written a few times about the so-called "anomeric effect", which relates to stereoelectronic…

1 month ago

Data Citation – a snapshot of the chemical landscape.

The recent release of the DataCite Data Citation corpus, which has the stated aim of…

2 months ago

Mechanistic templates computed for the Grubbs alkene-metathesis reaction.

Following on from my template exploration of the Wilkinson hydrogenation catalyst, I now repeat this…

2 months ago

3D Molecular model visualisation: 3 Million atoms +

In the late 1980s, as I recollected here the equipment needed for real time molecular…

3 months ago

The Macintosh computer at 40.

On 24th January 1984, the Macintosh computer was released, as all the media are informing…

3 months ago