Abstract:The observation of complex structural transitions in biological and abiological molecular objects within time scales amenable to molecular dynamics (MD) simulations is often hampered by significant free energy barriers associated with entangled movements. Importance-sampling algorithms, a powerful class of numerical schemes for the investigation of rare events, have been widely used to extend simulations beyond the time scale common to MD. However, probing processes spanning milliseconds through microsecond molecular simulations still constitutes in practice a daunting challenge because of the difficulty of taming the ruggedness of multidimensional free energy surfaces by means of naive transition coordinates. To address this limitation, in recent years we have elaborated importance-sampling methods relying on an adaptive biasing force (ABF). In this Account, we review recent developments of algorithms aimed at mapping rugged free energy landscapes that correspond to complex processes of physical, chemical, and biological relevance. Through these developments, we have broadened the spectrum of applications of the popular ABF algorithm while improving its computational efficiency, notably for multidimensional free energy calculations. One major algorithmic advance, coined meta-eABF, merges the key features of metadynamics and an extended Lagrangian variant of ABF (eABF) by simultaneously shaving the barriers and flooding the valleys of the free energy landscape, and it possesses a convergence rate up to 5-fold greater than those of other importance-sampling algorithms. Through faster convergence and enhanced ergodic properties, meta-eABF represents a significant step forward in the simulation of millisecond-time-scale events. Here we introduce extensions of the algorithm, notably its well-tempered and replica-exchange variants, which further boost the sampling efficiency while gaining in numerical stability, thus allowing quantum-mechanical/molecular-mechanical free energy calculations to be performed at a lower cost. As a paradigm to bridge microsecond simulations to millisecond events by means of free energy calculations, we have applied the ABF family of algorithms to decompose complex movements in molecular objects of biological and abiological nature. We show here how water lubricates the shuttling of an amide-based rotaxane by altering the mechanism that underlies the concerted translation and isomerization of the macrocycle. Introducing novel collective variables in a computational workflow for the rigorous determination of standard binding free energies, we predict with utmost accuracy the thermodynamics of protein-ligand reversible association. Because of their simplicity, versatility, and robust mathematical foundations, the algorithms of the ABF family represent an appealing option for the theoretical investigation of a broad range of problems relevant to physics, chemistry, and biology.

Overcoming Free-Energy Barriers with a Seamless Combination of a Biasing Force and a Collective Variable-Independent Boost Potential.

The Extended Generalized Adaptive Biasing Force Algorithm for Multidimensional Free-Energy Calculations.

Taming Rugged Free Energy Landscapes Using an Average Force.

Combining the Biased and Unbiased Sampling Strategy into One Convenient Free Energy Calculation Method.

Zooming Across the Free-Energy Landscape: Shaving Barriers, and Flooding Valleys.

Extended Adaptive Biasing Force Algorithm. An On-the-Fly Implementation for Accurate Free-Energy Calculations.

Improving the Exploration of High Dimensional Free Energy Landscape by a Combination of Temperature Accelerated Sliced Sampling and Parallel Biasing

Free energy calculations: An efficient adaptive biasing potential method

Overcoming Sampling Issues and Improving Computational Efficiency in Collective-Variable-Based Enhanced-Sampling Simulations: A Tutorial

Finding an Optimal Pathway on a Multidimensional Free-Energy Landscape.

BAR-based Optimum Adaptive Steered MD for Configurational Sampling.

The Adaptive Biasing Force Method: Everything You Always Wanted To Know but Were Afraid To Ask

Multireference Generalization of the Weighted Thermodynamic Perturbation Method.

Avoiding Non-Equilibrium Effects in Adaptive Biasing Force Calculations

Unrestrained Computation of Free Energy along A Path

Enhanced Sampling of Multidimensional Free-Energy Landscapes Using Adaptive Biasing Forces

Enhanced Sampling with Sub-optimal Collective Variables: Reconciling Accuracy and Convergence Speed

Sampling free energy surfaces as slices by combining umbrella sampling and metadynamics

Exploration of high dimensional free energy landscapes by a combination of temperature‐accelerated sliced sampling and parallel biasing

Exploring Valleys without Climbing Every Peak: More Efficient and Forgiving Metabasin Metadynamics via Robust On-the-Fly Bias Domain Restriction.

Improved Sampling of Adaptive Path Collective Variables by Stabilized Extended-System Dynamics