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 models).[4]
Additionally, the stereoelectronic principles that underlie nucleophiles adopting a proscribed range of Bürgi–Dunitz angles may contribute to the conformational stability of proteins[5] [6] and are invoked to explain the stability of particular conformations of molecules in one hypothesis of a chemical origin of life.[7]
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.[8]
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-formadehyde 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]
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 (H2C=O), gave a BD angle value of 107°.[2]
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.
In the structure of -methadone (above, left), note the tertiary amine projecting to the lower right, and the carbonyl (CO) 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 CO 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.
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.[10] [11] [12])
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).
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).[13]
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.
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.[14] [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]