Using quantum chemistry to estimate chemical shifts in biomolecules.
Ontology highlight
ABSTRACT:
SUBMITTER: Case DA
PROVIDER: S-EPMC7686263 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:We evaluate the performance of the automated fragmentation quantum mechanics/molecular mechanics approach (AF-QM/MM) on the calculation of protein and nucleic acid NMR chemical shifts. The AF-QM/MM approach models solvent effects implicitly through a set of surface charges computed using the Poisson-Boltzmann equation, and it can also be combined with an explicit solvent model through the placement of water molecules in the first solvation shell around the solute; the latter substantially improves the accuracy of chemical shift prediction of protons involved in hydrogen bonding with solvent. We also compare the performance of AF-QM/MM on proteins and nucleic acids with two leading empirical chemical shift prediction programs SHIFTS and SHIFTX2. Although the empirical programs outperform AF-QM/MM in predicting chemical shifts, the differences are in some cases small, and the latter can be applied to chemical shifts on biomolecules which are outside the training set employed by the empirical programs, such as structures containing ligands, metal centers, and non-standard residues. The AF-QM/MM described here is implemented in version 5 of the SHIFTS software, and is fully automated, so that only a structure in PDB format is required as input.
| S-EPMC6556433 | biostudies-literature
Project description:We present the ProCS method for the rapid and accurate prediction of protein backbone amide proton chemical shifts--sensitive probes of the geometry of key hydrogen bonds that determine protein structure. ProCS is parameterized against quantum mechanical (QM) calculations and reproduces high level QM results obtained for a small protein with an RMSD of 0.25 ppm (r = 0.94). ProCS is interfaced with the PHAISTOS protein simulation program and is used to infer statistical protein ensembles that reflect experimentally measured amide proton chemical shift values. Such chemical shift-based structural refinements, starting from high-resolution X-ray structures of Protein G, ubiquitin, and SMN Tudor Domain, result in average chemical shifts, hydrogen bond geometries, and trans-hydrogen bond ((h3)J(NC')) spin-spin coupling constants that are in excellent agreement with experiment. We show that the structural sensitivity of the QM-based amide proton chemical shift predictions is needed to obtain this agreement. The ProCS method thus offers a powerful new tool for refining the structures of hydrogen bonding networks to high accuracy with many potential applications such as protein flexibility in ligand binding.
| S-EPMC3877219 | biostudies-literature
Project description:Diamagnetic chemical exchange saturation transfer (CEST) contrast agents offer an alternative to Gd(3+) -based contrast agents for MRI. They are characterized by containing protons that can rapidly exchange with water and it is advantageous to have these protons resonate in a spectral window that is far removed from water. Herein, we report the first results of DFT calculations of the (1) H nuclear magnetic shieldings in 41 CEST agents, finding that the experimental shifts can be well predicted (R(2) =0.882). We tested a subset of compounds with the best MRI properties for toxicity and for activity as uncouplers, then obtained mice kidney CEST MRI images for three of the most promising leads finding 16 (2,4-dihydroxybenzoic acid) to be one of the most promising CEST MRI contrast agents to date. Overall, the results are of interest since they show that (1) H NMR shifts for CEST agents-charged species-can be well predicted, and that several leads have low toxicity and yield good in vivo MR images.
| S-EPMC4715718 | biostudies-literature
Project description:We present the results of the first quantum chemical investigations of 1H NMR hyperfine shifts in the blue copper proteins (BCPs): amicyanin, azurin, pseudoazurin, plastocyanin, stellacyanin, and rusticyanin. We find that very large structural models that incorporate extensive hydrogen bond networks, as well as geometry optimization, are required to reproduce the experimental NMR hyperfine shift results, the best theory vs experiment predictions having R2 = 0.94, a slope = 1.01, and a SD = 40.5 ppm (or approximately 4.7% of the overall approximately 860 ppm shift range). We also find interesting correlations between the hyperfine shifts and the bond and ring critical point properties computed using atoms-in-molecules theory, in addition to finding that hyperfine shifts can be well-predicted by using an empirical model, based on the geometry-optimized structures, which in the future should be of use in structure refinement.
| S-EPMC2761745 | biostudies-literature
Project description: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%.
| S-EPMC9499888 | biostudies-literature
Project description:Interaction of biomolecules typically proceeds in a highly selective and reversible manner, for which covalent bond formation has been largely avoided due to the potential difficulty of dissociation. However, employing reversible covalent warheads in drug design has given rise to covalent enzyme inhibitors that serve as powerful therapeutics, as well as molecular probes with exquisite target selectivity. This review article summarizes the recent advances in the development of reversible covalent chemistry for biological and medicinal applications. Specifically, we document the chemical strategies that allow for reversible modification of the three major classes of nucleophiles in biology: thiols, alcohols and amines. Emphasis is given to the chemical mechanisms that underlie the development of these reversible covalent reactions and their utilization in biology.
| S-EPMC5107367 | biostudies-literature
Project description:We developed an open-source chemical reaction equilibrium solver in Python (CASpy, https://github.com/omoultosEthTuDelft/CASpy) to compute the concentration of species in any reactive liquid-phase absorption system. We derived an expression for a mole fraction-based equilibrium constant as a function of excess chemical potential, standard ideal gas chemical potential, temperature, and volume. As a case study, we computed the CO2 absorption isotherm and speciation in a 23 wt % N-methyldiethanolamine (MDEA)/water solution at 313.15 K, and compared the results with available data from the literature. The results show that the computed CO2 isotherms and speciations are in excellent agreement with experimental data, demonstrating the accuracy and the precision of our solver. The binary absorptions of CO2 and H2S in 50 wt % MDEA/water solutions at 323.15 K were computed and compared with available data from the literature. The computed CO2 isotherms showed good agreement with other modeling studies from the literature while the computed H2S isotherms did not agree well with experimental data. The experimental equilibrium constants used as an input were not adjusted for H2S/CO2/MDEA/water systems and need to be adjusted for this system. Using free energy calculations with two different force fields (GAFF and OPLS-AA) and quantum chemistry calculations, we computed the equilibrium constant (K) of the protonated MDEA dissociation reaction. Despite the good agreement of the OPLS-AA force field (ln[K] = -24.91) with the experiments (ln[K] = -23.04), the computed CO2 pressures were significantly underestimated. We systematically investigated the limitations of computing CO2 absorption isotherms using free energy and quantum chemistry calculations and showed that the computed values of μiex are very sensitive to the point charges used in the simulations, which limits the predictive power of this method.
| S-EPMC10173452 | biostudies-literature
Project description:A recently popularized approach for the calculation of pseudocontact shifts (PCSs) based on first-principles quantum chemistry (QC) leads to different results than the classic "semiempirical" equation involving the susceptibility tensor. Studies that attempted a comparison of theory and experiment led to conflicting conclusions with respect to the preferred theoretical approach. In this Letter, we show that after inclusion of previously neglected terms in the full Hamiltonian, one can deduce the semiempirical equations from a rigorous QC-based treatment. It also turns out that in the long-distance limit, one can approximate the complete A tensor in terms of the g tensor. By means of Kohn-Sham density functional theory calculations, we numerically confirm the long-distance expression for the A tensor and the theoretically predicted scaling behavior of the different terms. Our derivation suggests a computational strategy in which one calculates the susceptibility tensor and inserts it into the classic equation for the PCS.
| S-EPMC7584370 | biostudies-literature
Project description:Reaction times, activation energies, branching ratios, yields, and many other quantitative attributes are important for precise organic syntheses and generating detailed reaction mechanisms. Often, it would be useful to be able to classify proposed reactions as fast or slow. However, quantitative chemical reaction data, especially for atom-mapped reactions, are difficult to find in existing databases. Therefore, we used automated potential energy surface exploration to generate 12,000 organic reactions involving H, C, N, and O atoms calculated at the ωB97X-D3/def2-TZVP quantum chemistry level. We report the results of geometry optimizations and frequency calculations for reactants, products, and transition states of all reactions. Additionally, we extracted atom-mapped reaction SMILES, activation energies, and enthalpies of reaction. We believe that this data will accelerate progress in automated methods for organic synthesis and reaction mechanism generation-for example, by enabling the development of novel machine learning models for quantitative reaction prediction.
| S-EPMC7210263 | biostudies-literature
Project description:Amyloid fibrils are the pathological hallmark of a large variety of neurodegenerative disorders. The structural characterization of amyloid fibrils, however, is challenging due to their non-crystalline, heterogeneous, and often dynamic nature. Thus, the structure of amyloid fibrils of many proteins is still unknown. We here show that the structure calculation program CS-Rosetta can be used to obtain insight into the core structure of amyloid fibrils. Driven by experimental solid-state NMR chemical shifts and taking into account the polymeric nature of fibrils CS-Rosetta allows modeling of the core of amyloid fibrils. Application to the Y145X stop mutant of the human prion protein reveals a left-handed ?-helix.
| S-EPMC3575924 | biostudies-literature