Joseph D. Larkin

A computational investigation of the nitrogen-boron interaction in o-(N,N-dialkylaminomethyl)arylboronate systems.

o-(N,N-Dialkylaminomethyl)arylboronate systems are an important class of compounds in diol-sensor development. We report results from a computational investigation of fourteen o-(N,N-dialkylaminomethyl)arylboronates using second-order Møller-Plesset (MP2) perturbation theory. Geometry optimizations were performed at the MP2/cc-pVDZ level and followed by single-point calculations at the MP2/aug-cc-pVDZ(cc-pVTZ) levels. These results are compared to those from density functional theory (DFT) at the PBE1PBE(PBE1PBE-D)/6-311++G(d,p)(aug-cc-pVDZ) levels, as well as to experiment. Results from continuum PCM and CPCM solvation models were employed to assess the effects of a bulk aqueous environment. Although the behavior of o-(N,N-dialkylaminomethyl) free acid and ester proved to be complicated, we were able to extract some important trends from our calculations: (1) for the free acids the intramolecular hydrogen-bonded B-O-H···N seven-membered ring conformers 12 and 16 are found to be slightly lower in energy than the dative-bonded N→B five-membered ring conformers 10 and 14 while conformers 13 and 17, with no direct boron-nitrogen interaction, are significantly higher in energy than 12 and 16; (2) for the esters where no intramolecular B-O-H···N bonded form is possible, the N→B conformers 18 and 21 are significantly lower in energy than the no-interaction forms 20 and 23; (3) H(2)O insertion reactions into the N→B structures 10, 14, 18, and 21 leading to the seven-membered intermolecular hydrogen-bonded B···OH(2)···N ring structures 11, 15, 19, and 22 are all energetically favorable.

The thymine radicals and their respective anions: molecular structures and electron affinities

The adiabatic electron affinities (AEA), vertical electron affinities (VEA) and vertical detachment energies (VDE) are predicted for six different radicals derived from thymine by the removal of one hydrogen atom. Geometry optimizations were carried out utilizing the DFT functionals B3LYP, BLYP and BP86 with double-c quality basis sets plus polarization and diffuse functions (DZP++). These methods have been carefully calibrated for the prediction of electron affinities [RIENSTRA-KIRACOFE, J. C., TSCHUMPER, G. S., SCHAEFER, H. F., NANDI, S., and ELLISON, G. B., 2002, Chem. Rev., 102, 231]. All optimized structures were confirmed to be minima via vibrational frequency analyses. Both the neutrals and the anion radicals were observed to possess Cs, symmetry, conserving the parent molecules qualitative conformation. The electron affinities ranged between 1.04 and 3.74eV for the different radicals, contrasting to the small electron affinities associated with the closed-shell thymine species. The radicals with a hydrogen atom removed from one of the nitrogens present the largest electron affinities of all six radicals investigated (3.22eV for NI and 3.74eV for N3).

Thermodynamics of boroxine formation from the aliphatic boronic acid monomers R-B(OH)2 (R = H, H3C, H2N, HO, and F): a computational investigation.

Boroxines are the six-membered cyclotrimeric dehydration products of organoboronic acids, 3R–B(OH)2 → R3B3O3 + 3H2O, and in recent years have emerged as a useful class of organoboron molecules with applications in organic synthesis both as reagents and catalysts, as structural components in boronic-acid-derived pharmaceutical agents, and as anion acceptors and electrolyte additives for battery materials [Korich, A. L.; Iovine, P. M. Dalton Trans. 2010, 39, 1423−1431]. Second-order Møller–Plesset perturbation theory, in conjunction with the Dunning–Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets, was used to investigate the structures and relative energies of the endo–exo, anti, and syn conformers of the aliphatic boronic acids R–B(OH)2 (R = H, H3C, H2N, HO, and F), as well as the thermodynamics of their boroxine formation; single-point calculations at the MP2/aug-cc-pVQZ, MP2/aug-cc-pV5Z, and CCSD(T)/aug-cc-pVTZ levels using the MP2/aug-cc-pVTZ optimized geometries were also performed in selected cases. The endo–exo conformer was generally lowest in energy in vacuo, as well as in PCM and CPCM models of aqueous and carbon tetrachloride media. The values of ΔH(298)(0) for boroxine formation via dehydration from the endo–exo conformers of these aliphatic boronic acids ranged from −2.9 for (H2N)3B3O3 to +12.2 kcal/mol for H3B3O3 at the MP2/aug-cc-pVTZ level in vacuo; for H3B3O3, the corresponding values in PCM/UFF implicit carbon tetrachloride and aqueous media were +11.2 and +9.8 kcal/mol, respectively. On the basis of our calculations, we recommend that ΔHf(298K) for boroxine listed in the JANAF compilation needs to be revised from −290.0 to approximately −277.0 kcal/mol.

Boronic acid based photoinduced electron transfer (PET) fluorescence sensors for saccharides

A simple three step synthesis was developed to provide six novel modular sensors, consisting of three para sensors, and three meta sensors with naphthalene, anthracene and pyrene fluorophores. The interaction of the six sensors with the saccharides: D-glucose, D-fructose, D-galactose, and D-mannose, were evaluated. All sensors displayed increasing fluorescence intensity upon the addition of these saccharides, with all of the sensors showing enhanced selectivity for D-glucose over D-galactose, D-fructose and D-mannose. High affinity (Kobs) was also observed for the meta sensors with respect to the para sensors. The naphthalene and anthracene meta sensors showed particularly high affinity (Kobs) for D-galactose. Circular dichroism spectroscopy was used to probe the structures of the complexes formed. Cyclic complexes were formed between all six sensors and D-glucose. Whilst naphthalene and anthracene meta sensors which displayed high affinity for D-galactose also formed cyclic complexes with that saccharide.

Efficient Calculation of QM/MM Frequencies with the Mobile Block Hessian

The calculation of the analytical second derivative matrix (Hessian) is the bottleneck for vibrational analysis in QM/MM systems when an electrostatic embedding scheme is employed. Even with a small number of QM atoms in the system, the presence of MM atoms increases the computational cost dramatically: the long-range Coulomb interactions require that additional coupled perturbed self-consistent field (CPSCF) equations need to be solved for each MM atom displacement. This paper presents an extension to the Mobile Block Hessian (MBH) formalism for QM/MM calculations with blocks in the MM region and its implementation in a parallel version of the Q-Chem/CHARMM interface. MBH reduces both the CPU time and the memory requirements compared to the standard full Hessian QM/MM analysis, without the need to use a cutoff distance for the electrostatic interactions. Special attention is given to the treatment of link atoms which are usually present when the QM/MM border cuts through a covalent bond. Computational efficiency improvements are highlighted using a reduced chorismate mutase enzyme system, consisting of 24 QM atoms and 306 MM atoms, as a test example. In addition, the drug bortezomib, used for cancer treatment of myeloma, has been studied as a test case with multiple MBH block choices and both a QM and QM/MM description. The accuracy of the calculated Hessians is quantified by imposing Eckart constraints, which allows for the assessment of numerical errors in second derivative procedures. The results show that MBH within the QM/MM description not only is a computationally attractive method but also produces accurate results.

Disaggregation is a Mechanism for Emission Turn-On of ortho-Aminomethylphenylboronic Acid-Based Saccharide Sensors

ortho-Aminomethylphenylboronic acid-based receptors with appended fluorophores are commonly used as molecular sensors for saccharides in aqueous media. The mechanism for fluorescence modulation in these sensors has been attributed to some form of photoinduced electron transfer (PET) quenching, which is diminished in the presence of saccharides. Using a well-known boronic acid-based saccharide sensor (3), this work reveals a new mechanism for fluorescence turn-on in these types of sensors. Compound 3 exhibits an excimer, and the associated ground-state aggregation is responsible for fluorescence modulation under certain conditions. When fructose was titrated into a solution of 3 in 2:1 water/methanol with NaCl, the fluorescence intensity increased. Yet, when the same titration was repeated in pure methanol, a solvent in which the sensor does not aggregate, no fluorescence response to fructose was observed. This reveals that the fluorescence increase is not fully associated with fructose binding, but instead disaggregation of the sensor in the presence of fructose. Further, an analogue of the sensor that does not contain a boronic acid (4) responded nearly identically to 3 in the presence of fructose, despite having no functional group with which to bind the saccharide. This further supports the claim that fluorescence modulation is not primarily a result of binding, but of disaggregation. Using an indicator displacement assay and isothermal titration calorimetry, it was confirmed that fructose does indeed bind to the sensor. Thus, our evidence reveals that while binding occurs with fructose in the aqueous solvent system used, it is not related to the majority of the fluorescence modulation. Instead, disaggregation dominates the signal turn-on, and is thus a mechanism that should be investigated in other ortho-aminomethylphenylboronic acid-based sensors.

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).

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).

The ability of silylenes to bind excess electrons: Electron affinities of SiX2, and SiXY species (X, Y=H, CH3, SiH3, F, Cl, Br)

Recently, Ishida and co-workers have isolated silylene radical anions via the one-electron reduction of isolable cyclic dialkylsilylenes, discovering these corresponding radical anions to be relatively stable at low temperatures. Herein we report theoretical predictions of the adiabatic electron affinities (AEA), vertical electron affinities, and vertical detachment energies of a series of methyl, silyl, and halosubstituted silylene compounds. This research utilizes the carefully calibrated DZP++ basis with the combination of the popular nonhybrid and hybrid DFT functionals, BLYP, B3LYP, and BHHLYP. The level of theory employed and the ensemble of species under study confirm the ability of silylenes to bind excess electrons with Si(SiH(3))(2) being the most effective, having a predicted AEA of 1.95 eV. While it is known that methyl substituents have a diminishing effect on the computed electron affinities (EAs), it is shown that fluorine shows an analogous negative effect. Similarly, previous suggestions that Si(CH(3))(2) will not bind an electron appear incorrect, with EA[Si(CH(3))(2)] predicted here to be 0.46 eV.

Reaction Path Optimization and Sampling Methods and Their Applications for Rare Events

Reaction mechanisms are an important tool for chemists in the determination of thermodynamic and kinetic properties of chemical reactions.(Hanggi & Borkovec 1990; Heidrich 1995; March 1992; Tolman 1925) The mechanisms are integral in the understanding of detailed molecular or chemical transitions from one equilibrium state (reactant) to another equilibrium state (product). In computational chemistry, the reaction mechanism is often represented as a reaction path on the Born-Oppenheimer potential energy surface (PES) of the system of interest through construction of a potential energy function of the nuclear coordinates.(Bader & Gangi 1975; Lewars 2011; Mezey 1987; Truhlar 2001; Wales 2003) The PES serves as an important theoretical construct to provide a framework to describe the transition between different states in detail. The equilibrium states correspond to local minimum on the PES with zero first order derivatives (gradient) in all directions and all positive eigenvalues of the second order derivative (Hessian) matrix, excluding rotation and translation degrees of freedom. The transition states (TSs), based on the transition state theory (TST),(Doll 2005; Eyring 1935; Laidler & King 1983; Pechukas 1981; Truhlar et al. 1983; Wigner 1938; Yamamoto 1960) are the first order saddle points with zero gradient and only one negative eigenvalue of the Hessian matrix. The equilibrium states are often easy to identify through experimental or computational studies. Understanding the detailed transition process between equilibrium states is of more interest in research, but unfortunately is very difficult to study experimentally. On a given PES, one can imagine that there could exist an infinite number of possible routes connecting two predefined states on that surface. However, not every route has the same weight in elucidating of reaction mechanisms. In the static point of view, the minimum energy path (MEP) is the route that needs the least amount of potential energy for the system to undergo the transition. The MEP connecting two local minima must go through one or more TSs, and is identified as a representative reaction path. In the statistical point of view, the minimum free energy path (MFEP) is the most probable transition path connects two metastable states. The simulation of the systems either through molecular dynamics (MD) or Monte Carlo (MC) sampling on the PES could generate an ensemble of transition paths, from which the MFEP can be identified. Both an MEP and an MFEP can be used to predict important properties, such as a reaction’s kinetic isotope effect.