MovableType Software for Fast Free Energy-Based Virtual Screening: Protocol Development, Deployment, Validation, and Assessment.

For decades, the complicated energy surfaces found in macromolecular protein:ligand structures, which require large amounts of computational time and resources for energy state sampling, have been an inherent obstacle to fast, routine free energy estimation in industrial drug discovery efforts. Beginning in 2013, the Merz research group addressed this cost with the introduction of a novel sampling methodology termed "Movable Type" (MT). Using numerical integration methods, the MT method reduces the computational expense for energy state sampling by independently calculating each atomic partition function from an initial molecular conformation in order to estimate the molecular free energy using ensembles of the atomic partition functions. In this work, we report a software package, the DivCon Discovery Suite with the MovableType module from QuantumBio Inc., that performs this MT free energy estimation protocol in a fast, fully encapsulated manner. We discuss the computational procedures and improvements to the original work, and we detail the corresponding settings for this software package. Finally, we introduce two validation benchmarks to evaluate the overall robustness of the method against a broad range of protein:ligand structural cases. With these publicly available benchmarks, we show that the method can use a variety of input types and parameters and exhibits comparable predictability whether the method is presented with "expensive" X-ray structures or "inexpensively docked" theoretical models. We also explore some next steps for the method. The MovableType software is available at http://www.quantumbioinc.com/.

Receptor-Ligand Binding Free Energies from a Consecutive Histograms Monte Carlo Sampling Method.

To obtain accurate and converged free energy calculations for ligand binding to biomolecular systems requires validated force fields and extensive sampling of the energy landscape, which requires exhaustive and effective conformational searching methods. Herein, we introduce the consecutive histograms Monte Carlo (CHMC) sampling protocol that generates receptor-ligand binding modes within a series of continuously distributed sampling units ranging from placement near the geometric center of the receptor's binding site to fully unbound states. This protocol employs independent energy-state sampling for calculating the ensemble energy within every predefined location along the receptor-ligand dissociation pathway, without the need to traverse the energy barriers as in molecular dynamic simulations during the dissociation procedure. We applied this method to a set of selected receptor targets with their corresponding ligands providing detailed studies of molecular binding free energy predictions. The results show that the CHMC gives an excellent accounting of the free energy surfaces and binding free energies at a reasonable computational cost.