1
|
Zhang J, Kriebel CN, Wan Z, Shi M, Glaubitz C, He X. Automated Fragmentation Quantum Mechanical Calculation of 15N and 13C Chemical Shifts in a Membrane Protein. J Chem Theory Comput 2023; 19:7405-7422. [PMID: 37788419 DOI: 10.1021/acs.jctc.3c00621] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
In this work, we developed an accurate and cost-effective automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) method to calculate the chemical shifts of 15N and 13C of membrane proteins. The convergence of the AF-QM/MM method was tested using Krokinobacter eikastus rhodopsin 2 as a test case. When the distance threshold of the QM region is equal to or larger than 4.0 Å, the results of the AF-QM/MM calculations are close to convergence. In addition, the effects of selected density functionals, basis sets, and local chemical environment of target atoms on the chemical shift calculations were systematically investigated. Our results demonstrate that the predicted chemical shifts are more accurate when important environmental factors including cross-protomer interactions, lipid molecules, and solvent molecules are taken into consideration, especially for the 15N chemical shift prediction. Furthermore, with the presence of sodium ions in the environment, the chemical shift of residues, retinal, and retinal Schiff base are affected, which is consistent with the results of the solid-state nuclear magnetic resonance (NMR) experiment. Upon comparing the performance of various density functionals (namely, B3LYP, B3PW91, M06-2X, M06-L, mPW1PW91, OB95, and OPBE), the results show that mPW1PW91 is a suitable functional for the 15N and 13C chemical shift prediction of the membrane proteins. Meanwhile, we find that the improved accuracy of the 13Cβ chemical shift calculations can be achieved by the employment of the triple-ζ basis set. However, the employment of the triple-ζ basis set does not improve the accuracy of the 15N and 13Cα chemical shift calculations nor does the addition of a diffuse function improve the overall prediction accuracy of the chemical shifts. Our study also underscores that the AF-QM/MM method has significant advantages in predicting the chemical shifts of key ligands and nonstandard residues in membrane proteins than most widely used empirical models; therefore, it could be an accurate computational tool for chemical shift calculations on various types of biological systems.
Collapse
Affiliation(s)
- Jinhuan Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Clara Nassrin Kriebel
- Institute of Biophysical Chemistry and Center for Biomolecular Magnetic Resonance, Goethe University Frankfurt, 60438 Frankfurt am Main, Germany
| | - Zheng Wan
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Man Shi
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Clemens Glaubitz
- Institute of Biophysical Chemistry and Center for Biomolecular Magnetic Resonance, Goethe University Frankfurt, 60438 Frankfurt am Main, Germany
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200062, China
| |
Collapse
|
2
|
Krivdin LB. Computational 1 H and 13 C NMR in structural and stereochemical studies. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2022; 60:733-828. [PMID: 35182410 DOI: 10.1002/mrc.5260] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Revised: 02/14/2022] [Accepted: 02/16/2022] [Indexed: 06/14/2023]
Abstract
Present review outlines the advances and perspectives of computational 1 H and 13 C NMR applied to the stereochemical studies of inorganic, organic, and bioorganic compounds, involving in particular natural products, carbohydrates, and carbonium ions. The first part of the review briefly outlines theoretical background of the modern computational methods applied to the calculation of chemical shifts and spin-spin coupling constants at the DFT and the non-empirical levels. The second part of the review deals with the achievements of the computational 1 H and 13 C NMR in the stereochemical investigation of a variety of inorganic, organic, and bioorganic compounds, providing in an abridged form the material partly discussed by the author in a series of parent reviews. Major attention is focused herewith on the publications of the recent years, which were not reviewed elsewhere.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Irkutsk, Russia
| |
Collapse
|
3
|
Payard PA, Perego LA, Grimaud L, Ciofini I. A DFT Protocol for the Prediction of 31P NMR Chemical Shifts of Phosphine Ligands in First-Row Transition-Metal Complexes. Organometallics 2020. [DOI: 10.1021/acs.organomet.0c00309] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Pierre-Adrien Payard
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Luca Alessandro Perego
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
- PSL University, Institute of Chemistry for Health and Life Sciences, I-CLeHS, CNRS-Chimie ParisTech, 11 rue P. et M. Curie, F-75005 Paris 05, France
| | - Laurence Grimaud
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Ilaria Ciofini
- PSL University, Institute of Chemistry for Health and Life Sciences, I-CLeHS, CNRS-Chimie ParisTech, 11 rue P. et M. Curie, F-75005 Paris 05, France
| |
Collapse
|
4
|
Unzueta PA, Beran GJO. Polarizable continuum models provide an effective electrostatic embedding model for fragment-based chemical shift prediction in challenging systems. J Comput Chem 2020; 41:2251-2265. [PMID: 32748418 DOI: 10.1002/jcc.26388] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Revised: 06/04/2020] [Accepted: 07/04/2020] [Indexed: 12/25/2022]
Abstract
Ab initio nuclear magnetic resonance chemical shift prediction provides an important tool for interpreting and assigning experimental spectra, but it becomes computationally prohibitive in large systems. The computational costs can be reduced considerably by fragmentation of the large system into a series of contributions from many smaller subsystems. However, the presence of charged functional groups and the need to partition the system across covalent bonds create complications in biomolecules that typically require the use of large fragments and careful descriptions of the electrostatic environment. The present work shows how a model that combines chemical shielding contributions from non-overlapping monomer and dimer fragments embedded in a polarizable continuum model provides a simple, easy-to-implement, and computationally inexpensive approach for predicting chemical shifts in complex systems. The model's performance proves rather insensitive to the continuum dielectric constant, making the selection of the optimal embedding dielectric less critical. The PCM-embedded fragment model is demonstrated to perform well across systems ranging from molecular crystals to proteins.
Collapse
Affiliation(s)
- Pablo A Unzueta
- Department of Chemistry, Univeristy of California, Riverside, California, USA
| | - Gregory J O Beran
- Department of Chemistry, Univeristy of California, Riverside, California, USA
| |
Collapse
|
5
|
Casabianca LB. Calculating nuclear magnetic resonance chemical shifts in solvated systems. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2020; 58:611-624. [PMID: 31916612 DOI: 10.1002/mrc.4994] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/04/2019] [Revised: 12/30/2019] [Accepted: 01/02/2020] [Indexed: 06/10/2023]
Abstract
The nuclear magnetic resonance (NMR) chemical shift is extremely sensitive to molecular geometry, hydrogen bonding, solvent, temperature, pH, and concentration. Calculated magnetic shielding constants, converted to chemical shifts, can be valuable aids in NMR peak assignment and can also give detailed information about molecular geometry and intermolecular effects. Calculating chemical shifts in solution is complicated by the need to include solvent effects and conformational averaging. Here, we review the current state of NMR chemical shift calculations in solution, beginning with an introduction to the theory of calculating magnetic shielding in general, then covering methods for inclusion of solvent effects and conformational averaging, and finally discussing examples of applications using calculated chemical shifts to gain detailed structural information.
Collapse
Affiliation(s)
- Leah B Casabianca
- Department of Chemistry, Clemson University, Clemson, South Carolina
| |
Collapse
|
6
|
Krivdin LB. Computational 1 H NMR: Part 3. Biochemical studies. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2020; 58:15-30. [PMID: 31286566 DOI: 10.1002/mrc.4895] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/20/2019] [Revised: 05/14/2019] [Accepted: 05/18/2019] [Indexed: 06/09/2023]
Abstract
This is the third and the last part of three closely interrelated reviews dealing with computation of 1 H nuclear magnetic resonance chemical shifts and 1 H-1 H spin-spin coupling constants. Present review deals with the computation of these parameters in biologically active natural products, carbohydrates, and other molecules of biological origin focusing on stereochemical applications of computational 1 H nuclear magnetic resonance to these objects.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Irkutsk, Russia
- Department of Chemistry, Angarsk State Technical University, Angarsk, Russia
| |
Collapse
|
7
|
Abstract
Since the introduction of the fragment molecular orbital method 20 years ago, fragment-based approaches have occupied a small but growing niche in quantum chemistry. These methods decompose a large molecular system into subsystems small enough to be amenable to electronic structure calculations, following which the subsystem information is reassembled in order to approximate an otherwise intractable supersystem calculation. Fragmentation sidesteps the steep rise (with respect to system size) in the cost of ab initio calculations, replacing it with a distributed cost across numerous computer processors. Such methods are attractive, in part, because they are easily parallelizable and therefore readily amenable to exascale computing. As such, there has been hope that distributed computing might offer the proverbial "free lunch" in quantum chemistry, with the entrée being high-level calculations on very large systems. While fragment-based quantum chemistry can count many success stories, there also exists a seedy underbelly of rarely acknowledged problems. As these methods begin to mature, it is time to have a serious conversation about what they can and cannot be expected to accomplish in the near future. Both successes and challenges are highlighted in this Perspective.
Collapse
Affiliation(s)
- John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
8
|
Krivdin LB. Computational protocols for calculating 13C NMR chemical shifts. PROGRESS IN NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY 2019; 112-113:103-156. [PMID: 31481156 DOI: 10.1016/j.pnmrs.2019.05.004] [Citation(s) in RCA: 58] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/20/2019] [Revised: 05/19/2019] [Accepted: 05/20/2019] [Indexed: 06/10/2023]
Abstract
The most recent results dealing with the computation of 13C NMR chemical shifts in chemistry (small molecules, saturated, unsaturated and aromatic compounds, heterocycles, functional derivatives, coordination complexes, carbocations, and natural products) are reviewed, paying special attention to theoretical background and accuracy, the latter involving solvent effects, vibrational corrections, and relativistic effects.
Collapse
Affiliation(s)
- Leonid B Krivdin
- A. E. Favorsky Irkutsk Institute of Chemistry, Siberian Branch of the Russian Academy of Sciences, Favorsky St. 1, 664033 Irkutsk, Russia.
| |
Collapse
|
9
|
Kobayashi R, Amos RD, Reid DM, Collins MA. Application of the Systematic Molecular Fragmentation by Annihilation Method to ab Initio NMR Chemical Shift Calculations. J Phys Chem A 2018; 122:9135-9141. [DOI: 10.1021/acs.jpca.8b09565] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Rika Kobayashi
- ANU Supercomputer Facility, Leonard Huxley Building 56, Mills Road, Canberra, ACT 2601, Australia
| | - Roger D. Amos
- ANU Supercomputer Facility, Leonard Huxley Building 56, Mills Road, Canberra, ACT 2601, Australia
| | - David M. Reid
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| | - Michael A. Collins
- Research School of Chemistry, Australian National University, Canberra, ACT 2601, Australia
| |
Collapse
|
10
|
Iron MA. Evaluation of the Factors Impacting the Accuracy of 13C NMR Chemical Shift Predictions using Density Functional Theory-The Advantage of Long-Range Corrected Functionals. J Chem Theory Comput 2017; 13:5798-5819. [PMID: 29016125 DOI: 10.1021/acs.jctc.7b00772] [Citation(s) in RCA: 58] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The various factors influencing the accuracy of 13C NMR calculations using density functional theory (DFT), including the basis set, exchange-correlation (XC) functional, and isotropic shielding calculation method, are evaluated. A wide selection of XC functionals (over 70) were considered, and it was found that long-range corrected functionals offer a significant improvement over the other classes of functionals. Based on a thorough study, it is recommended that for calculating NMR chemical shifts (δ) one should use the CSGT method, the COSMO solvation model, and the LC-TPSSTPSS exchange-correlation functional in conjunction with the cc-pVTZ basis set. A selection of problems in natural product identification are considered in light of the newly recommended level of theory.
Collapse
Affiliation(s)
- Mark A Iron
- Computational Chemistry Unit, Department of Chemical Research Support, Weizmann Institute of Science , Rehovot 7610001, Israel
| |
Collapse
|
11
|
Zhao D, Song R, Li W, Ma J, Dong H, Li S. Accurate Prediction of NMR Chemical Shifts in Macromolecular and Condensed-Phase Systems with the Generalized Energy-Based Fragmentation Method. J Chem Theory Comput 2017; 13:5231-5239. [DOI: 10.1021/acs.jctc.7b00380] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Dongbo Zhao
- Key
Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute
of Theoretical and Computational Chemistry, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
- Kuang
Yaming Honors School, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Ruiheng Song
- Kuang
Yaming Honors School, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Wei Li
- Key
Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute
of Theoretical and Computational Chemistry, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Jing Ma
- Key
Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute
of Theoretical and Computational Chemistry, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Hao Dong
- Kuang
Yaming Honors School, Nanjing University, Nanjing 210023, People’s Republic of China
| | - Shuhua Li
- Key
Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute
of Theoretical and Computational Chemistry, School of Chemistry and
Chemical Engineering, Nanjing University, Nanjing 210023, People’s Republic of China
| |
Collapse
|