Enhanced Sampling for Free Energy Profiles with Post-Transition-State Bifurcations

Abstract

We present a method to explore the free energy landscapes of chemical reactions with post-transition-state bifurcations using an enhanced sampling method based on well-tempered metadynamics. Obviating the need for computationally expensive DFT-level ab initio molecular dynamics simulations, we obtain accurate energetics by utilizing a free energy perturbation scheme and deep learning estimator for the single-point energies of substrate configurations. Using a pair of easily interpretable collective variables, we present a quantitative free energy surface that is compatible with harmonic transition state theory calculations and in which the bifurcations are clearly visible. We demonstrate our approach with the example of the SpnF-catalyzed Diels-Alder reaction, a cycloaddition reaction in which post-transition-state bifurcation leads to the [4+2] as well as the [6+4] cycloadduct. We obtain the free energy landscapes for different stereochemical reaction pathways and characterize the mechanistic continuum between relevant reaction channels without explicitly searching for the pertinent transition state structures.

arXiv comment27 pages, 5 figures

References in corpus (16)

• Canonical sampling through velocity-rescaling
• Escaping free-energy minima
• CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
• PLUMED 2: New feathers for an old bird
• Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method
• From Metadynamics to Dynamics
• Rethinking Metadynamics: from bias potentials to probability distributions
• Metadynamics with adaptive Gaussians
• Data-Driven Collective Variables for Enhanced Sampling
• Fast and Uncertainty-Aware Directional Message Passing for Non-Equilibrium Molecules
• Accelerated Computation of Free Energy Profile at ab Initio Quantum Mechanical/Molecular Mechanics Accuracy via a Semi-Empirical Reference Potential. I. Weighted Thermodynamics Perturbation
• From enhanced sampling to reaction profiles
• Molecular Mechanics-Driven Graph Neural Network with Multiplex Graph for Molecular Structures
• Variationally Optimized Free Energy Flooding for Rate Calculation
• Blind search for complex chemical pathways using Harmonic Linear Discriminant Analysis
• Exploration vs Convergence Speed in Adaptive-bias Enhanced Sampling