Relationship between NMR shielding and heme binding strength for a series of 7-substituted quinolines

Chemical shielding tensors are calculated for the carbons in a series of 4-aminoquinolines with different substituents at the 7-position. The sigma(11) component is used as a measure of the relative pi-electron density at each carbon. By comparing the pi-electron density at each carbon with the log K of binding to heme (Kaschula et al. J. Med. Chem. 2002, 45, 3531), the drug-heme association is found to increase with increasing pi-electron density at the carbons meta to the substituent and with decreasing pi-electron density at the carbons ortho and para to the substituent. The greatest change in pi-electron density is at the ortho carbons, and log K increases with a decrease in pi-electron density on the ring containing the substituent, which corresponds to an increase in the pi-dipole between the two rings. An examination of the solution structures of the pi-pi complexes formed by amodiaquine and quinine with heme (Leed et al. Biochemistry 2002, 41, 10245. de Dios et al. Inorg. Chem. 2004, 43, 8078) shows that the pi-dipoles in each drug and in the porphyrin ring of heme may be paired. The chloro-substituted compound has an association constant that is an order of magnitude higher than the other compounds in the series, but the pi-electron density at the ring containing the substituent is not correspondingly low. This lack of correlation indicates that the Cl-substituted compound may be binding to heme in a manner that differs from the other compounds in the series.     

Interactions between pairs of antimalarial drugs studied by experimental and ab initio (13)C NMR chemical shifts

The interactions between the antimalarial drugs chloroquine (CQ) and amodiaquine (AQ), chloroquine and quinine (QN), and amodiaquine and quinine are studied by (13)C NMR. Experimental changes in chemical shift are compared to nucleus-independent chemical shifts to determine the best structure of the complex formed by each drug pair in solution. Structures of the CQ-AQ and CQ-QN complexes are found to be similar to those found previously for the drug dimers. On the other hand, the best solution structure for the AQ-QN complex suggests that the quinoline rings of the two drugs are at an angle with respect to each other.     

Calculations of the contribution of ring currents to the chemical shielding anisotropy

Ring currents can exert a large effect upon the chemical shielding of NMR resonances. The analytical expression developed by Waugh and Fessenden (Waugh, J. S.; Fessenden, R. W. J. Am. Chem. Soc. 1957, 79, 846) and Johnson and Bovey (Johnson, C. E.; Bovey, F. A. J. Chem. Phys. 1957, 29, 1012) only quantifies the contribution of ring currents to the isotropic component of the shielding tensor. In the work described here an additional analytical expression is developed so that the contribution of ring currents to the full shielding tensor can be calculated, allowing an estimate of the influence of ring currents upon the chemical shielding anisotropy (CSA, Deltasigma). To test that this pair of analytical expressions can provide a reasonable estimate of the contribution of ring currents to the full shielding tensor a series of density functional calculations (DFT) were carried out. A shielding tensor in a model compound was calculated in two distinct ways. For the first series, DFT shielding calculations of the model compound were carried out in the presence of a benzene ring. For the second series a ring current contribution to the shielding tensor was calculated via the new expressions, and this was added to the result of a DFT shielding calculation which used in place of benzene the nonaromatic analogue 1,3 cyclohexadiene. The two series of results proved to be in excellent agreement. The pair of analytical expressions are used to calculate ring current contributions to the CSA (Deltasigma) of 1H(N) backbone amide resonances in a structure of the second type 2 module from the protein fibronectin. Significant CSA variations are predicted in particular for the 1H(N) of G42 which is most likely involved in a N-H...tpi aromatic hydrogen bond.     

Structural and solvent effects on the 13C and 15N NMR chemical shifts of indoloquinoline alkaloids: experimental and DFT study

Indoloquinoline alkaloids represent an important class of antimalarial, antibacterial and antiviral compounds. They have been shown to bind to DNA via intercalation preferentially at GC-rich sequences containing nonalternating CC sites. The stability of complexes formed with biological macromolecules depends on noncovalent binding. In the present study, the ability of indoloquinolines to form intermolecular interactions with solvents was investigated by using NMR spectroscopy and density functional theory (DFT) (B3LYP/6-31G**) calculations. NMR data measured for indoloquinoline bases and the corresponding hydrochlorides are discussed in relation to the structure. DFT calculations of shielding constants in vacuo and in solution allowed the investigation of the influence of the environment on the NMR parameters. Calculations incorporating solvent effects indicated significant changes in the anisotropy of the electron distribution, reflected in the span of the chemical shielding tensor (Omega = sigma11 - sigma33). Solvent effects on the span of the 13C and 15N shielding tensor depended on the type of atom and the data indicated a significant influence of solute-solvent interactions.     

Proton magnetic shielding tensor in liquid water

The nuclear magnetic shielding tensor is a sensitive probe of the local electronic environment, providing information about molecular structure and intermolecular interactions. The magnetic shielding tensor of the water proton has been determined in hexagonal ice, but in liquid water, where the tensor is isotropically averaged by rapid molecular tumbling, only the trace of the tensor has been measured. We report here the first determination of the proton shielding anisotropy in liquid water, which, when combined with chemical shift data, yields the principal shielding components parallel (sigma(parallel)) and perpendicular (sigma(perpendicular)) to the O-H bond. We obtained the shielding anisotropy sigma(parallel)-sigma(perpendicular) by measuring the proton spin relaxation rate as a function of magnetic induction field in a water sample where dipole-dipole couplings are suppressed by H/D isotope dilution. The temperature dependence of the shielding components, determined from 0 to 80 degrees C, reflects vibrational averaging over a distribution of instantaneous hydrogen-bond geometries in the liquid and thus contains unique information about the temperature-dependent structure of liquid water. The temperature dependence of the shielding anisotropy is found to be 4 times stronger than that of the isotropic shielding. We analyze the liquid water shielding components in the light of previous NMR and theoretical results for vapor and ice. We show that a simple two-state model of water structure fails to give a consistent interpretation of the shielding data and we argue that a more detailed analysis is needed that quantitatively relates the shielding components to hydrogen bond geometry.     

First-principles calculations within periodic boundary conditions of the NMR shielding tensor for a transition metal nucleus in a solid state system: The example of 51V in AlVO4

We present the first density functional theory based calculations of NMR shielding parameters for a transition metal nucleus using periodic boundary conditions. These calculations employ the gauge-including projected augmented-wave pseudopotential approach. The quality of this method is discussed by comparing experimental and calculated chemical shift tensor eigenvalues for the quadrupolar 51V nucleus in the diamagnetic solid-state compound AlVO4. Furthermore, the combination of shielding tensor with fast and accurate projector augmented-wave electric field gradient tensor calculations allows us to determine the relative orientation of these two tensors.     

Carbon-13 NMR shielding in the twenty common amino acids: comparisons with experimental results in proteins

We have used ab initio quantum chemical techniques to compute the (13)C(alpha) and (13)C(beta) shielding surfaces for the 14 amino acids not previously investigated (R. H. Havlin et al., J. Am. Chem. Soc. 1997, 119, 11951-11958) in their most popular conformations. The spans (Omega = sigma(33) - sigma(11)) of all the tensors reported here are large ( approximately 34 ppm) and there are only very minor differences between helical and sheet residues. This is in contrast to the previous report in which Val, Ile and Thr were reported to have large ( approximately 12 ppm) differences in Omega between helical and sheet geometries. Apparently, only the beta-branched (beta-disubstituted) amino acids have such large CSA span (Omega) differences; however, there are uniformly large differences in the solution-NMR-determined CSA (Deltasigma = sigma(orth) - sigma(par)) between helices and sheets in all amino acids considered. This effect is overwhelmingly due to a change in shielding tensor orientation. With the aid of such shielding tensor orientation information, we computed Deltasigma values for all of the amino acids in calmodulin/M13 and ubiquitin. For ubiquitin, we find only a 2.7 ppm rmsd between theory and experiment for Deltasigma over an approximately 45 ppm range, a 0.96 slope, and an R(2) = 0.94 value when using an average solution NMR structure. We also report C(beta) shielding tensor results for these same amino acids, which reflect the small isotropic chemical shift differences seen experimentally, together with similar C(beta) shielding tensor magnitudes and orientations. In addition, we describe the results of calculations of C(alpha), C(beta), C(gamma)1, C(gamma)2, and C(delta) shifts in the two isoleucine residues in bovine pancreatic trypsin inhibitor and the four isoleucines in a cytochrome c and demonstrate that the side chain chemical shifts are strongly influenced by chi(2) torsion angle effects. There is very good agreement between theory and experiment using either X-ray or average solution NMR structures. Overall, these results show that both C(alpha) backbone chemical shift anisotropy results as well as backbone and side chain (13)C isotropic shifts can now be predicted with good accuracy by using quantum chemical methods, which should facilitate solution structure determination/refinement using such shielding tensor surface information.     

Experimental and computational characterization of the 17O quadrupole coupling and magnetic shielding tensors for p-nitrobenzaldehyde and formaldehyde

We have used solid-state 17O NMR experiments to determine the 17O quadrupole coupling (QC) tensor and chemical shift (CS) tensor for the carbonyl oxygen in p-nitro-[1-(17)O]benzaldehyde. Analyses of solid-state 17O NMR spectra obtained at 11.75 and 21.15 T under both magic-angle spinning (MAS) and stationary conditions yield the magnitude and relative orientation of these two tensors: CQ = 10.7 +/- 0.2 MHz, etaQ = 0.45 +/- 0.10, delta11 = 1050 +/- 10, delta22 = 620 +/- 10, delta33 = -35 +/- 10, alpha = 90 +/- 10, beta = 90 +/- 2, gamma = 90 +/- 10 degrees. The principal component of the 17O CS tensor with the most shielding, delta33, is perpendicular to the H-C=O plane, and the tensor component with the least shielding, delta11, lies along the C=O bond. For the 17O QC tensor, the largest (chi(zz)) and smallest (chi(xx)) components are both in the H-C=O plane being perpendicular and parallel to the C=O bond, respectively. This study represents the first time that these two fundamental 17O NMR tensors have been simultaneously determined for the carbonyl oxygen of an aldehyde functional group by solid-state 17O NMR. The reported experimental solid-state 17O NMR results provide the first set of reliable data to allow evaluation of the effect of electron correlation on individual CS tensor components. We found that the electron correlation effect exhibits significant influence on 17O chemical shielding in directions within the H-C=O plane. We have also carefully re-examined the existing experimental data on the 17O spin-rotation tensor for formaldehyde and proposed a new set of best "experimental" 17O chemical shielding tensor components: sigma11 = -1139 +/- 80, sigma22 = -533 +/- 80, sigma33 = 431 +/- 5, and sigma(iso) = -414 +/- 60 ppm. Using this new set of data, we have evaluated the accuracy of quantum chemical calculations of the 17O CS tensors for formaldehyde at the Hartree-Fock (HF), density-functional theory (DFT), M?ller-Plesset second-order perturbation (MP2), and coupled-cluster singles and doubles (CCSD) levels of theory. The conclusion is that, while results from HF and DFT tend to underestimate the electron correlation effect, the MP2 method overestimates its contribution. The CCSD results are in good agreement with the experimental data.     

Intermolecular shielding from molecular magnetic susceptibility. A new view of intermolecular ring current effects

This paper presents calculations of the NICS (nuclear independent chemical shieldings) in a rectangular grid surrounding the molecules of benzene, naphthalene and coronene. Using the relationship between calculated NICS and the induced magnetic field, the calculated NICS are used to predict intermolecular effects due to molecular magnetic susceptibility or ring current effects. As expected from approximate ring current models, these intermolecular shielding effects are concentrated along the direction perpendicular to the molecular plane and they approach asymptotically to a dipolar functional dependence, i.e. (1-3 cos(2)theta)/r(3)). The deviations from the dipolar functional form require that the calculations of these intermolecular effects be done using a suitable interpolation scheme of the NICS calculated on the grid. The analysis of the NICS tensor components shows that these intermolecular shielding effects should be primarily expected on shielding components of the neighboring molecules nuclei, which are perpendicular to the molecular plane of the aromatic compound generating the induced field. The analysis of the calculated NICS along the series benzene, naphthalene and coronene shows that these intermolecular effects increase monotonically with the number of aromatic rings.     

Solid-state nitrogen-15 NMR and quantum chemical study of N,N-dimethylaniline derivatives

In this study the components of the nitrogen chemical shift (CS) tensor are examined for a series of para substituted N,N-dimethylaniline derivatives. This is done through measurement of the 15N NMR spectra of powder samples and through quantum chemical calculations on the isolated molecules. Experiments and calculations show that the isotropic CS, delta(iso), decreases with increasing electron donating ability of the para substituent, in agreement with previous solution studies. More importantly, this study shows that this decrease in the isotropic (solution) CS is due to decreasing values of the CS tensor component delta(11) and component delta(33). The component delta(22) is essentially invariant to the electron donating/withdrawing ability of the para substituent. Through Ramsey's theory of nuclear magnetic shielding, it can be seen that the variation in delta(11) and delta(33), and hence delta(iso), is due to changes in the n-pi* and the sigma-pi* energy gaps in N,N-dimethylaniline. This, in turn, is a result of the change in the energy of the pi* molecular orbital with change in the pi-electron donating ability of the para substituent. The effects of nitrogen inversion on the components of the nitrogen CS tensor components are also discussed. This study also shows the feasibility of performing 15N cross-polarization experiments on nonspinning powder samples at natural isotopic abundance.     

Calculation of binary magnetic properties and potential energy curve in xenon dimer: second virial coefficient of (129)Xe nuclear shielding

Quantum chemical calculations of the nuclear shielding tensor, the nuclear quadrupole coupling tensor, and the spin-rotation tensor are reported for the Xe dimer using ab initio quantum chemical methods. The binary chemical shift delta, the anisotropy of the shielding tensor Delta sigma, the nuclear quadrupole coupling tensor component along the internuclear axis chi( parallel ), and the spin-rotation constant C( perpendicular ) are presented as a function of internuclear distance. The basis set superposition error is approximately corrected for by using the counterpoise correction (CP) method. Electron correlation effects are systematically studied via the Hartree-Fock, complete active space self-consistent field, second-order M?ller-Plesset many-body perturbation, and coupled-cluster singles and doubles (CCSD) theories, the last one without and with noniterative triples, at the nonrelativistic all-electron level. We also report a high-quality theoretical interatomic potential for the Xe dimer, gained using the relativistic effective potential/core polarization potential scheme. These calculations used valence basis set of cc-pVQZ quality supplemented with a set of midbond functions. The second virial coefficient of Xe nuclear shielding, which is probably the experimentally best-characterized intermolecular interaction effect in nuclear magnetic resonance spectroscopy, is computed as a function of temperature, and compared to experiment and earlier theoretical results. The best results for the second virial coefficient, obtained using the CCSD(CP) binary chemical shift curve and either our best theoretical potential or the empirical potentials from the literature, are in good agreement with experiment. Zero-point vibrational corrections of delta, Delta sigma, chi (parallel), and C (perpendicular) in the nu=0, J=0 rovibrational ground state of the xenon dimer are also reported.     

A correlation study of quinoline derivatives and their pharmaceutical behavior by ab initio calculated NQR parameters

In this paper, ab initio calculated NQR parameters for some quinoline-containing derivatives are presented. The calculations are carried out in a search for the relationships between the charge distribution of these compounds and their ability to interact with haematin. On the basis of NQR parameters, pi-electron density on the nitrogen atom of the quinoline ring plays a dominant role in determining the ability of quinolines to interact with haematin. This point was confirmed with investigation of Fe+3 cation-pi quinoline ring interactions in 2- and 4-aminoquinoline. However, our results do not show any preference for those carbon atoms of the quinoline ring which previous reports have noted. In order to calculate the NQR parameters, the electric field gradient (EFG) should be evaluated at the site of a quadrupolar nucleus in each compound. EFGs are calculated by the Gaussian 98 program using the B3LYP/6-31 G* level of theory.     

A comparison of two-component and four-component approaches for calculations of spin-spin coupling constants and NMR shielding constants of transition metal cyanides

Relativistic density functional theory (DFT) calculations of nuclear spin-spin coupling constants and shielding constants have been performed for selected transition metal (11th and 12th group of periodic table) and thallium cyanides. The calculations have been carried out using zeroth-order regular approximation (ZORA) Hamiltonian and four-component Dirac-Kohn-Sham (DKS) theory with different nonrelativistic exchange-correlation functionals. Two recent approaches for representing the magnetic balance (MB) between the large and small components of four-component spinors, namely, mDKS-RMB and sMB, have been employed for shielding tensor calculations and their results have been compared. Relativistic effects have also been analysed in terms of scalar and spin-orbit contributions at the two-component level of theory, including discussion of heavy-atom-on-light-atom effects for (1)J(CN), σ(C), and σ(N). The results for molecules containing metals from 4th row of periodic table show that relativistic effects for them are small (especially for spin-spin coupling constants). The biggest effects are observed for the 6th row where nonrelativistic theory reproduces only about 50%-70% of the two-component ZORA results for (1)J(MeC) and about 75% for heavy metal shielding constants. It is important to employ a full Dirac picture for calculations of heavy metal shielding constants, since ZORA reproduces only 75%-90% of the DKS results. Smaller discrepancies between ZORA-DFT and DKS are observed for nuclear spin-spin coupling constants. No significant differences are observed between the results obtained using mDKS-RMB and sMB approaches for magnetic balance in four-component calculations of the shielding constants.     

13C and (15)N chemical shift tensors in adenosine,guanosine dihydrate, 2'-deoxythymidine,and cytidine

The (13)C and (15)N chemical shift tensor principal values for adenosine, guanosine dihydrate, 2'-deoxythymidine, and cytidine are measured on natural abundance samples. Additionally, the (13)C and (15)N chemical shielding tensor principal values in these four nucleosides are calculated utilizing various theoretical approaches. Embedded ion method (EIM) calculations improve significantly the precision with which the experimental principal values are reproduced over calculations on the corresponding isolated molecules with proton-optimized geometries. The (13)C and (15)N chemical shift tensor orientations are reliably assigned in the molecular frames of the nucleosides based upon chemical shielding tensor calculations employing the EIM. The differences between principal values obtained in EIM calculations and in calculations on isolated molecules with proton positions optimized inside a point charge array are used to estimate the contributions to chemical shielding arising from intermolecular interactions. Moreover, the (13)C and (15)N chemical shift tensor orientations and principal values correlate with the molecular structure and the crystallographic environment for the nucleosides and agree with data obtained previously for related compounds. The effects of variations in certain EIM parameters on the accuracy of the shielding tensor calculations are investigated.     

PHOTOSENSITIZATION BY ANTIMALARIAL DRUGS

Abstract The antimalarial drugs, chloroquine, hydroxychloroquine, quinine, quinacrine, amodiaquine and primaquine and the local anaesthetic, dibucaine, were tested for in vitro photosensitizing capability by irradiation with 365 nm UV light in aqueous solutions. The ability of these compounds to photosensitize the oxidation of 2,5-dimethylfuran, histidine, tryptophan or xanthine, and to initiate the free radical polymerization of acrylamide was examined in the pH range 2-12. Chloroquine and hydroxychloroquine show maximal photooxidative behaviour when in the monocation form at pH 9, in contrast to quinine which is extremely efficient as the dication below pH 4. This pattern appears to relate to the fluorescence yield as a function of pH. Chloroquine in the monocation or neutral form was found to undergo dechlorination upon irradiation, and this correlates directly with its ability to initiate photo-polymerization of acrylamide. Quinine also gives rise to small polymerization rates, attributed to photo-ionization in the quinoline ring, yielding a cation radical. Amodiaquine, primaquine and quinacrine do not have significant photochemical activity in aqueous solution. Dibucaine exhibits a strong photosensitizing capability at low pH, similar to quinine.     

Probing lead(II) bonding environments in 4-substituted pyridine adducts of (2,6-Me2C6H3S)2Pb: an X-ray structural and solid-state 207Pb NMR study

The effect of subtle changes in the sigma-electron donor ability of 4-substituted pyridine ligands on the lead(II) coordination environment of (2,6-Me(2)C(6)H(3)S)(2)Pb (1) adducts has been examined. The reaction of 1 with a series of 4-substituted pyridines in toluene or dichloromethane results in the formation of 1:1 complexes [(2,6-Me(2)C(6)H(3)S)(2)Pb(pyCOH)](2) (3), [(2,6-Me(2)C(6)H(3)S)(2)Pb(pyOMe)](2) (4), and (2,6-Me(2)C(6)H(3)S)(2)Pb(pyNMe(2)) (5) (pyCOH = 4-pyridinecarboxaldehyde; pyOMe = 4-methoxypyridine; pyNMe2 = 4-dimethylaminopyridine), all of which have been structurally characterized by X-ray crystallography. The structures of 3 and 4 are dimeric and have psi-trigonal bipyramidal S(3)N bonding environments, with the 4-substituted pyridine nitrogen and bridging sulfur atoms in axial positions and two thiolate sulfur atoms in equatorial sites. Conversely, compound 5 is monomeric and exhibits a psi-trigonal pyramidal S(2)N bonding environment at lead(II). The observed structures may be rationalized in terms of a simple valence bond model and the sigma-electron donor ability of the 4-pyridine ligands as derived from the analysis of proton affinity values. Solid-state (207)Pb NMR experiments are applied in combination with density functional theory (DFT) calculations to provide further insight into the nature of bonding in 4, 5, and (2,6-Me(2)C(6)H(3)S)(2)Pb(py)(2) (2). The lead chemical shielding (CS) tensor parameters of 2, 4, and 5 reveal some of the largest chemical shielding anisotropies (CSA) observed in lead coordination complexes to date. DFT calculations using the Amsterdam Density Functional (ADF) program, which take into account relativistic effects using the zeroth-order regular approximation (ZORA), yield lead CS tensor components and orientations. Paramagnetic contributions to the lead CS tensor from individual pairs of occupied and virtual molecular orbitals (MOs) are examined to gain insight into the origin of the large CSA. The CS tensor is primarily influenced by mixing of the occupied MOs localized on the sulfur and lead atoms with virtual MOs largely comprised of lead 6p orbitals.     

An ab initio study of amide proton shift tensor dependence on local protein structure.

Ab initio shielding tensor calculations were carried out on residues in human ubiquitin. Reported experimental data on isotropic and anisotropic components of the amide proton chemical shifts were used as benchmarks to test the validity of the chosen basis sets as well as methods in structure optimization and shielding calculations. The best agreement with the experimental values was observed when the 6-311**G and 6-311++G(2d,2p) basis sets were used to optimize the structure and to calculate the shielding tensor, respectively. The same method was employed in subsequent model calculations to characterize the dependence of amide proton shielding to the local structure. Both the isotropic and the anisotropic components of the symmetric tensor were found to depend very strongly on the hydrogen bond length. A weaker dependence can also be observed for the hydrogen bond angle. Antisymmetric tensor elements were found to be relatively small. This study permits separation of various local structure contributions to the amide proton shielding tensor that complements scarce experimental data.     

A stable neutral stannaaromatic compound: synthesis, structure and complexation of a kinetically stabilized 2-stannanaphthalene

A kinetically stabilized 2-stannanaphthalene, the first example of a stable, neutral aromatic compound containing a tin atom, has been synthesized and fully characterized. The results of spectroscopic and crystallographic structural analyses of the compound and theoretical calculations using model compounds strongly suggest that it has a delocalized 10pi-electron ring system as does naphthalene. In addition, it was found to behave as an eta6-arene ligand in the ligand exchange reaction with [Cr(CH3CN)3(CO)3] giving the first stable eta6-2-stannanaphthalene chromium complex.     

Experimental and theoretical investigation of the 13C and 15N chemical shift tensors in melanostatin-exploring the chemical shift tensor as a structural probe

The determination of backbone conformations in powdered peptides using 13C and 15N shift tensor information is explored. The 13C and 15N principal shift values in natural abundance 13C and 15N melanostatin (L-Pro-L-Leu-Gly amide) are measured using the FIREMAT technique. Furthermore, the orientation of the C-N bond in the 13C shift principal axis system for the backbone carbons is obtained from the presence of the 13C-14N dipolar coupling. The Ramachandran angles for the title compound are obtained from solid-state NMR data by comparing the experimentally determined shift tensor information to systematic theoretical shielding calculations on N-formyl-L-amino acid-amide models. The effects of geometry optimization and neglect of intermolecular interactions on the theoretical shielding values in the model compounds are investigated. The sets of NMR derived Ramachandran angles are assembled in a set of test structures that are compared to the available single-crystal X-ray structure. Shift tensor calculations on the test structures and the X-ray structure are used to further assess the importance of intermolecular interactions when the shift tensor is used as a structural probe in powdered peptides.