4.7 Article

Generation of Quantum Configurational Ensembles Using Approximate Potentials

Journal

JOURNAL OF CHEMICAL THEORY AND COMPUTATION
Volume 17, Issue 11, Pages 7021-7042

Publisher

AMER CHEMICAL SOC
DOI: 10.1021/acs.jctc.1c00532

Keywords

-

Funding

  1. EPSRC [EP/P020194/1, EP/T022213/1]
  2. AstraZeneca
  3. EPSRC Centre for Doctoral Training, Theory and Modelling in Chemical Sciences [EP/L015722/1]

Ask authors/readers for more resources

Conformational analysis is crucial in drug design, and molecular mechanics simulation methods are used to generate ensembles of structures to provide reliable structural information. Reparameterizing the force field can generate FFs that closely reproduce QM results, and the MC acceptance rate is strongly correlated with various phase space overlap measurements, serving as a robust metric to evaluate the similarity between MM and QM levels of theory.
Conformational analysis is of paramount importance in drug design: it is crucial to determine pharmacological properties, understand molecular recognition processes, and characterize the conformations of ligands when unbound. Molecular Mechanics (MM) simulation methods, such as Monte Carlo (MC) and molecular dynamics (MD), are usually employed to generate ensembles of structures due to their ability to extensively sample the conformational space of molecules. The accuracy of these MM-based schemes strongly depends on the functional form of the force field (FF) and its parametrization, components that often hinder their performance. High-level methods, such as ab initio MD, provide reliable structural information but are still too computationally expensive to allow for extensive sampling. Therefore, to overcome these limitations, we present a multilevel MC method that is capable of generating quantum configurational ensembles while keeping the computational cost at a minimum. We show that FF reparametrization is an efficient route to generate FFs that reproduce QM results more closely, which, in turn, can be used as low-cost models to achieve the gold standard QM accuracy. We demonstrate that the MC acceptance rate is strongly correlated with various phase space overlap measurements and that it constitutes a robust metric to evaluate the similarity between the MM and QM levels of theory. As a more advanced application, we present a self-parametrizing version of the algorithm, which combines sampling and FF parametrization in one scheme, and apply the methodology to generate the QM/MM distribution of a ligand in aqueous solution.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.7
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

No Data Available
No Data Available