220 likes | 244 Views
Explore methods for predicting reaction selectivity, including forward prediction, reduction, chemoselectivity, regioselectivity, and enantioselectivity. Learn about late-stage functionalization, C-H activation, electrophilic aromatic bromination, radical bromination, Pd-catalyzed processes, and calculating reaction selectivity using Boltzmann summation. Discover modeling catalyst stereoselectivity with QSAR, QSSR, QM energy screening, and Q2MM molecular mechanics. Develop virtual screening workflows, force fields, and asymmetric hydrogenation tests. Collaborate with experts for reliable chemical selectivity predictions.
E N D
Methods for predicting reaction selectivity Per-Ola Norrby BigChem Training School 9 May 2019
Chemical reactionselectivity, forward prediction Reduction Chemoselectivity Regioselectivity Enantioselectivity
Chemical reactionselectivity, cont. Late StageFunctionalization C-H activation ElectrophilicAromaticbromination Radicalbromination Pd-catalyzed C-H activation
G‡ Calculating reaction selectivity Boltzmann summation 1-10 paths 1000’s possible!
Modeling catalyst stereoselectivity (I): QSAR QSSR: QuantitativeStructureSelectivity Relationship Norrby, Lipkowits, Kozlowski, Bo, Harvey, Fey, ..., Sigman, Doyle, …
Modeling catalyst stereoselectivity (II): QM (R)- vs. (S)-energy, DFT-based screening Maeda, Morokuma, Zimmerman, ..., Wheeler 27% ee → 97% ee MUE 4 kJ/mol
Modeling catalyst stereoselectivity (III): Q2MM MolecularMechanics @ TS, conformationalsearch
Molecular Mechanics for Reactivity Westheimer 1946 Houk Transitionstate force fields, TSFF De Tar, Paterson, … 1950 1960 1970 1980 1990 Sprecher 1965 Menger 1989 ”Fudgit” Using 200 parameters to fit 20 data points is not a goodidea Can the overfitting problem be solved?
Parameterization from QM data • QMcalculationsgive abundant data • Determining force fields from quantum mechanical data: • Hagler, Kollman, Goddard, … • Bond lengths, angles from optimizedstructures • Charges from electrostatic potentials • Torsions from calculated rotation profiles • Force constants from calculated vibrations (Hessian data) • Problem: wrongcurvature
Reaction coordinate modification Diagonalization: Eigenvalue replacement: Forming the new Hessian:
Q2MM results MUE = 2.8 kJ/mol
Virtual screening ofcatalysts: CatVS • Buildstructures • Conformationalsearch • Calculatee.r.
Rh-catalyzedenamidehydrogenation Collaboration with Paul Helquist and Olaf Wiest, Notre Dame, IN Q2MM force field for TS Olaf Wiest Paul Helquist
Asymmetrichydrogenation – Test I ΔΔG* = 10 kJ/mol e.r. = 50 98:2 e.e. = 96% Underlying program changed - Energies no longer the same Current version: no model solvent Developed for: Tested on: Enrichment, 4 of 5!
Asymmetrichydrogenation – Test 2 ΔΔG* = 10 kJ/mol e.r. = 50 98:2 e.e. = 96% Developed for: Tested on: • 1 good, 1 fair, 1 false negative • Acceptable
CatVS - Summary Combinewithligandslibrary Draw substrate & Click ”submit” Email results for all ligands Determine all lowenergy TS • Needs Q2MM force fields; significant development time • Development tools on www.github.com/Q2MM/ • Significant enrichment • Makes experimental testing more efficient • Faster routes to chiral compounds Chemical selectivity can now be reliably predicted!
Ongoing Christian Sköld Maria Vergou • Asymmetric Suzuki-Miyaura Liam Byrne Peter Smith Rachel Munday
Hybrid ML-mechanistic model Machinelearning Mechanisticmodel Hybrid model Electrophilicaromatic substitution: Tomberg, Johansson, Norrby, J. Org. Chem. 2019,4695
Acknowledgment AstraZeneca Tobias Rein Per Ryberg Simone Tomasi Elaine Limé David Buttar Rachel Munday Maria Vergou Peter Smith Liam Byrne Anna Tomberg Notre Dame Paul Helquist Olaf Wiest Patrick Donoghue Elsa Kieken Vincenzo Verdolino Aaron Forbes Eric Hansen Brandon Tutkowski Anthony Rosales Taylor Quinn Jessica Wahlers Xin Zhang (Shenzhen) Uppsala Christian Sköld Denmark Tommy Liljefors Torben Rasmussen Peter Brandt Peter Fristrup Lund Ulf Ryde Patrik Rydberg Stockholm Björn Åkermark $$$ AstraZeneca, The University of Gothenburg, COST, FP7/SYNFLOW, The Swedish Research Council, C3SE/Gothenburg