Bürgi–Dunitz angle

Bürgi–Dunitz and Flippin–Lodge nucleophile approach trajectory angles, $\alpha _{BD}$ and $\alpha _{FL}$ , analogized to altitude- and azimuth-types of parameters in the celestial (horizontal) coordinate system. Both in this celestial coordinate system, and in the system to describe a nucleophile approaching a planar electrophile, the problem is to uniquely describe the location of a point off of a plane, relative to specific point on the plane. Hence, in both cases, the problem can be addressed using two angles, an altitude-type angle and an azimuth-type angle. Note, while elevation in celestial applications is most easily measured as the specific altitude shown, the elevation of the nucleophile, $\alpha _{BD}$ , is most easily measured as the supplementary angle, Nu-C-O, hence its values are most often >90° (see text and next image).

The Bürgi–Dunitz angle (BD angle) is one of two angles that fully define the geometry of "attack" (approach via collision) of a nucleophile on a trigonal unsaturated center in a molecule, originally the carbonyl center in an organic ketone, but now extending to aldehyde, ester, and amide carbonyls, and to alkenes (olefins) as well.^[1]^[2]^[3] The angle was named after crystallographers Hans-Beat Bürgi and Jack D. Dunitz, its first senior investigators.

Practically speaking, the Bürgi–Dunitz and Flippin–Lodge angles were central to the development of understanding of chiral chemical synthesis, and specifically of the phenomenon of asymmetric induction during nucleophilic attack at hindered carbonyl centers (see the Cram–Felkin–Anh and Nguyen^{[clarification needed]} models).^[4]^[5]

Additionally, the stereoelectronic principles that underlie nucleophiles adopting a proscribed range of Bürgi–Dunitz angles may contribute to the conformational stability of proteins^[6]^[7] and are invoked to explain the stability of particular conformations of molecules in one hypothesis of a chemical origin of life.^[8]

Definition

In the addition of a nucleophile (Nu) attack to a carbonyl, the BD angle is defined as the Nu-C-O bond angle. The BD angle adopted during an approach by a nucleophile to a trigonal unsaturated electrophile depends primarily on the molecular orbital (MO) shapes and occupancies of the unsaturated carbon center (e.g., carbonyl center), and only secondarily on the molecular orbitals of the nucleophile.^[1]

Of the two angles which define the geometry of nucleophilic "attack", the second describes the "offset" of the nucleophile's approach toward one of the two substituents attached to the carbonyl carbon or other electrophilic center, and was named the Flippin–Lodge angle (FL angle) by Clayton Heathcock after his contributing collaborators Lee A. Flippin and Eric P. Lodge.^[4]

These angles are generally construed to mean the angle measured or calculated for a given system, and not the historically observed value range for the original Bürgi–Dunitz aminoketones, or an idealized value computed for a particular system (such as hydride addition to formaldehyde, image at left). That is, the BD and FL angles of the hydride-formaldehyde system produce a given pair of values, while the angles observed for other systems may vary relative to this simplest of chemical systems.^[1]^[3]^[9]

Measurement

The original Bürgi-Dunitz measurements were of a series of intramolecular amine-ketone carbonyl interactions, in crystals of compounds bearing both functionalities—e.g., methadone and protopine. These gave a narrow range of BD angle values (105 ± 5°); corresponding computations—molecular orbital calculations of the SCF-LCAO-type—describing the approach of the s-orbital of a hydride anion (H⁻) to the pi-system of the simplest aldehyde, formaldehyde (H₂C=O), gave a BD angle value of 107°.^[2]^{[non-primary source needed]}

Hence, Bürgi, Dunitz, and thereafter many others noted that the crystallographic measurements of the aminoketones and the computational estimate for the simplest nucleophile-electrophile system were quite close to a theoretical ideal, the tetrahedral angle (internal angles of a tetrahedron, 109.5°), and so consistent with a geometry understood to be important to developing transition states in nucleophilic attacks at trigonal centers.^{[citation needed]}

Structure of L-methadone, a key molecule whose crystals were studied in the development of the BD angle concept.
Structure of protopine, an alkaloid, whose crystals were also studied in the development of the BD angle concept.
The amine-to-carbonyl n→π* interaction in protopine with an unusually short N···C distance of 2.555 Å and a Bürgi–Dunitz angle of 102°.^[10]

In the structure of L-methadone (above, left), note the tertiary amine projecting to the lower right, and the carbonyl (C=O) group at the center, which engage in an intramolecular interaction in the crystal structure (after rotation around the single bonds connecting them, during the crystallization process).

Similarly, in the structure of protopine (above, center), note the tertiary amine at the center of the molecule, part of a ten-membered ring, and the C=O group opposite it on the ring; these engage in an intramolecular interaction allowed by changes in the torsion angles of the atoms of the ring.

Theory

Hydride addition to formaldehyde

Estimated Bürgi–Dunitz (

\alpha _{BD}

) angle for this simplest H^(–) → H₂C=O nucleo-philic addition.

HOMO-LUMO interaction underlying

\alpha _{BD}

angle in simple systems.

(Left) Shown is nucleophilic attack by the charged nucleophile (Nu), hydride anion, on the unsaturated trigonal center of the aldehyde electrophile, formaldehyde (R,R'=H). The value computed as optimal for this system, 107°, is indicated, and is representative of the obtuse values observed in most experimental chemical systems.
(Right) Cartoon of the approach of a p-type highest occupied molecular orbital (HOMO) of a nucleophile such as chloride ion (green sphere), and the lowest unoccupied molecular orbital (LUMO) of the trigonal center of an electrophilic carbonyl of formaldehyde (black sphere carbon, red sphere oxygen, white spheres hydrogen). View is nearly side on, and the developing out-of-plane distortion of the carbonyl carbon atom is omitted for simplicity.

The convergence of observed BD angles can be viewed as arising from the need to maximize overlap between the highest occupied molecular orbital (HOMO) of the nucleophile, and the lowest unoccupied molecular orbital (LUMO) of the unsaturated, trigonal center of the electrophile.^[1] (See, in comparison, the related inorganic chemistry concept of the angular overlap model.^[11]^[12]^[13]^{[page needed]})

In the case of addition to a carbonyl, the HOMO is often a p-type orbital (e.g., on an amine nitrogen or halide anion), and the LUMO is generally understood to be the antibonding π* molecular orbital perpendicular to the plane containing the ketone C=O bond and its substituents (see figure at right above). The BD angle observed for nucleophilic attack is believed to approach the angle that would produce optimal overlap between HOMO and LUMO (based on the principle of the lowering of resulting new molecular orbital energies after such mixing of orbitals of similar energy and symmetry from the participating reactants). At the same time, the nucleophile avoids overlap with other orbitals of the electrophilic group that are unfavorable for bond formation (not apparent in image at right, above, because of the simplicity of the R=R'=H in formaldehyde).^{[citation needed]}

Complications

Electrostatic and Van der Waals interactions

To understand cases of real chemical reactions, the HOMO-LUMO-centered view is modified by understanding of further complex, electrophile-specific repulsive and attractive electrostatic and Van der Waals interactions that alter the altitudinal BD angle, and bias the azimuthal Flippin-Lodge angle toward one substituent or the other (see graphic above).^[14]^{[non-primary source needed]}

Linear and rotational dynamics

BD angle theory was developed based on "frozen" interactions in crystals where the impacts of dynamics at play in the system (e.g., easily changed torsional angles) may be negligible. However, much of the chemistry of general interest and use takes place via collisions of molecules tumbling in solution; accordingly, dynamics are taken into account in such cases.

Constrained environments in enzymes and nanomaterials

Moreover, in constrained reaction environments such as in enzyme and nanomaterial binding sites, early evidence suggests that BD angles for reactivity can be quite distinct, since reactivity concepts assuming orbital overlaps during random collision are not directly applicable.^[15]^[9]

For instance, the BD value determined for enzymatic cleavage of an amide by a serine protease (subtilisin) was 88°, quite distinct from the hydride-formaldehyde value of 107°; moreover, compilation of literature crystallographic BD angle values for the same reaction mediated by different protein catalysts clustered at 89 ± 7° (i.e., only slightly offset from directly above or below the carbonyl carbon). At the same time, the subtilisin FL value was 8°, and FL angle values from the careful compilation clustered at 4 ± 6° (i.e., only slightly offset from directly behind the carbonyl; see the Flippin–Lodge angle article).^[9]^{[non-primary source needed]}

References

^ ^a ^b ^c ^d Fleming, I. (2010) Molecular Orbitals and Organic Chemical Reactions: Reference Edition, John Wiley & Sons, pp. 214–215.
^ ^a ^b Bürgi, H.-B.; Dunitz, J. D.; Lehn. J.-M.; Wipff, G. (1974). "Stereochemistry of reaction paths at carbonyl centres". Tetrahedron. 30 (12): 1563–1572. doi:10.1016/S0040-4020(01)90678-7.
^ ^a ^b Cieplak, A.S. (2008) Organic addition and elimination reactions: Transformation paths of carbonyl derivatives In Structure Correlation, Vol. 1 (H.-B. Bürgi & J. D. Dunitz, eds.), New York:John Wiley & Sons, pp. 205–302, esp. 216-218. [doi:10.1002/9783527616091.ch06; ISBN 9783527616091 ]
^ ^a ^b Heathcock, C.H. (1990) Understanding and controlling diastereofacial selectivity in carbon-carbon bond-forming reactions, Aldrichimica Acta 23(4):94-111, esp. p. 101, see [1], accessed 9 June 2014.
^ Gawley, R.E. & Aube, J. 1996, Principles of Asymmetric Synthesis (Tetrahedron Organic Chemistry Series, Vo. 14), pp. 121-130, esp. pp. 127f.
^ Bartlett, G.J.; Choudhary, A.; Raines, R.T.; Woolfson, D.N. (2010). "n→π* interactions in proteins". Nat. Chem. Biol. 6 (8): 615–620. doi:10.1038/nchembio.406. PMC 2921280. PMID 20622857.
^ Fufezan, C. (2010). "The role of Buergi‐Dunitz interactions in the structural stability of proteins". Proteins. 78 (13): 2831–2838. doi:10.1002/prot.22800. PMID 20635415. S2CID 41838636.
^ Choudhary, A.; Kamer, K.J.; Powner, M.W.; Sutherland, J.D.; Raines, R.T. (2010). "A stereoelectronic effect in prebiotic nucleotide synthesis". ACS Chem. Biol. 5 (7): 655–657. doi:10.1021/cb100093g. PMC 2912435. PMID 20499895.
^ ^a ^b ^c Radisky, E.S. & Koshland, D.E. (2002), A clogged gutter mechanism for protease inhibitors, Proc. Natl. Acad. Sci. U.S.A., 99(16):10316-10321, see [2], accessed 28 November 2014.
^ Hall, S. R.; Ahmed, F. R. (1968). "The crystal structure of protopine, C₂₀H₁₉O₅N". Acta Crystallogr. B. 24 (3): 337–346. doi:10.1107/S0567740868002347.
^ Hoggard, P.E. (2004) Angular overlap model parameters, Struct. Bond. 106, 37.
^ Burdett, J.K. (1978) A new look at structure and bonding in transition metal complexes, Adv. Inorg. Chem. 21, 113.
^ Purcell, K.F. & Kotz, J.C. (1979) Inorganic Chemistry, Philadelphia, PA:Saunders Company.^{[page needed]}
^ Lodge, E.P. & Heathcock, C.H. (1987) Steric effects, as well as sigma*-orbital energies, are important in diastereoface differentiation in additions to chiral aldehydes, J. Am. Chem. Soc., 109:3353-3361.
^ See for instance, Light, S.H.; Minasov, G.; Duban, M.-E. & Anderson, W.F. (2014) Adherence to Bürgi-Dunitz stereochemical principles requires significant structural rearrangements in Schiff-base formation: Insights from transaldolase complexes, Acta Crystallogr. D 70(Pt 2):544-52, DOI: 10.1107/S1399004713030666, see [3], accessed 10 June 2014.