Molecular Dynamics of the Norbornyl Cation in Solution and Its Generation in Winstein-Trifan Solvolysis: The Timing of Sigma Bridging.

Molecular dynamics simulations were performed on the solvolyses of exo- and endo-norbornyl brosylate and for the "nonclassical" σ-bridged norbornyl cation in an acetic acid solution. This computational modeling of the original Winstein-Trifan experiment confirms that exo-solvolysis is accompanied by σ-bridging in the transition state, while endo-solvolysis is not; σ-bridging eventually occurs in a dynamically stepwise fashion. Simulations of the norbornyl cation in solution show typical vibrations due to zero-point and thermal vibrations but no tendency to sample localized "classical cation" geometries.

Cycloaddition Reactivities Analyzed by Energy Decomposition Analyses and the Frontier Molecular Orbital Model.

This Account describes our quest to understand and predict organic reactivity, a principal goal of physical and theoretical organic chemistry. The focus is on the development and testing of models for the prediction of cycloaddition reactivities and selectivities. We describe the involvement of the Houk group, and other groups, in the evolution of theoretical models that can achieve ever greater accuracy as well as provide practical heuristic models for understanding and prediction.Is the venerable frontier molecular orbital (FMO) model, the basis of Kenichi Fukui's 1981 Nobel Prize, still useful, or must it be replaced with more advanced models? In particular, models such as Conceptual Density Functional, the Pauli Exclusion Model, and the recent popularity of Electrostatic Potential Plots and Dispersion Energies have not only added to our understanding, but they have also created uncertainty about whether the simple FMO heuristic model has a place in 21st century discussions. This Account addresses this issue and asserts the value of the FMO model.Beginning with brief descriptions of selected models for cycloaddition reactivity starting with early donor-acceptor (nucleophile-electrophile) charge-transfer concepts, this Account reviews Fukui's frontier molecular orbital model, Salem and Klopman's orbital, electrostatic and Pauli repulsion model, the conceptual DFT model by Parr and later by Domingo and others, the recent Houk and Bickelhaupt Distortion/Interaction Activation Strain model, and the Bickelhaupt-Hamlin's Pauli-repulsion lowering model.Computations and analyses of four well-studied Diels-Alder cycloadditions, both normal and inverse electron-demand types, are presented. Most were studied earlier in our published work but are presented here with new insights from calculations with modern methods. Depending on the types of substrates (cycloaddends), the dominant factors controlling reactivity can be orbital interactions, electrostatics and polarization, or Pauli repulsion and dispersion effects, or a combination of all of these.By comparing orbital interactions, especially the frontier molecular orbital interactions, with the other factors that influence reactivity, we show why the FMO model is such a powerful─and theoretically meaningful─heuristic for understanding and predicting reactivity. We also present a method to understand Pauli repulsion effects on activation barriers, ρ(1.1). The use of a new reaction coordinate, the extent of Pauli repulsion along the reaction path, is advocated to emphasize the role of repulsive occupied orbital interactions on reactivity.Fukui's frontier molecular orbital model is effective because FMO interactions parallel all the quantities that influence reactivity. The FMO model continues to provide a practical model to understand and guide experiments.

Origins of Periselectivity and Regioselectivity in Ambimodal Tripericyclic [8+6]/[6+4]/[4+2] Intramolecular Cycloadditions of a Heptafulvenyl-Fulvene.

Quantum mechanical calculations and molecular dynamics simulations have elucidated the reaction mechanism for intramolecular cycloadditions of a heptafulvenyl-fulvene tethered by a trimethylene chain. Prior experiments by Liu and Houk reported the formation of only an endo-[8+6] cycloadduct at 185 °C. Liu et al. later reported an exo-[4+2] Diels-Alder cycloadduct as the major product at 140 °C (Tetrahedron, 1999, 55, 9171). Cycloadditions involve Diels-Alder and an ambimodal intramolecular tripericyclic [8+6]/[6+4]/[4+2] cycloaddition. The mechanistic details explain the experimental reports of temperature dependence on the periselectivity of intramolecular cycloadditions. Additional calculations with multireference-based methods CASSCF and NEVPT2 highlight the artifacts of DFT methods and single-reference wavefunction-based CCSD(T) in the description of complete potential energy surface involving various cycloadditions of the heptafulvenyl-fulvene.

An Enzymatic Platform for Primary Amination of 1-Aryl-2-alkyl Alkynes.

Propargyl amines are versatile synthetic intermediates with numerous applications in the pharmaceutical industry. An attractive strategy for efficient preparation of these compounds is nitrene propargylic C(sp3)-H insertion. However, achieving this reaction with good chemo-, regio-, and enantioselective control has proven to be challenging. Here, we report an enzymatic platform for the enantioselective propargylic amination of alkynes using a hydroxylamine derivative as the nitrene precursor. Cytochrome P450 variant PA-G8 catalyzing this transformation was identified after eight rounds of directed evolution. A variety of 1-aryl-2-alkyl alkynes are accepted by PA-G8, including those bearing heteroaromatic rings. This biocatalytic process is efficient and selective (up to 2610 total turnover number (TTN) and 96% ee) and can be performed on preparative scale.

Probative Evidence and Quantitative Energetics for the Unimolecular Mechanism of the 1,5-Sigmatropic Hydrogen Shift of Cyclopentadiene.

While the 1,5-sigmatropic hydrogen shift in cyclopentadiene is generally thought to be a unimolecular pericyclic reaction, Yamabe proposed a more complex bimolecular mechanism proceeding through the exo dimer of cyclopentadiene. DFT computations by Yamabe were claimed to show that the bimolecular mechanism was kinetically more favorable than the unimolecular mechanism. Reinvestigation of the unimolecular concerted mechanism and Yamabe's bimolecular mechanism with ωB97X-D and DLPNO-CCSD(T) calculations demonstrates a 25 kcal/mol preference for the unimolecular mechanism relative to the bimolecular mechanism. While Yamabe's calculations were performed with the less accurate B3LYP functional, the incorrect conclusion was the result of a different error discovered here. We have also computed corrections for tunneling that result in computed activation barriers within 1.5 kcal/mol of the experimental values.

Cycloadditions of Cyclopentadiene and Cycloheptatriene with Tropones: All Endo-[6+4] Cycloadditions Are Ambimodal.

The cycloadditions of cyclopentadiene and cycloheptatriene with tropone are some of the earliest published examples of [6+4] cycloaddition reactions. We report quantum mechanical studies (ωB97X-D and DLPNO-CCSD(T)) of transition structures and products of these reactions, as well as quasi-classical molecular dynamics simulations of reaction trajectories. The study reveals that these cycloadditions involve ambimodal transition states resulting in a web of products by pericyclic interconversion pathways. Combined with these studies, calculations of simple parent systems and a thorough meta-analysis of literature examples reveal the general concept that all endo-[6+4] cycloadditions are ambimodal.

Parameterization of Monovalent Ions for the OPC3, OPC, TIP3P-FB, and TIP4P-FB Water Models.

Monovalent ions play significant roles in various biological and material systems. Recently, four new water models (OPC3, OPC, TIP3P-FB, and TIP4P-FB), with significantly improved descriptions of condensed phase water, have been developed. The pairwise interaction between the metal ion and water necessitates the development of ion parameters specifically for these water models. Herein, we parameterized the 12-6 and the 12-6-4 nonbonded models for 12 monovalent ions with the respective four new water models. These monovalent ions contain eight cations including alkali metal ions (Li+, Na+, K+, Rb+, Cs+), transition-metal ions (Cu+ and Ag+), and Tl+ from the boron family, along with four halide anions (F-, Cl-, Br-, I-). Our parameters were designed to reproduce the target hydration free energies (the 12-6 hydration free energy (HFE) set), the ion-oxygen distances (the 12-6 ion-oxygen distance (IOD) set), or both of them (the 12-6-4 set). The 12-6-4 parameter set provides highly accurate structural features overcoming the limitations of the routinely used 12-6 nonbonded model for ions. Specifically, we note that the 12-6-4 parameter set is able to reproduce experimental hydration free energies within 1 kcal/mol and experimental ion-oxygen distances within 0.01 Å simultaneously. We further reproduced the experimentally determined activity derivatives for salt solutions, validating the ion parameters for simulations of ion pairs. The improved performance of the present water models over our previous parameter sets for the TIP3P, TIP4P, and SPC/E water models (Li, P. et al J. Chem. Theory Comput. 2015 11 1645 1657) highlights the importance of the choice of water model in conjunction with the metal ion parameter set.

Thermodynamics of Transition Metal Ion Binding to Proteins.

Modeling the thermodynamics of a transition metal (TM) ion assembly be it in proteins or in coordination complexes affords us a better understanding of the assembly and function of metalloclusters in diverse application areas including metal organic framework design, TM-based catalyst design, the trafficking of TM ions in biological systems, and drug design in metalloprotein platforms. While the structural details of TM ions bound to metalloproteins are generally well understood via experimental and computational approaches, accurate studies describing the thermodynamics of TM ion binding are rare. Herein, we demonstrate that we can obtain accurate structural and absolute binding free energies of Co2+ and Ni2+ to the enzyme glyoxalase I using an optimized 12-6-4 (m12-6-4) potential. Critically, this model simultaneously reproduces the solvation free energy of the individual TM ions and reproduces the thermodynamics of TM ion-ligand coordination as well as the thermodynamics of TM ion binding to a protein active site unlike extant models. We find the incorporation of the thermodynamics associated with protonation state changes for the TM ion (un)binding to be crucial. The high accuracy of m12-6-4 potential in this study presents an accurate route to explore more complicated processes associated with TM cluster assembly and TM ion transport.

Eliminating Systematic Errors in DFT via Connectivity-Based Hierarchy: Accurate Bond Dissociation Energies of Biodiesel Methyl Esters.

We present a computational study focusing on the determination of accurate bond dissociation energies (BDEs) involved in the combustion of biodiesel methyl esters. We have adapted our previously developed efficient error-cancellation protocols, based on the systematic "connectivity-based hierarchy" (CBH), to derive accurate BDEs of biodiesel molecules at a modest computational cost. Using DFT energies on the full biodiesel molecule in conjunction with accurate G4 energies on the small fragments involved in the CBH reaction schemes, systematic errors in the DFT methods can be cancelled efficiently. Herein, we apply our G4-corrected ΔCBH-2 and ΔCBH-3 schemes in conjunction with several popular DFT methods to determine accurate bond dissociation energies of different C-C, C-H, and C-O bonds in biodiesel surrogate molecules. We first evaluate the performance of different DFT methods using a test set of 21 reactions involving various bond dissociations in small to medium biodiesel surrogates (up to methyl decanoate, a C10-methyl ester) by calibration against accurate values calculated with multireference methods (MRACPF2), reported by Carter and co-workers. The CBH-2 corrections for all tested dispersion-corrected functionals yield mean absolute deviations (MADs) in a narrow range of 1.3-1.5 kcal/mol, the best performance being obtained for B97-D3 and ωB97X-D functionals (MAD = 1.3 kcal/mol). Further, significant improvement yielding a MAD of only 0.9 kcal/mol is obtained using the G4-corrected CBH-3 scheme. Finally, the protocol has been applied to derive accurate BDEs of eight different bonds in the larger biodiesel molecule, methyl linolenate, yielding a MAD of only 1.13 kcal/mol using the ΔCBH-3 error correction scheme. The results suggest that our protocol in conjunction with different DFT methods should be broadly applicable to yield accurate BDEs for a variety of large biodiesel molecules.

Simulating the Chelate Effect.

Despite the rich history of experimental studies focusing on the thermochemistry and kinetics associated with the chelate effect, molecular-level computational studies on the chelate ring opening/ring closure are scarce. The challenge lies in an accurate description of both the metal ion and its aqueous environment. Herein, we demonstrate that an optimized 12-6-4 Lennard-Jones (LJ) model can capture the thermodynamics and provide detailed structural and mechanistic insights into the formation of ethylenediamine (en) complexes with metal ions. The water molecules in the first solvation shell of the metal ion are found to facilitate the chelate ring formation. The optimized parameters further simulate the formation of bis and tris(en) complexes representing the wide applicability of the model to simulate coordination chemistry and self-assembly processes.

Anion-Binding Macrocycles Operate Beyond the Electrostatic Regime: Interaction Distances Matter.

Anion recognition impacts many areas of chemistry and often relies on receptors with multiple hydrogen-bond donors. Previous studies of these donors in small molecules have long promoted the idea that electrostatic interactions alone correlate with association strength, yet this correlation has not been critically evaluated in the framework of larger, macrocyclic receptors. Here, we provide that assessment by evaluating how much electrostatics contributes to the gas-phase binding energy of macrocyclic receptors with various anions. Whereas small-molecule complexes behave as expected, we find that electrostatic interactions fail to accurately describe total binding energies of many common macrocyclic receptors: calix[4]pyrroles, dipyrrolyldiketones, indolocarbazoles, amido-pyrroles, triazolophanes, and cyanostars. This deviation arises from the fact that most macrocycles have multiple points of contact with the anion. Whereas the hydrogen-bond donors collectively stabilize the anion, the interaction distances are typically larger than equilibrium values seen with small molecules. This leads to increases in the relative contributions of the attractive components such as induction (e.g., induced dipoles) and dispersion, which are found to be as high as 32 % for CH-donor based tricarbazole triazolophane complex with large polarizable ClO4- . This study augments previous observations of the importance of dispersion and induction towards anion binding of macrocyclic receptors in solution.

Accurate Thermochemistry for Organic Cations via Error Cancellation using Connectivity-Based Hierarchy.

Connectivity-Based Hierarchy (CBH) is an effective error-cancellation scheme for the determination of chemically accurate thermochemical properties of a variety of organic and biomolecules. Neutral molecules and open-shell radicals have already been treated successfully by this approach utilizing inexpensive computational methods such as density functional theory. Herein, we present an extension of the method to a new class of molecules, specifically, organic cations. Because of the presence of structural rearrangements involving hydrogen migrations as well as unusual structures such as bridged cations, the application of the standard CBH protocol to a test set of 25 cations leads to significant errors due to ineffective bond-type matching. We propose an adjusted protocol to overcome such limitations to achieve highly effective error cancellation. The modified CBH methods, in conjunction with a wide range of density functionals, reproduce G4 energies for the test set of organic cations accurately within 1-2 kcal/mol at a reduced computational cost.

Flexibility Coexists with Shape-Persistence in Cyanostar Macrocycles.

Shape-persistent macrocycles are attractive functional targets for synthesis, molecular recognition, and hierarchical self-assembly. Such macrocycles are noncollapsible and geometrically well-defined, and they are traditionally characterized by having repeat units and low conformational flexibility. Here, we find it necessary to refine these ideas in the face of highly flexible yet shape-persistent macrocycles. A molecule is shape-persistent if it has a small change in shape when perturbed by external stimuli (e.g., heat, light, and redox chemistry). In support of this idea, we provide the first examination of the relationships between a macrocycle's shape persistence, its conformational space, and the resulting functions. We do this with a star-shaped macrocycle called cyanostar that is flexible as well as being shape-persistent. We employed molecular dynamics (MD), density functional theory (DFT), and NMR experiments. Considering a thermal bath as a stimulus, we found a single macrocycle has 332 accessible conformers with olefins undergoing rapid interconversion by up-down and in-out motions on short time scales (0.2 ns). These many interconverting conformations classify single cyanostars as flexible. To determine and confirm that cyanostars are shape-persistent, we show that they have a high 87% shape similarity across these conformations. To further test the idea, we use the binding of diglyme to the single macrocycle as guest-induced stimulation. This guest has almost no effect on the conformational space. However, formation of a 2:1 sandwich complex involving two macrocycles enhances rigidity and dramatically shifts the conformer distribution toward perfect bowls. Overall, the present study expands the scope of shape-persistent macrocycles to include flexible macrocycles if, and only if, their conformers have similar shapes.

Breaking a bottleneck: Accurate extrapolation to "gold standard" CCSD(T) energies for large open shell organic radicals at reduced computational cost.

Open Shell organic radicals are principal species involved in many diverse areas such as combustion, photochemistry, and polymer chemistry. Computational studies of such species with an accurate method like coupled-cluster with single and double and perturbative triple (CCSD(T)) may be restricted to systems of modest size due to the steep computational scaling of the method. Herein, we assess the accuracy of extrapolated CCSD(T) energies determined using the connectivity-based hierarchy (CBH) method on medium to large sized radicals. In our method, an MP2 calculation on the target radical is coupled with CCSD(T) energies of fragments determined uniquely by our hierarchy to perform accurate extrapolations. A careful assessment is done with a robust CBH-rad49 test set comprising of 49 diverse cyclic and acyclic radicals with a variety of functional groups. We demonstrate that the extrapolation method with CBH-2 or CBH-3 is sufficient to obtain sub-kcal accuracy. ROMP2 and PMP2 calculations with both Pople-style and Dunning-style basis-sets resulted in mean absolute errors for CCSD(T) extrapolation (full CCSD(T)-extrapolated CCSD(T)) within 0.5 kcal/mol. Further speedup for such CCSD(T) extrapolations are obtained with ROHF-based RI-MP2 calculations. Challenging systems with (a) high ring strain, (b) delocalized character, and (c) spin contamination are identified and analyzed in detail. Finally, we apply our extrapolation method on 10 larger radicals containing 10-15 heavy atoms, where accurate CCSD(T) energies are obtained at a fractional cost of full CCSD(T) calculations.

Anions Stabilize Each Other inside Macrocyclic Hosts.

Contrary to the simple expectations from Coulomb's law, Weinhold proposed that anions can stabilize each other as metastable dimers, yet experimental evidence for these species and their mutual stabilization is missing. We show that two bisulfate anions can form such dimers, which stabilize each other with self-complementary hydrogen bonds, by encapsulation inside a pair of cyanostar macrocycles. The resulting 2:2 complex of the bisulfate homodimer persists across all states of matter, including in solution. The bisulfate dimer's OH⋅⋅⋅O hydrogen bonding is seen in a 1 H NMR peak at 13.75 ppm, which is consistent with borderline-strong hydrogen bonds.

Electrostatic and Allosteric Cooperativity in Ion-Pair Binding: A Quantitative and Coupled Experiment-Theory Study with Aryl-Triazole-Ether Macrocycles.

Cooperative binding of ion pairs to receptors is crucial for the manipulation of salts, but a comprehensive understanding of cooperativity has been elusive. To this end, we combine experiment and theory to quantify ion-pair binding and to separate allostery from electrostatics to understand their relative contributions. We designed aryl-triazole-ether macrocycles (MC) to be semiflexible, which allows ion pairs (NaX; X = anion) to make contact, and to be monocyclic to simplify analyses. A multiequilibrium model allows us to quantify, for the first time, the experimental cooperativity, α, for the equilibrium MC·Na(+) + MC·X(-) ⇌ MC·NaX + MC, which is associated with contact ion-pair binding of NaI (α = 1300, ΔGα = -18 kJ mol(-1)) and NaClO4 (α = 400, ΔGα = -15 kJ mol(-1)) in 4:1 dichloromethane-acetonitrile. We used accurate energies from density functional theory to deconvolute how the electrostatic effects and the allosteric changes in receptor geometry individually contribute to cooperativity. Computations, using a continuum solvation model (dichloromethane), show that allostery contributes ∼30% to overall positive cooperativity. The calculated trend of electrostatic cooperativity using pairs of spherical ions (NaCl > NaBr > NaI) correlates to experimental observations (NaI > NaClO4). We show that intrinsic ionic size, which dictates charge separation distance in contact ion pairs, controls electrostatic cooperativity. This finding supports the design principle that semiflexible receptors can facilitate optimal electrostatic cooperativity. While Coulomb's law predicts the size-dependent trend, it overestimates electrostatic cooperativity; we suggest that binding of the individual anion and cation to their respective binding sites dilutes their effective charge. This comprehensive understanding is critical for rational designs of ion-pair receptors for the manipulation of salts.

Stereoselective Construction of ß-, γ-, and δ-Lactam Rings via Enzymatic C-H Amidation.

Lactam rings are found in many biologically active natural products and pharmaceuticals, including important classes of antibiotics. Methods for the asymmetric synthesis of these molecules are therefore highly desirable, particularly through the selective functionalization of unreactive aliphatic C-H bonds. Here we show the development of a strategy for the asymmetric synthesis of ß-, γ-, and δ-lactams via hemoprotein-catalysed intramolecular C-H amidation reaction with readily available dioxazolone reagents. Engineered myoglobin variants serve as excellent biocatalysts for this transformation yielding the desired lactam products in high yields, high enantioselectivity, and on preparative scale. Mechanistic and computational studies elucidate the nature of the C-H amination and enantiodetermining steps and provide insights into protein-mediated control of regioselectivity and stereoselectivity. Additionally, an alkaloid natural product and a drug molecule were synthesized chemoenzymatically in much fewer steps (7-8 vs. 11-12) than previously reported, further demonstrating the power of biosynthetic strategy for the preparation of complex bioactive molecules.

Biocatalytic strategy for the construction of sp3-rich polycyclic compounds from directed evolution and computational modelling.

Catalysis with engineered enzymes has provided more efficient routes for the production of active pharmaceutical agents. However, the potential of biocatalysis to assist in early-stage drug discovery campaigns remains largely untapped. In this study, we have developed a biocatalytic strategy for the construction of sp3-rich polycyclic compounds via the intramolecular cyclopropanation of benzothiophenes and related heterocycles. Two carbene transferases with complementary regioisomer selectivity were evolved to catalyse the stereoselective cyclization of benzothiophene substrates bearing diazo ester groups at the C2 or C3 position of the heterocycle. The detailed mechanisms of these reactions were elucidated by a combination of crystallographic and computational analyses. Leveraging these insights, the substrate scope of one of the biocatalysts could be expanded to include previously unreactive substrates, highlighting the value of integrating evolutionary and rational strategies to develop enzymes for new-to-nature transformations. The molecular scaffolds accessed here feature a combination of three-dimensional and stereochemical complexity with 'rule-of-three' properties, which should make them highly valuable for fragment-based drug discovery campaigns.