A detailed description of the events ruling ligand/protein interaction and an accurate estimation of the drug affinity to its target is of great help in speeding drug discovery strategies. We have developed a metadynamics-based approach, named funnel metadynamics, that allows the ligand to enhance the sampling of the target binding sites and its solvated states. This method leads to an efficient characterization of the binding free-energy surface and an accurate calculation of the absolute protein-ligand binding free energy. We illustrate our protocol in two systems, benzamidine/ trypsin and SC-558/cyclooxygenase 2. In both cases, the X-ray conformation has been found as the lowest free-energy pose, and the computed protein-ligand binding free energy in good agreement with experiments. Furthermore, funnel metadynamics unveils important information about the binding process, such as the presence of alternative binding modes and the role of waters. The results achieved at an affordable computational cost make funnel metadynamics a valuable method for drug discovery and for dealing with a variety of problems in chemistry, physics, and material science.
Studying the molecular interactions between a drug and its target helps in understanding the target functional mechanism and offers the possibility for exogenous control of its physiological activity. In recent years, a vast experimental and computational effort has revealed in ever-more-precise detail the ligand/target recognition mechanism (1, 2). In this context, an accurate estimation of the ligand-binding affinity is in great demand because it would facilitate many steps of the drug discovery pipeline, such as structure-based drug design and lead optimization; this is not, however, a simple task. In fact, an accurate estimation of the binding affinity or, equivalently, the absolute protein-ligand binding free energy, requires an accurate description of the ligand/ protein interactions, their flexibility, and the solvation process. Many methods have been proposed to tackle this problem. For instance, docking protocols are widely used to generate and rank candidate poses based on empirical scoring functions, either physically or statistically based (3-5). These techniques have been proven to be highly efficient in screening a large number of compounds in a short time (6); this, however, at the price of limited accuracy in estimating affinities (7).Alternatively, a variety of methods to describe ligand/protein interactions in a more accurate way at higher computational cost have been proposed. These techniques can be grouped in two categories: (i) endpoint and (ii) pathway methods. The former group is composed of those techniques that sample ligand and protein in unbound and bound states and compute the proteinligand binding free energy by taking the difference between the absolute free energy of these two states. Examples include microscopic linear response approximation (8), linear interaction energy (9, 10), protein dipoles Langevin dipoles (11), as well as mol...