MIT Libraries logoDSpace@MIT

MIT
View Item 
  • DSpace@MIT Home
  • MIT Open Access Articles
  • MIT Open Access Articles
  • View Item
  • DSpace@MIT Home
  • MIT Open Access Articles
  • MIT Open Access Articles
  • View Item
JavaScript is disabled for your browser. Some features of this site may not work without it.

Systematic Quantum Mechanical Region Determination in QM/MM Simulation

Author(s)
Karelina, Maria; Kulik, Heather Janine
Thumbnail
DownloadKarelina_OA.pdf (2.723Mb)
PUBLISHER_POLICY

Publisher Policy

Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use.

Terms of use
Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use.
Metadata
Show full item record
Abstract
Hybrid quantum mechanical-molecular mechanical (QM/MM) simulations are widely used in enzyme simulation. Over ten convergence studies of QM/MM methods have revealed over the past several years that key energetic and structural properties approach asymptotic limits with only very large (ca. 500–1000 atom) QM regions. This slow convergence has been observed to be due in part to significant charge transfer between the core active site and the surrounding protein environment, which cannot be addressed by improvement of MM force fields or the embedding method employed within QM/MM. Given this slow convergence, it becomes essential to identify strategies for the most atom-economical determination of optimal QM regions and to gain insight into the crucial interactions captured only in large QM regions. Here, we extend and develop two methods for quantitative determination of QM regions. First, in the charge shift analysis (CSA) method, we probe the reorganization of electron density when core active site residues are removed completely, as determined by large-QM region QM/MM calculations. Second, we introduce the highly parallelizable Fukui shift analysis (FSA), which identifies how core/substrate frontier states are altered by the presence of an additional QM residue in smaller initial QM regions. We demonstrate that the FSA and CSA approaches are complementary and consistent on three test case enzymes: catechol O-methyltransferase, cytochrome P450cam, and hen eggwhite lysozyme. We also introduce validation strategies and test the sensitivities of the two methods to geometric structure, basis set size, and electronic structure methodology. Both methods represent promising approaches for the systematic, unbiased determination of quantum mechanical effects in enzymes and large systems that necessitate multiscale modeling.
Date issued
2017-01
URI
http://hdl.handle.net/1721.1/113023
Department
Massachusetts Institute of Technology. Department of Chemical Engineering; Massachusetts Institute of Technology. Department of Electrical Engineering and Computer Science
Journal
Journal of Chemical Theory and Computation
Publisher
American Chemical Society (ACS)
Citation
Karelina, Maria, and Kulik, Heather J. “Systematic Quantum Mechanical Region Determination in QM/MM Simulation.” Journal of Chemical Theory and Computation 13, 2 (January 2017): 563–576 © 2017 American Chemical Society
Version: Author's final manuscript
ISSN
1549-9618
1549-9626

Collections
  • MIT Open Access Articles

Browse

All of DSpaceCommunities & CollectionsBy Issue DateAuthorsTitlesSubjectsThis CollectionBy Issue DateAuthorsTitlesSubjects

My Account

Login

Statistics

OA StatisticsStatistics by CountryStatistics by Department
MIT Libraries
PrivacyPermissionsAccessibilityContact us
MIT
Content created by the MIT Libraries, CC BY-NC unless otherwise noted. Notify us about copyright concerns.