We present a new concept of free energy calculations of chemical reactions by means of extended sampling molecular dynamics simulations. Biasing potentials are applied on partial atomic charges, which may be combined with atomic coordinates either in a single collective variable or in multi-dimensional biasing simulations. The necessary additional gradients are obtained by solving coupled-perturbed equations within the approximative density-functional tight-binding method. The new computational scheme was implemented in a combination of Gromacs and Plumed. As a prospective application, proton-coupled electron transfer in a model molecular system is studied. Two collective variables are introduced naturally, one for the proton transfer and the other for the electron transfer. The results are in qualitative agreement with the extended free simulations performed for reference. Free energy minima as well as the mechanism of the process are identified correctly, while the topology of the transition region and the height of the energy barrier are only reproduced qualitatively. The application also illustrates possible difficulties with the new methodology. These may be inefficient sampling of spatial coordinates when atomic charges are biased exclusively and a decreased stability of the simulations. Still, the new approach represents a viable alternative for free energy calculations of a certain class of chemical reactions, for instance a proton-coupled electron transfer in proteins.
Skip Nav Destination
Article navigation
21 August 2018
Research Article|
June 27 2018
Coupled-perturbed DFTB-QM/MM metadynamics: Application to proton-coupled electron transfer
Special Collection:
Enhanced Sampling for Molecular Systems
Natacha Gillet;
Natacha Gillet
1
Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT)
, 76131 Karlsruhe, Germany
Search for other works by this author on:
Marcus Elstner;
Marcus Elstner
1
Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT)
, 76131 Karlsruhe, Germany
2
Institute for Biological Interfaces (IBG-2), KIT
, 76131 Karlsruhe, Germany
Search for other works by this author on:
Tomáš Kubař
Tomáš Kubař
a)
1
Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT)
, 76131 Karlsruhe, Germany
3
Center for Functional Nanostructures, KIT
, 76131 Karlsruhe, Germany
Search for other works by this author on:
a)
Electronic mail: tomas.kubar@kit.edu
J. Chem. Phys. 149, 072328 (2018)
Article history
Received:
February 27 2018
Accepted:
June 07 2018
Citation
Natacha Gillet, Marcus Elstner, Tomáš Kubař; Coupled-perturbed DFTB-QM/MM metadynamics: Application to proton-coupled electron transfer. J. Chem. Phys. 21 August 2018; 149 (7): 072328. https://doi.org/10.1063/1.5027100
Download citation file:
Sign in
Don't already have an account? Register
Sign In
You could not be signed in. Please check your credentials and make sure you have an active account and try again.
Pay-Per-View Access
$40.00
Citing articles via
DeePMD-kit v2: A software package for deep potential models
Jinzhe Zeng, Duo Zhang, et al.