Solvent Boundary Potentials for Hybrid QM/MM Computations Using Classical Drude Oscillators: A Fully Polarizable Model
Accurate quantum mechanical/molecular mechanical (QM/MM) treatments should account for MM polarization and properly include long-range electrostatic interactions. We report on a development that covers both these aspects. Our approach combines the classical Drude oscillator (DO) model for the electr...
Saved in:
| Published in: | Journal of chemical theory and computation Vol. 8; no. 11; p. 4527 |
|---|---|
| Main Authors: | , |
| Format: | Journal Article |
| Language: | English |
| Published: |
United States
13.11.2012
|
| ISSN: | 1549-9618 |
| Online Access: | Get more information |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Abstract | Accurate quantum mechanical/molecular mechanical (QM/MM) treatments should account for MM polarization and properly include long-range electrostatic interactions. We report on a development that covers both these aspects. Our approach combines the classical Drude oscillator (DO) model for the electronic polarizability of the MM atoms with the generalized solvent boundary Potential (GSBP) and the solvated macromolecule boundary potential (SMBP). These boundary potentials (BP) are designed to capture the long-range effects of the outer region of a large system on its interior. They employ a finite difference approximation to the Poisson-Boltzmann equation for computing electrostatic interactions and take into account outer-region bulk solvent through a polarizable dielectric continuum (PDC). This approach thus leads to fully polarizable three-layer QM/MM-DO/BP methods. As the mutual responses of each of the subsystems have to be taken into account, we propose efficient schemes to converge the polarization of each layer simultaneously. For molecular dynamics (MD) simulations using GSBP, this is achieved by considering the MM polarizable model as a dynamical degree of freedom, and hence contributions from the boundary potential can be evaluated for a frozen state of polarization at every time step. For geometry optimizations using SMBP, we propose a dual self-consistent field approach for relaxing the Drude oscillators to their ideal positions and converging the QM wave function with the proper boundary potential. The chosen coupling schemes are evaluated with a test system consisting of a glycine molecule in a water ball. Both boundary potentials are capable of properly reproducing the gradients at the inner-region atoms and the Drude oscillators. We show that the effect of the Drude oscillators must be included in all terms of the boundary potentials to obtain accurate results and that the use of a high dielectric constant for the PDC does not lead to a polarization catastrophe of the DO models. Optimum values for some key parameters are discussed. We also address the efficiency of these approaches compared to standard QM/MM-DO calculations without BP. In the SMBP case, computation times can be reduced by around 40% for each step of a geometry optimization, with some variation depending on the chosen QM method. In the GSBP case, the computational advantages of using the boundary potential increase with system size and with the number of MD steps. |
|---|---|
| AbstractList | Accurate quantum mechanical/molecular mechanical (QM/MM) treatments should account for MM polarization and properly include long-range electrostatic interactions. We report on a development that covers both these aspects. Our approach combines the classical Drude oscillator (DO) model for the electronic polarizability of the MM atoms with the generalized solvent boundary Potential (GSBP) and the solvated macromolecule boundary potential (SMBP). These boundary potentials (BP) are designed to capture the long-range effects of the outer region of a large system on its interior. They employ a finite difference approximation to the Poisson-Boltzmann equation for computing electrostatic interactions and take into account outer-region bulk solvent through a polarizable dielectric continuum (PDC). This approach thus leads to fully polarizable three-layer QM/MM-DO/BP methods. As the mutual responses of each of the subsystems have to be taken into account, we propose efficient schemes to converge the polarization of each layer simultaneously. For molecular dynamics (MD) simulations using GSBP, this is achieved by considering the MM polarizable model as a dynamical degree of freedom, and hence contributions from the boundary potential can be evaluated for a frozen state of polarization at every time step. For geometry optimizations using SMBP, we propose a dual self-consistent field approach for relaxing the Drude oscillators to their ideal positions and converging the QM wave function with the proper boundary potential. The chosen coupling schemes are evaluated with a test system consisting of a glycine molecule in a water ball. Both boundary potentials are capable of properly reproducing the gradients at the inner-region atoms and the Drude oscillators. We show that the effect of the Drude oscillators must be included in all terms of the boundary potentials to obtain accurate results and that the use of a high dielectric constant for the PDC does not lead to a polarization catastrophe of the DO models. Optimum values for some key parameters are discussed. We also address the efficiency of these approaches compared to standard QM/MM-DO calculations without BP. In the SMBP case, computation times can be reduced by around 40% for each step of a geometry optimization, with some variation depending on the chosen QM method. In the GSBP case, the computational advantages of using the boundary potential increase with system size and with the number of MD steps. Accurate quantum mechanical/molecular mechanical (QM/MM) treatments should account for MM polarization and properly include long-range electrostatic interactions. We report on a development that covers both these aspects. Our approach combines the classical Drude oscillator (DO) model for the electronic polarizability of the MM atoms with the generalized solvent boundary Potential (GSBP) and the solvated macromolecule boundary potential (SMBP). These boundary potentials (BP) are designed to capture the long-range effects of the outer region of a large system on its interior. They employ a finite difference approximation to the Poisson-Boltzmann equation for computing electrostatic interactions and take into account outer-region bulk solvent through a polarizable dielectric continuum (PDC). This approach thus leads to fully polarizable three-layer QM/MM-DO/BP methods. As the mutual responses of each of the subsystems have to be taken into account, we propose efficient schemes to converge the polarization of each layer simultaneously. For molecular dynamics (MD) simulations using GSBP, this is achieved by considering the MM polarizable model as a dynamical degree of freedom, and hence contributions from the boundary potential can be evaluated for a frozen state of polarization at every time step. For geometry optimizations using SMBP, we propose a dual self-consistent field approach for relaxing the Drude oscillators to their ideal positions and converging the QM wave function with the proper boundary potential. The chosen coupling schemes are evaluated with a test system consisting of a glycine molecule in a water ball. Both boundary potentials are capable of properly reproducing the gradients at the inner-region atoms and the Drude oscillators. We show that the effect of the Drude oscillators must be included in all terms of the boundary potentials to obtain accurate results and that the use of a high dielectric constant for the PDC does not lead to a polarization catastrophe of the DO models. Optimum values for some key parameters are discussed. We also address the efficiency of these approaches compared to standard QM/MM-DO calculations without BP. In the SMBP case, computation times can be reduced by around 40% for each step of a geometry optimization, with some variation depending on the chosen QM method. In the GSBP case, the computational advantages of using the boundary potential increase with system size and with the number of MD steps.Accurate quantum mechanical/molecular mechanical (QM/MM) treatments should account for MM polarization and properly include long-range electrostatic interactions. We report on a development that covers both these aspects. Our approach combines the classical Drude oscillator (DO) model for the electronic polarizability of the MM atoms with the generalized solvent boundary Potential (GSBP) and the solvated macromolecule boundary potential (SMBP). These boundary potentials (BP) are designed to capture the long-range effects of the outer region of a large system on its interior. They employ a finite difference approximation to the Poisson-Boltzmann equation for computing electrostatic interactions and take into account outer-region bulk solvent through a polarizable dielectric continuum (PDC). This approach thus leads to fully polarizable three-layer QM/MM-DO/BP methods. As the mutual responses of each of the subsystems have to be taken into account, we propose efficient schemes to converge the polarization of each layer simultaneously. For molecular dynamics (MD) simulations using GSBP, this is achieved by considering the MM polarizable model as a dynamical degree of freedom, and hence contributions from the boundary potential can be evaluated for a frozen state of polarization at every time step. For geometry optimizations using SMBP, we propose a dual self-consistent field approach for relaxing the Drude oscillators to their ideal positions and converging the QM wave function with the proper boundary potential. The chosen coupling schemes are evaluated with a test system consisting of a glycine molecule in a water ball. Both boundary potentials are capable of properly reproducing the gradients at the inner-region atoms and the Drude oscillators. We show that the effect of the Drude oscillators must be included in all terms of the boundary potentials to obtain accurate results and that the use of a high dielectric constant for the PDC does not lead to a polarization catastrophe of the DO models. Optimum values for some key parameters are discussed. We also address the efficiency of these approaches compared to standard QM/MM-DO calculations without BP. In the SMBP case, computation times can be reduced by around 40% for each step of a geometry optimization, with some variation depending on the chosen QM method. In the GSBP case, the computational advantages of using the boundary potential increase with system size and with the number of MD steps. |
| Author | Boulanger, Eliot Thiel, Walter |
| Author_xml | – sequence: 1 givenname: Eliot surname: Boulanger fullname: Boulanger, Eliot organization: Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany – sequence: 2 givenname: Walter surname: Thiel fullname: Thiel, Walter organization: Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany |
| BackLink | https://www.ncbi.nlm.nih.gov/pubmed/26605612$$D View this record in MEDLINE/PubMed |
| BookMark | eNo1kMlOwzAYhH0oogsceAHkI5dQL1kcbqVQitSoIMo5sp0_yMiJS-xUKk9PEOU0Gmn0zWimaNS6FhC6ouSWEkbnOnBCMsZghCY0ifMoT6kYo6n3n4RwHjN-jsYsTUmSUjZBhzdnD9AGfO_6tpLdEb-4MHgjrce16_D6qDpT4ddiXhR46Zp9H2QwrvX43Zv2Ay-t9N5oafFD11eAt14ba2Vwnb_DC7zqrf1lWtmZb6ks4MJVYC_QWT00wOVJZ2i3etwt19Fm-_S8XGwiyWkSIgaqTlRCgSoOItNES-BQSZJLktKMC8VILmrJMwEVo4wLohkRXAidKxmzGbr5w-4799WDD2VjvIZhXwuu9yXN4jwllMZ0iF6for1qoCr3nWmGO8r_q9gPzIZqng |
| CitedBy_id | crossref_primary_10_1002_qua_24889 crossref_primary_10_1039_C8CP00487K crossref_primary_10_3390_molecules23112830 crossref_primary_10_1002_ange_201402118 crossref_primary_10_1016_j_saa_2023_123485 crossref_primary_10_3390_molecules29184378 crossref_primary_10_1002_qua_25711 crossref_primary_10_1002_anie_201403217 crossref_primary_10_1016_j_jmb_2023_168358 crossref_primary_10_1016_j_molliq_2023_122861 crossref_primary_10_1016_j_sbi_2018_01_003 crossref_primary_10_3389_fchem_2018_00275 crossref_primary_10_1016_j_comptc_2014_03_022 crossref_primary_10_1016_j_jmgm_2014_09_003 crossref_primary_10_1063_5_0086275 crossref_primary_10_1080_08927022_2014_911870 crossref_primary_10_1002_jcc_24295 crossref_primary_10_1088_0953_8984_29_1_013002 crossref_primary_10_1063_5_0203247 crossref_primary_10_1002_wcms_1163 crossref_primary_10_1002_qua_25719 crossref_primary_10_1039_C5CP00986C crossref_primary_10_1146_annurev_physchem_040412_110040 crossref_primary_10_3390_nano11020299 crossref_primary_10_1146_annurev_physchem_090419_104031 crossref_primary_10_1002_jcc_27550 crossref_primary_10_1002_qua_25684 crossref_primary_10_1002_anie_201402118 crossref_primary_10_1002_cbic_201900363 crossref_primary_10_1039_c3cp44417a crossref_primary_10_1002_ange_201403217 crossref_primary_10_1007_s00214_018_2247_7 crossref_primary_10_3390_molecules25245853 crossref_primary_10_1002_wcms_1674 crossref_primary_10_1002_jcc_27018 crossref_primary_10_1016_j_dyepig_2023_111227 crossref_primary_10_1002_wcms_1731 crossref_primary_10_1002_wcms_1515 crossref_primary_10_1063_5_0064075 crossref_primary_10_1038_s41570_019_0092_4 crossref_primary_10_1007_s00214_018_2220_5 crossref_primary_10_1002_jcc_23716 crossref_primary_10_1039_C5CP00905G |
| ContentType | Journal Article |
| DBID | NPM 7X8 |
| DOI | 10.1021/ct300722e |
| DatabaseName | PubMed MEDLINE - Academic |
| DatabaseTitle | PubMed MEDLINE - Academic |
| DatabaseTitleList | PubMed MEDLINE - Academic |
| Database_xml | – sequence: 1 dbid: NPM name: PubMed url: http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed sourceTypes: Index Database – sequence: 2 dbid: 7X8 name: MEDLINE - Academic url: https://search.proquest.com/medline sourceTypes: Aggregation Database |
| DeliveryMethod | no_fulltext_linktorsrc |
| Discipline | Chemistry |
| ExternalDocumentID | 26605612 |
| Genre | Journal Article |
| GroupedDBID | 4.4 53G 55A 5GY 5VS 7~N AABXI ABJNI ABMVS ABQRX ABUCX ACGFS ACIWK ACS ADHLV AEESW AENEX AFEFF AHGAQ ALMA_UNASSIGNED_HOLDINGS AQSVZ BAANH CS3 CUPRZ D0L DU5 EBS ED~ EJD F5P GGK GNL IH9 J9A JG~ LG6 NPM P2P RNS ROL UI2 VF5 VG9 W1F 7X8 ABBLG ABLBI |
| ID | FETCH-LOGICAL-a315t-2ebf5b51e1b3e87c0cae3eda09a061738b2098fa378ed212380c208388c9ba42 |
| IEDL.DBID | 7X8 |
| ISICitedReferencesCount | 117 |
| ISICitedReferencesURI | http://www.webofscience.com/api/gateway?GWVersion=2&SrcApp=Summon&SrcAuth=ProQuest&DestLinkType=CitingArticles&DestApp=WOS_CPL&KeyUT=000311191900054&url=https%3A%2F%2Fcvtisr.summon.serialssolutions.com%2F%23%21%2Fsearch%3Fho%3Df%26include.ft.matches%3Dt%26l%3Dnull%26q%3D |
| ISSN | 1549-9618 |
| IngestDate | Fri Jul 11 16:25:06 EDT 2025 Thu Jan 02 22:26:02 EST 2025 |
| IsPeerReviewed | true |
| IsScholarly | true |
| Issue | 11 |
| Language | English |
| LinkModel | DirectLink |
| MergedId | FETCHMERGED-LOGICAL-a315t-2ebf5b51e1b3e87c0cae3eda09a061738b2098fa378ed212380c208388c9ba42 |
| Notes | ObjectType-Article-1 SourceType-Scholarly Journals-1 ObjectType-Feature-2 content type line 23 |
| PMID | 26605612 |
| PQID | 1749601141 |
| PQPubID | 23479 |
| ParticipantIDs | proquest_miscellaneous_1749601141 pubmed_primary_26605612 |
| PublicationCentury | 2000 |
| PublicationDate | 2012-11-13 |
| PublicationDateYYYYMMDD | 2012-11-13 |
| PublicationDate_xml | – month: 11 year: 2012 text: 2012-11-13 day: 13 |
| PublicationDecade | 2010 |
| PublicationPlace | United States |
| PublicationPlace_xml | – name: United States |
| PublicationTitle | Journal of chemical theory and computation |
| PublicationTitleAlternate | J Chem Theory Comput |
| PublicationYear | 2012 |
| SSID | ssj0033423 |
| Score | 2.4089072 |
| Snippet | Accurate quantum mechanical/molecular mechanical (QM/MM) treatments should account for MM polarization and properly include long-range electrostatic... |
| SourceID | proquest pubmed |
| SourceType | Aggregation Database Index Database |
| StartPage | 4527 |
| Title | Solvent Boundary Potentials for Hybrid QM/MM Computations Using Classical Drude Oscillators: A Fully Polarizable Model |
| URI | https://www.ncbi.nlm.nih.gov/pubmed/26605612 https://www.proquest.com/docview/1749601141 |
| Volume | 8 |
| WOSCitedRecordID | wos000311191900054&url=https%3A%2F%2Fcvtisr.summon.serialssolutions.com%2F%23%21%2Fsearch%3Fho%3Df%26include.ft.matches%3Dt%26l%3Dnull%26q%3D |
| hasFullText | |
| inHoldings | 1 |
| isFullTextHit | |
| isPrint | |
| link | http://cvtisr.summon.serialssolutions.com/2.0.0/link/0/eLvHCXMwpV1LS8NAEF7UCnrx_agvVvAamuwmzcaL1EfpwdSKRXor2c0EhJLUNi303zuTh54EwUtOyRJmh51vHvt9jN0YSEQSQGyJONaW6-nAUm3bt1yNsVAjJFdF9_z92e_31WgUDKqC27waq6zPxOKgjjNDNfIWImcE24jenbvpp0WqUdRdrSQ01llDIpShkS5_9N1FkMRuV_ClusRC6aiaWUg4LZNLIs0W8DuyLCJMd_e__7bHdipsyTulM-yzNUgP2NZDLel2yJZv2YTmG_l9IaY0W_FBltO4EPogR_TKeyu6wMVfw1YY8lLwoazo8WK0gBcSmrSt_HG2iIG_YABFPyLFnlve4ZTP0pqYLdOw2AQ4Sa1Njtiw-zR86FmV8IIVScfLLQE68bTngKMlKN_YJgIJcWQHESEeqbSwA5VE0lcQU-xTthGI5ZQygY5cccw20iyFU8btJDaulMbYHrgQAYIx0PiKF3hEa9NusuvaomO0BTUrohSyxXz8Y9MmOym3ZTwtCTjGCCoo8RFnf_j6nG0jxhF0fdCRF6yRoEXhkm2aZf4xn10VDoPP_iD8AmRvzGY |
| linkProvider | ProQuest |
| openUrl | ctx_ver=Z39.88-2004&ctx_enc=info%3Aofi%2Fenc%3AUTF-8&rfr_id=info%3Asid%2Fsummon.serialssolutions.com&rft_val_fmt=info%3Aofi%2Ffmt%3Akev%3Amtx%3Ajournal&rft.genre=article&rft.atitle=Solvent+Boundary+Potentials+for+Hybrid+QM%2FMM+Computations+Using+Classical+Drude+Oscillators%3A+A+Fully+Polarizable+Model&rft.jtitle=Journal+of+chemical+theory+and+computation&rft.au=Boulanger%2C+Eliot&rft.au=Thiel%2C+Walter&rft.date=2012-11-13&rft.issn=1549-9618&rft.volume=8&rft.issue=11&rft.spage=4527&rft_id=info:doi/10.1021%2Fct300722e&rft.externalDBID=NO_FULL_TEXT |
| thumbnail_l | http://covers-cdn.summon.serialssolutions.com/index.aspx?isbn=/lc.gif&issn=1549-9618&client=summon |
| thumbnail_m | http://covers-cdn.summon.serialssolutions.com/index.aspx?isbn=/mc.gif&issn=1549-9618&client=summon |
| thumbnail_s | http://covers-cdn.summon.serialssolutions.com/index.aspx?isbn=/sc.gif&issn=1549-9618&client=summon |