David A. Bardwell

As part of a programme of work to extend central-group coverage in the Cambridge Crystallographic Data Centre's (CCDC) IsoStar knowledge base of intermolecular interactions, we have studied the hydrogen-bonding abilities of squaric acid (H2SQ) and its mono- and dianions (HSQ(-) and SQ(2-)) using the Cambridge Structural Database (CSD) along with dispersion-corrected density functional theory (DFT-D) calculations for a range of hydrogen-bonded dimers. The -OH and -C=O groups of H2SQ, HSQ(-) and SQ(2-) are potent donors and acceptors, as indicated by their hydrogen-bond geometries in available crystal structures in the CSD, and by the attractive energies calculated for their dimers with acetone and methanol, which were used as model acceptors and donors. The two anions have sufficient examples in the CSD for their addition as new central groups in IsoStar. It is also shown that charge- and resonance-assisted hydrogen bonds involving H2SQ and HSQ(-) are similar in strength to those made by carboxylate COO(-) acceptors, while hydrogen bonds made by the dianion SQ(2-) are somewhat stronger. The study reinforces the value of squaric acid and its anions as cocrystal formers and their actual and potential importance as isosteric replacements for carboxylic acid and carboxylate functions.

Bioisosterism involving replacement of a carboxylic acid substituent by 1H-tetrazole, yielding deprotonated carboxylate and tetrazolate under physiological conditions, is a well-known synthetic strategy in medicinal chemistry. To improve our overall understanding of bioisosterism, we have used this example to study the geometrical and energetic aspects of the functional group replacement. Specifically, we use crystal structure informatics and high-level ab initio calculations to study the hydrogen bond (H-bond) energy landscapes of the protonated and deprotonated bioisosteric pairs. Each pair exhibits very similar H-bond environments in crystal structures retrieved from the CSD, and the attractive energies of these H-bonds are also very similar. However, by comparison with −COOH and −COO^–, the H-bond environments around 1H-tetrazole and tetrazolate substituents extend further, by about 1.2 Å, from the core of the connected molecule. Analysis of pairs of PDB structures containing ligands which differ only in having a tetrazole or a carboxyl substituent and which are bound to the same protein indicates that the protein binding site must flex sufficiently to form strong H-bonds to either substituent. A survey of DrugBank shows a rather small number of tetrazole-containing drugs in the ‘approved’ and ‘experimental’ drug sections of that database.

Following on from the success of the previous crystal structure prediction blind tests (CSP1999, CSP2001, CSP2004 and CSP2007), a fifth such collaborative project (CSP2010) was organized at the Cambridge Crystallographic Data Centre. A range of methodologies was used by the participating groups in order to evaluate the ability of the current computational methods to predict the crystal structures of the six organic molecules chosen as targets for this blind test. The first four targets, two rigid molecules, one semi-flexible molecule and a 1:1 salt, matched the criteria for the targets from CSP2007, while the last two targets belonged to two new challenging categories – a larger, much more flexible molecule and a hydrate with more than one polymorph. Each group submitted three predictions for each target it attempted. There was at least one successful prediction for each target, and two groups were able to successfully predict the structure of the large flexible molecule as their first place submission. The results show that while not as many groups successfully predicted the structures of the three smallest molecules as in CSP2007, there is now evidence that methodologies such as dispersion-corrected density functional theory (DFT-D) are able to reliably do so. The results also highlight the many challenges posed by more complex systems and show that there are still issues to be overcome.

The magnetic properties of two series of dinuclear complexes, and one trinuclear complex, have been examined as a function of the bridging pathway between the metal centers. The first series of dinuclear complexes is [{Mo(V)(O)(Tp)Cl}(2)(&mgr;-OO)], where "OO" is [1,4-O(C(6)H(4))(n)O](2)(-) (n = 1, 1; n = 2, 3), [4,4'-O(C(6)H(3)-2-Me)(2)O](2)(-) (4), or [1,3-OC(6)H(4)O](2)(-) (2) [Tp = tris(3,5-dimethylpyrazolyl)hydroborate]. The second series of dinuclear complexes is [{Mo(I)(NO)(Tp)Cl}(2)(&mgr;-NN)], where "NN" is 4,4'-bipyridyl (5), 3,3'-dimethyl-4,4'-bipyridine (6), 3,8-phenanthroline (7), or 2,7-diazapyrene (8). The trinuclear complex is [{Mo(V)(O)(Tp)Cl}(3)(1,3,5-C(6)H(3)O(3))] (9), whose crystal structure was determined [9.5CH(2)Cl(2): C(56)H(81)B(3)Cl(13)Mo(3)N(18)O(6); monoclinic, P2(1)/n; a = 13.443, b = 41.46(2), c = 14.314(6) Å; beta = 93.21(3) degrees; V = 7995(5) Å(3); Z = 4; R(1) = 0.106]. In these complexes, the sign and magnitude of the exchange coupling constant J is clearly related to both the topology and the conformation of the bridging ligand [where J is derived from H = -JS(1)().S(2)() for 1-8 and H = -J(S(1)().S(2)() + S(2)().S(3)() + S(1)().S(3)()) for 9]. The values are as follows: 1, -80 cm(-)(1); 2, +9.8 cm(-)(1); 3, -13.2 cm(-)(1); 4, -2.8 cm(-)(1); 5, -33 cm(-)(1); 6, -3.5 cm(-)(1); 7, -35.6 cm(-)(1); 8, -35.0 cm(-)(1); 9, +14.4 cm(-)(1). In particular the following holds: (1) J is negative (antiferromagnetic exchange) across the para-substituted bridges ligands of 1 and 3-8 but positive (ferromagnetic exchange) across the meta-substituted bridging ligands of 2 and 9. (2) J decreases in magnitude dramatically as the bridging ligand conformation changes from planar to twisted (compare 3 and 4, or 6 and 8). These observations are consistent with a spin-polarization mechanism for the exchange interaction, propagated across the pi-system of the bridging ligand by via overlap of bridging ligand p(pi) orbitals with the d(pi) magnetic orbitals of the metals. The EPR spectrum of 9 is characteristic of a quartet species and shows weak Deltam(s) = 2 and Deltam(s) = 3 transitions at one-half and one-third, respectively, of the field strength of the principal Deltam(s) = 1 component.

Lanthanide complexes of two polydentate N-donor ligands containing a mixture of pyridyl and pyrazolyl donors have been prepared. Dihydrobis[3-(2-pyridyl)pyrazolyl]borate (L1)- is a tetradentate ligand with two bidentate chelating pyridyl/pyrazolyl arms linked by an apical BH2 group; 2,6-bis(1H-pyrazol-3-yl)pyridin