Matthew N. Grayson

Mechanistic Insights into the BINOL-Derived Phosphoric Acid-Catalyzed Asymmetric Allylboration of Aldehydes

BINOL-derived phosphoric acids catalyze the asymmetric allylboration of aldehydes. DFT and QM/MM hybrid calculations showed that the reaction proceeds via a transition state involving both a hydrogen-bonding interaction from the catalyst hydroxyl group to the pseudoaxial oxygen of the cyclic boronate and a stabilizing interaction from the phosphoryl oxygen of the catalyst to the formyl hydrogen of the aldehyde. These interactions lower the energy of the transition structure and provide extra rigidity to the system. This mechanistic pathway is consistent with the experimentally observed enantioselectivity except in one case. We have used our models predictions to guide our own experimental work. The conflict is resolved in favor of our calculations.

Theory and Modeling of Asymmetric Catalytic Reactions

Modern density functional theory and powerful contemporary computers have made it possible to explore complex reactions of value in organic synthesis. We describe recent explorations of mechanisms and origins of stereoselectivities with density functional theory calculations. The specific functionals and basis sets that are routinely used in computational studies of stereoselectivities of organic and organometallic reactions in our group are described, followed by our recent studies that uncovered the origins of stereocontrol in reactions catalyzed by (1) vicinal diamines, including cinchona alkaloid-derived primary amines, (2) vicinal amidophosphines, and (3) organo-transition-metal complexes. Two common cyclic models account for the stereoselectivity of aldol reactions of metal enolates (Zimmerman-Traxler) or those catalyzed by the organocatalyst proline (Houk-List). Three other models were derived from computational studies described in this Account. Cinchona alkaloid-derived primary amines and other vicinal diamines are venerable asymmetric organocatalysts. For α-fluorinations and a variety of aldol reactions, vicinal diamines form enamines at one terminal amine and activate electrophilically with NH(+) or NF(+) at the other. We found that the stereocontrolling transition states are cyclic and that their conformational preferences are responsible for the observed stereoselectivity. In fluorinations, the chair seven-membered cyclic transition states is highly favored, just as the Zimmerman-Traxler chair six-membered aldol transition state controls stereoselectivity. In aldol reactions with vicinal diamine catalysts, the crown transition states are favored, both in the prototype and in an experimental example, shown in the graphic. We found that low-energy conformations of cyclic transition states occur and control stereoselectivities in these reactions. Another class of bifunctional organocatalysts, the vicinal amidophosphines, catalyzes the (3 + 2) annulation reaction of allenes with activated olefins. Stereocontrol here is due to an intermolecular hydrogen bond that activates the electrophilic partner in this reaction. We have also studied complex organometallic catalysts. Krisches ruthenium-catalyzed asymmetric hydrohydroxyalkylation of butadiene involves two chiral ligands at Ru, a chiral diphosphine and a chiral phosphate. The size of this combination strains the limits of modern computations with over 160 atoms, multiple significant steps, and a variety of ligand coordinations and conformations possible. We found that carbon-carbon bond formation occurs via a chair Zimmerman-Traxler-type transition structure and that a formyl CH···O hydrogen bond from aldehyde CH to phosphate oxygen, as well as steric interactions of the two chiral ligands, control the stereoselectivity.

Cinchona Alkaloid-Catalyzed Asymmetric Conjugate Additions: The Bifunctional Brønsted Acid–Hydrogen Bonding Model

Wynbergs report from 1977 that natural cinchona alkaloids catalyze the asymmetric conjugate addition of aromatic thiols to cycloalkenones is a landmark discovery in hydrogen bonding organocatalysis. Wynberg proposed that this reaction proceeded via the formation of a thiolate-alkylammonium tight ion pair and activation of the enone electrophile by a hydrogen bond from the catalysts hydroxyl group. This reaction model provided the mechanistic basis for understanding Wynbergs reaction and many other asymmetric transformations since. Our quantum mechanical calculations reveal a different model should be used to explain the results: the alkylammonium ion activates the enone by Brønsted acid catalysis, and the catalysts hydroxyl group orients the thiolate nucleophile. The new model rationalizes the stereoselective outcome of Wynbergs reaction and provides a new, general model for asymmetric cinchona organocatalysis.

Cinchona Urea-Catalyzed Asymmetric Sulfa-Michael Reactions: The Brønsted Acid−Hydrogen Bonding Model

The cinchona alkaloid-derived urea-catalyzed asymmetric conjugate addition of aromatic thiols to cycloalkenones was studied using density functional theory (DFT). Deprotonation of the thiol gives a protonated amine that activates the electrophile by Brønsted acid catalysis, while the urea group binds the nucleophilic thiolate by hydrogen bonding. These results demonstrate the generality of the Brønsted acid-hydrogen bonding transition state (TS) model for cinchona alkaloid catalysis that we recently showed to be favored over Wynbergs widely accepted ion pair-hydrogen bonding model and represent the first detailed mechanistic study of a cinchona urea-catalyzed reaction. The conformation of the catalyst methoxy group has a strong effect on the TS, an effect overlooked in previous mechanistic studies of reactions catalyzed by cinchona alkaloids.

Ruthenium-Catalyzed Asymmetric Hydrohydroxyalkylation of Butadiene: The Role of the Formyl Hydrogen Bond in Stereochemical Control

The catalyst generated in situ from RuH2(CO)(PPh3)3, (S)-SEGPHOS, and a chiral phosphoric acid promotes asymmetric hydrohydroxyalkylation of butadiene and affords enantioenriched α-methyl homoallylic alcohols. The observed diastereo- and enantioselectivities are determined by both the chiral phosphine and chiral phosphate ligands. Density functional theory calculations (M06/SDD-6-311G(d,p)-IEFPCM(acetone)//B3LYP/SDD-6-31G(d)) predict that the product distribution is controlled by the kinetics of carbon-carbon bond formation, and this process occurs via a closed-chair Zimmerman-Traxler-type transition structure (TS). Chiral-phosphate-dependent stereoselectivity arising from this TS is enabled through a hydrogen bond between the phosphoryl oxygen and the aldehyde formyl proton present in TADDOL-derived catalysts. This interaction is absent in the corresponding BINOL-derived systems, and the opposite diastereo- and enantioselectivity is observed. Additional factors influencing the stereochemical control are determined.

Asymmetric Boronate Addition to o-Quinone Methides: Ligand Exchange, Solvent Effects, and Lewis Acid Catalysis

Density functional theory calculations suggest that asymmetric boronate addition to o-quinone methides proceeds via a Lewis acid catalyzed process through a closed six-membered transition structure. The BINOL-derived catalyst undergoes an exchange process with the original ethoxide boronate ligands. This activation mode successfully accounts for the sense and level of enantioselectivity observed experimentally. A qualitative model which accurately predicts the observed enantioselectivity has been developed and is consistent with results from our study of ketone propargylation demonstrating the reaction models generality. The effects of replacing the BINOL framework with H8-BINOL have been rationalized.

Efficient Biosynthesis of Fungal Polyketides Containing the Dioxabicyclo-octane Ring System

Aurovertins are fungal polyketides that exhibit potent inhibition of adenosine triphosphate synthase. Aurovertins contain a 2,6-dioxabicyclo[3.2.1]octane ring that is proposed to be derived from a polyene precursor through regioselective oxidations and epoxide openings. In this study, we identified only four enzymes required to produce aurovertin E. The core polyketide synthase produces a polyene α-pyrone. Following pyrone O-methylation by a methyltransferase, a flavin-dependent mono-oxygenase and an epoxide hydrolase can iteratively transform the terminal triene portion of the precursor into the dioxabicyclo[3.2.1]octane scaffold. We demonstrate that a tetrahydrofuranyl polyene is the first stable intermediate in the transformation, which can undergo epoxidation and anti-Baldwin 6-endo-tet ring opening to yield the cyclic ether product. Our results further demonstrate the highly concise and efficient ways in which fungal biosynthetic pathways can generate complex natural product scaffolds.

Chronology of CH···O Hydrogen Bonding from Molecular Dynamics Studies of the Phosphoric Acid-Catalyzed Allylboration of Benzaldehyde

CH···O hydrogen bonds involving formyl groups have been invoked as a crucial factor controlling many asymmetric transformations. We conducted quasi-classical direct molecular dynamics simulations on the phosphoric acid-catalyzed allylboration of benzaldehyde to understand the synergy between the phosphoric acid OH···O hydrogen bond and the secondary CH···O formyl hydrogen bond as the reaction occurs. In the gas phase, both the CH···O and OH···O hydrogen bonds are enhanced from reactants to transition states. In toluene, the trend of H-bond enhancement is observed with a smaller magnitude because of solvent caging. The strength of the formyl hydrogen bond in the TS, a second CH···O interaction between the P═O oxygen and ortho-hydrogen of the phenyl ring and the OH···O hydrogen bond were determined using quantum mechanical calculations (4.6, 1.0, and 14.5 kcal mol-1, respectively).

Base-Mediated Cascade Rearrangements of Aryl-Substituted Diallyl Ethers

Two base-mediated cascade rearrangement reactions of diallyl ethers were developed leading to selective [2,3]-Wittig-oxy-Cope and isomerization-Claisen rearrangements. Both diaryl and arylsilyl-substituted 1,3-substituted propenyl substrates were examined, and each exhibits unique reactivity and different reaction pathways. Detailed mechanistic and computational analysis was conducted, which demonstrated that the role of the base and solvent was key to the reactivity and selectivity observed. Crossover experiments also suggest that these reactions proceed with a certain degree of dissociation, and the mechanistic pathway is highly complex with multiple competing routes.

Using Transition State Modeling To Predict Mutagenicity for Michael Acceptors.

The Ames mutagenicity assay is a long established in vitro test to measure the mutagenicity potential of a new chemical used in regulatory testing globally. One of the key computational approaches to modeling of the Ames assay relies on the formation of chemical categories based on the different electrophilic compounds that are able to react directly with DNA and form a covalent bond. Such approaches sometimes predict false positives, as not all Michael acceptors are found to be Ames-positive. The formation of such covalent bonds can be explored computationally using density functional theory transition state modeling. We have applied this approach to mutagenicity, allowing us to calculate the activation energy required for α,β-unsaturated carbonyls to react with a model system for the guanine nucleobase of DNA. These calculations have allowed us to identify that chemical compounds with activation energies greater than or equal to 25.7 kcal/mol are not able to bind directly to DNA. This allows us to reduce the false positive rate for computationally predicted mutagenicity assays. This methodology can be used to investigate other covalent-bond-forming reactions that can lead to toxicological outcomes and learn more about experimental results.