Phosphate chemistry is involved in many key biological processes, yet the underlying mechanism often remains unclear. For theoretical analysis to effectively complement experimental mechanistic analysis, it is essential to develop computational methods that can capture the complexity of the underlying potential energy surface and allow for sufficient sampling of the configurational space. To this end, we report the parametrization of an approximate density functional theory, the Self-Consistent-Charge Density-Functional Tight-Binding (SCC-DFTB) method for systems containing phosphorus. Compared to high-level density functional theory and ab initio (IVIP2 and G3B3) results, the standard second-order parametrization is shown to give reliable structures for a diverse set of phosphate compounds but inaccurate energetics. With the on-site third-order terms included, referred to as SCC-DFTBPA, calculated proton affinities of phosphate compounds are substantially improved, although it remains difficult to obtain reliable proton affinity for both phosphates and compounds that do not contain phosphorus, indicating that further improvement in the formulation of SCC-DFTB is still a challenge to meet. To make SCC-DFTB applicable to phosphate reactions in the current (on-site-third-order-only) formulation, a "reaction-specific" parametrization, referred to as SCC-DFTBPR, is developed based on hydrolysis reactions of model phosphate species. Benchmark calculations in both the gas phase and solution phase indicate that SCC-DFTBPR gives reliable structural properties and semiquantitative energetics for phosphate hydrolysis reactions. Since the number of reaction-specific parameters is small, it is likely that SCC-DFTBPR is applicable to a broad set of phosphate species. Indeed, for 56 reaction exothermicities and 47 energy barriers related to RNA catalysis model reactions collected from the QCRNA database, which involve molecules rather different from those used to parametrize SCC-DFTBPR, the corresponding root-mean-square difference between SCC-DFTBPR and high-level DFT results is only 5.3 kcal/mol. We hope that the parametrized SCC-DFTB models will complement NDDO based reaction-specific models (e.g., AM1-d/PhoT) and high-level ab initio QM/MM methods in better understanding the mechanism of phosphate chemistry in condensed phase, particularly biological systems.
Description of Phosphate Hydrolysis Reactions with the Self-Consistent-Charge Density-Functional-Tight-Binding (SCC-DFTB) Theory. 1. Parameterization
Creators
Yang Yang - Drexel University, Chemistry
Haibo Yu - University of Wisconsin–Madison
Darrin York - University of Wisconsin–Madison
Marcus Elstner - University of Wisconsin–Madison
Qiang Cui - University of Wisconsin–Madison
Publication Details
Journal of chemical theory and computation, v 4(12), pp 2067-2084
Publisher
Amer Chemical Soc
Number of pages
18
Grant note
Alfred P. Sloan Foundation
National Center for Supercomputing Applications at the University of Illinois
Minnesota Supercomputing Institute (MSI)
R01-GM071428; GM62248 / National Institutes of Health; United States Department of Health & Human Services; National Institutes of Health (NIH) - USA
R01GM071428 / NATIONAL INSTITUTE OF GENERAL MEDICAL SCIENCES; United States Department of Health & Human Services; National Institutes of Health (NIH) - USA; NIH National Institute of General Medical Sciences (NIGMS)
Resource Type
Journal article
Language
English
Academic Unit
Chemistry
Web of Science ID
WOS:000261613800009
Scopus ID
2-s2.0-58149235077
Other Identifier
991021874215504721
InCites Highlights
Data related to this publication, from InCites Benchmarking & Analytics tool:
Collaboration types
Domestic collaboration
International collaboration
Web of Science research areas
Chemistry, Physical
Physics, Atomic, Molecular & Chemical
Research Home Page
Browse by research and academic units
Learn about the ETD submission process at Drexel
Learn about the Libraries’ research data management services