Posts Tagged ‘Conformational isomerism’

Organocatalytic cyclopropanation of an enal: (computational) assignment of absolute configurations.

Saturday, September 1st, 2018

I am exploring the fascinating diverse facets of a recently published laboratory experiment for undergraduate students.[1] Previously I looked at a possible mechanistic route for the reaction between an enal (a conjugated aldehyde-alkene) and benzyl chloride catalysed by base and a chiral amine, followed by the use of NMR coupling constants to assign relative stereochemistries. Here I take a look at some chiroptical techniques which can be used to assign absolute stereochemistries (configurations).

I will focus on the compound 4a, the major stereochemical product of this student laboratory reaction, with the stereochemistry as represented in e.g. the abstract of the main article[1] and shown below with added CIP (Cahn-Ingold-Prelog) notation as (1S,2R,3R);

Its enantiomer (not shown in the article) is of course;

In the article supporting information[1]), the major diasteromer of 4a deriving from use of the S stereoisomer of the prolinol catalyst is reported as having an optical rotation (ORP) [α]D25 of -62.4°, p6 or -58.1°, p5), but the stereo-labels are not added there. On  p1 (“based on a student report”) 4a was however labelled as (1R,2S,3S) and the chirality (S) of the catalyst used was also noted in the adjacent experimental procedure. One might then reasonably match (1R,2S,3S)-4a to the S-catalyst and hence (1S,2R,3R)-4a to the R-catalyst.  However, in a laboratory environment where both S and R catalysts are in circulation, it is always useful to have procedures available for independent checks.

There are two methods of assigning absolute chirality, crystallography and chiroptical spectroscopy. The former does require crystalline samples; the latter can use solutions. To cut to the chase, the former method was used for a related compound where the n-heptyl group above is replaced by a p-chlorophenyl substituent (perhaps because the latter imparts suitable crystallinity). On p S123 of the SI of an earlier article[2] the assignment for the p-chlorophenyl derivative was as (1R,2S,3S)-4a for S-catalyst (see DOI: 10.5517/ccdc.csd.cc1mcqg5 OZAXEU). But this procedure is not entirely foolproof; the stereochemistry is decided by interactions between often bulky substituents at the transition state and it might be that here the p-chlorophenyl derivative has different properties from n-heptyl. Moreover bulk solutions may be different in their composition from single crystals. So it is useful to obtain independent proof.

An absolute assignment procedure based on chiroptical methods was first famously used by Kirkwood in 1951 (the Fischer convention is confirmed as a structurally correct representation of absolute configuration).[3] Such calculations need to take into account e.g. rotational conformers about the two bonds labelled in red above. In the previous post, I had noted variation of up to 2Hz in the calculated 3JHH coupling constants as a result of this mobility. This variation is probably too small to really influence any relative stereochemical interpretations, but is the same true for chiroptical assignments?

In Table 1 we can see whether this is still true for the predicted optical rotation of compound 4a, using two different functionals for the calculation (B3LYP and M062X respectively). The results rather surprised me; a simple bond rotation of an aryl or carbonyl group can invert the sign of the rotation. Clearly the observed optical rotation of -62.4° arises from a suitable combination of different Boltzmann populations of the individual bond rotamers, but to combine these accurately you would need to know the solution populations themselves very accurately and that is quite a challenge. So at this stage, we do not really have totally convincing independent evidence of whether the observed negative optical rotation corresponds to (1S,2R,3R)-4a or to its enantiomer (1R,2S,3S).

Table 1. Calculated Optical rotations for (1S,2R,3R)-4a. 

FAIR Data DOI: 10.14469/hpc/4678

Conformer

ORP [α]D, B3LYP+GD3BJ/Def2-TZVPP/SCRF=chloroform

ORP [α]D, M062X/Def2-TZVPP/SCRF=chloroform

4 +376 +238
3 -335 -301
2 -247 -223
1 +710 +522

Next, another chiroptical technique, electronic circular dichroism, or ECD. Here, the sign of the difference in absorption of polarized light (Δε), and known at the Cotton effect, characterises the specific enantiomer. The experimental Cotton effect for compound 4a obtained from S-catalyst (known as 3d in the SI, p S142[2]) can be simply summarised as +ve@315nm and -ve@275nm. Comparison with calculated spectra (Figure S17, p S146-7[2])  was performed using a Boltzmann-averaging (albeit based on enthalpies rather than the formally correct free energies), for three significant populations and this procedure matched to (1R,2S,3S).  Since the reported calculations were apparently for gas phase (and replacing n-heptyl with methyl) here I have repeated them in the actual solvent used (acetonitrile) and with the heptyl present. Although the ECD responses can still be severely dependent on the conformation, three of the spectra qualitatively agree that the responses at ~300nm and 260 nm are respectively -ve and +ve. This confirms that (1S,2R,3R)-4a is the wrong enantiomer for S-catalyst and that the correct assignment is therefore (1R,2S,3S), as was indeed reported.[2]

Table 2. Calculated electronic circular dichroism for

 (1S,2R,3R)-4a. FAIR Data DOI: 10.14469/hpc/4678

Conformer

ECD calculation, ωB97XD/Def2-TZVPP

4
3
2
1

It is still true that the overall the fit between chiroptical experiment and theory can be sensitive to the Boltzmann population, as obtained from e.g. ΔΔG = -RT ln [1]/[2]), where 1 and 2 are two different conformers. ΔΔG is a difficult energy difference to compute accurately. Here is a suggested exercise in the statistics of error propagation. How does an error in ΔΔG propagate to the ratio of concentrations of two conformers [1]/[2]? Or, how accurately must ΔΔG be calculated in order to predict conformer populations to say better than 5%.

One more go at chiroptics, this time Vibrational Circular Dichroism, or VCD. The nature of the chromophore is different, but the principle is the same as ECD. I have deliberately truncated the spectrum to cut off all vibrations below 1000 cm-1 (these being the modes associated with group rotations) but to no avail, the four conformations all still look too different to avoid doing a Boltzmann averaging.

Table 3. Calculated VCD spectra for (1S,2R,3R)-4a. 

Conformer Spectrum
4
3
2
1

A modern VCD instrument does have one trick up its sleeve for coping with the conformer problem. The sample (as a thin-film) can be annealed down to very low temperatures before the spectrum is recorded. This effectively removes all higher energy forms, leaving just the most stable conformation as the only species present. However, that is an expensive experiment (and instrument!) to use.

There are perhaps some 2 million scalemic molecules (substances where one chiral form is in excess over the mirror image) for which chiroptical properties have been reported, but probably <50,000 crystal structures where absolute configurations have been assigned. Thus the vast majority of absolute configuration assignments have been done either chiroptically or by synthetic correlations (chemical transformations from molecules of known absolute configuration, with the assumption that you know how each transformation affects the chiral centres present). Given some of the difficulties and challenges noted above, it is tempting to conclude that a significant proportion of those 2 million molecules may have been mis-assigned (I once estimated up to 20%). However, we may conclude that the molecules discussed here are safely assigned correctly! 


No CIP-stereolabels appear in the article itself.[1] Perhaps this assignment is omitted in order to provide a student exercise? There are many errors in stereochemical assignments in the literature. A good many of them may be the result of simple sample mis-labelling.[4] The caption to Figure S17 states All the simulations are for the 1R,2R,3S absolute configuration. This is probably an error and should read 1R,2S,3SA correction of ~+15nm is sometimes applied to these values, but not done here.

 

References

  1. M. Meazza, A. Kowalczuk, S. Watkins, S. Holland, T.A. Logothetis, and R. Rios, "Organocatalytic Cyclopropanation of (<i>E</i>)-Dec-2-enal: Synthesis, Spectral Analysis and Mechanistic Understanding", Journal of Chemical Education, vol. 95, pp. 1832-1839, 2018. https://doi.org/10.1021/acs.jchemed.7b00566
  2. M. Meazza, M. Ashe, H.Y. Shin, H.S. Yang, A. Mazzanti, J.W. Yang, and R. Rios, "Enantioselective Organocatalytic Cyclopropanation of Enals Using Benzyl Chlorides", The Journal of Organic Chemistry, vol. 81, pp. 3488-3500, 2016. https://doi.org/10.1021/acs.joc.5b02801
  3. W.W. Wood, W. Fickett, and J.G. Kirkwood, "The Absolute Configuration of Optically Active Molecules", The Journal of Chemical Physics, vol. 20, pp. 561-568, 1952. https://doi.org/10.1063/1.1700491
  4. H.S. Rzepa, "The Chiro-optical Properties of a Lemniscular Octaphyrin", Organic Letters, vol. 11, 2009. https://doi.org/10.1021/ol901172g

Tautomeric polymorphism.

Thursday, June 1st, 2017

Conformational polymorphism occurs when a compound crystallises in two polymorphs differing only in the relative orientations of flexible groups (e.g. Ritonavir). At the Beilstein conference, Ian Bruno mentioned another type;  tautomeric polymorphism, where a compound can crystallise in two forms differing in the position of acidic protons. Here I explore three such examples.

The term occurs in the title of this article,[1] for a compound known as Omeprazole.

When the bottom structure (the 6-methoxy) is used to search the CSD, two separate series are found. The first of these is UDAVIF (DOI:  10.5517/ccp82qq,  6-Methoxy-2-((4-methoxy-3,5-dimethyl-2-pyridinyl)methylsulfinyl)-1H-benzimidazole). There is no information regarding the absolute configuration of the chiral S-centre. Although the downloaded coordinates show it as R it is probably a racemic mixture. A note added to the structure declares disorder: “Omeprazole exists as solid solutions of the two tautomers. The structure is mixed 5-methoxy/6-methoxy with occupancies 0.078:0.922“, which indicates 7.8% is present as in the upper structure above. 

The second hit is VAYXOI (DOI: 10.5517/ccp82pp, rac-6-Methoxy-2-(((4-methoxy-3,5-dimethyl-2-pyridinyl)methyl)sulfinyl)-1H-benzimidazole) which now contains no disorder; the contaminating 5-methoxy tautomer is no longer present. Perhaps not quite a true tautomeric polymorph, since the 5-methoxy tautomer is never observed in pure form.

This does occur with a second example. DEBFAR[2] represents the keto form on the right which crystallises from methanol, whilst YUYDOL as the enol form on the left crystallises from n-hexane. 

Calculations shed some light on this behaviour. DEBFAR has a computed (DOI: 10.14469/hpc/2591)  dipole moment of 11D, whereas YUYDOL (DOI: 10.14469/hpc/2590) is 2.5D. In chloroform solutions (~half way between the two solvent polarities), the keto form is ~6.1 kcal/mol lower in ΔG than the enol. The crystal packing for the two forms is very different and the differences in this packing must clearly amount to >6.1 kcal/mol to over-ride the lesser stability of DEBFAR in solution.


The final example [3] is illustrated using scheme 2 from that article, one entitled tautomeric species of 4-hydroxynicotinic acid:

The original diagram has two unfortunate bond errors which are NOT reproduced above (and which perhaps are a good topic for discussion in tutorials with students), along with an unusual interpretation of the term tautomerism. The blue arrows above are mine and I suggest the isomerism between the connected species is resonance isomerism, and not tautomerism. So three possible different true tautomers then. Five crystal structures are reported which I list below.

  1. 10.5517/cctswjz (KUXPUP, 4-oxo-1,4-dihydropyridine-3-carboxylic acid, no H2O),  10.5517/ccdc.csd.cc1kfyxv (KUXPUP01 no H2O) and 10.5517/ccdc.csd.cc1kfyzx (KUXPUP02 no H2O)
  2. 10.5517/ccx59s4 (AVEMUK, 4-Oxo-1,4-dihydropyridine-3-carboxylic acid hemihydrate) and  10.5517/ccdc.csd.cc1kfz21 (AVEMUK01)
  3. 10.5517/ccdc.csd.cc1kfz54 (AKIHIN, 4-hydroxypyridin-1-ium-3-carboxylate monohydrate) 
  4. 10.5517/ccdc.csd.cc1kfz10 (AKIHAF, 4-hydroxypyridin-1-ium-3-carboxylate)

KUXPUP and AVEMUK differ only in the presence of one solvent water molecule and both represent tautomer 2 above. AKIHIN and AKIHAF similarly represent tautomer 3 above; both are represented as 3a in the CSD and not as 3b. There are no examples of tautomer 1 in the crystal structure database; it may only exist in the gas phase. So the equilibrium 2 ⇌ 3 is another genuine example of tautomeric polymorphism, with the keto form favoured by more polar solvents, as was noted for the previous example.

With this last article,[3] comprehensive calculations at a good level were reported, including modelling the periodic cell using the Crystal program and including corrections such as BSSE (basis set superposition error) and dispersion terms. I was hopeful that this might lead me to something as simple as the computed dipole moments of the (isolated) species (as I reported above for the previous system), but these were not mentioned in the text of the article. Unfortunately, the supporting information also had no details of any such calculations, which left me frustrated again at how difficult it can be in (it has to be said) the vast majority of articles which report calculations to get details of such calculations. 

Tautomeric polymorphism remains a very rare phenomenon. SciFinder for example only has 19 references citing it (2 of which are to conference talks). Perhaps the most intriguing[4] claims that 2-thiobarbituric acid has the richest collection of tautomeric polymorphs with five. Since no calculations are reported there, I might try these out and report back here.

Postscript:  Here is some analysis of 2-thiobarbituric.

  1. THBARB (DOI 10.5517/cctbxcd10.5517/cctbxfg  and 10.5517/cctbxgh) are three polymorphs of  the keto tautomer, the isolated molecule having a small calculated dipole moment (DOI: 10.14469/hpc/2632).
  2. PABNAJ (DOI: 10.5517/cctbxbc) is a polymorph in the enol form, with a much larger calculated dipole moment (DOI: 10.14469/hpc/2633)
  3. PABNIR (DOI: 10.5517/cctbxdf) is a mixed polymorph with one enol paired with one keto form. 

The relative free-energies of the isolated molecules are 0.0 (keto) and 9.0 (enol). The keto-enol pair is 0.4 kcal/mol more stable than the isolated components. This again shows the effect that crystal packing can have on the relative energies and also shows that a  simple inspection of the dipole moment may cast light on the polymorphism.

 

References

  1. P.M. Bhatt, and G.R. Desiraju, "Tautomeric polymorphism in omeprazole", Chemical Communications, pp. 2057, 2007. https://doi.org/10.1039/b700506g
  2. Y. Akama, M. Shiro, T. Ueda, and M. Kajitani, "Keto and Enol Tautomers of 4-Benzoyl-3-methyl-1-phenyl-5(2H)-pyrazolone", Acta Crystallographica Section C Crystal Structure Communications, vol. 51, pp. 1310-1314, 1995. https://doi.org/10.1107/s0108270194007389
  3. S. Long, M. Zhang, P. Zhou, F. Yu, S. Parkin, and T. Li, "Tautomeric Polymorphism of 4-Hydroxynicotinic Acid", Crystal Growth & Design, vol. 16, pp. 2573-2580, 2016. https://doi.org/10.1021/acs.cgd.5b01639
  4. M. Chierotti, L. Ferrero, N. Garino, R. Gobetto, L. Pellegrino, D. Braga, F. Grepioni, and L. Maini, "The Richest Collection of Tautomeric Polymorphs: The Case of 2‐Thiobarbituric Acid", Chemistry – A European Journal, vol. 16, pp. 4347-4358, 2010. https://doi.org/10.1002/chem.200902485

The conformation of enols: revealed and explained.

Thursday, April 6th, 2017

Enols are simple compounds with an OH group as a substituent on a C=C double bond and with a very distinct conformational preference for the OH group. Here I take a look at this preference as revealed by crystal structures, with the theoretical explanation.

First, a search of the Cambridge structure database (CDS), using the search query shown below (DOI: 10.14469/hpc/2429)


The first search (no errors, no disorder, R < 0.05) is unconstrained in the sense that the HO group is free to hydrogen bond itself. The syn conformer has the torsion of 0° and it has a distinct preponderance over the anti isomer (180°). There is the first hint that the most probable C=C distance for the syn isomer may be longer than that for the anti, but this is not yet entirely convincing.
To try to make it so, a constrained search is now performed, in which only structures where the HO group has no contact (hydrogen bonding) interaction are included. This is achieved using a “Boolean” search;

The number of hits approximately halves, but the proportion of syn examples increases considerably. There is an interesting double “hot-spot” distribution, which amplifies the lengthening of the C=C bond compared to the anti orientation.

The next constraint added is that the data collection must be <100K (to reduce thermal noise) which reduces the hits considerably but now shows the lengthening of the C=C bond for the syn isomer very clearly.

A final plot is of the C=C length vs the C-O length (no temperature, but HO interaction constraint). If there were no correlation, the distribution would be ~circular. In fact it clearly shows that as the C=C bond lengthens, the C-O bond contracts.

Now for some calculations (ωB97XD/Def2-TZVPP, DOI: 10.14469/hpc/2429) which reveal the following:

  1. The free energy of the syn isomer is 1.2 kcal/mol lower than that of the syn. The effect is small, and hence easily masked by other interactions such as hydrogen bonding to the OH group. Hence the reason why removing such interactions from the search above increased the syn population compared to anti.
  2. The syn C=C bond length (1.325Å) is longer than the anti (1.322Å). 
  3. The syn isomer has one unique σO-Lp*C-C NBO orbital interaction (below) with a value of E(2) 7.7 kcal/mol, which is absent in the anti form. As it happens, a πO*C=C interaction is present in both forms but is also stronger in the syn isomer (E(2)= 46.8 vs 44.2 kcal/mol).
    unoccupied NBO, σ*C-C
    Occupied NBO, σO-Lp
  4. The overlap of the filled σO-Lp with the empty σ*C-C orbital is shown below (blue overlaps with purple, red overlaps with orange).

    To view the overlap in rotatable 3D, click on any of the colour diagrams above.

It is nice to see how experiment (crystal structures) and theory (the calculation of geometries and orbital interactions) can quickly and simply be reconciled. Both these searches and the calculations can be done in just one day of “laboratory time” and I think it would make for an interesting undergraduate chemistry lab experiment.


This visualisation uses Java. Increasingly this browser plugin is becoming more onerous to activate (because of increased security) and some browsers do not support it at all. The macOS Safari browser is one that still does, but you do have to allow it via the security permissions.