Dimers of boroglycine and methylamine boronic acid: a computational comparison of the relative importance of dative versus hydrogen bonding

Physical binding energies using the electron localization function in 4-hydroxyphenylboronic acid co-crystals with aza donors

Binding energies are traditionally simulated using cluster models by computation of each synthon for each individual co-crystal former. However, our investigation of the binding strengths using the electron localization function (ELF) reveals that these can be determined directly from the crystal supercell computations. We propose a new modeling protocol for the computation of physical binding energies directly from bulk simulations using ELF analysis. In this work, we establish a correlation between ELF values and binding energies calculated for co-crystals of 4-hydroxyphenylboronic acid (4HPBA) with four different aza donors using density functional theory with varying descriptions of dispersion. Boronic acids are gaining significant interest in the field of crystal engineering, but theoretical studies on their use in materials are still very limited. Here, we present a systematic investigation of the non-covalent interactions in experimentally realized co-crystals. Prior diffraction studies on these complexes have shown the competitive nature between the boronic acid functional group and the para-substituted phenolic group forming heteromeric interactions with aza donors. We determine the stability of the co-crystals by simulating their lattice energies, and the different dispersion descriptions show similar trends in lattice energies and lattice parameters. Our study bolsters the experimental observation of the boronic acid group as a competitive co-crystal former in addition to the well-studied phenolic group. Further research on correlating ELF values for physical binding could potentially transform this approach to a viable alternative for the computation of binding energies.

Heats of Formation for the Boronic Acids R-B(OH)2 and Boroxines R3B3O3 (R=H, Li, HBe, H2B, H3C, H2N, HO, F, and Cl) Calculated at the G2, G3, and G4 Levels of Theory

Boronic acids (R-B(OH)₂) and their boroxine (R₃B₃O₃) dehydration products have emerged as important classes of compounds with a multitude of diverse applications. However, the available heats of formation for these compounds are not always as accurate as would be required for further use. In this study the heats of formation at 298.15 K of R-B(OH)₂ and R₃B₃O₃ (R = H, Li, HBe, H₂B, H₃C, H₂N, HO, F, and Cl) have been calculated at the G2, G3[G3B3], and G4 levels of theory and used to determine the enthalpy changes for the dehydration reactions: 3 R-B(OH)₂ → R₃B₃O₃ + 3 H₂O; comparisons are made with other rigorous levels of theory, e.g. CBS-Q[CBS-QB3] and W1U, as well as with experimental values wherever possible. Enthalpy changes for the dehydration reactions have also been calculated using second-order Møller-Plesset perturbation theory (MP2) with the Dunning-Woon correlation-consistent aug-cc-pVDZ and aug-cc-pVTZ basis sets, and B3LYP density functional theory with the 6-311++G(2df,2pd) basis set. With the exception of H₂N-B(OH)₂, the dehydration reactions are consistently predicted to be exothermic. Our results provide a cautionary note for the use of the B3LYP functional in the calculation of structures and energies of boronic acids and boroxines. Where comparisons could be made, the G4 and W1U predictions for the heats of formation of these boron compounds differ significantly.

Computational investigation of the oxidative deboronation of boroglycine, H2N-CH2-B(OH)2, Using H2O and H2O2

We report results from a computational investigation of the oxidative deboronation of boroglycine, H2N-CH2-B(OH)2, using H2O and H2O2 as the reactive oxygen species (ROS) to yield aminomethanol, H2N-CH2-OH; these results complement our study on the protodeboronation of boroglycine to produce methylamine, H2N-CH3 (Larkin et al. J. Phys. Chem. A 2007, 111, 6489-6500). Second-order Møller-Plesset (MP2) perturbation theory with Dunning-Woon correlation-consistent (cc) basis sets were used for the calculations with comparisons made to results from density functional theory (DFT) at the PBE1PBE/6-311++G(d,p)(cc-pVDZ) levels. The effects of a bulk aqueous environment were also incorporated into the calculations employing PCM and CPCM methodology. Using H2O as the ROS, the reaction H2O + H2N-CH2-B(OH)2 --> H2N-CH2-OH + H-B(OH)2 was calculated to be endothermic; the value of DeltaH(298)(0) was +12.0 kcal/mol at the MP2(FC)/cc-pVTZ computational level in vacuo and +13.7 kcal/mol in PCM aqueous media; the corresponding value for the activation barrier, DeltaH(double dagger), was +94.3 kcal/mol relative to the separated reactants in vacuo and +89.9 kcal/mol in PCM aqueous media. In contrast, the reaction H2O2 + H2N-CH2-B(OH)2 --> H2N-CH2-OH + B(OH)3 was calculated to be highly exothermic with an DeltaH(298)(0) value of -100.9 kcal/mol at the MP2(FC)/cc-pVTZ computational level in vacuo and -99.6 kcal/mol in CPCM aqueous media; the highest-energy transition state for the multistep process associated with this reaction involved the rearrangement of H2N-CH2-B(OH)(OOH) to H2N-CH2-O-B(OH)2 with a DeltaH(double dagger) value of +23.2 kcal/mol in vacuo relative to the separated reactants. These computational results for boroglycine are in accord with the experimental observations for the deboronation of the FDA approved anticancer drug bortezomib (Velcade, PS-341), where it was found to be the principle deactivation pathway (Labutti et al. Chem. Res. Toxicol. 2006, 19, 539-546).

Conformational analysis of hydroxymatairesinol in aqueous solution with molecular dynamics

Molecular dynamics simulations were performed on the naturally occuring lignan hydroxymatairesinol (HMR) using the GROMACS software. The aim of this study was to explore the conformational behavior of HMR in aqueous solution adopting the TIP4P model. The topology of HMR was constructed by hand and HMR was modeled with the OPLS-AA force field implemented in GROMACS. The five torsional angles in HMR were properly analyzed during the simulations. Correlations through certain patterns were observed between the angles. The determining property for the conformation preferred in aqueous solution was found to be the dipole moment and not the lowest energy in gas phase. The solvation effects on HMR was also studied by quantum chemical calculations applying the COnductorlike Screening MOdel (COSMO), the results of which were compared with results from a previous study using the Polarized Continuum Model (PCM). In the present work, COSMO was found to give more credible relative energies than PCM.

A computational characterization of boron-oxygen multiple bonding in HN=CH-CH=CH-NH-BO

Structures, relative energies, and bonding characteristics for various conformers of 3-imino-N-(oxoboryl)prop-1-en-1-amine, HN=CH-CH=CH-NH-BO, and the corresponding borocycle (-HN=CH-CH=CH-NH-B-)O are discussed using results from second-order Møller-Plesset (MP2) perturbation theory with the Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, and cc-pVTZ basis sets. These MP2 results are compared to those from computationally efficient density functional theory (DFT) calculations using the LDA, PBE, TPSS, BLYP, B3LYP, BVP86, OLYP, O3LYP, and PBE1PBE functionals in conjunction with the economical Pople-type 6-311++G(d,p) basis set to evaluate the suitability of these DFT/6-311++G(d,p) levels for use with larger boron-containing systems. The effects of an aqueous environment were incorporated into the calculations using COSMO methodology. The calculated boron-oxygen bond lengths, orbital compositions, and bond orders in all the (acyclic) HN=CH-CH=CH-NH-BO conformers were consistent with the presence of a boron-oxygen triple bond, similar to that found in H-BO and H2N-BO. The (-HN=CH-CH=CH-NH-B-)O borocycle is predicted to be planar (C2v symmetry), and it is approximately 30 kcal/mol lower in energy than any of the (acyclic) HN=CH-CH=CH-NH-BO conformers; the boron-oxygen bond in this borocycle has significant double bond character, a bonding scheme for which there has been only one experimental structure reported in the literature (Vidovic, D. ; et al. J. Am. Chem. Soc. 2005, 127, 4566- 4569).