Mechanistic Analysis of Alkyne Haloboration: A DFT, MP2, and DLPNO-CCSD(T) Study.

Stereocontrol of the alkyne haloboration reaction has received attention in many experimental but few theoretical studies. Here we present a detailed quantum-chemical study of mechanisms leading to Z versus E isomers of haloboration products, considering acetylene and propyne combined with BCl3, BBr3, and BI3. Calculations using B3LYP-D3, MP2, and DLPNO-CCSD(T) methods are used to study polar reactions between the alkyne and BX3 in the absence and presence of an additional halide anion whose content in the reaction mixture can be controlled experimentally. The formation of anti-haloboration products via radical mechanisms is also explored, namely, by adding BX3 to (Z)-halovinyl radical. For the anti-haloboration of propyne, the radical route is prohibited by the regiochemistry of the initiating halopropenyl radical, while the polar route is unlikely due to a competitive allene generation. In contrast, energetically accessible routes exist for both syn- and anti-bromoboration of acetylene; hence, careful control of reaction conditions is necessary to steer the stereochemical outcome. Methodologically, MP2 results correspond better to the DLPNO-CCSD(T) energies than the B3LYP-D3 results in terms of both reaction barrier heights and relative ordering of energetically close stationary points.

Improving IDP theoretical chemical shift accuracy and efficiency through a combined MD/ADMA/DFT and machine learning approach.

This work extends the multi-scale computational scheme for the quantum mechanics (QM) calculations of Nuclear Magnetic Resonance (NMR) chemical shifts (CSs) in proteins that lack a well-defined 3D structure. The scheme couples the sampling of an intrinsically disordered protein (IDP) by classical molecular dynamics (MD) with protein fragmentation using the adjustable density matrix assembler (ADMA) and density functional theory (DFT) calculations. In contrast to our early investigation on IDPs (Pavlíková Precechtelová et al., J. Chem. Theory Comput., 2019, 15, 5642-5658) and the state-of-the art NMR calculations for structured proteins, a partial re-optimization was implemented on the raw MD geometries in vibrational normal mode coordinates to enhance the accuracy of the MD/ADMA/DFT computational scheme. In addition, machine-learning based cluster analysis was performed on the scheme to explore its potential in producing protein structure ensembles (CLUSTER ensembles) that yield accurate CSs at a reduced computational cost. The performance of the cluster-based calculations is validated against results obtained with conventional structural ensembles consisting of MD snapshots extracted from the MD trajectory at regular time intervals (REGULAR ensembles). CS calculations performed with the refined MD/ADMA/DFT framework employed the 6-311++G(d,p) basis set that outperformed IGLO-III calculations with the same density functional approximation (B3LYP) and both explicit and implicit solvation. The partial geometry optimization did not universally improve the agreement of computed CSs with the experiment but substantially decreased errors associated with the ensemble averaging. A CLUSTER ensemble with 50 structures yielded ensemble averages close to those obtained with a REGULAR ensemble consisting of 500 MD frames. The cluster based calculations thus required only a fraction of the computational time.

Free Radical Isomerizations in Acetylene Bromoboration Reaction.

The experimentally motivated question of the acetylene bromoboration mechanism was addressed in order to suggest possible radical isomerization pathways for the syn-adduct. Addition-elimination mechanisms starting with a bromine radical attack at the "bromine end" or the "boron end" of the C=C bond were considered. Dispersion-corrected DFT and MP2 methods with the SMD solvation model were employed using three all-electron bases as well as the ECP28MWB ansatz. The rate-determining, elimination step had a higher activation energy (12 kcal mol-1) in case of the "bromine end" attack due to intermediate stabilization at both the MP2 and DFT levels. In case of the "boron end" attack, two modes of C-C bond rotation were followed and striking differences in MP2 vs. DFT potential energy surfaces were observed. Employing MP2, addition was followed by either a 180° rotation through an eclipsed conformation of vicinal bromine atoms or by an opposite rotation avoiding that conformation, with 5 kcal mol-1 of elimination activation energy. Within B3LYP, the addition and rotation proceeded simultaneously, with a 9 (7) kcal mol-1 barrier for rotation involving (avoiding) eclipsed conformation of vicinal bromines. For weakly bound complexes, ZPE corrections with MP2 revealed significant artifacts when diffuse bases were included, which must be considered in the Gibbs free energy profile interpretation.

Stereoselective Bromoboration of Acetylene with Boron Tribromide: Preparation and Cross-Coupling Reactions of (Z)-Bromovinylboronates.

The mechanism of acetylene bromoboration in neat boron tribromide was studied carefully by means of experiment and theory. Besides the syn-addition mechanism through a four-center transition state, radical and polar anti-addition mechanisms are postulated, both triggered by HBr, which is evidenced also to take part in the Z/E isomerization of the product. The proposed mechanism is well supported by ab initio calculations at the MP2/6-31+G* level with Ahlrichs' SVP all-electron basis for Br. Implicit solvation in CH2Cl2 has been included using the PCM and/or SMD continuum solvent models. Comparative case studies have been performed involving the B3LYP/6-31+G* with Ahlrichs' SVP for Br and MP2/Def2TZVPP levels. The mechanistic studies resulted in development of a procedure for stereoselective bromoboration of acetylene yielding E/Z mixtures of dibromo(bromovinyl)borane with the Z-isomer as a major product (up to 85%). Transformation to the corresponding pinacol and neopentyl glycol boronates and stereoselective decomposition of their E-isomer provided pure (Z)-(2-bromovinyl)boronates in 57-60% overall yield. Their reactivity in a Negishi cross-coupling reaction was tested. An example of the one-pot reaction sequence of Negishi and Suzuki-Miyaura cross-couplings for synthesis of combretastatin A4 is also presented.