Abstract
The use of small organic molecules has become one of the most popular strategies in computer-aided drug design (CADD) to facilitate the identification of potential drug-like compounds in the early stages of drug development. In this scenario, novel computational approaches such as the use of the fragment dissolved molecular dynamics (fdMD) methodology emerged as a new framework for the modeling of ligand-receptor interactions. Consisting of molecular dynamics (MD) simulations of the target protein solvated with multiple copies of the same fragment, the original approach is able to identify the most favorable binding site for the system studied in a reasonable simulation time scale (0.2-1 μs). In the present work, we have introduced the use of Gaussian accelerated molecular dynamics (GaMD) to facilitate system exploration, accelerate binding site identification and additionally enhance binding mode elucidation. For this purpose, up to 12 different systems with crystallographic information available have been employed for validation.