The 1,3-dipolar cycloaddition of cyclooctynes with azides, also called "copper-free click chemistry", is a bioorthogonal reaction with widespread applications in biological discovery. The kinetics of this reaction are of paramount importance for studies of dynamic processes, particularly in living subjects. Here we performed a systematic analysis of the effects of strain and electronics on the reactivity of cyclooctynes with azides through both experimental measurements and computational studies using a density functional theory (DFT) distortion/interaction transition state model. In particular, we focused on biarylazacyclooctynone (BARAC) because it reacts with azides faster than any other reported cyclooctyne and its modular synthesis facilitated rapid access to analogs. We found that substituents on BARAC's aryl rings can alter the calculated transition state interaction energy through electronic effects or the calculated distortion energy through steric effects. Experimental data confirmed that electronic perturbation of BARAC's aryl rings has a modest effect on reaction rate, whereas steric hindrance in the transition state can significantly retard the reaction. Drawing on these results, we analyzed the relationship between alkyne bond angles, which we determined using X-ray crystallography, and reactivity, quantified by experimental second-order rate constants, for a range of cyclooctynes. Our results suggest a correlation between decreased alkyne bond angle and increased cyclooctyne reactivity. Finally, we obtained structural and computational data that revealed the relationship between the conformation of BARAC's central lactam and compound reactivity. Collectively, these results indicate that the distortion/interaction model combined with bond angle analysis will enable predictions of cyclooctyne reactivity and the rational design of new reagents for copper-free click chemistry.

The mechanism of cyclohexyne insertion into a C(O)-Cα bond of cyclic ketones, explored experimentally by the Carreira group, has been investigated using density functional theory. B3LYP and M06-2X calculations were performed in both gas phase and THF (CPCM, UAKS radii). The reaction proceeds through a stepwise [2+2] cycloaddition of cyclohexyne to the enolate, followed by three disparate ring opening possibilities of the cyclobutene alkoxide to give the product: 1) thermally-allowed conrotatory electrocyclic ring opening, 2) thermally-forbidden disrotatory electrocyclic ring opening, or 3) non- pericyclic C-C bond cleavage. Our computational results for the model alkoxide and potassium alkoxide systems show that the thermally-allowed electrocyclic ring opening pathway is favored by less than 1 kcal/mol. In more complex systems containing a potassium alkoxide (e-f), the barrier of the allowed conrotatory ring opening is disfavored by 4-8 kcal/mol. This suggests that the thermodynamically more stable disrotatory product can be formed directly through a "forbidden" pathway. Analysis of geometrical parameters and atomic charges throughout the ring opening pathways provides evidence for a non-pericyclic C-C bond cleavage, rather than a thermally-forbidden disrotatory ring opening. A true forbidden disrotatory ring opening transition structure was computed for the cyclobutene alcohol; however, it was 19 kcal/mol higher in energy than the allowed conrotatory transition structure. An alternate mechanism in which the disrotatory product forms via isomerization of the conrotatory product was also explored for the alkoxide and potassium alkoxide systems.

The decomposition of a Z-selective ruthenium metathesis catalyst and structurally similar analogs has been investigated utilizing X-ray crystallography and density functional theory. Isolated X-ray crystal structures suggest that recently reported C-H activated catalysts undergo decomposition via insertion of the alkylidene moiety into the chelating ruthenium-carbon bond followed by hydride elimination, which is supported by theoretical calculations. The resulting ruthenium hydride intermediates have been implicated in previously observed olefin migration, and thus lead to unwanted byproducts in cross metathesis reactions. Preventing these decomposition modes will be essential in the design of more active and selective Z-selective catalysts.

Derivatives of 4-hydroxyproline with a series of hydrophobic groups in well-defined orientations have been tested as catalysts for the aldol reactions. All of the modified proline catalysts carry out the intermolecular aldol reaction in water and provide high diastereoselectivity and enantioselectivity. Modified prolines with aromatic groups syn to the carboxylic acid are better catalysts than those with small hydrophobic groups (1a is 43.5 times faster than 1f). Quantum mechanical calculations provide transition structures, TS-1awater and TS-1fwater that support the hypothesis that a stabilizing hydrophobic interaction occurs with 1a.

A software suite, SABER (Selection of Active/Binding sites for Enzyme Redesign), has been developed for the analysis of atomic geometries in protein structures, using a geometric hashing algorithm (Barker and Thornton, Bioinformatics 2003;19:1644-1649). SABER is used to explore the Protein Data Bank (PDB) to locate proteins with a specific 3D arrangement of catalytic groups to identify active sites that might be redesigned to catalyze new reactions. As a proof-of-principle test, SABER was used to identify enzymes that have the same catalytic group arrangement present in o-succinyl benzoate synthase (OSBS). Among the highest-scoring scaffolds identified by the SABER search for enzymes with the same catalytic group arrangement as OSBS were L-Ala D/L-Glu epimerase (AEE) and muconate lactonizing enzyme II (MLE), both of which have been redesigned to become effective OSBS catalysts, demonstrated by experiments. Next, we used SABER to search for naturally existing active sites in the PDB with catalytic groups similar to those present in the designed Kemp elimination enzyme KE07. From over 2000 geometric matches to the KE07 active site, SABER identified 23 matches that corresponded to residues from known active sites. The best of these matches, with a 0.28 Å catalytic atom RMSD to KE07, was then redesigned to be compatible with the Kemp elimination using RosettaDesign. We also used SABER to search for potential Kemp eliminases using a theozyme predicted to provide a greater rate acceleration than the active site of KE07, and used Rosetta to create a design based on the proteins identified.

Quantum mechanical calculations reveal the origin of diastereo- and enantioselectivities of aldol reactions between aldehydes catalyzed by histidine, and differences between related reactions catalyzed by proline. A stereochemical model that explains both the sense and the high levels of the experimentally observed stereoselectivity is proposed. The computations suggest that both the imidazolium and the carboxylic acid functionalities of histidine are viable hydrogen-bond donors that can stabilize the cyclic aldolization transition state. The stereoselectivity is proposed to arise from minimization of gauche interactions around the forming C-C bond.

Polycyclic polyether natural products have fascinated chemists and biologists alike owing to their useful biological activity, highly complex structure and intriguing biosynthetic mechanisms. Following the original proposal for the polyepoxide origin of lasalocid and isolasalocid and the experimental determination of the origins of the oxygen and carbon atoms of both lasalocid and monensin, a unified stereochemical model for the biosynthesis of polyether ionophore antibiotics was proposed. The model was based on a cascade of nucleophilic ring closures of postulated polyepoxide substrates generated by stereospecific oxidation of all-trans polyene polyketide intermediates. Shortly thereafter, a related model was proposed for the biogenesis of marine ladder toxins, involving a series of nominally disfavoured anti-Baldwin, endo-tet epoxide-ring-opening reactions. Recently, we identified Lsd19 from the Streptomyces lasaliensis gene cluster as the epoxide hydrolase responsible for the epoxide-opening cyclization of bisepoxyprelasalocid A to form lasalocid A. Here we report the X-ray crystal structure of Lsd19 in complex with its substrate and product analogue to provide the first atomic structure-to our knowledge-of a natural enzyme capable of catalysing the disfavoured epoxide-opening cyclic ether formation. On the basis of our structural and computational studies, we propose a general mechanism for the enzymatic catalysis of polyether natural product biosynthesis.

The alcoholysis of cyclic meso-anhydrides catalyzed by β-amino alcohols has been investigated with DFT quantum mechanics to determine the mechanism of this reaction. Both nucleophilic catalysis and general base catalysis pathways are explored for methanol-induced ring opening of an anhydride catalyzed by a chiral amino alcohol. The nucleophilic pathway involves a late transition state with a high energy barrier. In this mechanism, methanolysis is expected to take place following the amine-induced ring opening of the anhydride. In the base-catalyzed mechanism, methanol attack on one carbonyl group of the meso-anhydride is assisted by the β-amino alcohol; the amine functionality abstracts the methanol proton. The chiral amino alcohol also catalyzes the reaction by stabilizing the oxyanion that forms upon ring opening of the anhydride by hydrogen bonding with its alcoholic moiety. Both stepwise and concerted pathways have been studied for the general base catalysis route. Transition structures for both are found to be lower in energy than in the nucleophilic mechanism. Overall this study has shed light on the mechanism of the β-amino alcohol-catalyzed alcoholysis of cyclic meso-anhydrides, showing that the nucleophilic pathway is approximately 100 kJ mol(-1) higher in energy than the general base pathway.

A general approach for the computational design of enzymes to catalyze arbitrary reactions is a goal at the forefront of the field of protein design. Recently, computationally designed enzymes have been produced for three chemical reactions through the synthesis and screening of a large number of variants. Here, we present an iterative approach that has led to the development of the most catalytically efficient computationally designed enzyme for the Kemp elimination to date. Previously established computational techniques were used to generate an initial design, HG-1, which was catalytically inactive. Analysis of HG-1 with molecular dynamics simulations (MD) and X-ray crystallography indicated that the inactivity might be due to bound waters and high flexibility of residues within the active site. This analysis guided changes to our design procedure, moved the design deeper into the interior of the protein, and resulted in an active Kemp eliminase, HG-2. The cocrystal structure of this enzyme with a transition state analog (TSA) revealed that the TSA was bound in the active site, interacted with the intended catalytic base in a catalytically relevant manner, but was flipped relative to the design model. MD analysis of HG-2 led to an additional point mutation, HG-3, that produced a further threefold improvement in activity. This iterative approach to computational enzyme design, including detailed MD and structural analysis of both active and inactive designs, promises a more complete understanding of the underlying principles of enzymatic catalysis and furthers progress toward reliably producing active enzymes.

Think before you act: a computational approach is reported for evaluating the synthetic potential of heterocyclic arynes. Routine and rapid calculations of arene dehydrogenation energies and aryne angle distortion predict the likelihood that a given hetaryne can be generated, as well as the degree of regioselectivity expected in a reaction between a given hetaryne and a nucleophilic trapping agent.