Ab-inito Quantum Mechanical Calculations of NMR Chemical Shifts in Nucleic Acids Constituents II Conformational Dependence of the lH and 13C Chemical Shifts in the Ribose

Abstract

The magnetic shielding constant of the different ¹³C and ¹³H nuclei of a deoxyribose are calculated for the C2′ endo and C3′ endo puckerings of the furanose ring as a function of the conformation about the C4′C5′ bond. For the carbons the calculated variations are of several ppm, the C3′ endo puckering corresponding in most cases to a larger shielding than the C2′ endo one. For the protons the calculated variations of chemical shifts are all smaller than 1.3 ppm, that is of the order of magnitude of the variation of the geometrical shielding produced on these protons by the other units of a DNA double helix, with a change of the overall structure of the helix. The computations carried out on the deoxyribose ?3′ and 5′ phosphates for several conformations of the phosphate group tend to show that the changes of conformation of the charged group of atoms produce chemical shift variations smaller than the two conformational parameters of the deoxyribose itself. The calculations carried out for a ribose do give the general features of the differences between the carbon and proton spectra of deoxynucleosides and nucleosides.

The comparison of the measured and calculated phosphorylation shifts tend to show that the counterion contributes significantly, for some nuclei of the deoxyribose, to the shifts measured. The calculated magnitude of this polarization effect on carbon shifts suggests a tentative qualitative interpretation of carbon spectra of the ribose part of DNA double helices.     

Regression formulae for ab initio and density functional calculated chemical shifts

Linear regression formulae are given for converting ¹H and ¹³C magnetic shielding constants calculated at common ab initio and density functional theory levels of calculation into chemical shifts relative to tetramethylsilane. Accuracies of roughly ±2.2 ppm (¹³C) and ±0.15 ppm (¹H) or better are found for the training set for most levels. The highest level calculations do not always give better results than economical standard calculations.Electronic Supplementary Material Supplementary material is available for this article at     

First-principles electronic structure study of rhizoferrin and its Fe(III) complexes

We have determined the structure and coordination chemistry of rhizoferrin (Rf), which is a particular type of siderophore, and its Fe(III) complexes using density functional theory calculations. Our results show that the Fe(III) ion binds in an octahedral coordination, with a low-spin (S = 1/2) charge-neutral chiral complex having the largest binding energy of the investigated complexes. We have also calculated nuclear magnetic resonance parameters, such as chemical shifts for ¹H and ¹³C, and indirect nuclear spin–spin couplings for ¹H–¹H and ¹³C–¹H in free Rf and in a low-spin neutral Rf metal complex, as well as nuclear quadrupole interaction parameters, such as asymmetry parameter and nuclear quadrupole coupling constants for ¹⁴N. Our calculated values for the chemical shifts for free Rf are in excellent agreement with experimental data while the calculated NMR parameters for Fe(III) complexes are predictions for future experimental work.     

Spectroscopic aspects of the coordination modes of 2,4-dithiohydantoins: Experimental and theoretical study on copper and nickel complexes of cyclohexanespiro-5-(2,4-dithiohydantoin)

The complexation of cyclohexanespiro-5-(2,4-dithiohydantoin), L, with copper and nickel was studied by means of experimental and theoretical methods. The Cu(I) and Ni(II) complexes were synthesized and characterized using ¹³C CPMAS NMR, IR and FAB-MS. Reduction of Cu(II) ions and the formation of Cu(I) complexes with dithiohydantoin was proved. Various coordination modes were investigated on the basis of calculated (DFT-GIAO) shielding constants of the free ligand and model structures of the complexes. General trends in the changes of spectroscopic parameters (NMR chemical shifts, vibrational modes) upon different types of coordination were outlined. Dimeric structures for the Cu(I) and Ni(II) complexes were proposed in which the ligands were coordinated in N3^S4- and N3^S2-bridging ways, respectively, acting as monoanions. The results demonstrate that the combined experimental (¹³C CPMAS NMR, IR) and theoretical (DFT) approach can be used to characterize the molecular structure of solid complexes for which crystallographic data are not available.     

C, N CP MAS and high resolution multinuclear NMR study of methyl

The ¹³C NMR spectra of some isoquinoline and tetrahydroisoquinoline alkaloids and their corresponding N-methosalts and of the bisbenzylisoquinoline alkaloid isochondodendrine were recorded and the signals assigned. The substituent shielding effects and the ¹³C¹H long range couplings were analysed and utilized in the spectral interpretation.     

13C NMR spectral analysis of some isoquinoline alkaloids

The ¹³C NMR spectra of some isoquinoline and tetrahydroisoquinoline alkaloids and their corresponding N-methosalts and of the bisbenzylisoquinoline alkaloid isochondodendrine were recorded and the signals assigned. The substituent shielding effects and the ¹³C¹H long range couplings were analysed and utilized in the spectral interpretation.     

1H, 13C and 199Hg NMR characteristics of new amine-mercuric chloride complexes

Reduction of some N-alkylimines has been achieved with NaBH₄ to give the corresponding secondary amines with high yields (85–95%). These amines were characterized on the bases of their ¹H and ¹³C NMR spectra. The reaction of these amines with mercuric chloride to afford the corresponding complexes was found to occur through a weak dative bond between the nitrogen lone pair of electrons and the mercury atom to form HgCl₂L₂ complexes. The ¹H, ¹³C and ¹⁹⁹Hg NMR chemical shifts have been obtained as well as ¹J(¹³CH) and ²J(¹³CH) coupling constants. Labelling with nitrogen-15 revealed that there is a weak coupling between the nitrogen and the ¹⁹⁹Hg.     

Backbone Structure Confirmation and Side Chain Conformation Refinement of a Bradykinin Mimic BKM-824 by Comparing Calculated 1H, 13C and 19F Chemical Shifts with Experiment

Abstract

Calculated and experimental ¹H, ¹³C and ¹⁹F chemical shifts were compared in BKM-824, a cyclic bradykinin antagonist mimic, c[Ava¹-Igl²-Ser³-DF5F⁴-Oic⁵-Arg⁶] (Ava=5-amino- valeric acid, Igl=α-(2-indanyl)glycine, DF5F=pentafluorophenylalanine, Oic=(2S,3aS,7aS)- octahydroindole-2-carboxylic acid). The conformation of BKM-824 has been studied earlier by NMR spectroscopy (M. Miskolzie et al., J. Biomolec. Struct. Dyn. 17, 947–955 (2000)). All NMR structures have qualitatively the same backbone structure but there is considerable variation in the side chain conformations. We have carried out quantum mechanical optimization for three representative NMR structures at the B3LYP/6–31G* level, constraining the backbone dihedral angles at their NMR structure values, followed by NMR chemical shift calculations at the optimized structures with the 6–311G** basis set. There is an intramolecular hydrogen bond at Ser³ in the optimized structures.

The experimental ¹³C chemical shifts at five C^α positions as well as at the C^β, C^γ and Cδ position of Ava¹, which forms part of the backbone, are well reproduced by the calculations, confirming the NMR backbone structure. A comparison between the calculated and experimental H^β chemical shifts in Igl² shows that the dominant conformation at this residue is gauche. Changes of proton chemical shifts with the scan of the χ1 angle in DF5F⁴ suggest that χ1 ≈180°. The calculated ¹H and ¹³C chemical shifts are in good agreement with experiment at the rigid residue Oic⁵. None of the models gives accurate results for Arg⁶, presumably because of its positive charge. Our study indicates that calculated NMR shifts can be used as additional constraints in conjunction with NMR data to determine protein conformations. However, to be computationally effective, a database of chemical shifts in small peptide fragments should be precalculated.     

Synthesis and conformational analysis of carbasugar bioisosteres of α-l-iduronic acid and its methyl glycoside

The synthesis of two novel carbasugar analogues of α-l-iduronic acid is described in which the ring-oxygen is replaced by a methylene group. In analogy with the conformational equilibrium described for α-l-IdopA, the conformation of the carbasugars was investigated by ¹H and ¹³C NMR spectroscopy. Hadamard transform NMR experiments were utilised for rapid acquisition of ¹H,¹³C-HSQC spectra and efficient measurements of heteronuclear long-range coupling constants. Analysis of ¹H NMR chemical shifts and J_H,H coupling constants extracted by a total-lineshape fitting procedure in conjunction with J_H,C coupling constants obtained by three different 2D NMR experiments, viz., ¹H,¹³C-HSQC-HECADE, J-HMBC and IPAP-HSQC-TOCSY-HT, as well as effective proton-proton distances from 1D ¹H,¹H T-ROE and NOE experiments showed that the conformational equilibrium ⁴C₁?²S_5a?¹C₄ is shifted towards ⁴C₁ as the predominant or exclusive conformation. These carbasugar bioisosteres of α-l-iduronic acid do not as monomers show the inherent flexibility that is anticipated to be necessary for biological activity.     

Conformational analysis of the 2'-deoxyribofuranose ring from proton-proton coupling constants: analysis of a nucleoside-carcinogen adduct formed from 2-acetylaminofluorene utilizing a three-state model

The conformation of the sugar moiety of 8-(N-fluoren-2-ylamino)-2′-deoxyguanosine in solution has been examined as a function of temperature by ¹H-nmr spectroscopy. Analysis of coupling constants shows that lowering the temperature to ?50°C in methanol shifts the conformational equilibrium of the sugar ring resulting in a C2′-endo conformation at a mole fraction of 0.97. The computed phase angle of pseudorotation and amplitude of pucker are 154° and 36°, respectively, with very little discrepancy between the five calculated coupling constants and coupling constants extrapolated from the temperature profiles. A computer program has been written enabling a three-state best-fit analysis. The three-state analysis indicates an equilibrium between C2′-endo, C3′-endo, and 04′-endo conformations. In aqueous solution, the computed mole fraction of the 04′-endo form is 0.18 at 30°C. The conformation associated with the sugar ring and the C4′? C5′ bond is compared to that of 2′-deoxyguanosine.     

Spectroscopic investigation on glutamic acid by Coulomb-attenuating and double hybrid density functional theory methods

This study deals with the identification of glutamic acid by means of quantum chemical approach. FT-IR, FT-Raman and UV–vis spectra were recorded in the region 4000–400, 4000–50 cm^? 1 and 200–600 nm, respectively. CAM-B3LYP/6-31G(d,p) and B2PLYP/6-31G(d,p) calculations were performed to obtain the optimised molecular structures, vibrational frequencies and corresponding vibrational assignment, thermodynamic properties and natural bonding orbital (NBO) analysis. The results show that the obtained optimised geometric parameters (bond lengths, bond angles and bond dihedrals) and vibrational frequencies were found to be in good agreement with the experimental results. The calculations of the electronic spectra were compared with the experimental ones. Furthermore, highest occupied molecular orbital and lowest unoccupied molecular orbital analyses and UV–vis spectral analysis were also performed to determine the energy band gaps and transition states. NBO analysis, calculated using density functional theory methods (CAM-B3LYP/6-31G(d,p) and B2PLYP/6-31G(d,p)), was induced to find inter-molecular atoms. ¹³C and ¹H NMR isotropic chemical shifts were calculated and the assignments made were compared with the ChemDraw Ultra values.     

Molecular structure and spectroscopic (FT-IR,FT-Raman, 13C, 1H NMR and UV) studies of 3,4-dihydroxy-l-phenylalanine using density functional theory

Fourier transform infrared (FT-IR) and Fourier transform Raman (FT-Raman) spectra of 3,4-dihydroxy-l-phenylalanine (3,4-DPA) in solid phase were recorded and analysed in this research. Along with this, the IR spectra in CHCl₃ and the use of acetone as solvents of 3,4-DPA were also recorded. The equilibrium geometry, bonding features and harmonic vibrational frequencies were investigated with the help of density functional theory (DFT) method. The ¹H and ¹³C nuclear magnetic resonance chemical shifts of the molecule were calculated using the gauge including atomic orbital method and compared with experimental results. Stability of the molecule arising from hyperconjugative interactions and charge delocalisation was analysed using natural bond orbital analysis. The results show that charge in electron density (E _D) in the σ* and π* antibonding orbitals and second-order delocalisation energies E(2) confirms the occurrence of intramolecular charge transfer within the molecule. UV–vis spectrum of the compound was recorded and the electronic properties, such as HOMO and LUMO energies, were analysed using the time-dependent (TD)-DFT approach. Finally, the calculation results were applied to simulate infrared and Raman spectra of the title compound, which showed good agreement with the observed spectra.     

DFT study of (17)O, (1)H and (13)C NMR chemical shifts in two forms of native cellulose, I(α) and I(β)

This computational study is intended to shed light on the crystalline and molecular structure, together with the hydrogen bonding (H-bonding) differences between two forms of native cellulose. DFT calculations were carried out to characterize the ¹⁷O, ¹H and ¹³C nuclear magnetic resonance (NMR) parameters in cellulose I_α and I_β with the B3LYP functional employing the 6–311++G7 and 6–31+G1 basis sets. Geometry optimization revealed that the average HB length is shortened by 0.01–0.08 Å when the chains are aligned, whereas the average bond angle increases by about 4–8° exhibiting the enhancement of HB strength. For the isolated cellotetramer chains, the isotropic ¹⁷O–H chemical shifts were plotted as a function of HB length. Our results indicated that as the HB length in cellotetramer I_α increases, the ¹⁷O–H chemical shift isotropy increases, but this parameter changes in the opposite direction for the other structure. Moreover, B3LYP/6–311++G7 calculations reveal that there is an acceptable correlation between the calculated ¹³C chemical shifts of the two structures and their experimental values.     

13C nuclear magnetic resonance study of the complexation of calcium by taurine

¹³C Nuclear magnetic resonance chemical shifts, ¹J_C-C scalar coupling constants, spin-lattice relaxation times, and nuclear Overhauser effects were determined for taurine-[1, 2 ¹³C] and a taurine-[1 ¹³C] and taurine-[2 ¹³C] mixture in the presence and absence of calcium. Ionization constants for taurine amino and sulfonic acid groups and chemical shifts of N-methylene and S-methylene carbons of the taurine cation, zwitterion, and anion were obtained from simultaneous least squares analysis of ¹³C titration curves of both taurine carbons. Comparison of taurine titration shifts to values for related compounds reveals some unusual electronic properties of the taurine molecule. Stability constants of 1:1 calcium complexes with taurine zwitterions and anions, as well as their ¹³C chemical shifts, were obtained by least squares analysis of titration curves measured in the presence of calcium. The stability constants of calcium-taurine complexes were significantly lower than previous values and led to estimates that only approximately one percent of intracellular calcium of mammalian myocardial cells would exist in a taurine complex. The implications of these results with respect to the effect of taurine on calcium ion flux are discussed.     

1H and 13C NMR spectra,protonation, deprotonation,and host-guest complexation-induced shifts of some fluorescence dyes

¹H and ¹³C NMR signal assignments for 8-anilinonaphthalenesulfonic acid (ANS) and dansyl amide (DNSA) are achieved using high-field spectra, decoupling, and two-dimensional NMR techniques as well as shift differences between conjugate acid and bases. Complexation of ANS and DNSA with a macrocyclic azoniacyclophane is measured by fluorescence and by NMR shift titration, furnishing an independent check for the equilibrium constant determination. The complexation-induced shifts (CIS) for ANS and DNSA are analyzed on the basis of aromatic ring current and linear electric field effect models. Comparison of equilibrium constants of the cyclophane and different substrates shows that, e.g., for ANS, lipophilic/hydrophobic binding dominates over electrostatic effects despite the presence of charges and the absence of a lipophilic cap or bottom on the receptor molecule.     

Rapid and accurate calculation of protein 1H, 13C and 15N chemical shifts

A computer program (SHIFTX) is described which rapidly and accurately calculates the diamagnetic ¹H, ¹³C and ¹⁵N chemical shifts of both backbone and sidechain atoms in proteins. The program uses a hybrid predictive approach that employs pre-calculated, empirically derived chemical shift hypersurfaces in combination with classical or semi-classical equations (for ring current, electric field, hydrogen bond and solvent effects) to calculate ¹H, ¹³C and ¹⁵N chemical shifts from atomic coordinates. The chemical shift hypersurfaces capture dihedral angle, sidechain orientation, secondary structure and nearest neighbor effects that cannot easily be translated to analytical formulae or predicted via classical means. The chemical shift hypersurfaces were generated using a database of IUPAC-referenced protein chemical shifts – RefDB (Zhang et al., 2003), and a corresponding set of high resolution (<2.1 Å) X-ray structures. Data mining techniques were used to extract the largest pairwise contributors (from a list of 20 derived geometric, sequential and structural parameters) to generate the necessary hypersurfaces. SHIFTX is rapid (< 1 CPU second for a complete shift calculation of 100 residues) and accurate. Overall, the program was able to attain a correlation coefficient (r) between observed and calculated shifts of 0.911 (¹H), 0.980 (¹³C), 0.996 (¹³C), 0.863 (¹³CO), 0.909 (¹⁵N), 0.741 (¹HN), and 0.907 (sidechain ¹H) with RMS errors of 0.23, 0.98, 1.10, 1.16, 2.43, 0.49, and 0.30 ppm, respectively on test data sets. We further show that the agreement between observed and SHIFTX calculated chemical shifts can be an extremely sensitive measure of the quality of protein structures. Our results suggest that if NMR-derived structures could be refined using heteronuclear chemical shifts calculated by SHIFTX, their precision could approach that of the highest resolution X-ray structures. SHIFTX is freely available as a web server at http://redpoll.pharmacy.ualberta.ca.     

Automated and assisted RNA resonance assignment using NMR chemical shift statistics

The three-dimensional structure determination of RNAs by NMR spectroscopy relies on chemical shift assignment, which still constitutes a bottleneck. In order to develop more efficient assignment strategies, we analysed relationships between sequence and ¹H and ¹³C chemical shifts. Statistics of resonances from regularly Watson–Crick base-paired RNA revealed highly characteristic chemical shift clusters. We developed two approaches using these statistics for chemical shift assignment of double-stranded RNA (dsRNA): a manual approach that yields starting points for resonance assignment and simplifies decision trees and an automated approach based on the recently introduced automated resonance assignment algorithm FLYA. Both strategies require only unlabeled RNAs and three 2D spectra for assigning the H2/C2, H5/C5, H6/C6, H8/C8 and H1′/C1′ chemical shifts. The manual approach proved to be efficient and robust when applied to the experimental data of RNAs with a size between 20 nt and 42 nt. The more advanced automated assignment approach was successfully applied to four stem-loop RNAs and a 42 nt siRNA, assigning 92–100% of the resonances from dsRNA regions correctly. This is the first automated approach for chemical shift assignment of non-exchangeable protons of RNA and their corresponding ¹³C resonances, which provides an important step toward automated structure determination of RNAs.     

Relationship between1H and13C NMR chemical shifts and the secondary and tertiary structure of a zinc finger peptide

Summary Essentially complete assignments have been obtained for the¹H and protonated¹³C NMR spectra of the zinc finger peptide Xfin-31 in the presence and absence of zinc. The patterns observed for the¹H and¹³C chemical shifts of the peptide in the presence of zinc, relative to the shifts in the absence of zinc, reflect the zinc-mediated folding of the unstructured peptide into a compact globular structure with distinct elements of secondary structure. Chemical shifts calculated from the 3D solution structure of the peptide in the presence of zinc and the observed shifts agree to within ca. 0.2 and 0.6 ppm for the backbone C^aH and NH protons, respectively. In addition, homologous zinc finger proteins exhibit similar correlations between their¹H chemical shifts and secondary structure.     

Spectroscopic (FTIR,FT-Raman,NMR and UV) and molecular structure investigations of 1,5-diphenylpenta-2,4-dien-1-one: a combined experimental and theoretical approach

The title molecule 1,5-diphenylpenta-2,4-dien-1-one (cinnamylideneacetophenone, CA) has been synthesised and characterised by FTIR, FT-Raman, NMR and UV–vis spectral analyses. The possible stable conformers of the CA molecule were searched by potential energy surface scan at B3LYP level of theory. The molecular geometry from X-ray determination of the CA molecule in the ground state has been compared using the density functional theory (DFT) with 6-31G(d,p) basis set. The harmonic vibrational modes, the corresponding wavenumbers and IR and Raman intensities of most stable conformer were calculated by the DFT method. The assignments of the fundamentals were proposed on the basis of total energy distribution calculations. The calculated ¹³C and ¹H NMR chemical shifts using gauge including atomic orbitals approach are in good agreement with the observed chemical shifts. The molecular stability and bond strength have been investigated by applying natural bond orbital analysis. Using the time-dependent DFT method, the electronic absorption spectrum of the title compound has been predicted and the electronic transitions within the molecule have been interpreted. The molecular electrostatic potential map was used for predicting possible hydrogen and oxygen bonding sites in the CA molecule.     

Complete H and C NMR chemical shift assignments of mono-, di-, and trisaccharides as basis for NMR chemical shift predictions of polysaccharides using the computer program casper

The computer program casper uses ¹H and ¹³C NMR chemical shift data of mono- to trisaccharides for the prediction of chemical shifts of oligo- and polysaccharides. In order to improve the quality of these predictions the ¹H and ¹³C, as well as ³¹P when applicable, NMR chemical shifts of 30 mono-, di-, and trisaccharides were assigned. The reducing sugars gave two distinct sets of NMR resonances due to the α- and β-anomeric forms. In total 35 ¹H and ¹³C NMR chemical shift data sets were obtained from the oligosaccharides. One- and two-dimensional NMR experiments were used for the chemical shift assignments and special techniques were employed in some cases such as 2D ¹H,¹³C-HSQC Hadamard Transform methodology which was acquired approximately 45 times faster than a regular t₁ incremented ¹H,¹³C-HSQC experiment and a 1D ¹H,¹H-CSSF-TOCSY experiment which was able to distinguish spin-systems in which the target protons were only 3.3 Hz apart. The ¹H NMR chemical shifts were subsequently refined using total line-shape analysis with the PERCH NMR software. The acquired NMR data were then utilized in the casper program (http://www.casper.organ.su.se/casper/) for NMR chemical shift predictions of the O-antigen polysaccharides from Klebsiella O5, Shigella flexneri serotype X, and Salmonella arizonae O62. The data were compared to experimental data of the polysaccharides from the two former strains and the lipopolysaccharide of the latter strain showing excellent agreement between predicted and experimental ¹H and ¹³C NMR chemical shifts.