Posts Tagged ‘search results’

Deviations from tetrahedral four-coordinate carbon: a statistical exploration.

Sunday, September 6th, 2015

An article entitled “Four Decades of the Chemistry of Planar Hypercoordinate Compounds[1] was recently reviewed by Steve Bacharach on his blog, where you can also see comments. Given the recent crystallographic themes here, I thought I might try a search of the CSD (Cambridge structure database) to see whether anything interesting might emerge for tetracoordinate carbon.

The search definition is shown below using a  simple carbon with four ligands, the ligands themselves also being tetracoordinate carbon. The search is restricted to data collected below temperatures of 140K, as well as R-factor <5%, no errors and no disorder. Cyclic species are allowed and a statistically reasonable 2773 hits emerged from the search.

Scheme

Recollect that the idealised angle subtended at the centre is 109.47°. I show below three separate heat plots of the search results. Why three? The way the search software (Conquest) works is that one could define four C-C distances and six angles, and then plot any combination of one distance and one angle. I show just three combinations here, but could have included many more.

There appear to be four distinct clusters of values for this angle that emerge from the three plots shown below (the “bin size” is 100, and the frequency colour code indicates how many hits there are in each bin).

  1. The hotspot is unsurprisingly ~109° with a corresponding C-C distance of ~1.54Å.
  2. There may be two clusters at angles of ~60° (cyclopropane), with C-C values ranging from ~1.47 to ~1.55Å.
  3. A collection at ~90° (mostly cyclobutane?), with C-C values up to 1.6Å.
  4. A collection at ~140° (again small rings), now with much shorter C-C values of ~1.46Å. This reminds of the approximation that the hybridisation in e.g. cyclopropane is a combination of sp5 and sp3.

Scheme

Scheme

Scheme

Ideally, what one might want to plot would be sums of four angles; for a pure tetrahedral carbon the sum would always be 438° (4*109.47°) but for a pure planar carbon it could be as low as 360° (4*90°). One could then see how closely the distribution approaches to the latter and hence reveal whether there are any true planar tetracoordinate carbon species known. Although the Conquest software cannot analyse in such terms, a Python-based API has recently been released that should allow this to be done, although I should state that this requires a commercial license and it is not open access code. If we manage to get it working, I will report!


As a teaser I also include a plot of six-coordinate carbon, in which the ligands can be any non-metal. Note the clusters at angles of 60, ~112 and ~120-130°. It is worth pointing out that the definition of the connection between a carbon and a ligand as a “bond” becomes increasingly arbitrary as the coordination becomes “hyper”. Because crystallography does not measure electron densities in “bonds”, we know nothing of its topology in this region. It is therefore quite possible that the appearance of the heat plot below might be related just as much to whatever convention is being used in creating the entry in the CSD as it would be to a quantum analysis of the bonding.

Scheme

References

  1. L. Yang, E. Ganz, Z. Chen, Z. Wang, and P.V.R. Schleyer, "Four Decades of the Chemistry of Planar Hypercoordinate Compounds", Angewandte Chemie International Edition, vol. 54, pp. 9468-9501, 2015. https://doi.org/10.1002/anie.201410407

The conformational preference of s-cis amides.

Sunday, February 10th, 2013

Amides with an H-N group are a component of the peptide linkage (O=C-NH). Here I ask what the conformation (it could also be called a configuration) about the C-N bond is. A search of the following type can be defined:

cis-amide

The dihedral shown is for H-N-C=O (but this is equivalent to the C-C-N-C dihedral, which is also often called the dihedral angle associated with the peptide group). I have also added a distance, from a C-H to the carbonyl oxygen. Other search constraints include T ≤ 175K, R < 0.05, no disorder, no errors, that neither N-C bonds are part of a ring and that the two carbons marked T4 both have four connected bonds. The search results in 619 hits (January 2013 version of the CCDC database), and these are displayed below.

cis-amide-search-heat

The horizontal axis reveals the highest concentration (red) at ~2.4Å due to a syn-co-planar alignment of the C-H bond with the plane of the C=O bond in the s-cis conformer (the significantly smaller hot-spot at ~3.9A may be due to an anti-co-planar alignment of this C-H bond).

s-cis-amide

The vertical axis shows a clear preference for a dihedral of 179° (in fact no hits with a dihedral of less than 14o° were found) and this can only arise from the s-cis conformation in which the H-N bond is oriented antiperiplanar to the axis of the C=O bond. This preference can be rationalised by filled/empty NBO-orbital interactions, which include:

  1. Antiperiplanar interaction between the N-H as donor and the C=O as a σ-acceptor (E(2) = 4.1 kcal/mol)
  2. Antiperiplanar interaction between the N-H as acceptor and C-H as donor (E(2) = 4.7 kcal/mol)
Click for 3D

H-N/C=O. Click for 3D

 

Click for 3D.

Click for 3D.

This latter overlap conspires to bring the C-H hydrogen close to the oxygen (~2.35Å, DIST1 in the diagram above). So one might be entitled to ask: is this a hydrogen bond? There are (at least) two ways of testing this.

  1. The NBO E(2) interaction energy between the oxygen in-plane lone pair and the H-C as acceptor is 0.8 kcal/mol. For hydrogen bonds, such E(2) energies more or less resemble the actual H-bond strengths, i.e. a strong H-bond has an E(2) energy of ~ 8 kcal/mol; and a medium O…H-C hydrogen bond weighs in at around 3 kcal/mol.  So this one is very weak. This is due to poor overlap resulting from the small ring size (5).
  2. The NCI (non-covalent-interaction) surface does reveal a feature in the CH…O region, but the colour coding (which indicates how attractive/repulsive this is) is both pale blue (attractive) and yellow (repulsive). Again this is only consistent with a very weak overall H-bond.
NCI surface. Click for 3D.

NCI surface. Click for 3D.

I end by reminding that the s-cis H-N-C=O conformation is a very common feature in peptides (the CCDC database comprises mostly small molecules, not larger peptides and proteins) arising from really quite subtle orbital interactions.