Development of an Automatic Pipeline for Participation in the CELPP Challenge

The prediction of how a ligand binds to its target is an essential step for Structure-Based Drug Design (SBDD) methods. Molecular docking is a standard tool to predict the binding mode of a ligand to its macromolecular receptor and to quantify their mutual complementarity, with multiple applications in drug design. However, docking programs do not always find correct solutions, either because they are not sampled or due to inaccuracies in the scoring functions. Quantifying the docking performance in real scenarios is essential to understanding their limitations, managing expectations and guiding future developments. Here, we present a fully automated pipeline for pose prediction validated by participating in the Continuous Evaluation of Ligand Pose Prediction (CELPP) Challenge. Acknowledging the intrinsic limitations of the docking method, we devised a strategy to automatically mine and exploit pre-existing data, defining—whenever possible—empirical restraints to guide the docking process. We prove that the pipeline is able to generate predictions for most of the proposed targets as well as obtain poses with low RMSD values when compared to the crystal structure. All things considered, our pipeline highlights some major challenges in the automatic prediction of protein–ligand complexes, which will be addressed in future versions of the pipeline.

Extracting Atomic Contributions to Binding Free Energy Using Molecular Dynamics Simulations with Mixed Solvents (MDmix)

BACKGROUND

Mixed solvents MD (MDmix) simulations have proved to be a useful and increasingly accepted technique with several applications in structure-based drug discovery. One of the assumptions behind the methodology is the transferability of free energy values from the simulated cosolvent molecules to larger drug-like molecules. However, the binding free energy maps (ΔGbind) calculated for the different moieties of the cosolvent molecules (e.g. a hydroxyl map for the ethanol) are largely influenced by the rest of the solvent molecule and do not reflect the intrinsic affinity of the moiety in question. As such, they are hardly transferable to different molecules.

METHOD

To achieve transferable energies, we present here a method for decomposing the molecular binding free energy into accurate atomic contributions.

RESULT

We demonstrate with two qualitative visual examples how the corrected energy maps better match known binding hotspots and how they can reveal hidden hotspots with actual drug design potential.

CONCLUSION

Atomic decomposition of binding free energies derived from MDmix simulations provides transferable and quantitative binding free energy maps.