Computational protocols for prediction of solute NMR relative chemical shifts. A case study of L‐tryptophan in aqueous solution

In this study, we have applied two different spanning protocols for obtaining the molecular conformations of L‐tryptophan in aqueous solution, namely a molecular dynamics simulation and a molecular mechanics conformational search with subsequent geometry re‐optimization of the stable conformers using a quantum mechanically based method. These spanning protocols represent standard ways of obtaining a set of conformations on which NMR calculations may be performed. The results stemming from the solute–solvent configurations extracted from the MD simulation at 300 K are found to be inferior to the results stemming from the conformations extracted from the MM conformational search in terms of replicating an experimental reference as well as in achieving the correct sequence of the NMR relative chemical shifts of L‐tryptophan in aqueous solution. We find this to be due to missing conformations visited during the molecular dynamics run as well as inaccuracies in geometrical parameters generated from the classical molecular dynamics simulations. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011

[1]  Raymond J. Abraham,et al.  1H chemical shifts in NMR. Part 18. Ring currents and π-electron effects in hetero-aromatics , 2002 .

[2]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[3]  Kurt V. Mikkelsen,et al.  Coupled Cluster/Molecular Mechanics Method: Implementation and Application to Liquid Water , 2003 .

[4]  J. Olsen,et al.  Linear and nonlinear response functions for an exact state and for an MCSCF state , 1985 .

[5]  R. Buckingham,et al.  The Classical Equation of State of Gaseous Helium, Neon and Argon , 1938 .

[6]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[7]  Trygve Helgaker,et al.  Gaussian basis sets for high-quality ab initio calculations , 1988 .

[8]  Anders Wallqvist,et al.  A molecular dynamics study of polarizable water , 1989 .

[9]  M. Monduzzi,et al.  Proton and carbon-13 nuclear magnetic resonance spectra and INDO study of serotonine, tryptamine and L-tryptophan , 1983 .

[10]  Frank Herman,et al.  Improved Statistical Exchange Approximation for Inhomogeneous Many-Electron Systems , 1969 .

[11]  Kurt V. Mikkelsen,et al.  The combined multiconfigurational self-consistent-field/molecular mechanics wave function approach , 2001 .

[12]  Henry Margenau,et al.  Theory of intermolecular forces , 1969 .

[13]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[14]  Jacob Kongsted,et al.  Solvatochromic Shifts in Uracil: A Combined MD-QM/MM Study. , 2010, Journal of chemical theory and computation.

[15]  Antonio Rizzo,et al.  Computational study of the one- and two-photon absorption and circular dichroism of (L)-tryptophan. , 2010, The journal of physical chemistry. B.

[16]  Kenneth Levenberg.  A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

[17]  Peter A. Kollman,et al.  AMBER: Assisted model building with energy refinement. A general program for modeling molecules and their interactions , 1981 .

[18]  Jacob Kongsted,et al.  Excited States in Solution through Polarizable Embedding , 2010 .

[19]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[20]  Alessandro Bagno,et al.  Prediction of the 1H and 13C NMR spectra of alpha-D-glucose in water by DFT methods and MD simulations. , 2007, The Journal of organic chemistry.

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

[22]  Dieter Cremer,et al.  Sum‐over‐states density functional perturbation theory: Prediction of reliable 13C, 15N, and 17O nuclear magnetic resonance chemical shifts , 1996 .

[23]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[24]  Jacob Kongsted,et al.  On the existence of the H3 tautomer of adenine in aqueous solution. Rationalizations based on hybrid quantum mechanics/molecular mechanics predictions. , 2010, Physical chemistry chemical physics : PCCP.

[25]  Kenneth Ruud,et al.  Nuclear magnetic shielding constants of liquid water: insights from hybrid quantum mechanics/molecular mechanics models. , 2007, The Journal of chemical physics.

[26]  J. Tomasi,et al.  Quantum mechanical continuum solvation models. , 2005, Chemical reviews.

[27]  Frank Jensen,et al.  The basis set convergence of the Hartree–Fock energy for H3+, Li2 and N2 , 2000 .

[28]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[29]  Kenneth Ruud,et al.  Solvent effects on NMR isotropic shielding constants. a comparison between explicit polarizable discrete and continuum approaches. , 2007, The journal of physical chemistry. A.

[30]  V. Barone,et al.  Toward reliable density functional methods without adjustable parameters: The PBE0 model , 1999 .

[31]  P. Dirac.  Note on Exchange Phenomena in the Thomas Atom , 1930, Mathematical Proceedings of the Cambridge Philosophical Society.

[32]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[33]  P. Hohenberg,et al.  Inhomogeneous Electron Gas , 1964 .

[34]  Norman L. Allinger,et al.  The MMP2 calculational method , 1987 .

[35]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[36]  Norman L. Allinger,et al.  Molecular mechanics. The MM3 force field for hydrocarbons. 1 , 1989 .

[37]  H. Friebolin,et al.  Basic one- and two-dimensional NMR spectroscopy , 1991 .

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

[39]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[40]  N. Handy,et al.  Left-right correlation energy , 2001 .

[41]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

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

[43]  Paweł Sałek,et al.  Linear-scaling formation of Kohn-Sham Hamiltonian: application to the calculation of excitation energies and polarizabilities of large molecular systems. , 2004, The Journal of chemical physics.

[44]  Thomas W Keal,et al.  A semiempirical generalized gradient approximation exchange-correlation functional. , 2004, The Journal of chemical physics.

[45]  Kurt V. Mikkelsen,et al.  The QM/MM approach for wavefunctions, energies and response functions within self-consistent field and coupled cluster theories , 2002 .

[46]  G. Scuseria,et al.  Assessment of the Perdew–Burke–Ernzerhof exchange-correlation functional , 1999 .

[47]  A. Schäfer,et al.  Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr , 1994 .

[48]  Elfi Kraka,et al.  Homotropenylium Cation: Structure, Stability, and Magnetic Properties , 1991 .

[49]  S. Lifson,et al.  Consistent force field studies of intermolecular forces in hydrogen-bonded crystals. 1. Carboxylic acids, amides, and the C:O.cntdot..cntdot..cntdot.H- hydrogen bonds , 1979 .

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

[51]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[52]  Cynthia J. Jameson.  Gas-phase NMR spectroscopy , 1991 .

[53]  Vincenzo Barone,et al.  Spectroscopic properties in the liquid phase: combining high-level ab initio calculations and classical molecular dynamics. , 2006, Chemphyschem : a European journal of chemical physics and physical chemistry.

[54]  Graham Webb,et al.  Theory of NMR parameters , 1983 .

[55]  W. L. Jorgensen.  Quantum and statistical mechanical studies of liquids. 10. Transferable intermolecular potential functions for water, alcohols, and ethers. Application to liquid water , 2002 .

[56]  S. Goedecker.  Linear scaling electronic structure methods , 1999 .

[57]  Trygve Helgaker,et al.  GIAO shielding constants and indirect spin–spin coupling constants: performance of density functional methods , 2004 .

[58]  Jacob Kongsted,et al.  Density functional self-consistent quantum mechanics/molecular mechanics theory for linear and nonlinear molecular properties: Applications to solvated water and formaldehyde. , 2007, The Journal of chemical physics.

[59]  M. Bühl,et al.  Calculation of NMR and EPR parameters : theory and applications , 2004 .

[60]  Paul von Ragué Schleyer,et al.  Application and evaluation of ab initio chemical shift calculations for boranes and carboranes. How reliable are "accurate" experimental structures? , 1992 .

[61]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[62]  Dennis R. Salahub,et al.  NUCLEAR MAGNETIC RESONANCE SHIELDING TENSORS CALCULATED WITH A SUM-OVER-STATES DENSITY FUNCTIONAL PERTURBATION THEORY , 1994 .

[63]  Jacob Kongsted,et al.  On the Accuracy of Density Functional Theory to Predict Shifts in Nuclear Magnetic Resonance Shielding Constants due to Hydrogen Bonding. , 2008, Journal of chemical theory and computation.

[64]  Norman L. Allinger,et al.  Molecular mechanics. The MM3 force field for hydrocarbons. 2. Vibrational frequencies and thermodynamics , 1989 .

[65]  Paolo Lazzeretti,et al.  Gauge invariant calculations of nuclear magnetic shielding constants using the continuous transformation of the origin of the current density approach. II. Density functional and coupled cluster theory. , 2007, The Journal of chemical physics.

[66]  Paweł Sałek,et al.  Density-functional theory of linear and nonlinear time-dependent molecular properties , 2002 .

[67]  Stephan P. A. Sauer,et al.  Molecular Electromagnetism: A Computational Chemistry Approach , 2011 .

[68]  Jenn-Huei Lii,et al.  The MM3 force field for amides, polypeptides and proteins , 1991 .

[69]  Norman L. Allinger,et al.  Molecular mechanics. The MM3 force field for hydrocarbons. 3. The van der Waals' potentials and crystal data for aliphatic and aromatic hydrocarbons , 1989 .

[70]  Frank Jensen,et al.  The basis set convergence of the density functional energy for H2 , 2000 .

[71]  David J. Tozer,et al.  Improved NMR chemical shifts in density functional theory , 2003 .

[72]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[73]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[74]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[75]  Hans W. Horn,et al.  Fully optimized contracted Gaussian basis sets for atoms Li to Kr , 1992 .