Acceleration of an aromatic Claisen rearrangement via a designed spiroligozyme catalyst that mimics the ketosteroid isomerase catalytic dyad

Non-Covalent Interactions Enforce Conformation in Switchable and Water-Soluble Diketopiperazine-Pyridine Foldamers

To reach their potential as mimics of the dynamic molecules present in biological systems, foldamers must be designed to display stimulus-responsive behavior. Here we report such a foldamer architecture based on alternating pyridine-diketopiperazine linkers. Epimerization is conveniently prevented through a copper-catalyzed coupling protocol. The compounds' native unswitched conformation is first discovered in the solid and solution state. The foldamers can be solubilized in DMSO and pH 9.5 buffer, retaining conformational control to a large degree. Lastly, dynamic switching is demonstrated through treatment with acid, leading to behaviour we describe as stimulus-responsive sidechain reconfiguration.

A spirocyclic backbone accesses new conformational space in an extended, dipole-stabilized foldamer

Most aromatic foldamers adopt uniform secondary structures, offering limited potential for the exploration of conformational space and the formation of tertiary structures. Here we report the incorporation of spiro bis-lactams to allow controlled rotation of the backbone of an iteratively synthesised foldamer. This enables precise control of foldamer shape along two orthogonal directions, likened to the aeronautical yaw and roll axes. XRD, NMR and computational data suggest that homo-oligomers adopt an extended right-handed helix with a pitch of over 30 Å, approximately that of B-DNA. Compatibility with extant foldamers to form hetero-oligomers is demonstrated, allowing greater structural complexity and function in future hybrid foldamer designs.

DyKAT by DiCat: Stereoconvergent Dienamine-Catalyzed Claisen Rearrangements

This Claisen rearrangement establishes the feasibility of DyKAT of γ-epimeric enals via dienamine formation to afford enantioenriched products. γ-Aryl and -alkyl enals, and exocyclic enals that introduce quaternary centers, are all amenable substrates. Products are readily converted into pyrrolidines or cyclopentenols. Notably, a reactive dienamine intermediate has been isolated from a catalytic reaction, fully characterized, and converted to product upon reexposure to reaction conditions. Product configuration arises from a directing C-H···π interaction in the transition state.

Development of Fmoc-Protected Bis-Amino Acids toward Automated Synthesis of Highly Functionalized Spiroligomers

We report the fluorenylmethoxycarbonyl (Fmoc) protection of functionalized bis-amino acid building blocks using a temporary Cu²⁺ complexation strategy, together with an efficient multikilogram-scale synthesis of bis-amino acid precursors. This allows the synthesis of stereochemically and functionally diverse spiroligomers utilizing solid-phase Fmoc/tBu chemistry to facilitate the development of applications. Four tetramers were assembled on a semiautomated microwave peptide synthesizer. We determined their secondary structures with two-dimensional nuclear magnetic resonance spectroscopy.

Trends in computational molecular catalyst design

Computational methods have emerged as a powerful tool to augment traditional experimental molecular catalyst design by providing useful predictions of catalyst performance and decreasing the time needed for catalyst screening. In this perspective, we discuss three approaches for computational molecular catalyst design: (i) the reaction mechanism-based approach that calculates all relevant elementary steps, finds the rate and selectivity determining steps, and ultimately makes predictions on catalyst performance based on kinetic analysis, (ii) the descriptor-based approach where physical/chemical considerations are used to find molecular properties as predictors of catalyst performance, and (iii) the data-driven approach where statistical analysis as well as machine learning (ML) methods are used to obtain relationships between available data/features and catalyst performance. Following an introduction to these approaches, we cover their strengths and weaknesses and highlight some recent key applications. Furthermore, we present an outlook on how the currently applied approaches may evolve in the near future by addressing how recent developments in building automated computational workflows and implementing advanced ML models hold promise for reducing human workload, eliminating human bias, and speeding up computational catalyst design at the same time. Finally, we provide our viewpoint on how some of the challenges associated with the up-and-coming approaches driven by automation and ML may be resolved.

Metal Cation-Binding Mechanisms of Q-Proline Peptoid Macrocycles in Solution

The rational design of foldable and functionalizable peptidomimetic scaffolds requires the concerted application of both computational and experimental methods. Recently, a new class of designed peptoid macrocycle incorporating spiroligomer proline mimics (Q-prolines) has been found to preorganize when bound by monovalent metal cations. To determine the solution-state structure of these cation-bound macrocycles, we employ a Bayesian inference method (BICePs) to reconcile enhanced-sampling molecular simulations with sparse ROESY correlations from experimental NMR studies to predict and design conformational and binding properties of macrocycles as functional scaffolds for peptidomimetics. Conformations predicted to be most populated in solution were then simulated in the presence of explicit cations to yield trajectories with observed binding events, revealing a highly preorganized all-trans amide conformation, whose formation is likely limited by the slow rate of cis/trans isomerization. Interestingly, this conformation differs from a racemic crystal structure solved in the absence of cation. Free energies of cation binding computed from distance-dependent potentials of mean force suggest Na⁺ has a higher affinity to the macrocycle than K⁺, with both cations binding much more strongly in acetonitrile than water. The simulated affinities are able to correctly rank the extent to which different macrocycle sequences exhibit preorganization in the presence of different metal cations and solvents, suggesting our approach is suitable for solution-state computational design.

Identification of Thermal Conduits That Link the Protein-Water Interface to the Active Site Loop and Catalytic Base in Enolase

We report here on the salient role of protein mobility in accessing conformational landscapes that enable efficient enzyme catalysis. We are focused on yeast enolase, a highly conserved lyase with a TIM barrel domain and catalytic loop, as part of a larger study of the relationship of site selective protein motions to chemical reactivity within superfamilies. Enthalpically hindered variants were developed by replacement of a conserved hydrophobic side chain (Leu 343) with smaller side chains. Leu343 is proximal to the active site base in enolase, and comparative pH rate profiles for the valine and alanine variants indicate a role for side chain hydrophobicity in tuning the pK_a of the catalytic base. However, the magnitude of a substrate deuterium isotope effect is almost identical for wild-type (WT) and Leu343Ala, supporting an unchanged rate-determining proton abstraction step. The introduced hydrophobic side chains at position 343 lead to a discontinuous break in both activity and activation energy as a function of side chain volume. Hydrogen-deuterium exchange mass spectrometry (HDX-MS) experiments were performed as a function of time and temperature for WT and Leu343Ala, and provide a spatially resolved map of changes in protein flexibility following mutation. Impacts on protein flexibility are localized to specific networks that arise at the protein-solvent interface and terminate in a loop that has been shown by X-ray crystallography to close over the active site. These interrelated effects are discussed in the context of long-range, solvent-accessible and thermally activated networks that play key roles in tuning the precise distances and interactions among reactants.

Model Selection Using BICePs: A Bayesian Approach for Force Field Validation and Parameterization

The Bayesian Inference of Conformational Populations (BICePs) algorithm reconciles theoretical predictions of conformational state populations with sparse and/or noisy experimental measurements. Among its key advantages is its ability to perform objective model selection through a quantity we call the BICePs score, which reflects the integrated posterior evidence in favor of a given model, computed through free energy estimation methods. Here, we explore how the BICePs score can be used for force field validation and parametrization. Using a 2D lattice protein as a toy model, we demonstrate that BICePs is able to select the correct value of an interaction energy parameter given ensemble-averaged experimental distance measurements. We show that if conformational states are sufficiently fine-grained, the results are robust to experimental noise and measurement sparsity. Using these insights, we apply BICePs to perform force field evaluations for all-atom simulations of designed β-hairpin peptides against experimental NMR chemical shift measurements. These tests suggest that BICePs scores can be used for model selection in the context of all-atom simulations. We expect this approach to be particularly useful for the computational foldamer design as a tool for improving general-purpose force fields given sparse experimental measurements.

One-Pot Synthesis of Chiral, Spirocyclic 4-Hydantoin-Proline Derivatives for Incorporation into Spiroligomer-Based Macromolecules

Derivatives of 4-hydantoin-proline have been synthesized via a direct two-step alkylation method. This method is valuable in the development of applications of N,N'-disubstituted hydantoin bearing α-amino acids by improving yields, reducing the time and number of steps required to synthesize these substituted molecules, and enabling late stage functionalization of spiroligomer termini. Over 20 unique electrophiles have been tested, highlighting the inherent versatility of this chemistry.

Claisen rearrangements of benzyl vinyl ethers: theoretical investigation of mechanism, substituent effects, and regioselectivity

Theoretical calculations are reported which examine the mechanisms of Claisen rearrangements of benzyl vinyl ethers and the ways in which substituents influence reactivity and regioselectivity.

Mechanism of Silver-Mediated Geminal Difluorination of Styrenes with a Fluoroiodane Reagent: Insights into Lewis-Acid-Activation Model

Fluorination mediated by the cyclic hypervalent fluoroiodane reagent (1) often requires an exogenous Lewis acid. The widely accepted Lewis-acid-activation model is that a given Lewis acid binds to the oxygen atom of 1 (O-coordination) to polarize the I-O bond. Computational studies of silver-mediated geminal difluorination of styrenes with 1 reveal a new "F-coordination" model that is energetically much preferred over the commonly accepted "O-coordination" model. The calculations rationalize the regioselective formation of the geminal difluorination product.

Synthesis of Spiroligomer-Containing Macrocycles

We demonstrate the synthesis and characterization of the solution conformations of a collection of functionalized spiroligomer-based macrocycles. These macrocycles contain 14 independently controllable stereocenters and four independently controllable functional groups on a highly preorganized scaffold. These molecules are being developed to display complex, preorganized surfaces for binding proteins and to create enzyme-like active sites. In this work, we demonstrate the convergent synthetic approach to this new class of macrocycles and demonstrate that the conformational properties of these molecules can be changed by altering the configuration stereocenters within the backbone.

Combining a Clostridial enzyme exhibiting unusual active site plasticity with a remarkably facile sigmatropic rearrangement: rapid, stereocontrolled entry into densely functionalized fluorinated phosphonates for chemical biology

Described is an efficient stereocontrolled route into valuable, densely functionalized fluorinated phosphonates that takes advantage of (i) a Clostridial enzyme to set the absolute stereochemistry and (ii) a new [3,3]-sigmatropic rearrangement of the thiono-Claisen variety that is among the fastest sigmatropic rearrangements yet reported. Here, a pronounced rate enhancement is achieved by distal fluorination. This rearrangement is completely stereoretentive, parlaying the enzymatically established β-C-O stereochemistry in the substrate into the δ-C-S stereochemistry in the product. The final products are of interest to chemical biology, with a platform for Zn-aminopeptidase A inhibitors being constructed here. The enzyme, Clostridium acetobutylicum (CaADH), recently expressed by our group, reduces a spectrum of γ,δ-unsaturated β-keto-α,α-difluorophosphonate esters (93-99% ee; 10 examples). The resultant β-hydroxy-α,α-difluorophosphonates possess the "L"-stereochemistry, opposite to that previously observed for the CaADH-reduction of ω-keto carboxylate esters ("D"), indicating an unusual active site plasticity. For the thiono-Claisen rearrangement, a notable structure-reactivity relationship is observed. Measured rate constants vary by over 3 orders of magnitude, depending upon thiono-ester structure. Temperature-dependent kinetics reveal an unusually favorable entropy of activation (ΔS(‡) = 14.5 ± 0.6 e.u.). Most notably, a 400-fold rate enhancement is seen upon fluorination of the distal arene ring, arising from favorable enthalpic (ΔΔH(‡) = -2.3 kcal/mol) and entropic (ΔΔS(‡) = 4 e.u., i.e. 1.2 kcal/mol at rt) contributions. The unusual active site plasticity seen here is expected to drive structural biology studies on CaADH, while the exceptionally facile sigmatropic rearrangement is expected to drive computational studies to elucidate its underlying entropic and enthalpic basis.

Quantification of electrophilic activation by hydrogen-bonding organocatalysts

A spectrophotometric sensor is described that provides a useful assessment of the LUMO-lowering provided by catalysts in Diels–Alder and Friedel–Crafts reactions. A broad range of 33 hydrogen-bonding catalysts was assessed with the sensor, and the relative rates in the above reactions spanned 5 orders of magnitude as determined via ¹H- and ²H NMR spectroscopic measurements, respectively. The differences between the maximum wavelength shift of the sensor with and without catalyst (Δλ_max^–1) were found to correlate linearly with ln(k_rel) values for both reactions, even though the substrate feature that interacts with the catalyst differs significantly (ketone vs nitro). The sensor provides an assessment of both the inherent reactivity of a catalyst architecture as well as the sensitivity of the reaction to changes within an architecture. In contrast, catalyst pK_a values are a poor measure of reactivity, although correlations have been identified within catalyst classes.