1
|
Yesiltepe Y, Govind N, Metz TO, Renslow RS. An initial investigation of accuracy required for the identification of small molecules in complex samples using quantum chemical calculated NMR chemical shifts. J Cheminform 2022; 14:64. [PMID: 36138446 PMCID: PMC9499888 DOI: 10.1186/s13321-022-00587-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2020] [Accepted: 02/06/2022] [Indexed: 11/24/2022] Open
Abstract
The majority of primary and secondary metabolites in nature have yet to be identified, representing a major challenge for metabolomics studies that currently require reference libraries from analyses of authentic compounds. Using currently available analytical methods, complete chemical characterization of metabolomes is infeasible for both technical and economic reasons. For example, unambiguous identification of metabolites is limited by the availability of authentic chemical standards, which, for the majority of molecules, do not exist. Computationally predicted or calculated data are a viable solution to expand the currently limited metabolite reference libraries, if such methods are shown to be sufficiently accurate. For example, determining nuclear magnetic resonance (NMR) spectroscopy spectra in silico has shown promise in the identification and delineation of metabolite structures. Many researchers have been taking advantage of density functional theory (DFT), a computationally inexpensive yet reputable method for the prediction of carbon and proton NMR spectra of metabolites. However, such methods are expected to have some error in predicted 13C and 1H NMR spectra with respect to experimentally measured values. This leads us to the question–what accuracy is required in predicted 13C and 1H NMR chemical shifts for confident metabolite identification? Using the set of 11,716 small molecules found in the Human Metabolome Database (HMDB), we simulated both experimental and theoretical NMR chemical shift databases. We investigated the level of accuracy required for identification of metabolites in simulated pure and impure samples by matching predicted chemical shifts to experimental data. We found 90% or more of molecules in simulated pure samples can be successfully identified when errors of 1H and 13C chemical shifts in water are below 0.6 and 7.1 ppm, respectively, and below 0.5 and 4.6 ppm in chloroform solvation, respectively. In simulated complex mixtures, as the complexity of the mixture increased, greater accuracy of the calculated chemical shifts was required, as expected. However, if the number of molecules in the mixture is known, e.g., when NMR is combined with MS and sample complexity is low, the likelihood of confident molecular identification increased by 90%.
Collapse
Affiliation(s)
- Yasemin Yesiltepe
- The Gene and Linda Voiland School of Chemical Engineering and Bioengineering, Washington State University, Pullman, WA, USA.,Earth and Biological Sciences Directorate, Pacific Northwest National Laboratory, Richland, WA, USA
| | - Niranjan Govind
- Earth and Biological Sciences Directorate, Pacific Northwest National Laboratory, Richland, WA, USA
| | - Thomas O Metz
- The Gene and Linda Voiland School of Chemical Engineering and Bioengineering, Washington State University, Pullman, WA, USA
| | - Ryan S Renslow
- The Gene and Linda Voiland School of Chemical Engineering and Bioengineering, Washington State University, Pullman, WA, USA. .,Earth and Biological Sciences Directorate, Pacific Northwest National Laboratory, Richland, WA, USA.
| |
Collapse
|
2
|
Nazarski RB. Summary of DFT calculations coupled with current statistical and/or artificial neural network (ANN) methods to assist experimental NMR data in identifying diastereomeric structures. Tetrahedron Lett 2021. [DOI: 10.1016/j.tetlet.2020.152548] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
3
|
Hernandes IS, Da Silva HC, Dos Santos HF, De Almeida WB. Conformational Analysis of 5,4'-Dihydroxy-7,5',3'-trimethoxyisoflavone in Solution Using 1H NMR: A Density Functional Theory Approach. J Phys Chem A 2020; 124:5182-5193. [PMID: 32466649 DOI: 10.1021/acs.jpca.0c02996] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Among 20 compounds isolated from the extracts of Ouratea ferruginea the 5,4'-dihydroxy-7,5',3'-trimethoxyisoflavone (9) showed the best inhibitory effect on glutathione S-transferase (GST) and so deserves our attention. In this work we investigated the preferred molecular structure of 9 in chloroform solution using the density functional theory (DFT) and molecular dynamics simulation. Comparison between experimental 1H NMR data in CDCl3 solution and calculated chemical shifts enabled us to precisely determine the conformation adopted by 9 in solution, which can be used in further theoretical studies involving interaction with biological targets. Moreover, the experimental NMR data were used as reference to assess the ability of DFT based methods to predict 1H NMR spectrum in solution for organic compounds. Among various DFT functionals the hybrid B3LYP was the most adequate for the calculation of chemical shifts in what CHn protons are concerned. Regarding the OH hydrogen, inclusion of explicit CHCl3 solvent molecules adequately placed around the solute led to good agreement with the experimental chemical shifts (in CDCl3). It is a well-known fact that theoretical prediction of chemical shifts for OH hydrogens poses as a challenge and also revealed that the way the solvent effects are included in the DFT calculations is crucial for the right prediction of the whole 1H NMR spectrum. It was found in this work that a supermolecule solute-solvent calculation with a minimum of four CHCl3 molecules is enough to correctly reproduce the 1H NMR experimental profile observed in solution, revealing that the calculated solvated structure used to reproduce the NMR chemical shifts is not unique.
Collapse
Affiliation(s)
- Isabel S Hernandes
- Laboratório de Química Computacional e Modelagem Molecular (LQC-MM), Departamento de Química Inorgânica, Instituto de Química, Universidade Federal Fluminense (UFF), Outeiro de São João Batista s/n, Campus do Valonguinho, Centro, Niterói, RJ 24020-141, Brazil
| | - Haroldo C Da Silva
- Laboratório de Química Computacional e Modelagem Molecular (LQC-MM), Departamento de Química Inorgânica, Instituto de Química, Universidade Federal Fluminense (UFF), Outeiro de São João Batista s/n, Campus do Valonguinho, Centro, Niterói, RJ 24020-141, Brazil
| | - Hélio F Dos Santos
- Núcleo de Estudos em Química Computacional (NEQC), Departamento de Química, ICE, Universidade Federal de Juiz de Fora (UFJF), Campus Universitário, Martelos, Juiz de Fora, MG 36036-330, Brazil
| | - Wagner B De Almeida
- Laboratório de Química Computacional e Modelagem Molecular (LQC-MM), Departamento de Química Inorgânica, Instituto de Química, Universidade Federal Fluminense (UFF), Outeiro de São João Batista s/n, Campus do Valonguinho, Centro, Niterói, RJ 24020-141, Brazil
| |
Collapse
|
4
|
Oveisi AR, Karimi P, Delarami HS, Daliran S, Khorramabadi-Zad A, Khajeh M, Sanchooli E, Ghaffari-Moghaddam M. New porphyrins: synthesis, characterization, and computational studies. Mol Divers 2019; 24:335-344. [PMID: 31062142 DOI: 10.1007/s11030-019-09955-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2019] [Accepted: 04/29/2019] [Indexed: 10/26/2022]
Abstract
New trans-A2B2-porphyrins substituted at phenyl positions were synthesized from 4-methylphthalic acid as a starting material through sequential multistep reactions. These macrocycles were characterized by 1H NMR, 13C NMR, 19F NMR, 1H-1H COSY NMR, and MALDI-TOF mass spectrometry. Computational studies were performed on the porphyrins to investigate various factors such as structural features, electronic energy, energy gaps, and aromaticity. Energy band gap values of these compounds especially N-hydroxyphthalimide-functionalized porphyrins were small that makes them as good candidates for solar cell systems and photocatalysis. Relationships between electronic energies and aromaticity of the compounds were then investigated. The data indicated that the aromaticity features at the center of two series of these compounds (fluorinated and non-fluorinated porphyrins) were in the opposite manner.
Collapse
Affiliation(s)
- Ali Reza Oveisi
- Department of Chemistry, Faculty of Science, University of Zabol, Zabol, Iran.
| | - Pouya Karimi
- Department of Chemistry, Faculty of Science, University of Zabol, Zabol, Iran
| | | | - Saba Daliran
- Faculty of Chemistry, Bu-Ali Sina University, Hamedan, 6517838695, Iran
| | | | - Mostafa Khajeh
- Department of Chemistry, Faculty of Science, University of Zabol, Zabol, Iran
| | - Esmael Sanchooli
- Department of Chemistry, Faculty of Science, University of Zabol, Zabol, Iran
| | | |
Collapse
|
5
|
Kapusta K, Voronkov E, Okovytyy S, Korobov V, Leszczynski J. Reconstruction of STO-3G Family Basis Set for the Accurate Calculation of Magnetic Properties. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2019. [DOI: 10.1134/s0036024418130174] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
6
|
An automated framework for NMR chemical shift calculations of small organic molecules. J Cheminform 2018; 10:52. [PMID: 30367288 PMCID: PMC6755567 DOI: 10.1186/s13321-018-0305-8] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2018] [Accepted: 10/09/2018] [Indexed: 12/04/2022] Open
Abstract
When using nuclear magnetic resonance (NMR) to assist in chemical identification in complex samples, researchers commonly rely on databases for chemical shift spectra. However, authentic standards are typically depended upon to build libraries experimentally. Considering complex biological samples, such as blood and soil, the entirety of NMR spectra required for all possible compounds would be infeasible to ascertain due to limitations of available standards and experimental processing time. As an alternative, we introduce the in silico Chemical Library Engine (ISiCLE) NMR chemical shift module to accurately and automatically calculate NMR chemical shifts of small organic molecules through use of quantum chemical calculations. ISiCLE performs density functional theory (DFT)-based calculations for predicting chemical properties—specifically NMR chemical shifts in this manuscript—via the open source, high-performance computational chemistry software, NWChem. ISiCLE calculates the NMR chemical shifts of sets of molecules using any available combination of DFT method, solvent, and NMR-active nuclei, using both user-selected reference compounds and/or linear regression methods. Calculated NMR chemical shifts are provided to the user for each molecule, along with comparisons with respect to a number of metrics commonly used in the literature. Here, we demonstrate ISiCLE using a set of 312 molecules, ranging in size up to 90 carbon atoms. For each, calculation of NMR chemical shifts have been performed with 8 different levels of DFT theory, and with solvation effects using the implicit solvent Conductor-like Screening Model. The DFT method dependence of the calculated chemical shifts have been systematically investigated through benchmarking and subsequently compared to experimental data available in the literature. Furthermore, ISiCLE has been applied to a set of 80 methylcyclohexane conformers, combined via Boltzmann weighting and compared to experimental values. We demonstrate that our protocol shows promise in the automation of chemical shift calculations and, ultimately, the expansion of chemical shift libraries.
![]()
Collapse
|
7
|
Abstract
We investigate by explicit parameter optimization to what extent basis sets of polarized double-ζ quality can introduce compensating errors in five different density functional methods. It is shown that minor changes in the contraction coefficients of the valence functions in the basis sets can have a significant impact and allow different density functional methods to achieve very similar performances. This holds for nuclear magnetic shielding constants and for isomerization energies, barrier heights, and noncovalent interactions. It is furthermore shown that errors due to neglect of vibrational and solvent effects can be absorbed in the combined method and basis set errors. These findings hold for data sets consisting of 50-150 data points. This raises the question of whether the common practice of identifying combinations of density functional methods and basis sets that have a good performance against a selected set of reference data should be considered as data fitting in the combined parameter space spanned by the method and basis set.
Collapse
Affiliation(s)
- Frank Jensen
- Department of Chemistry , Aarhus University , DK-8000 Aarhus , Denmark
| |
Collapse
|
8
|
Di(N-heterocyclic carbene) gold(III) imidate complexes obtained by oxidative addition of N-halosuccinimides. J Organomet Chem 2018. [DOI: 10.1016/j.jorganchem.2018.03.014] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
|
9
|
Guzzo RN, Rezende MJC, Kartnaller V, Carneiro JWDM, Stoyanov SR, Costa LMD. Experimental and DFT evaluation of the 1H and 13C NMR chemical shifts for calix[4]arenes. J Mol Struct 2018. [DOI: 10.1016/j.molstruc.2017.12.038] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
10
|
Synthesis and structure elucidation of allyl Pd(II) complexes of NHC ligands derived from substituted imidazo[1,5-a]quinolin-1(2H)-ylidene. J Organomet Chem 2017. [DOI: 10.1016/j.jorganchem.2017.09.036] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
|
11
|
Mota AJ, Neuhold J, Drescher M, Lemouzy S, González L, Maulide N. Intramolecular hydrogen bonding in conformationally semi-rigid α-acylmethane derivatives: a theoretical NMR study. Org Biomol Chem 2017; 15:7572-7579. [PMID: 28858370 DOI: 10.1039/c7ob01834g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Conformational mobility is a core property of organic compounds, and conformational analysis has become a pervasive tool for synthetic design. In this work, we present experimental and computational (employing Density Functional Theory) evidence for unusual intramolecular hydrogen bonding interactions in a series of α-acylmethane derivatives, as well as a discussion of the consequences thereof for their NMR spectroscopic properties.
Collapse
Affiliation(s)
- Antonio J Mota
- Department of Inorganic Chemistry, Faculty of Sciences, Avda. Fuentenueva s/n, 18002 Granada, Spain.
| | | | | | | | | | | |
Collapse
|
12
|
Fritz M, Quinn CM, Wang M, Hou G, Lu X, Koharudin LMI, Polenova T, Gronenborn AM. Toward Closing the Gap: Quantum Mechanical Calculations and Experimentally Measured Chemical Shifts of a Microcrystalline Lectin. J Phys Chem B 2017; 121:3574-3585. [PMID: 28001418 PMCID: PMC5465307 DOI: 10.1021/acs.jpcb.6b09479] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
NMR chemical shifts are exquisitely sensitive probes for conformation and dynamics in molecules and supramolecular assemblies. Although isotropic chemical shifts are easily measured with high accuracy and precision in conventional NMR experiments, they remain challenging to calculate quantum mechanically, particularly in inherently dynamic biological systems. Using a model benchmark protein, the 133-residue agglutinin from Oscillatoria agardhii (OAA), which has been extensively characterized by us previously, we have explored the integration of X-ray crystallography, solution NMR, MAS NMR, and quantum mechanics/molecular mechanics (QM/MM) calculations for analysis of 13Cα and 15NH isotropic chemical shifts. The influence of local interactions, quaternary contacts, and dynamics on the accuracy of calculated chemical shifts is analyzed. Our approach is broadly applicable and expected to be beneficial in chemical shift analysis and chemical-shift-based structure refinement for proteins and protein assemblies.
Collapse
Affiliation(s)
- Matthew Fritz
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
- Pittsburgh center for HIV Protein Interactions, University of Pittsburgh School of Medicine, 1051 Biomedical Science Tower 3, 3501 Fifth Ave., Pittsburgh, PA 15261, United States
| | - Caitlin M. Quinn
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
- Pittsburgh center for HIV Protein Interactions, University of Pittsburgh School of Medicine, 1051 Biomedical Science Tower 3, 3501 Fifth Ave., Pittsburgh, PA 15261, United States
| | - Mingzhang Wang
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
- Pittsburgh center for HIV Protein Interactions, University of Pittsburgh School of Medicine, 1051 Biomedical Science Tower 3, 3501 Fifth Ave., Pittsburgh, PA 15261, United States
| | - Guangjin Hou
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| | - Xingyu Lu
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
- Pittsburgh center for HIV Protein Interactions, University of Pittsburgh School of Medicine, 1051 Biomedical Science Tower 3, 3501 Fifth Ave., Pittsburgh, PA 15261, United States
| | - Leonardus M. I. Koharudin
- Pittsburgh center for HIV Protein Interactions, University of Pittsburgh School of Medicine, 1051 Biomedical Science Tower 3, 3501 Fifth Ave., Pittsburgh, PA 15261, United States
- Department of Structural Biology, University of Pittsburgh School of Medicine,3501 Fifth Ave., Pittsburgh, PA 15261, United States
| | - Tatyana Polenova
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
- Pittsburgh center for HIV Protein Interactions, University of Pittsburgh School of Medicine, 1051 Biomedical Science Tower 3, 3501 Fifth Ave., Pittsburgh, PA 15261, United States
| | - Angela M. Gronenborn
- Pittsburgh center for HIV Protein Interactions, University of Pittsburgh School of Medicine, 1051 Biomedical Science Tower 3, 3501 Fifth Ave., Pittsburgh, PA 15261, United States
- Department of Structural Biology, University of Pittsburgh School of Medicine,3501 Fifth Ave., Pittsburgh, PA 15261, United States
| |
Collapse
|
13
|
Abstract
We propose a method to quantify positional uncertainties in crystal structures determined by chemical-shift-based NMR crystallography. The method combines molecular dynamics simulations and density functional theory calculations with experimental and computational chemical shift uncertainties. In this manner we find the average positional accuracy as well as the isotropic and anisotropic positional accuracy associated with each atom in a crystal structure determined by NMR crystallography. The approach is demonstrated on the crystal structures of cocaine, flutamide, flufenamic acid, the K salt of penicillin G, and form 4 of the drug 4-[4-(2-adamantylcarbamoyl)-5-tert-butylpyrazol-1-yl]benzoic acid (AZD8329). We find that, for the crystal structure of cocaine, the uncertainty corresponds to a positional RMSD of 0.17 Å. This is a factor of 2.5 less than for single-crystal X-ray-diffraction-based structure determination.
Collapse
Affiliation(s)
- Albert Hofstetter
- Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL) , CH-1015 Lausanne, Switzerland
| | - Lyndon Emsley
- Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL) , CH-1015 Lausanne, Switzerland
| |
Collapse
|
14
|
Dangalov M, Petrov P, Vassilev NG. Fluxional allyl Pd(II) and Pt(II) complexes of NHC ligands derived from substituted 1,8-naphthalimides - Synthesis and structure elucidation. J Organomet Chem 2016. [DOI: 10.1016/j.jorganchem.2016.10.010] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
|
15
|
Holmes ST, Iuliucci RJ, Mueller KT, Dybowski C. Critical Analysis of Cluster Models and Exchange-Correlation Functionals for Calculating Magnetic Shielding in Molecular Solids. J Chem Theory Comput 2015; 11:5229-41. [DOI: 10.1021/acs.jctc.5b00752] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sean T. Holmes
- Department
of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| | - Robbie J. Iuliucci
- Department
of Chemistry, Washington and Jefferson College, Washington, Pennsylvania 15301, United States
| | - Karl T. Mueller
- Department
of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Cecil Dybowski
- Department
of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|