Benchmark Study of the SCC-DFTB Approach for a Biomolecular Proton Channel

Overview

Journal J Chem Theory Comput

Publisher American Chemical Society

Specialties Biochemistry
Chemistry

Date 2014 Aug 9

PMID 25104919

Citations 16

Authors

Ruibin Liang

Jessica M J Swanson

Gregory A Voth

Affiliations

Soon will be listed here.

Abstract

The self-consistent charge density functional tight binding (SCC-DFTB) method has been increasingly applied to study proton transport (PT) in biological environments. However, recent studies revealing some significant limitations of SCC-DFTB for proton and hydroxide solvation and transport in bulk aqueous systems call into question its accuracy for simulating PT in biological systems. The current work benchmarks the SCC-DFTB/MM method against more accurate DFT/MM by simulating PT in a synthetic leucine-serine channel (LS2), which emulates the structure and function of biomolecular proton channels. It is observed that SCC-DFTB/MM produces over-coordinated and less structured pore water, an over-coordinated excess proton, weak hydrogen bonds around the excess proton charge defect and qualitatively different PT dynamics. Similar issues are demonstrated for PT in a carbon nanotube, indicating that the inaccuracies found for SCC-DFTB are not due to the point charge based QM/MM electrostatic coupling scheme, but rather to the approximations of the semiempirical method itself. The results presented in this work highlight the limitations of the present form of the SCC-DFTB/MM approach for simulating PT processes in biological protein or channel-like environments, while providing benchmark results that may lead to an improvement of the underlying method.

Citing Articles

Quantum chemical investigation of the effect of alkali metal ions on the dynamic structure of water in aqueous solutions.

Khanal R, Irle S RSC Adv. 2022; 12(39):25500-25510.

PMID: 36275866 PMC: 9480497. DOI: 10.1039/d2ra04563j.

Multiscale Simulation of an Influenza A M2 Channel Mutant Reveals Key Features of Its Markedly Different Proton Transport Behavior.

Watkins L, DeGrado W, Voth G J Am Chem Soc. 2022; 144(2):769-776.

PMID: 34985907 PMC: 8834648. DOI: 10.1021/jacs.1c09281.

Using Machine Learning to Greatly Accelerate Path Integral Molecular Dynamics.

Li C, Voth G J Chem Theory Comput. 2022; 18(2):599-604.

PMID: 34982562 PMC: 8864787. DOI: 10.1021/acs.jctc.1c01085.

Quantum Phonon Transport in Nanomaterials: Combining Atomistic with Non-Equilibrium Green's Function Techniques.

Medrano Sandonas L, Gutierrez R, Pecchia A, Croy A, Cuniberti G Entropy (Basel). 2020; 21(8).

PMID: 33267449 PMC: 7515264. DOI: 10.3390/e21080735.

Understanding and Tracking the Excess Proton in Ab Initio Simulations; Insights from IR Spectra.

Li C, Swanson J J Phys Chem B. 2020; 124(27):5696-5708.

PMID: 32515957 PMC: 7448536. DOI: 10.1021/acs.jpcb.0c03615.

References

Randa H, Forrest L, Voth G, Sansom M . Molecular dynamics of synthetic leucine-serine ion channels in a phospholipid membrane. Biophys J. 1999; 77(5):2400-10. PMC: 1300517. DOI: 10.1016/S0006-3495(99)77077-3. View

Goedecker , Teter , HUTTER . Separable dual-space Gaussian pseudopotentials. Phys Rev B Condens Matter. 1996; 54(3):1703-1710. DOI: 10.1103/physrevb.54.1703. View

Yang Y, Yu H, York D, Cui Q, Elstner M . Extension of the self-consistent-charge density-functional tight-binding method: third-order expansion of the density functional theory total energy and introduction of a modified effective coulomb interaction. J Phys Chem A. 2007; 111(42):10861-73. DOI: 10.1021/jp074167r. View

MacKerell A, Bashford D, Bellott M, Dunbrack R, Evanseck J, Field M . All-atom empirical potential for molecular modeling and dynamics studies of proteins. J Phys Chem B. 2014; 102(18):3586-616. DOI: 10.1021/jp973084f. View

Brooks B, Brooks 3rd C, MacKerell Jr A, Nilsson L, Petrella R, Roux B . CHARMM: the biomolecular simulation program. J Comput Chem. 2009; 30(10):1545-614. PMC: 2810661. DOI: 10.1002/jcc.21287. View

Wu Y, Voth G . A computer simulation study of the hydrated proton in a synthetic proton channel. Biophys J. 2003; 85(2):864-75. PMC: 1303208. DOI: 10.1016/S0006-3495(03)74526-3. View

Wu Y, Ilan B, Voth G . Charge delocalization in proton channels, II: the synthetic LS2 channel and proton selectivity. Biophys J. 2006; 92(1):61-9. PMC: 1697870. DOI: 10.1529/biophysj.106.091942. View

Choi T, Liang R, Maupin C, Voth G . Application of the SCC-DFTB method to hydroxide water clusters and aqueous hydroxide solutions. J Phys Chem B. 2013; 117(17):5165-79. DOI: 10.1021/jp400953a. View

Shepherd L, Morrison C . Simulating proton transport through a simplified model for trans-membrane proteins. J Phys Chem B. 2010; 114(20):7047-55. DOI: 10.1021/jp910262d. View

10.

Goyal P, Ghosh N, Phatak P, Clemens M, Gaus M, Elstner M . Proton storage site in bacteriorhodopsin: new insights from quantum mechanics/molecular mechanics simulations of microscopic pK(a) and infrared spectra. J Am Chem Soc. 2011; 133(38):14981-97. PMC: 3178665. DOI: 10.1021/ja201568s. View

11.

Markovitch O, Chen H, Izvekov S, Paesani F, Voth G, Agmon N . Special pair dance and partner selection: elementary steps in proton transport in liquid water. J Phys Chem B. 2008; 112(31):9456-66. DOI: 10.1021/jp804018y. View

12.

Riccardi D, Yang S, Cui Q . Proton transfer function of carbonic anhydrase: Insights from QM/MM simulations. Biochim Biophys Acta. 2009; 1804(2):342-51. PMC: 6787916. DOI: 10.1016/j.bbapap.2009.07.026. View

13.

Riccardi D, Schaefer P, Cui Q . pKa calculations in solution and proteins with QM/MM free energy perturbation simulations: a quantitative test of QM/MM protocols. J Phys Chem B. 2006; 109(37):17715-33. DOI: 10.1021/jp0517192. View

14.

Hou G, Zhu X, Elstner M, Cui Q . A modified QM/MM Hamiltonian with the Self-Consistent-Charge Density-Functional-Tight-Binding Theory for highly charged QM regions. J Chem Theory Comput. 2013; 8(11):4293-4304. PMC: 3529911. DOI: 10.1021/ct300649f. View

15.

MacKerell Jr A, Feig M, Brooks 3rd C . Extending the treatment of backbone energetics in protein force fields: limitations of gas-phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations. J Comput Chem. 2004; 25(11):1400-15. DOI: 10.1002/jcc.20065. View

16.

Laino T, Mohamed F, Laio A, Parrinello M . An Efficient Linear-Scaling Electrostatic Coupling for Treating Periodic Boundary Conditions in QM/MM Simulations. J Chem Theory Comput. 2015; 2(5):1370-8. DOI: 10.1021/ct6001169. View

17.

Riccardi D, Konig P, Prat-Resina X, Yu H, Elstner M, Frauenheim T . "Proton holes" in long-range proton transfer reactions in solution and enzymes: A theoretical analysis. J Am Chem Soc. 2006; 128(50):16302-11. PMC: 2561195. DOI: 10.1021/ja065451j. View

18.

Lear J, Wasserman Z, Degrado W . Synthetic amphiphilic peptide models for protein ion channels. Science. 1988; 240(4856):1177-81. DOI: 10.1126/science.2453923. View

19.

Voth G . The computer simulation of proton transport in biomolecular systems. Front Biosci. 2003; 8:s1384-79. DOI: 10.2741/1213. View

20.

Gaus M, Goez A, Elstner M . Parametrization and Benchmark of DFTB3 for Organic Molecules. J Chem Theory Comput. 2015; 9(1):338-54. DOI: 10.1021/ct300849w. View