Methods for computing nuclear magnetic resonance chemical shielding in large systems. Multiple cluster and charge field approaches

Efficient Calculation of NMR Shielding Constants Using Composite Method Approximations and Locally Dense Basis Sets

This paper presents a systematic study of applying composite method approximations with locally dense basis sets (LDBS) to efficiently calculate NMR shielding constants in small and medium-sized molecules. The pcSseg-n series of basis sets are shown to have similar accuracy to the pcS-n series when n ≥ 1 and can slightly reduce computational costs. We identify two different LDBS partition schemes that perform very effectively for density functional calculations. We select a large subset of the recent NS372 database containing 290 H, C, N, and O shielding values evaluated by reference methods on 106 molecules to carefully assess methods of the high, medium, and low computational costs to make practical recommendations. Our assessment covers conventional electronic structure methods (density functional theory and wave function) with global basis calculations, as well as their use in one of the satisfactory LDBS approaches, and a range of composite approaches, also with and without LDBS. Altogether 99 methods are evaluated. On this basis, we recommend different methods to reach three different levels of accuracy and time requirements across the four nuclei considered.

Automated fragmentation quantum mechanical calculation of 13C and 1H chemical shifts in molecular crystals

In this work, the automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach was applied to calculate the ¹³C and ¹H nuclear magnetic resonance (NMR) chemical shifts in molecular crystals. Two benchmark sets of molecular crystals were selected to calculate the NMR chemical shifts. Systematic investigation was conducted to examine the convergence of AF-QM/MM calculations and the impact of various density functionals with different basis sets on the NMR chemical shift prediction. The result demonstrates that the calculated NMR chemical shifts are close to convergence when the distance threshold for the QM region is larger than 3.5 Å. For ¹³C chemical shift calculations, the mPW1PW91 functional is the best density functional among the functionals chosen in this study (namely, B3LYP, B3PW91, M06-2X, M06-L, mPW1PW91, OB98, and OPBE), while the OB98 functional is more suitable for the ¹H NMR chemical shift prediction of molecular crystals. Moreover, with the B3LYP functional, at least a triple-ζ basis set should be utilized to accurately reproduce the experimental ¹³C and ¹H chemical shifts. The employment of diffuse basis functions will further improve the accuracy for ¹³C chemical shift calculations, but not for the ¹H chemical shift prediction. We further proposed a fragmentation scheme of dividing the central molecule into smaller fragments. By comparing with the results of the fragmentation scheme using the entire central molecule as the core region, the AF-QM/MM calculations with the fragmented central molecule can not only achieve accurate results but also reduce the computational cost. Therefore, the AF-QM/MM approach is capable of predicting the ¹³C and ¹H NMR chemical shifts for molecular crystals accurately and effectively, and could be utilized for dealing with more complex periodic systems such as macromolecular polymers and biomacromolecules. The AF-QM/MM program for molecular crystals is available at https://github.com/shiman1995/NMR.

Accurate and cost-effective NMR chemical shift predictions for proteins using a molecules-in-molecules fragmentation-based method

We have developed an efficient protocol using our two-layer Molecules-in-Molecules (MIM2) fragmentation-based quantum chemical method for the prediction of NMR chemical shifts of large biomolecules. To investigate the performance of our fragmentation approach and demonstrate its applicability, MIM-NMR calculations are first calibrated on a test set of six proteins. The MIM2-NMR method yields a mean absolute deviation (MAD) from unfragmented full molecule calculations of 0.01 ppm for 1H and 0.06 ppm for 13C chemical shifts. Thus, the errors from fragmentation are only about 3% of our target accuracy of ∼0.3 ppm for 1H and 2-3 ppm for 13C chemical shifts. To compare with experimental chemical shifts, a standard protocol is first derived using two smaller proteins 2LHY (176 atoms) and 2LI1 (146 atoms) for obtaining an appropriate protein structure for NMR chemical shift calculations. The effect of the solvent environment on the calculated NMR chemical shifts is incorporated through implicit, explicit, or explicit-implicit solvation models. The expensive first solvation shell calculations are replaced by a micro-solvation model in which only the immediate interaction between the protein and the explicit solvation environment is considered. A single explicit water molecule for each amine and amide proton is found to be sufficient to yield accurate results for 1H chemical shifts. The 1H and 13C NMR chemical shifts calculated using our protocol give excellent agreement with experiments for two larger proteins, 2MC5 (the helical part with 265 atoms) and 3UMK (33 residue slice with 547 atoms). Overall, our target accuracy of ∼0.3 ppm for 1H and ∼2-3 ppm for 13C has been achieved for the larger proteins. The proposed MIM-NMR method is accurate and computationally cost-effective and should be applicable to study a wide range of large proteins.

Calculating nuclear magnetic resonance chemical shifts in solvated systems

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.

Layer selection effect on solid state 13C and 15N chemical shifts calculation using ONIOM approach

Solid state (13)C and (15)N chemical shifts of uracil and imidazole have been calculated using a 2-layer ONIOM approach at 32 levels of theory. The effect of electron correlation between two layers has been investigated by choosing two different kinds of layer selection. Factorial design has been applied as a multivariate technique to analyze the effect of wave function and layer selection on solid state (13)C and (15)N chemical shifts calculations. PBEPBE/6-311+G(d,p) was recommended as an optimally selected level of theory for high layer in both models. It is illustrated that considering the electron correlation of two layers of ONIOM models is important factor to calculate solid state (15)N chemical shifts. The agreement between the calculated and experimental values of solid state (13)C and (15)N chemical shifts using ONIOM (PBEPBE/6-311+G(d,p):AM1) for both uracil and imidazole confirmed the reliability of the selected wave functions and layer selection.

Automated Fragmentation QM/MM Calculation of Amide Proton Chemical Shifts in Proteins with Explicit Solvent Model

We have performed a density functional theory (DFT) calculation of the amide proton NMR chemical shift in proteins using a recently developed automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach. Systematic investigation was carried out to examine the influence of explicit solvent molecules, cooperative hydrogen bonding effects, density functionals, size of the basis sets, and the local geometry of proteins on calculated chemical shifts. Our result demonstrates that the predicted amide proton ((1)HN) NMR chemical shift in explicit solvent shows remarkable improvement over that calculated with the implicit solvation model. The cooperative hydrogen bonding effect is also shown to improve the accuracy of (1)HN chemical shifts. Furthermore, we found that the OPBE exchange-correlation functional is the best density functional for the prediction of protein (1)HN chemical shifts among a selective set of DFT methods (namely, B3LYP, B3PW91, M062X, M06L, mPW1PW91, OB98, OPBE), and the locally dense basis set of 6-311++G**/4-31G* is shown to be sufficient for (1)HN chemical shift calculation. By taking ensemble averaging into account, (1)HN chemical shifts calculated by the AF-QM/MM approach can be used to validate the performance of various force fields. Our study underscores that the electronic polarization of protein is of critical importance to stabilizing hydrogen bonding, and the AF-QM/MM method is able to describe the local chemical environment in proteins more accurately than most widely used empirical models.

Characterization of solid polymer dispersions of active pharmaceutical ingredients by 19F MAS NMR and factor analysis

In this contribution the ability of (19)F MAS NMR spectroscopy to probe structural variability of poorly water-soluble drugs formulated as solid dispersions in polymer matrices is discussed. The application potentiality of the proposed approach is demonstrated on a moderately sized active pharmaceutical ingredient (API, Atorvastatin) exhibiting extensive polymorphism. In this respect, a range of model systems with the API incorporated in the matrix of polvinylpyrrolidone (PVP) was prepared. The extent of mixing of both components was determined by T(1)((1)H) and T(1ρ)((1)H) relaxation experiments, and it was found that the API forms nanosized domains. Subsequently it was found out that the polymer matrix induces two kinds of changes in (19)F MAS NMR spectra. At first, this is a high-frequency shift reaching 2-3 ppm which is independent on molecular structure of the API and which results from the long-range polarization of the electron cloud around (19)F nucleus induced by electrostatic fields of the polymer matrix. At second, this is broadening of the signals and formation of shoulders reflecting changes in molecular arrangement of the API. To avoid misleading in the interpretation of the recorded (19)F MAS NMR spectra, because both the contributions act simultaneously, we applied chemometric approach based on multivariate analysis. It is demonstrated that factor analysis of the recorded spectra can separate both these spectral contributions, and the subtle structural differences in the molecular arrangement of the API in the nanosized domains can be traced. In this way (19)F MAS NMR spectra of both pure APIs and APIs in solid dispersions can be directly compared. The proposed strategy thus provides a powerful tool for the analysis of new formulations of fluorinated pharmaceutical substances in polymer matrices.

Fragment density functional theory calculation of NMR chemical shifts for proteins with implicit solvation

Fragment density functional theory (DFT) calculation of NMR chemical shifts for several proteins (Trp-cage, Pin1 WW domain, the third IgG-binding domain of Protein G (GB3) and human ubiquitin) has been carried out. The present study is based on a recently developed automatic fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach but the solvent effects are included by using the PB (Poisson-Boltzmann) model. Our calculated chemical shifts of (1)H and (13)C for these four proteins are in excellent agreement with experimentally measured values and represent clear improvement over that from the gas phase calculation. However, although the inclusion of the solvent effect also improves the computed chemical shifts of (15)N, the results do not agree with experimental values as well as (1)H and (13)C. Our study also demonstrates that AF-QM/MM calculated results accurately reproduce the separation of α-helical and β-sheet chemical shifts for (13)C(α) atoms in proteins, and using the (1)H chemical shift to discriminate the native structure of proteins from decoys is quite remarkable.

Protein NMR chemical shift calculations based on the automated fragmentation QM/MM approach

An automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach has been developed to routinely calculate ab initio protein NMR chemical shielding constants. The AF-QM/MM method is linear-scaling and trivially parallel. A general fragmentation scheme is employed to generate each residue-centric region which is treated by quantum mechanics, and the environmental electrostatic field is described with molecular mechanics. The AF-QM/MM method shows good agreement with standard self-consistent field (SCF) calculations of the NMR chemical shieldings for the mini-protein Trp cage. The root-mean-square errors (RMSEs) for 1H, 13C, and 15N NMR chemical shieldings are equal to or less than 0.09, 0.32, and 0.78 ppm, respectively, for all Hartree-Fock (HF) and density functional theory (DFT) calculations reported in this work. The environmental electrostatic potential is necessary to accurately reproduce the NMR chemical shieldings using the AF-QM/MM approach. The point-charge models provided by AMBER, AM1/CM2, PM3/CM1, and PM3/CM2 all effectively model the electrostatic field. The latter three point-charge models are generated via semiempirical linear-scaling SCF calculations of the entire protein system. The correlations between experimental 1H NMR chemical shifts and theoretical predictions are >0.95 for AF-QM/MM calculations using B3LYP with the 6-31G**, 6-311G**, and 6-311++G** basis sets. Our study, not unexpectedly, finds that conformational changes within a protein structure play an important role in the accurate prediction of experimental NMR chemical shifts from theory.

Density functional calculations of chemical shielding of backbone 15N in helical residues of protein G

We performed density functional calculations of backbone (15)N chemical shielding tensors in selected helical residues of protein G. Here we describe a computationally efficient methodology to include most of the important effects in the calculation of chemical shieldings of backbone (15)N. We analyzed the role of long-range intra-protein electrostatic interactions by comparing models with different complexity in vacuum and in charge field. Our results show that the dipole moment of the alpha-helix can cause significant deshielding of (15)N; therefore, it needs to be considered when calculating (15)N chemical shielding. We found that it is important to include interactions with the side chains that are close in space when the charged form for ionizable side chains is adopted in the calculation. We also illustrate how the ionization state of these side chains can affect the chemical shielding tensor elements. Chemical shielding calculations using a 8-residue fragment model in vacuum and adopting the charged form of ionizable side chains yield a generally good agreement with experimental data.

Deciphering the NMR fingerprints of the disordered system with quantum chemical studies

Recent developments in solid-state NMR techniques helped acquire high-resolution NMR spectra for solid systems with structural disorder. But the structural origin of the observed chemical shift nonequivalence in these systems has not been revealed. We report a quantum chemical investigation of the solid-state NMR spectrum in N,N-bis(diphenylphosphino)-N-((S)-alpha-methylbenzyl)amine, where eight nonequivalent (31)P NMR chemical shifts were resolved with a range of 13.0 ppm. Results from using different quantum chemical methods, computational algorithms, intermolecular effects, and structures indicate that for the disordered system, geometry optimization gives the best accord with experimental NMR chemical shifts, which has a theory-versus-experiment correlation R(2) = 0.949 and SD = 1.1 ppm, or R(2) = 0.994 and SD = 0.4 ppm when the average of two unassigned NMR shifts for each molecule is used. In addition, these calculations indicate that the experimental chemical shift nonequivalence in this system is mainly a consequence of the different geometries around the phosphorus atoms due to disordered environments. The experimental (31)P NMR chemical shifts are well correlated (R(2) = 0.981) with two conformation angles and one bond length, each associated with one of the three bonding interactions around the phosphorus atoms. These results will facilitate the use of quantum chemical techniques in structural characterization of disordered solids and elucidation of NMR properties.

Chemical shift computations on a crystallographic basis: some reflections and comments

Computations for chemical shifts of molecular organic compounds using the gauge-including projector augmented wave method and the NMR-CASTEP code are reviewed. The methods are briefly introduced, and some general aspects involving the sources of uncertainty in the results are explored. The limitations are outlined. Successful applications of the computations to problems of interpretation of NMR results are discussed and the range of areas in which useful information is obtained is illustrated by examples. The particular value of the computations for comparing shifts between resonances where the same chemical site is involved is emphasised. Such cases arise for shifts between different crystallographically independent molecules of the same chemical species, between polymorphs and for shift anisotropies and asymmetries.

Structural characterization of an anhydrous polymorph of paclitaxel by solid-state NMR

The three-dimensional structure of a unique polymorph of the anticancer drug paclitaxel (Taxol) is established using solid state NMR (SSNMR) tensor ((13)C & (15)N) and heteronuclear correlation ((1)H-(13)C) data. The polymorph has two molecules per asymmetric unit (Z' = 2) and is thus the first conformational characterization with Z' > 1 established solely by SSNMR. Experimental data are correlated with structure through a series of computational models that extensively sample all conformations. For each computational model, corresponding tensor values are computed to supply comparisons with experimental information which, in turn, establishes paclitaxel's structure. Heteronuclear correlation data at thirteen key positions provide shift assignments to the asymmetric unit for each comparison. The two distinct molecules of the asymmetric unit possess nearly identical baccatin III moieties with matching conformations of the C10 acetyl moiety and, specifically, the torsion angle formed by C30-O-C10-C9. Additionally, both are found to exhibit an extended conformation of the phenylisoserine sidechain at C13 with notable differences in the dihedral angles centered around the rotation axes of O-C13, C2'-C1' and C3'-C2'.

A Solid State 13C NMR, Crystallographic, and Quantum Chemical Investigation of Phenylalanine and Tyrosine Residues in Dipeptides and Proteins

We report the results of a solid-state NMR and quantum chemical investigation of the 13C gamma NMR chemical shifts in phenylalanine and tyrosine in dipeptides and proteins. Accurate computation of the experimental shifts is shown to require a good description of local electrostatic field effects, and we find the best results (R2=0.94, rmsd=1.6 ppm, range = 17.1 ppm, N=14) by using a self-consistent reaction field continuum model. There are no obvious correlations with phi, psi, chi 1, or chi2 torsion angles, unlike the results seen with other amino acids. There is, however, a linear relation between computed C gamma atomic charges and shifts for the 14 peptide as well as 18 protein residues investigated. This result is similar to the correlation reported in the 1960s between pi-electron density and 13C shifts for classical 4n + 2 (n=0, 1, 2) pi-electron aromatic species, such as cyclopentadienide and the tropylium cation, and in fact, we found that the shielding/atomic charge correlation seen in the peptides and proteins is virtually identical to that seen with a broad range of aromatic carbocations/carbanions. These results suggest the dominance of an electrostatic field polarization model in which increasing pi electron density results in an increase in C gamma atomic charge and increased shielding (of sigma 11 and sigma 22, perpendicular to the pi orbital) in Phe and Tyr, as well as in the other aromatic species. These results are of general interest since they demonstrate the importance of electrostatic field effects on Phe and Tyr C gamma chemical shifts in peptides and proteins and imply that inclusion of these effects will be necessary in order to interpret the shifts of other aromatic species, such as drug molecules, bound to proteins.

Revisiting the calculation of (13)C chemical shift tensors in cadmium acetate dihydrate with EIM and EIM/cluster methods

The chemical shift tensors of the acetate anions in cadmium acetate dihydrate are calculated using a cluster approach, the embedded ion method (EIM), and a combination of the two in the EIM/cluster method. The results of these calculations are compared with those completed on the isolated acetate anion and show the need for the inclusion of intermolecular interactions. The RMS difference between experiment and theory improves from over 60 ppm when the calculation is completed on an isolated anion, to below 10 ppm when interactions to nearby atoms are included. The best cluster model includes three cadmium acetate dihydrate and gives an RMS result of 4.4 ppm. The EIM method, which uses point charges to account for the intermolecular effects, achieves an RMS of 7.7 ppm on individual anions alone. A combination of the two, the EIM/cluster method, shows that the only necessary atom to explicitly add is the nearest cadmium; this addition results in an RMS of 4.1 ppm. These results are also discussed in terms of the computational cost of the different calculations.

31P NMR Chemical Shifts in Hypervalent Oxyphosphoranes and Polymeric Orthophosphates

We report the first quantum chemical investigation of the solid- and solution-state 31P NMR chemical shifts in models for phosphoryl transfer enzyme reaction intermediates and in polymeric inorganic phosphates. The 31P NMR chemical shifts of five- and six-coordinate oxyphosphoranes containing a variety of substitutions at phosphorus, as well as four-coordinate polymeric orthophosphates and four-coordinate phosphonates, are predicted with a slope of 1.00 and an R2= 0.993 (N = 34), corresponding to a 3.8 ppm (or 2.1%) error over the entire 178.3 ppm experimental chemical shift range, using Hartree-Fock methods. For the oxyphosphoranes, we used either experimental crystallographic structures or, when these were not available, fully geometry optimized molecular structures. For the four-coordinate phosphonates we used X-ray structures together with charge field perturbation, to represent lattice interactions. For the three-dimensional orthophosphates (BPO4, AlPO4, GaPO4 we again used X-ray structures, but for these inorganic systems we employed a self-consistent charge field perturbation approach on large clusters, to deduce peripheral atom charges. For pentaoxyphosphoranes, the solvent effect on 31P NMR chemical shieldings was found to be very small (<0.5 ppm). The 31P NMR chemical shielding tensors in the pentaoxyphosphoranes were in most cases found to be close to axially symmetric and were dominated by changes in the shielding tensor components in the equatorial plane (sigma22 and sigma33). The isotropic shifts were highly correlated (R2= 0.923) with phosphorus natural bonding orbital charges, with the larger charges being associated with shorter axial P-O bond lengths and, hence, more shielding. Overall, these results should facilitate the use of 31P NMR techniques in investigating the structures of more complex systems, such as phosphoryl transfer enzymes, as well as in investigating other, complex oxide structures.

Quantum chemical studies of protein structure

Quantum chemical methods now permit the prediction of many spectroscopic observables in proteins and related model systems, in addition to electrostatic properties, which are found to be in excellent accord with those determined from experiment. I discuss the developments over the past decade in these areas, including predictions of nuclear magnetic resonance chemical shifts, chemical shielding tensors, scalar couplings and hyperfine (contact) shifts, the isomer shifts and quadrupole splittings in Mössbauer spectroscopy, molecular energies and conformations, as well as a range of electrostatic properties, such as charge densities, the curvatures, Laplacians and Hessians of the charge density, electrostatic potentials, electric field gradients and electrostatic field effects. The availability of structure/spectroscopic correlations from quantum chemistry provides a basis for using numerous spectroscopic observables in determining aspects of protein structure, in determining electrostatic properties which are not readily accessible from experiment, as well as giving additional confidence in the use of these techniques to investigate questions about chemical bonding and chemical reactions.

Quantum capping potentials with point charges: A simple QM/MM approach for the calculation of large-molecule NMR shielding tensors

A simple quantum-mechanics/molecular-mechanics (QM/MM) approach for calculating NMR shielding tensors (sigma) is presented. The method involves capping the QM region with quantum capping potentials (QCPs) and representing the MM region with point charges. Test calculations on simple systems without MM charges show that calculated sigma values improve relative to the full QM results with increasing distance between the capped bond and chromophore. Calculations on the histidine amino acid and cytosine monophosphate (CMP) nucleic acid show that the use of QCPs with point charges result in mean errors in the isotropic component of sigma that are less than 1.6 ppm. The results also reveal that, contrary to previous work, the explicit effect of point charges on sigma through coupling with gauge factors, as in the gauge including atomic orbital approach, is minimal for the CMP molecule. The present QM/MM approach for calculating sigma is easy to apply and requires no code modification.

Tryptophan Chemical Shift in Peptides and Proteins: A Solid State Carbon-13 Nuclear Magnetic Resonance Spectroscopic and Quantum Chemical Investigation

We have obtained the carbon-13 nuclear magnetic resonance spectra of a series of tryptophan-containing peptides and model systems, together with their X-ray crystallographic structures, and used quantum chemical methods to predict the (13)C NMR shifts (or shieldings) of all nonprotonated aromatic carbons (C(gamma), C(delta 2) and C(epsilon 2). Overall, there is generally good accord between theory and experiment. The chemical shifts of Trp C(gamma) in several proteins, hen egg white lysozyme, horse myoglobin, horse heart cytochrome c, and four carbonmonoxyhemoglobins, are also well predicted. The overall Trp C(gamma) shift range seen in the peptides and proteins is 11.4 ppm, and individual shifts (or shieldings) are predicted with an rms error of approximately 1.4 ppm (R value = 0.86). Unlike C(alpha) and N(H) chemical shifts, which are primarily a function of the backbone phi,psi torsion angles, the Trp C(gamma) shifts are shown to be correlated with the side-chain torsion angles chi(1) and chi(2) and appear to arise, at least in part, from gamma-gauche interactions with the backbone C' and N(H) atoms. This work helps solve the problem of the chemical shift nonequivalences of nonprotonated aromatic carbons in proteins first identified over 30 years ago and opens up the possibility of using aromatic carbon chemical shift information in structure determination.

An experimental and theoretical study of the 13C and 31P chemical shielding tensors in solid O-phosphorylated amino acids

A series of l and dl forms of O-phosphorylated amino acids (serine, threonine, tyrosine) have been studied by using solid-state multinuclear NMR spectroscopy and ab initio calculations. Principal elements of the (13)C and (31)P chemical shielding tensors have been measured and discussed in relation to zwitterionic structures and intermolecular contacts. DFT calculations have been compared with experimental data showing their ability to reproduce experimentally obtained tensor values in this challenging class of compounds. The changes of orientation of (31)P chemical shielding tensor with respect to the molecular frame in the presence of hydrogen bonds have been revealed and discussed on the ground of theoretical calculations. The measurements of internuclear P...P distances, based on Zeeman magnetization exchange between (31)P spins with differing chemical shielding tensor orientations, were exploited for a clear distinction between enantiomers and racemates.

13C NMR investigation of solid-state polymorphism in 10-deacetyl baccatin III

To investigate the origins of solid-state NMR shift differences in polymorphs, carbon NMR chemical shift tensors are measured for two forms of solid 10-deacetyl baccatin III: a dimethyl sulfoxide (DMSO) solvate and an unsolvated form. A comparison of ab initio computed tensors that includes and omits the DMSO molecules demonstrates that lattice interactions cannot fully account for the shift differences in the two forms. Instead, conformational differences in the cyclohexenyl, benzoyl, and acetyl moieties are postulated to create the differences observed. X-ray analysis of six baccatin III analogues supports the suggested changes in the cyclohexenyl and benzoyl systems. The close statistical match of the (13)C chemical shifts of both polymorphic forms with those calculated using the X-ray geometry of 10-deacetyl baccatin III supports the contention that the B, C, and D rings are fairly rigid. Therefore, the observed tensor differences appear to arise primarily from conformational variations in ring substituents and the cyclohexenyl ring.

Chemical shifts in amino acids, peptides, and proteins: from quantum chemistry to drug design

This chapter discusses recent progress in the investigation and use of (13)C, (15)N, and (19)F nuclear magnetic resonance (NMR) chemical shifts and chemical shift tensors in proteins and model systems primarily using quantum chemical (ab initio Hartree-Fock and density functional theory) techniques. Correlations between spectra and structure are made and the techniques applied to other spectroscopic and electrostatic properties as well, including hydrogen bonding, ligand binding to heme proteins, J-couplings, electric field gradients, and atoms-in-molecules theory, together with a brief review of the use of NMR chemical shifts in drug design.

Polarization of cinnamoyl-CoA substrates bound to enoyl-CoA hydratase: correlation of (13)C NMR with quantum mechanical calculations and calculation of electronic strain energy

When alpha,beta-unsaturated substrates bind to the active site of enoyl-CoA hydratase, large spectral changes can be observed [D'Ordine, R. L., et al. (1994) Biochemistry 33, 12635-12643]. The differences in the isotropic magnetic shieldings of the free and active site-bound forms of the carbonyl, alpha-, and beta-carbons of the substrates, hexadienoyl-CoA, cinnamoyl-CoA, and (N,N-dimethyl-p-amino)cinnamoyl-CoA have been experimentally determined. The carbonyl and beta-carbons are all deshielded, while the alpha-carbons show increased shielding. These chemical shift perturbations are interpreted to suggest that the pi-electrons of the enoyl thiolester are polarized when bound at the active site. Using the crystal structure of (N,N-dimethyl-p-amino)cinnamoyl-CoA bound at the enzyme active site, the shielding tensors were calculated at three different levels of theory, up to a density functional theory model that included all of the contiguous active site residues. These calculations successfully reproduced the observed spectral changes and permitted the electronic polarization of the substrate to be quantified as an electron density difference map. The calculated electron density difference confirms the loss of electrons at the electrophilic beta-carbon and carbonyl carbon, while a slight increase in electron density at the alpha-carbon where proton donation occurs during the hydration reaction and a larger increase in electron density at the carbonyl oxygen are predicted. The energy required to polarize the electrons to the observed extent was calculated to be 3.2 kcal/mol. The force that provides the requisite energy for the polarization is the interaction of the electric field generated by the protein at the enzyme active site with the polarizable electrons of the substrate. Because the induced electronic polarization is along the predicted reaction pathway, the extent of substrate activation by the induced electronic strain is catalytically relevant.

Characterization of stereochemistry and molecular conformation using solid-state NMR tensors

A solid-state NMR technique is described for establishing stereochemistry using the natural product terrein as a model compound. This method involves comparison of experimental (13)C tensor principal values with ab initio computed values for all possible computer-generated stereoisomers. In terrein the relative stereochemistry is confirmed by NMR to be 2R*,3S with high statistical probability (>99.5%). The proposed approach also simultaneously verifies the molecular conformation of the two hydroxy groups in terrein established by X-ray data. It is sufficient to use only shift tensor values at carbons 2 and 3, the stereocenters, to characterize both the stereochemistry and molecular conformations. The solid-state NMR method appears to be especially useful for determining relative stereochemistry of compounds or their derivatives that are difficult to crystallize.