A comparison of quantum chemical models for calculating NMR shielding parameters in peptides: Mixed basis set and ONIOM methods combined with a complete basis set extrapolation

This article compares several quantum mechanical approaches to the computation of chemical shielding tensors in peptide fragments. First, we describe the effects of basis set quality up to the complete basis set (CBS) limit and level of theory (HF, MP2, and DFT) for four different atoms in trans N‐methylacetamide. For both isotropic shielding and shielding anisotropy, the MP2 results in the CBS limit show the best agreement with experiment. The HF values show quite a different tendency to MP2, and even in the CBS limit they are far from experiment for not only the isotropic shielding of carbonyl carbon but also most shielding anisotropies. In most cases, the DFT values differ systematically from MP2, and small basis‐set (double‐ or triple‐zeta) results are often fortuitously in better agreement with the experiment than the CBS ones. Second, we compare the mixed basis set and ONIOM methods, combined with CBS extrapolation, for chemical shielding calculations at a DFT level using various model peptides. From the results, it is shown that the mixed basis set method provides better results than ONIOM, compared to CBS calculations using the nonpartitioned full systems. The information studied here will be useful in guiding the selection of proper quantum chemical models, which are in a tradeoff between accuracy and cost, for shielding studies of peptides and proteins. © 2006 Wiley Periodicals, Inc. J Comput Chem 27: 825–836, 2006

[1]  P. Markwick,et al.  Site-specific variations of carbonyl chemical shift anisotropies in proteins. , 2004, Journal of the American Chemical Society.

[2]  P. A. Christiansen,et al.  Quantum capping potentials with point charges: a simple QM/MM approach for the calculation of large-molecule NMR shielding tensors. , 2004, The Journal of chemical physics.

[3]  András Perczel,et al.  Toward direct determination of conformations of protein building units from multidimensional NMR experiments. V. NMR chemical shielding analysis of N‐formyl‐serinamide, a model for polar side‐chain containing peptides , 2003, J. Comput. Chem..

[4]  B. Ruscic,et al.  Hartree-Fock and density functional complete basis-set (CBS) predicted nuclear shielding anisotropy and shielding tensor components. , 2003, Solid state nuclear magnetic resonance.

[5]  David A. Case,et al.  Probing multiple effects on 15N, 13Cα, 13Cβ, and 13C′ chemical shifts in peptides using density functional theory , 2002 .

[6]  B. Ruscic,et al.  Toward Hartree-Fock and density functional complete basis set predicted NMR parameters. , 2002 .

[7]  Margaret M. Hurley,et al.  Simple one-electron quantum capping potentials for use in hybrid QM/MM studies of biological molecules , 2002 .

[8]  Eric Oldfield,et al.  Carbon-13 NMR shielding in the twenty common amino acids: comparisons with experimental results in proteins. , 2002, Journal of the American Chemical Society.

[9]  Ulrich Sternberg,et al.  Chemical shift driven geometry optimization , 2002, J. Comput. Chem..

[10]  B. Brooks,et al.  An ab initio study of amide proton shift tensor dependence on local protein structure. , 2002, Journal of the American Chemical Society.

[11]  D. Case,et al.  Use of chemical shifts in macromolecular structure determination. , 2002, Methods in enzymology.

[12]  D. Case,et al.  Automated prediction of 15N, 13Cα, 13Cβ and 13C′ chemical shifts in proteins using a density functional database , 2001, Journal of biomolecular NMR.

[13]  J H Prestegard,et al.  A dipolar coupling based strategy for simultaneous resonance assignment and structure determination of protein backbones. , 2001, Journal of the American Chemical Society.

[14]  D E Wemmer,et al.  An experimental and theoretical investigation of the chemical shielding tensors of (13)C(alpha) of alanine, valine, and leucine residues in solid peptides and in proteins in solution. , 2001, Journal of the American Chemical Society.

[15]  J H Prestegard,et al.  Partial alignment of biomolecules: an aid to NMR characterization. , 2001, Current opinion in chemical biology.

[16]  A. Császár,et al.  Toward direct determination of conformations of protein building units from multidimensional NMR experiments part II: a theoretical case study of formyl-L-valine amide. , 2001, Chemistry.

[17]  J. Brender,et al.  Orientation of amide-nitrogen-15 chemical shift tensors in peptides: a quantum chemical study. , 2001, Journal of the American Chemical Society.

[18]  Nathaniel R. Luman,et al.  Predicting 15N amide chemical shifts in proteins. I. An additive model for the backbone contribution , 2001, J. Comput. Chem..

[19]  A. Bax,et al.  Measurement of Proton, Nitrogen, and Carbonyl Chemical Shielding Anisotropies in a Protein Dissolved in a Dilute Liquid Crystalline Phase , 2000 .

[20]  S. Sugai,et al.  Molten globule structure of equine beta-lactoglobulin probed by hydrogen exchange. , 2000, Journal of molecular biology.

[21]  F. Castellino,et al.  Structure–function relationships of the NMDA receptor antagonist peptide, conantokin‐R , 2000, FEBS letters.

[22]  N. Handy,et al.  Density functional predictions for magnetizabilities and nuclear shielding constants , 1999 .

[23]  Kenneth B. Wiberg,et al.  Comparison of density functional theory models' ability to reproduce experimental 13C‐NMR shielding values , 1999, J. Comput. Chem..

[24]  Trygve Helgaker,et al.  Ab Initio Methods for the Calculation of NMR Shielding and Indirect Spin-Spin Coupling Constants , 1999 .

[25]  Fred A. Hamprecht,et al.  Development and assessment of new exchange-correlation functionals , 1998 .

[26]  Vincenzo Barone,et al.  Exchange functionals with improved long-range behavior and adiabatic connection methods without adjustable parameters: The mPW and mPW1PW models , 1998 .

[27]  David A. Case,et al.  Density Functional Calculations of Proton Chemical Shifts in Model Peptides , 1997 .

[28]  H. Dyson,et al.  Chemical shift dispersion and secondary structure prediction in unfolded and partly folded proteins , 1997, FEBS letters.

[29]  E. Oldfield,et al.  PREDICTING CHEMICAL SHIFTS IN PROTEINS : STRUCTURE REFINEMENT OF VALINE RESIDUES BY USING AB INITIO AND EMPIRICAL GEOMETRY OPTIMIZATIONS , 1997 .

[30]  A. Bax,et al.  MAGNETIC FIELD DEPENDENT AMIDE 15N CHEMICAL SHIFTS IN A PROTEIN-DNA COMPLEX RESULTING FROM MAGNETIC ORDERING IN SOLUTION , 1997 .

[31]  A. Dios,et al.  Chemical Shift Tensors in Peptides: A Quantum Mechanical Study , 1997 .

[32]  E. Oldfield,et al.  Determination of Dihedral Angles in Peptides through Experimental and Theoretical Studies of α-Carbon Chemical Shielding Tensors , 1997 .

[33]  K. Burke,et al.  Generalized Gradient Approximation Made Simple [Phys. Rev. Lett. 77, 3865 (1996)] , 1997 .

[34]  Wang,et al.  Generalized gradient approximation for the exchange-correlation hole of a many-electron system. , 1996, Physical review. B, Condensed matter.

[35]  K. Burke,et al.  Rationale for mixing exact exchange with density functional approximations , 1996 .

[36]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[37]  E. Oldfield,et al.  AB INITIO STUDIES OF AMIDE-15N CHEMICAL SHIFTS IN DIPEPTIDES : APPLICATIONS TO PROTEIN NMR SPECTROSCOPY , 1996 .

[38]  L Mayne,et al.  Mechanisms and uses of hydrogen exchange. , 1996, Current opinion in structural biology.

[39]  Feliu Maseras,et al.  IMOMM: A new integrated ab initio + molecular mechanics geometry optimization scheme of equilibrium structures and transition states , 1995, J. Comput. Chem..

[40]  E. Oldfield,et al.  A Basis Size Dependence Study of Carbon-13 Nuclear Magnetic Resonance Spectroscopic Shielding in Alanyl and Valyl Fragments: Toward Protein Shielding Hypersurfaces , 1995 .

[41]  E. Oldfield,et al.  Chemical Shifts of Carbonyl Carbons in Peptides and Proteins , 1994 .

[42]  E. Oldfield,et al.  Ab Initio Study of the Effects of Torsion Angles on Carbon-13 Nuclear Magnetic Resonance Chemical Shielding in N-Formyl-L-alanine Amide, N-Formyl-L-valine Amide, and Some Simple Model Compounds: Applications to Protein NMR Spectroscopy , 1994 .

[43]  M. Schindler,et al.  Carbonyl carbon and nitrogen chemical shift tensors of the amide fragment of acetanilide and N-methylacetanilide , 1994 .

[44]  D. B. Chesnut,et al.  Use of locally dense basis sets for nuclear magnetic resonance shielding calculations , 1993, J. Comput. Chem..

[45]  Eric Oldfield,et al.  NMR chemical shifts and structure refinement in proteins , 1993, Journal of biomolecular NMR.

[46]  E. Oldfield,et al.  Secondary and tertiary structural effects on protein NMR chemical shifts: an ab initio approach. , 1993, Science.

[47]  A. Becke Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[48]  Charles J. Pouchert,et al.  The Aldrich library of [13]C and [1]H FT NMR spectra , 1993 .

[49]  S. Kuroki,et al.  Hydrogen-bonding effect on carbon-13 NMR chemical shifts of L-alanine residue carbonyl carbons of peptides in the solid state , 1992 .

[50]  C. Jameson Gas‐Phase NMR Spectroscopy , 1991 .

[51]  T. M. Duncan,et al.  A compilation of chemical shift anisotropies , 1990 .

[52]  D. B. Chesnut,et al.  Locally dense basis sets for chemical shift calculations , 1989 .

[53]  David Neuhaus,et al.  The Nuclear Overhauser Effect in Structural and Conformational Analysis , 1989 .

[54]  T. H. Dunning Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[55]  K. Wüthrich Protein structure determination in solution by nuclear magnetic resonance spectroscopy. , 1989, Science.

[56]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[57]  D. Torchia,et al.  Determination of 15N chemical shift tensor via 15N-2H dipolar coupling in Boc-glycylglycyl[15N]glycine benzyl ester , 1988 .

[58]  F. Dahlquist,et al.  The Amide 15N Chemical Shift Tensors of Four Peptides Determined from 13C Dipole-Coupled Chemical Shift Powder Patterns , 1987 .

[59]  K. Wüthrich NMR of proteins and nucleic acids , 1988 .

[60]  U. Singh,et al.  A combined ab initio quantum mechanical and molecular mechanical method for carrying out simulations on complex molecular systems: Applications to the CH3Cl + Cl− exchange reaction and gas phase protonation of polyethers , 1986 .

[61]  J. Perdew,et al.  Density-functional approximation for the correlation energy of the inhomogeneous electron gas. , 1986, Physical review. B, Condensed matter.

[62]  Warren J. Hehre,et al.  AB INITIO Molecular Orbital Theory , 1986 .

[63]  M. Billeter,et al.  Calibration of the angular dependence of the amide proton-Cα proton coupling constants, 3JHNα, in a globular protein: Use of 3JHNα for identification of helical secondary structure , 1984 .

[64]  D. Torchia,et al.  15N chemical shift and 15N-13C dipolar tensors for the peptide bond in [1-13C]glycyl[15N]glycine hydrochloride monohydrate , 1984 .

[65]  R. Ditchfield,et al.  Self-consistent perturbation theory of diamagnetism , 1974 .

[66]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .