Journal
JOURNAL OF CHEMICAL THEORY AND COMPUTATION
Volume 14, Issue 6, Pages 3279-3288Publisher
AMER CHEMICAL SOC
DOI: 10.1021/acs.jctc.8b00031
Keywords
-
Funding
- National Institutes of Health [GM32136]
- Yale University Faculty of Arts and Sciences High Performance Computing Center
- NATIONAL INSTITUTE OF GENERAL MEDICAL SCIENCES [R01GM032136] Funding Source: NIH RePORTER
Ask authors/readers for more resources
The generation of a complete ensemble of geometrical configurations is required to obtain reliable estimations of absolute binding free energies by alchemical free energy methods. Molecular dynamics (MD) is the most popular sampling method, but the representation of large biomolecular systems may be incomplete owing to energetic barriers that impede efficient sampling of the configurational space. Monte Carlo (MC) methods can possibly overcome this issue by adapting the attempted movement sizes to facilitate transitions between alternative local-energy minima. In this study, we present an MC statistical mechanics algorithm to explore the protein-ligand conformational space with emphasis on the motions of the protein backbone and side chains. The parameters for each MC move type were optimized to better reproduce conformational distributions of 18 dipeptides and the well-studied T4-lysozyme L99A protein. Next, the performance of the improved MC algorithms was evaluated by computing absolute free energies of binding for L99A lysozyme with benzene and seven analogs. Results for benzene with L99A lysozyme from MD and the optimized MC protocol were found to agree within 0.6 kcal/mol, while a mean unsigned error of 1.2 kcal/mol between MC results and experiment was obtained for the seven benzene analogs. Significant advantages in computation speed are also reported with MC over MD for similar extents of configurational sampling.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available