CIFDock: A novel CHARMM-based flexible receptor-flexible ligand docking protocol.

Docking studies play a critical role in the current workflow of drug discovery. However, limitations may often arise through factors including inadequate ligand sampling, a lack of protein flexibility, scoring function inadequacies (e.g., due to metals, co-factors, etc.), and difficulty in retaining explicit water molecules. Herein, we present a novel CHARMM-based induced fit docking (CIFDock) workflow that can circumvent these limitations by employing all-atom force fields coupled to enhanced sampling molecular dynamics procedures. Self-guided Langevin dynamics simulations are used to effectively sample relevant ligand conformations, side chain orientations, crystal water positions, and active site residue motion. Protein flexibility is further enhanced by dynamic sampling of side chain orientations using an expandable rotamer library. Steps in the procedure consisting of fixing individual components (e.g., the ligand) while sampling the other components (e.g., the residues in the active site of the protein) allow for the complex to adapt to conformational changes. Ultimately, all components of the complex-the protein, ligand, and waters-are sampled simultaneously and unrestrained with SGLD to capture any induced fit effects. This modular flexible docking procedure is automated using CHARMM scripting, interfaced with SLURM array processing, and parallelized to use the desired number of processors. We validated the CIFDock procedure by performing cross-docking studies using a data set comprised of 21 pharmaceutically relevant proteins. Five variants of the CHARMM-based SWISSDOCK scoring functions were created to quantify the results of the final generated poses. Results obtained were comparable to, or in some cases improved upon, commercial docking program data.

Ameliorative Properties of Boronic Compounds in In Vitro and In Vivo Models of Alzheimer's Disease.

Alzheimer's disease (AD) is characterized by amyloid (Aß) aggregation, hyperphosphorylated tau, neuroinflammation, and severe memory deficits. Reports that certain boronic compounds can reduce amyloid accumulation and neuroinflammation prompted us to compare trans-2-phenyl-vinyl-boronic-acid-MIDA-ester (TPVA) and trans-beta-styryl-boronic-acid (TBSA) as treatments of deficits in in vitro and in vivo models of AD. We hypothesized that these compounds would reduce neuropathological deficits in cell-culture and animal models of AD. Using a dot-blot assay and cultured N2a cells, we observed that TBSA inhibited Aß42 aggregation and increased cell survival more effectively than did TPVA. These TBSA-induced benefits were extended to C. elegans expressing Aß42 and to the 5xFAD mouse model of AD. Oral administration of 0.5 mg/kg dose of TBSA or an equivalent amount of methylcellulose vehicle to groups of six- and 12-month-old 5xFAD or wild-type mice over a two-month period prevented recognition- and spatial-memory deficits in the novel-object recognition and Morris-water-maze memory tasks, respectively, and reduced the number of pyknotic and degenerated cells, Aß plaques, and GFAP and Iba-1 immunoreactivity in the hippocampus and cortex of these mice. These findings indicate that TBSA exerts neuroprotective properties by decreasing amyloid plaque burden and neuroinflammation, thereby preventing neuronal death and preserving memory function in the 5xFAD mice.

Modeling Boronic Acid Based Fluorescent Saccharide Sensors: Computational Investigation of d-Fructose Binding to Dimethylaminomethylphenylboronic Acid.

Designing organic saccharide sensors for use in aqueous solution is a nontrivial endeavor. Incorporation of hydrogen bonding groups on a sensor's receptor unit to target saccharides is an obvious strategy but not one that is likely to ensure analyte-receptor interactions over analyte-solvent or receptor-solvent interactions. Phenylboronic acids are known to reversibly and covalently bind saccharides (diols in general) with highly selective affinity in aqueous solution. Therefore, recent work has sought to design such sensors and understand their mechanism for allowing fluorescence with bound saccharides. In past work, binding orientations of several saccharides were determined to dimethylaminomethylphenylboronic acid (DMPBA) receptors with an anthracene fluorophore; however, the binding orientation of d-fructose to such a sensor could not be determined. In this work, we investigate the potential binding modes by generating 20 possible bidentate and six possible tridentate modes between fructose and DMPBA, a simplified receptor model. Gas phase and implicit solvent geometry optimizations, with a myriad functional/basis set pairs, were carried out to identify the lowest energy bidentate and tridentate binding modes of d-fructose to DMPBA. An interesting hydrogen transfer was observed during selected bidentate gas phase optimizations; this transfer suggests a strong sharing of the hydrogen atom between the boronate hydroxyl and amine nitrogen.

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.

Elucidating a chemical defense mechanism of Antarctic sponges: A computational study.

In 2000, a novel secondary metabolite (erebusinone, Ereb) was isolated from the Antarctic sea sponge, Isodictya erinacea. The bioactivity of Ereb was investigated, and it was found to inhibit molting when fed to the arthropod species Orchomene plebs. Xanthurenic acid (XA) is a known endogenous molt regulator present in arthropods. Experimental studies have confirmed that XA inhibits molting by binding to either (or both) of two P450 enzymes (CYP315a1 or CYP314a1) that are responsible for the final two hydroxylations in the production of the molt-inducing hormone, 20-hydroxyecdysone (20E). The lack of crystal structures and biochemical assays for CYP315a1 or CYP314a1, has prevented further experimental exploration of XA and Ereb's molt inhibition mechanisms. Herein, a wide array of computational techniques - homology modeling, molecular dynamics simulations, binding site bioinformatics, flexible receptor-flexible ligand docking, and molecular mechanics-generalized Born surface area calculations - have been employed to elucidate the structure-function relationships between the aforementioned P450s and the two described small molecule inhibitors (Ereb and XA). Results indicate that Ereb likely targets CYP315a1 by interacting with a network of aromatic residues in the binding site, while XA may inhibit both CYP315a1 and CYP314a1 because of its aromatic, as well as charged nature.

Monosubstituted Phenylboronic Acids, R-B(OH)2 (R = C6H5, C6H4CH3, C6H4NH2, C6H4OH, and C6H4F): A Computational Investigation.

Phenylboronic acids (PBAs) are an important class of compounds with diverse applications in synthetic, biological, medicinal, and materials chemistry. In this investigation we report structural and thermochemical parameters for several monosubstituted ortho, meta, and para PBAs, R-B(OH)2 (R = C6H5, C6H4CH3, C6H4NH2, C6H4OH, and C6H4F). Equilibrium geometries of all the PBAs discussed in this article were obtained using second-order Møller-Plesset perturbation theory (MP2) with the Dunning-Woon aug-cc-pVDZ basis set; heats of formation (HOF) were calculated at the Gaussian-3 (G3) level of theory. The endo-exo conformers of all the positional isomers of these PBAs were lowest in energy. Using HOF for the monosubstituted PBAs calculated at the G3 level of theory, in conjunction with the experimental HOF for benzene, toluene, aniline, phenol, and fluorobenzene, the values of [Formula: see text] for the transfer processes C6H6 + C6H4X-B(OH)2 → C6H5X + C6H5-B(OH)2 (X = CH3, NH2, OH, and F) are found to be in good agreement with values of [Formula: see text] calculated at the MP2(FC)/aug-cc-pVTZ//MP2(FC)/aug-cc-pVTZ computational level; the bonding in the reactants and products for these transfer reactions are well-matched and thermochemical calculations at this level are expected to be very accurate, providing checks on the G3 HOF calculations.

A Comparison of the Structure and Bonding in the Aliphatic Boronic R-B(OH)2 and Borinic R-BH(OH) acids (R=H; NH2, OH, and F): A Computational Investigation.

Boronic acids, R-B(OH)2, play an important role in synthetic, biological, medicinal, and materials chemistry. This investigation compares the structure and bonding surrounding the boron atoms in the simple aliphatic boronic acids, R-B(OH)2 (R = H; NH2, OH, and F) and the analogous borinic acids, R-BH(OH). Geometry optimizations were performed using second-order Møller-Plesset perturbation theory (MP2) with the Dunning-Woon aug-cc-pVTZ, aug-cc-pVQZ and aug-cc-pV5Z basis sets; single-point CCSD(FC)/aug-cc-pVTZ//MP2(FC)/aug-cc-pVTZ level calculations were used to generate a QCI density for Natural Bond Orbital analyses of the bonding. The optimized boron-oxygen bond lengths for the X-B-Ot-H trans-branch of the endo-exo form of the boronic acids and for the X-B-O-H cis-branch of the boronic and borinic acids (X = N, O, and F respectively) decrease as the electronegativity of X increases. The boron-oxygen bond lengths are generally longer in the endo-exo or anti forms of the boronic acids than in the corresponding borinic acids. NBO analyses suggest the boron-oxygen bond in H2BOH is a double bond; the boron-oxygen bonding in the remaining boronic and borinic acids in this study have a significant contribution from dative pπ-pπ bonding. Values for [Formula: see text] for the highly balanced reaction, R-B(OH)2 + R-BH2 → 2 R-BH(OH), suggest that the bonding surrounding the boron atom is stronger in the borinic acid than in the corresponding boronic acid.

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)2) and their boroxine (R3B3O3) 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)2 and R3B3O3 (R = H, Li, HBe, H2B, H3C, H2N, 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)2 → R3B3O3 + 3 H2O; 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 H2N-B(OH)2, 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.

Comparison of Three Chain-of-States Methods: Nudged Elastic Band and Replica Path with Restraints or Constraints.

Chain-of-state methods are becoming important tools in studying the chemical reaction mechanisms, especially for biomacromolecules. In this article, three chain-of-state methods, nudged elastic band (NEB) method and the replica path method with restraints or constraints, were tested and compared using three model systems with various sizes and at different levels of theory: alanine dipeptide isomerization, ß-alanine intramolecular condensation, and the matrix metalloproteinase 2 inhibition mechanism. The levels of theory used to describe the three model systems include molecular mechanics (MM), quantum mechanics (QM), and combined quantum mechanics and molecular mechanics (QM/MM). All three methods could correctly determine a reaction path with reasonable estimation of reaction barriers in most cases. The RMSD measurement with additional weighting schemes provides practically infinite choices of reaction coordinates to describe the reaction progress. These findings demonstrate that the chain-of-state methods are powerful tools when being used carefully to generate a plausible reaction mechanism with full pathway for complex systems at an affordable computational cost.

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.

MSCALE: A General Utility for Multiscale Modeling.

The combination of theoretical models of macromolecules that exist at different spatial and temporal scales has become increasingly important for addressing complex biochemical problems. This work describes the extension of concurrent multiscale approaches, introduces a general framework for carrying out calculations, and describes its implementation into the CHARMM macromolecular modeling package. This functionality, termed MSCALE, generalizes both the additive and subtractive multiscale scheme (e.g. QM/MM ONIOM-type), and extends its support to classical force fields, coarse grained modeling (e.g. ENM, GNM, etc.), and a mixture of them all. The MSCALE scheme is completely parallelized with each subsystem running as an independent, but connected calculation. One of the most attractive features of MSCALE is the relative ease of implementation using the standard MPI communication protocol. This allows external access to the framework and facilitates the combination of functionality previously isolated in separate programs. This new facility is fully integrated with free energy perturbation methods, Hessian based methods, and the use of periodicity and symmetry, which allows the calculation of accurate pressures. We demonstrate the utility of this new technique with four examples; (1) subtractive QM/MM and QM/QM calculations; (2) multi-force field alchemical free energy perturbation; (3) integration with the SANDER module of AMBER and the TINKER package to gain access to potentials not available in CHARMM; and (4) mixed resolution (i.e. coarse grain / all-atom) normal mode analysis. The potential of this new tool is clearly established and in conclusion an interesting mathematical problem is highlighted and future improvements are proposed.

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.

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.

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

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

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

Boronic acids are widely used in materials science, pharmacology, and the synthesis of biologically active compounds. In this Article, geometrical structures and relative energies of dimers of boroglycine, H2N-CH2-B(OH)2, and its constitutional isomer H3C-NH-B(OH)2, were computed using second-order Møller-Plesset perturbation theory and density functional theory; Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the MP2 calculations, and the Pople 6-311++G(d,p) basis set was employed for a majority of the DFT calculations. Effects of an aqueous environment were incorporated into the results using PCM and COSMO-RS methodology. The lowest-energy conformer of the H2N-CH2-B(OH)2 dimer was a six-membered ring structure (chair conformation; Ci symmetry) with two intermolecular B:N dative-bonds; it was 14.0 kcal/mol lower in energy at the MP2/aug-cc-pVDZ computational level than a conformer with the classic eight-centered ring structure (Ci symmetry) in which the boroglycine monomers are linked by a pair of H-O...H bonds. Compared to the results of MP2 calculations with correlation-consistent basis sets, DFT calculations using the PBE1PBE and TPSS functionals with the 6-311++G(d,p) basis set were significantly better at predicting relative conformational energies of the H2N-CH2-B(OH)2 and H3C-NH-B(OH)2 dimers than corresponding calculations using the BLYP, B3LYP, OLYP, and O3LYP functionals, particularly with respect to dative-bonded structures.

A computational investigation of the geometrical structure and protodeboronation of boroglycine, H2N-CH2-B(OH)2.

In this article the geometrical structure of the simple, achiral, alpha-amino boronic acid boroglycine, H2N-CH2-B(OH)2, was investigated using density functional theory (DFT), second-order Møller-Plesset (MP2) perturbation theory, and coupled cluster methodology with single- and double-excitations (CCSD); the effects of an aqueous environment were incorporated into the results by using a few explicit water molecules and/or self-consistent reaction field (SCRF) calculations with the IEF polarizable continuum model (PCM). Neutral reaction mechanisms were investigated for the direct protodeboronation (hydrolysis) of boroglycine (H2O+H2N-CH2-B(OH)2-->B(OH)3+H2N-CH3), for which DeltaH degrees 298 was -21.9 kcal/mol at the MP2(FC)/aug-cc-pVDZ level, and for the 1,2-carbon-to-nitrogen shift of the -B(OH)2 moiety (H2N-CH2-B(OH)2-->H3C-NH-B(OH)2), for which the corresponding value of DeltaH degrees 298 was -18.2 kcal/mol. A boron-oxygen double-bonded intermediate was found to play an important role in the 1,2-rearrangement mechanism.

Hydrogen bridging in the compounds X2H (X=Al,Si,P,S).

X2H hydrides (X=Al, Si, P, and S) have been investigated using coupled cluster theory with single, double, and triple excitations, the latter incorporated as a perturbative correction [CCSD(T)]. These were performed utilizing a series of correlation-consistent basis sets augmented with diffuse functions (aug-cc-pVXZ, X=D, T, and Q). Al2H and Si2H are determined to have H-bridged C2v structures in their ground states: the Al2H ground state is of 2B1 symmetry with an Al-H-Al angle of 87.6 degrees, and the Si2H ground state is of 2A1 symmetry with a Si-H-Si angle of 79.8 degrees. However, P2H and S2H have nonbridged, bent Cs structures: the P2H ground state is of 2A' symmetry with a P-P-H angle of 97.0 degrees, and the S2H ground state is of 2A' symmetry with a S-S-H angle of 93.2 degrees. Ground state geometries, vibrational frequencies, and electron affinities have been computed at all levels of theory. Our CCSD(T)/aug-cc-pVQZ adiabatic electron affinity of 2.34 eV for the Si2H radical is in excellent agreement with the photoelectron spectroscopy experiments of Xu et al. [J. Chem. Phys. 108, 7645 (1998)], where the electron affinity was determined to be 2.31+/-0.01 eV.

Structure of the boronic acid dimer and the relative stabilities of its conformers.

Despite the widespread use of boronic acids in materials science and as pharmaceutical agents, many aspects of their structure and reactivity are not well understood. In this research the boronic acid dimer, [HB(OH)(2)](2), was studied by second-order Møller-Plesset (MP2) perturbation theory and coupled cluster methodology with single and double excitations (CCSD). Pople split-valence 6-31+G*, 6-311G**, and 6-311++G** and Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the calculations. A doubly hydrogen-bonded conformer (1) of the dimer was consistently found to be lowest in energy; the structure of 1 was planar (C(2h)) at most computational levels employed but was significantly nonplanar (C(2)) at the MP2/6-311++G** and CCSD/6-311++G** levels, the result of an intrinsic problem with Pople-type sp-diffuse basis functions on heavy atoms. The dimerization energy, enthalpy, and free energy for the formation of (1) from the exo-endo conformer of the monomer were -10.8, -9.2, and +1.2 kcal/mol, respectively, at the MP2/aug-cc-pVTZ level. Several other hydrogen-bonded conformers of the dimer were local minima on the potential energy surface (PES) and ranged from 2 to 5 kcal/mol higher in energy than 1. Nine doubly OH-bridged conformers, in which the boron atoms were tetracoordinated, were also local minima on the PES, but they were all greater than 13 kcal/mol higher in energy than 1; doubly H-bridged structures proved to be transition states. MP2 and CCSD results were compared to those from the BLYP, B3LYP, OLYP, O3LYP, PBE1PBE, and TPSS functionals with the 6-311++G** and aug-cc-pVTZ basis sets; the PBE1PBE functional performed best relative to the MP2 and CCSD results. Self-consistent reaction field (SCRF) calculations predict that boronic acid dimerization is less favorable in solution than in vacuo.

Assessing alkyl-, silyl-, and halo-substituent effects on the electron affinities of silyl radicals.

Neutral anion energy differences for a large class of alpha-substituted silyl radicals have been computed to determine the effect of alkyl, silyl, and halo substituents on their electron affinities. In particular, we report theoretical predictions of the adiabatic electron affinities (AEAs), vertical electron affinities (VEAs), and vertical detachment energies (VDEs) for a series of methyl-, silyl-, and halo-substituted silyl radical compounds. This work utilizes the carefully calibrated DZP++ basis set, in conjunction with the pure BLYP and OLYP functionals, as well as with the hybrid B3LYP, BHLYP, PBE1PBE, MPW1K, and O3LYP functionals. Bromine has the largest effect in stabilizing the anions, and the BLYP/DZP++ AEA for SiBr(3) is 3.29 eV. The other predicted electron affinities are for SiH(3) (1.37 eV), SiH(2)CH(3) (1.09 eV), SiH(2)F (1.54 eV), SiH(2)Cl (1.94 eV), SiH(2)Br (2.05 eV), SiH(2)(SiH(3)) (1.77 eV), SiH(CH(3))(2) (0.92 eV), SiHF(2) (1.86 eV), SiHCl(2) (2.53 eV), SiHBr(2) (2.67 eV), Si(CH(3))(3) (0.86 eV), SiF(3) (2.66 eV), SiCl(3) (3.21 eV), Si(SiH(3))(3) (2.25 eV), and SiFClBr (3.13 eV). For the five silyl radicals where experimental data are available, the BLYP functional gives the most accurate determination of AEAs; the average absolute error is 0.04(1) eV, whereas the corresponding errors for the O3LYP, MPW1K, PBE1PBE, B3LYP, OLYP, and BHLYP functionals are 0.05(8), 0.06(0), 0.06(3), 0.08(5), 0.11(5), and 0.15(3) eV, respectively.