title: Periodic Boundary Conditions in QM/MM Calculations: Implementation and Tests.
abstract: Quantum mechanics/molecular mechanics (QM/MM) simulations of reactions in solutions and in solvated enzymes can be performed using the QM/MM-Ewald approach with periodic boundary conditions (PBC) or a nonperiodic treatment with a finite solvent shell (droplet model). To avoid the changes in QM codes that are required in standard QM/MM-Ewald implementations, we present a general method (Gen-Ew) for periodic QM/MM calculations that can be used with any QM method in the QM/MM framework. The Gen-Ew approach approximates the QM/MM-Ewald method by representing the PBC potential by virtual charges on a sphere and the QM density by electrostatic potential (ESP) charges. Test calculations show that the deviations between Gen-Ew and QM/MM-Ewald results are generally small enough to justify the application of the Gen-Ew method in the absence of a suitable QM/MM-Ewald implementation. We compare the results from periodic QM/MM calculations (QM/MM-Ewald, Gen-Ew) to their nonperiodic counterparts (droplet model) for five test reactions in water and for the Claisen rearrangement in chorismate mutase. The periodic and nonperiodic QM/MM treatments give similar free energy profiles for the reactions in solution (umbrella sampling, free energy deviations of the order of 1 kcal/mol) and essentially the same energy profile (constrained geometry optimizations) for the Claisen rearrangement in chorismate mutase. In all cases considered, long-range electrostatic interactions are thus well captured by nonperiodic QM/MM calculations in a water droplet of reasonable size (radius of 15-20 Å). This provides further justification for the widespread use of the computationally efficient droplet model in QM/MM studies of reactions in solution and in enzymes.
authors.affiliation: Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany.
authors.forename: Tatiana
authors.initials: T
authors.lastname: Vasilevskaya
authors.affiliation: Max-Planck-Institut für Kohlenforschung, 45470 Mülheim an der Ruhr, Germany.
authors.forename: Walter
authors.initials: W
authors.lastname: Thiel
mesh_terms: D000643:Ammonium Chloride / Q000737:chemistry; D008958:Models, Molecular; D008968:Molecular Conformation; D056004:Molecular Dynamics Simulation; D010710:Phosphates / Q000737:chemistry; D011789:Quantum Theory; D055672:Static Electricity; D013816:Thermodynamics; D014867:Water / Q000737:chemistry
chemical_list: D010710:Phosphates; C024788:ammonium phosphate; D000643:Ammonium Chloride; D014867:Water