Abstract:Conformational sampling of complex biomolecules is an emerging frontier in drug discovery. Indeed, advances in lab-based structural biology and related computational approaches like AlphaFold have made great strides in obtaining static protein structures. However, biology is in constant motion and many important biological processes rely on conformationally-driven events. Unrestrained molecular dynamics (MD) simulations require that the simulated time be comparable to the real time of the biological processes of interest, rendering pure MD impractical for many drug design projects, where conformationally-driven biological events can take microseconds to milliseconds or longer. An alternative approach is to accelerate the sampling of specific motions by applying restraints, guided by insights about the underlying biological process of interest. A plethora of restraints exist to limit the size of conformational search space, although each has drawbacks when simulating complex biological motions. In this work, we introduce a new kind of restraint for molecular dynamics simulations (MD) that is particularly well suited for complex conformationally-driven biological events, such as protein-ligand binding, allosteric modulations, conformational signalling, and membrane permeability. The new restraint, which relies on a barrier function (the scaled reciprocal function) is particularly beneficial to MD, where hard-wall restraints are needed with zero tolerance to restraint violation. We have implemented this restraint within a hybrid sampling framework that combines metadynamics and extended-Lagrangian adaptive biasing force (meta-eABF). We use two particular examples to demonstrate the value of this approach: (1) quantification of the approach of E3-loaded ubiquitin to a protein of interest as part of the Cullin ring ligase and (2) membrane permeability of heterobifunctional degrader molecules with a large degree of conformational flexibility. Future work will involve extension to additional systems and benchmarking of this approach compared with other methods.

Enhanced molecular dynamics sampling of drug target conformations

Sampling of conformational ensemble for virtual screening using molecular dynamics simulations and normal mode analysis

Steered Molecular Dynamics Simulation in Rational Drug Design

Predicting the bioactive conformations of macrocycles: a molecular dynamics-based docking procedure with DynaDock

Molecular docking as a popular tool in drug design, an in silico travel

Molecular Docking Improved With Human Spatial Perception Using Virtual Reality

Assessing small molecule conformational sampling methods in molecular docking

Molecular Docking: Principles, Advances, and its Applications in Drug Discovery

Advances in Molecular Dynamics Simulations and Enhanced Sampling Methods for the Study of Protein Systems

Docking Challenge: Protein Sampling and Molecular Docking Performance

Dynamic Docking: A Paradigm Shift in Computational Drug Discovery

Improving Protein-Ligand Docking Results with High-Throughput Molecular Dynamics Simulations

Dealing with induced fit, conformational selection and secondary poses in molecular dynamics simulations for reliable free energy predictions

Molecular Docking and Structure-Based Drug Design Strategies

Interactive molecular dynamics in virtual reality for accurate flexible protein-ligand docking

Enhancing Ligand Pose Sampling for Molecular Docking

Dynamic Docking of Macrocycles in Bound and Unbound Protein Structures with DynaDock

Molecular dynamics simulations and novel drug discovery

A Robust Induced Fit Docking Approach with the Combination of the Hybrid All-Atom/United-Atom/Coarse-Grained Model and Simulated Annealing

A Different Kind of Restraint Suitable for Molecular Dynamics Simulations

Roles of Accelerated Molecular Dynamics Simulations in Predictions of Binding Kinetic Parameters