Treating electrostatics with Wolf summation in combined quantum mechanical and molecular mechanical simulations

PyDFT-QMMM: A modular, extensible software framework for DFT-based QM/MM molecular dynamics

PyDFT-QMMM is a Python-based package for performing hybrid quantum mechanics/molecular mechanics (QM/MM) simulations at the density functional level of theory. The program is designed to treat short-range and long-range interactions through user-specified combinations of electrostatic and mechanical embedding procedures within periodic simulation domains, providing necessary interfaces to external quantum chemistry and molecular dynamics software. To enable direct embedding of long-range electrostatics in periodic systems, we have derived and implemented force terms for our previously described QM/MM/PME approach [Pederson and McDaniel, J. Chem. Phys. 156, 174105 (2022)]. Communication with external software packages Psi4 and OpenMM is facilitated through Python application programming interfaces (APIs). The core library contains basic utilities for running QM/MM molecular dynamics simulations, and plug-in entry-points are provided for users to implement custom energy/force calculation and integration routines, within an extensible architecture. The user interacts with PyDFT-QMMM primarily through its Python API, allowing for complex workflow development with Python scripting, for example, interfacing with PLUMED for free energy simulations. We provide benchmarks of forces and energy conservation for the QM/MM/PME and alternative QM/MM electrostatic embedding approaches. We further demonstrate a simple example use case for water solute in a water solvent system, for which radial distribution functions are computed from 100 ps QM/MM simulations; in this example, we highlight how the solvation structure is sensitive to different basis-set choices due to under- or over-polarization of the QM water molecule's electron density.

Mapping Free Energy Pathways for ATP Hydrolysis in the E. coli ABC Transporter HlyB by the String Method

HlyB functions as an adenosine triphosphate (ATP)-binding cassette (ABC) transporter that enables bacteria to secrete toxins at the expense of ATP hydrolysis. Our previous work, based on potential energy profiles from combined quantum mechanical and molecular mechanical (QM/MM) calculations, has suggested that the highly conserved H-loop His residue H662 in the nucleotide binding domain (NBD) of E. coli HlyB may catalyze the hydrolysis of ATP through proton relay. To further test this hypothesis when entropic contributions are taken into account, we obtained QM/MM minimum free energy paths (MFEPs) for the HlyB reaction, making use of the string method in collective variables. The free energy profiles along the MFEPs confirm the direct participation of H662 in catalysis. The MFEP simulations of HlyB also reveal an intimate coupling between the chemical steps and a local protein conformational change involving the signature-loop residue S607, which may serve a catalytic role similar to an Arg-finger motif in many ATPases and GTPases in stabilizing the phosphoryl-transfer transition state.

Representation of the QM Subsystem for Long-Range Electrostatic Interaction in Non-Periodic Ab Initio QM/MM Calculations

In QM/MM calculations, it is essential to handle electrostatic interactions between the QM and MM subsystems accurately and efficiently. To achieve maximal efficiency, it is convenient to adopt a hybrid scheme, where the QM electron density is used explicitly in the evaluation of short-range QM/MM electrostatic interactions, while a multipolar representation for the QM electron density is employed to account for the long-range QM/MM electrostatic interactions. In order to avoid energy discontinuity at the cutoffs, which separate the short- and long-range QM/MM electrostatic interactions, a switching function should be utilized to ensure a smooth potential energy surface. In this study, we benchmarked the accuracy of such hybrid embedding schemes for QM/MM electrostatic interactions using different multipolar representations, switching functions and cutoff distances. For test systems (neutral and anionic oxyluciferin in MM (aqueous and enzyme) environments), the best accuracy was acquired with a combination of QM electrostatic potential (ESP) charges and dipoles and two switching functions (long-range electrostatic corrections (LREC) and Switch) in the treatment of long-range QM/MM electrostatics. It allowed us to apply a 10Å distance cutoff and still obtain QM/MM electrostatics/polarization energies within 0.1 kcal/mol and time-dependent density functional theory (TDDFT)/MM vertical excitation energies within 10-3 eV from theoretical reference values.

Performance evaluation of the zero-multipole summation method in modern molecular dynamics software

The zero-multiple summation method (ZMM) is a cutoff-based method for calculating electrostatic interactions in molecular dynamics simulations, utilizing an electrostatic neutralization principle as a physical basis. Since the accuracies of the ZMM have been revealed to be sufficient in previous studies, it is highly desirable to clarify its practical performance. In this paper, the performance of the ZMM is compared with that of the smooth particle mesh Ewald method (SPME), where the both methods are implemented in molecular dynamics software package GROMACS. Extensive performance comparisons against a highly optimized, parameter-tuned SPME implementation are performed for various-sized water systems and two protein-water systems. We analyze in detail the dependence of the performance on the potential parameters and the number of CPU cores. Even though the ZMM uses a larger cutoff distance than the SPME does, the performance of the ZMM is comparable to or better than that of the SPME. This is because the ZMM does not require a time-consuming electrostatic convolution and because the ZMM gains short neighbor-list distances due to the smooth damping feature of the pairwise potential function near the cutoff length. We found, in particular, that the ZMM with quadrupole or octupole cancellation and no damping factor is an excellent candidate for the fast calculation of electrostatic interactions. © 2018 Wiley Periodicals, Inc.

Long-range electrostatic corrections in multipolar/polarizable QM/MM simulations

Taking long-range electrostatic effects into account in classical and hybrid quantum mechanics-molecular mechanics (QM/MM) simulations is necessary for an accurate description of the system under study. We have recently developed a method, termed long-range electrostatic corrections (LREC), for monopolar QM/MM calculations. Here, we present an extension of LREC for multipolar/polarizable QM/MM simulations within the LICHEM software package. Reaction barriers and QM-MM interaction energies converge with a LREC cutoff between 20 and 25 Å, in agreement with our previous results. Additionally, the LREC approach for the QM-MM interactions can be smoothly combined with standard shifting or Ewald summation methods in the MM calculations. We recommend the use of QM(LREC)/MM(PME), where the QM region is treated with LREC and the MM region is treated with particle mesh Ewald (PME) summation. This combination is an excellent compromise between simplicity, speed, and accuracy for large QM/MM simulations.