Computational evaluation of unsaturated carbonitriles as neutral receptor model for beryllium(II) recognition

Design of neutral receptor molecules (ionophores) for beryllium(II) using unsaturated carbonitrile models has been carried out via density functional theory, G3, and G4 calculations. The first part of this work focuses on gas phase binding energies between beryllium(II) and 2-cyano butadiene (2-CN BD), 3-cyano propene (3-CN P), and simpler models with two separate fragments; acrylonitrile and ethylene. Interactions between beryllium(II) and cyano nitrogen and terminal olefin in the models have been examined in terms of geometrical changes, distribution of charge over the entire π-system, and rehybridization of vinyl carbon orbitals. NMR shieldings and vibrational frequencies probed charge centers and strength of interactions. The six-membered cyclic complexes have planar structures with the rehybridized carbon slightly out of plane (16° in 2-CN BD). G3 results show that in 2-CN BD complex participation of vinyl carbon further stabilizes the cyclic adduct by 16.3 kcal mol(-1), whereas, in simpler models, interaction between beryllium(II) and acetonitrile is favorable by 46.4 kcal mol(-1) compared with that of ethylene. The terminal vinyl carbon in 2-CN BD rehybridizes to sp (3) with an increase of 7 % of s character to allow interaction with beryllium(II). G4 calculations show that the Be(II) and 2-CN BD complex is more strongly bound than those with Mg(II) and Ca(II) by 98.5 and 139.2 kcal mol(-1) (-1), respectively. QST2 method shows that the cyclic and acyclic forms of Be(II)-2-CN BD complexes are separated by 12.3 kcal mol(-1) barrier height. Overlap population analysis reveals that Ca(II) can be discriminated based on its tendency to form ionic interaction with the receptor models.