Temperature Dependence of Proton NMR Chemical Shift As a Criterion To Identify Low-Barrier Hydrogen Bonds

Discovering the Solid-State Secrets of Lorlatinib by NMR Crystallography: To Hydrogen Bond or not to Hydrogen Bond

Lorlatinib is an active pharmaceutical ingredient (API) used in the treatment of lung cancer. Here, an NMR crystallography analysis is presented whereby the single-crystal X-ray diffraction structure (CSD: 2205098) determination is complemented by multinuclear (1H, 13C, 14/15N, 19F) magic-angle spinning (MAS) solid-state NMR and gauge-including projector augmented wave (GIPAW) calculation of NMR chemical shifts. Lorlatinib crystallises in the P21 space group, with two distinct molecules in the asymmetric unit cell, Z' = 2. Three of the four NH2 hydrogen atoms form intermolecular hydrogen bonds, N30-H…N15 between the two distinct molecules and N30-H…O2 between two equivalent molecules. This is reflected in one of the NH21H chemical shifts being significantly lower, 4.0 ppm compared to 7.0 ppm. Two-dimensional 1H-13C, 14N-1H and 1H (double-quantum, DQ)-1H (single-quantum, SQ) MAS NMR spectra are presented. The 1H resonances are assigned and specific HH proximities corresponding to the observed DQ peaks are identified. The resolution enhancement at a 1H Larmor frequency of 1 GHz as compared to 500 or 600 MHz is demonstrated.

Direct Spectroscopic Evidence for Charge-Assisted Hydrogen-Bond Formation between Ionizable Organic Chemicals and Carbonaceous Materials

The direct evidence for the formation of charge-assisted hydrogen bond (CAHB) between the charged groups of ionizable organic chemicals (IOCs) and carbonaceous materials with similar proton affinity remains elusive. We therefore selected three pharmaceutical contaminants (PCs) as representative IOCs to provide the direct evidence of CAHB formation between IOCs and functionalized carbon nanotubes (CNTs) and its intensity/contribution to PC sorption on CNTs by NMR, FTIR, and DFT analyses. Sorption of PCs on functionalized CNTs resulted in the FTIR characteristic peak that appeared at a higher frequency (3787 cm^-1) and the ¹H NMR characteristic peak that emerged at an extremely low-field region (<18.0 ppm), which can be used as the direct spectroscopic evidence for CAHB formation. Both homonuclear CAHB (HM-CAHB, e.g., [O-H···O]^-) and heteronuclear CAHB (HT-CAHB, e.g., [N⁺-H···O^-]/[O-H···N]⁺) exhibited a much higher sorption energy (|E_ads| ≥ 56.24 kJ/mol) than ordinary hydrogen bond (OHB, |E_ads| ≤ 6.136 kJ/mol), leading to a greater sorption contribution (HM-/HT-CAHB ≥ 42.3%, OHB ≤ 36.5%) and irreversibility (hysteresis index: HM-/HT-CAHB ≥ 1.69, OHB ≤ 0.43) of PCs on CNTs. This work presents the direct evidence for CAHB formation between IOCs and CNTs, which is significant for understanding and predicting the environmental fate and risk of IOCs, thus providing new insights for controlling their pollution using specifically designed carbonaceous materials.

Importance of Nuclear Quantum Effects for Molecular Cocrystals with Short Hydrogen Bonds

Many efforts have been recently devoted to the design and investigation of multicomponent pharmaceutical solids, such as salts and cocrystals. The experimental distinction between these solid forms is often challenging. Here, we show that the transformation of a salt into a cocrystal with a short hydrogen bond does not occur as a sharp phase transition but rather a smooth shift of the positional probability of the hydrogen atoms. A combination of solid-state NMR spectroscopy, X-ray diffraction, and diffuse reflectance measurements with density functional theory calculations that include nuclear quantum effects (NQEs) provides evidence of temperature-induced hydrogen atom shift in cocrystals with short hydrogen bonds. We demonstrate that for the predictions of the salt/cocrystal solid forms with short H-bonds, the computations have to include NQEs (particularly hydrogen nuclei delocalization) and temperature effects.

Variable-Temperature NMR Analysis of the Thermodynamics of Polymer Partitioning between Aqueous and Drug-Rich Phases and Its Significance for Amorphous Formulations

We previously reported that the polymers used in amorphous solid dispersion (ASD) formulations, such as polyvinylpyrrolidone (PVP), polyvinylpyrrolidone/vinyl acetate (PVP-VA), and hypromellose (HPMC), distribute into the drug-rich phase of ibuprofen (IBP) formed by liquid-liquid phase separation, resulting in a reduction in the maximum drug supersaturation in the aqueous phase. Herein, the mechanism underlying the partitioning of the polymer into the drug-rich phase was investigated from a thermodynamic perspective. The dissolved IBP concentration in the aqueous phase and the amount of polymer distributed into the IBP-rich phase were quantitatively analyzed in IBP-supersaturated solutions containing different polymers using variable-temperature solution-state nuclear magnetic resonance (NMR) spectroscopy. The polymer weight ratio in the IBP-rich phase increased at higher temperatures, leading to a more notable reduction of IBP amorphous solubility. Among the polymers, the amorphous solubility reduction was the greatest for the PVP-VA solution at lower temperatures, while HPMC reduced the amorphous solubility to the greatest extent at higher temperatures. The change in the order of polymer impact on the amorphous solubility resulted from the differences in the temperature dependency of polymer partitioning. The van't Hoff plot of the polymer partition coefficient revealed that both enthalpy and entropy changes for polymer transfer into the IBP-rich phase from the aqueous phase (ΔH_{aqueous→IBP-rich} and ΔS_{aqueous→IBP-rich}) gave positive values for most of the measured temperature range, indicating that polymer partitioning into the IBP-rich phase was an endothermic but entropically favorable process. The polymer transfer into the IBP-rich phase was more endothermic for HPMC than for PVP and PVP-VA. The solid-state NMR analysis of the IBP/polymer ASD implied that the newly formed IBP/polymer interactions in the IBP-rich phase upon polymer incorporation were weaker for HPMC, providing a rationale for the larger positive transfer enthalpy for HPMC. The change in Gibbs free energy for polymer transfer (ΔG_{aqueous→IBP-rich}) showed negative values across the experimental temperature range, decreasing with an increase in temperature, indicating that the distribution of the polymer into the IBP-rich phase is favored at higher temperatures. Moreover, ΔG_{aqueous→IBP-rich} for HPMC showed the greatest decrease with the temperature, likely reflecting the temperature-induced dehydration of HPMC in the aqueous phase. This study contributes fundamental insights into the phenomenon of polymer partitioning into drug-rich phases, furthering the understanding of achievable supersaturation levels and ultimately providing information on polymer selection for ASD formulations.

Covalence and π-electron delocalization influence on hydrogen bonds in proton transfer process of o-hydroxy aryl Schiff bases: A combined NMR and QTAIM analysis

Within the framework of the density functional theory approach, we studied the relationship between the chemical nature of intramolecular hydrogen bonds (HBs) and nuclear magnetic resonance (NMR) parameters, J-couplings and ¹H-chemical shifts [δ(¹H)], of the atoms involved in such bonds in o-hydroxyaryl Schiff bases during the proton transfer process. For the first time, the shape of the dependence of the degree of covalence in HBs on ¹J(N-H), ¹J(O-H), ^2hJ(O-N), and δ(¹H) during the proton transfer process in o-hydroxyaryl Schiff bases was analyzed. Parameters obtained from Bader's theory of atoms in molecules were used to assess the dependence of covalent character in HBs with both the NMR properties. The influence of π-electronic delocalization on ^2hJ(N-O) under the proton transfer process was investigated. ^2hJ(O-N) in a Mannich base was also studied in order to compare the results with an unsaturated system. In addition, substituent effects on the phenolic ring were investigated. Our results indicate that the covalent character of HBs on both sides of the transition state undergoes a smooth exponential increase as the δ(¹H) moves downfield. The degree of covalence of the N⋯H (O⋯H) bond increases linearly as ¹J(N-H) (¹J(O-H)) becomes more negative, even after reaching the transition state. Non-vanishing values of spin dipolar (SD) and paramagnetic spin orbital terms of ^2hJ(O-N) show that π-electronic delocalization has a non-negligible effect on tautomeric equilibrium and gives evidence of the presence of the resonance assisted HB.Variation of the SD term of ^2hJ(O-N) follows a similar pattern as the change in the para-delocalization aromaticity index of the chelate ring.

A review of NMR analysis in polysaccharide structure and conformation: Progress, challenge and perspective

Nuclear magnetic resonance (NMR) has been widely used as an analytical chemistry technique to investigate the molecular structure and conformation of polysaccharides. Combined with 1D spectra, chemical shifts and coupling constants in both homo- and heteronuclear 2D NMR spectra are able to infer the linkage and sequence of sugar residues. Besides, NMR has also been applied in conformation, quantitative analysis, cell wall in situ, degradation, polysaccharide mixture interaction analysis, as well as carbohydrates impurities profiling. This review summarizes the principle and development of NMR in polysaccharides analysis, and provides NMR spectra data collections of some common polysaccharides. It will help to promote the application of NMR in complex polysaccharides of biochemical interest, and provide valuable information on commercial polysaccharide products.

High-resolution 2-D NMR spectroscopy based on the Radon transform and pure shift technique for studying chemical shifts perturbations

Chemical shift plays an important role in molecular analysis. However, chemical shifts are influenced by temperature, solvent concentration, pressure, and so forth. Therefore, measuring chemical shift perturbations caused by these factors is helpful to molecular studies. A new form of 2-D spectroscopy (projection spectroscopy) has been introduced whose indirect dimension is derived by implementing the Radon transform on a series of conventional 1-D proton spectra and indicates such perturbations. However, signal overlap may exist in the conventional 1-D spectra and hence in the resulting projection spectra, hampering clear multiplet analysis and accurate extraction of perturbations. Here, the pure shift decoupling technique is employed to obtain clearer projection spectrum with higher spectral resolution. The combination of pure shift technique and the Radon transform is helpful to accurately extract chemical shift perturbations. It is believed that this application will open up a vast prospect for molecular analysis.

Impact of nanoconfinement on acetylacetone Equilibria in Ordered Mesoporous Silicates

Nanoconfinement is one of the most intriguing nanoscale effects and affects several physical and chemical properties of molecules and materials, including viscosity, reaction kinetics, and glass transition temperature. In this work, liquid nuclear magnetic resonance (NMR) was used to analyze the behavior of 2,4-pentadienone in ordered mesoporous materials with a pore diameter of between 3 and 10 nm. The liquid NMR results showed meaningful changes in the hydrogen chemical shift and the keto-enol chemical equilibrium, which were associated with the pore diameter, allowing the authors to observe the effects of nanoconfinement. An interesting phenomenon was observed where the chemical equilibria of 2,4-pentadienone confined in a mesoporous material with a pore diameter of 3.5 nm was similar to that obtained with free (bulk) 2,4-pentadienone in larger pore materials. Another interesting result was observed for the enthalpy and entropy of the tautomeric equilibria of 2,4-pentadienone confined in mesoporous materials with a 5.5 nm pore diameter being -7.9 kJ mol^-1 and -15.9 J mol^-1.K. These values are similar to those obtained by dimethyl sulfoxide. This phenomenon indicates the possible use of ordered mesoporous materials as a reaction substitute in organic solvents. It was further observed that while the values of enthalpy (ΔH) and entropy (ΔS) had been modified by confinement, the Gibbs free energy (ΔG) value remained closer to that observed in free (bulk) 2,4-pentadienone. It is expected that this study will help in understanding the effects of nanoconfinement and provide a simple method to employ NMR techniques to analyze these phenomena.

MAS NMR Investigation of Molecular Order in an Ionic Liquid Crystal

The structure and molecular order in the thermotropic ionic liquid crystal (ILC), [choline][geranate(H)octanoate], an analogue of Choline And GEranate (CAGE), which has potential for use as a broad-spectrum antimicrobial and transdermal and oral delivery agent, were investigated by magic-angle spinning (MAS) nuclear magnetic resonance (NMR), polarizing optical microscopy, small-angle X-ray scattering (SAXS), and mass spectrometry. Mass spectrometry and the ¹H NMR chemical shift reveal that CAGE-oct is a dynamic system, with metathesis (the exchange of interacting ions) and hydrogen exchange occurring between hydrogen-bonded/ionic complexes such as [(choline)(geranate)(H)(octanoate)], [(choline)(octanoate)₂(H)], and [(choline)(geranate)₂(H)]. These clusters, which are shown by mass spectrometry to be significantly more stable than expected for typical electrostatic ion clusters, involve hydrogen bonding between the carboxylic acid, carboxylate, and hydroxyl groups, with rapid hydrogen bond breaking and re-formation observed to average the ¹H chemical shifts. The formation of a partial bilayer liquid crystal (LC) phase was identified by SAXS and polarizing optical microscopy at temperatures below ∼293 K. The occurrence of this transition close to room temperature could be utilized as a potential temperature-induced "switch" of the anisotropic properties for particular applications. The presence of an isotropic component of approximately 23% was observed to coexist with the LC phase, as detected by polarizing optical microscopy and quantified by both ¹H-¹³C dipolar-chemical shift correlation (DIPSHIFT) and ¹H double-quantum (DQ) MAS NMR experiments. At temperatures above the LC-to-isotropic transition, intermediate-range order (clustering of polar and nonpolar domains), a feature of many ILs, persists. Site-specific order parameters for the LC phase of CAGE-oct were obtained from the MAS NMR measurement of the partially averaged ¹³C-¹H dipolar couplings (D_CH) by cross-polarization (CP) build-up curves and DIPSHIFT experiments, and ¹H-¹H dipolar couplings (D_HH) by double-quantum (DQ) build-up curves. The corresponding order parameters, S_CH and S_HH, are in the range 0-0.2 and are lower compared to those for smectic (i.e., layered) phases of conventional nonionic liquid crystals, resembling those of lamellar phases formed by lyotropic surfactant-solvent systems.

Spin System Modeling of Nuclear Magnetic Resonance Spectra for Applications in Metabolomics and Small Molecule Screening

The exceptionally rich information content of nuclear magnetic resonance (NMR) spectra is routinely used to identify and characterize molecules and molecular interactions in a wide range of applications, including clinical biomarker discovery, drug discovery, environmental chemistry, and metabolomics. The set of peak positions and intensities from a reference NMR spectrum generally serves as the identifying signature for a compound. Reference spectra normally are collected under specific conditions of pH, temperature, and magnetic field strength, because changes in conditions can distort the identifying signatures of compounds. A spin system matrix that parametrizes chemical shifts and coupling constants among spins provides a much richer feature set for a compound than a spectral signature based on peak positions and intensities. Spin system matrices expand the applicability of NMR spectral libraries beyond the specific conditions under which data were collected. In addition to being able to simulate spectra at any field strength, spin parameters can be adjusted to systematically explore alterations in chemical shift patterns due to variations in other experimental conditions, such as compound concentration, pH, or temperature. We present methodology and software for efficient interactive optimization of spin parameters against experimental 1D-¹H NMR spectra of small molecules. We have used the software to generate spin system matrices for a set of key mammalian metabolites and are also using the software to parametrize spectra of small molecules used in NMR-based ligand screening. The software, along with optimized spin system matrix data for a growing number of compounds, is available from http://gissmo.nmrfam.wisc.edu/ .

Hydrogen Peroxide Complex of Zinc

Metal(H2O2) complexes have been implicated in kinetic and computational studies but have never been observed. Accordingly, H2O2 has been described as a very weak ligand. We report the first metal(H2O2) adduct, which is made possible by incorporating intramolecular hydrogen-bonding interactions with bound H2O2. This Zn(II)(H2O2) complex decays in solution by a second-order process that is slow enough to enable characterization of this species by X-ray crystallography. This report speaks to the intermediacy of metal(H2O2) adducts in chemistry and biology and opens the door to exploration of these species in oxidation catalysis.

(15)N and (1)H Solid-State NMR Investigation of a Canonical Low-Barrier Hydrogen-Bond Compound: 1,8-Bis(dimethylamino)naphthalene

Strong or low-barrier hydrogen bonds have often been proposed in proteins to explain enzyme catalysis and proton-transfer reactions. So far (1)H chemical shifts and scalar couplings have been used as the main NMR spectroscopic signatures for strong H-bonds. In this work, we report simultaneous measurements of (15)N and (1)H chemical shifts and N-H bond lengths by solid-state NMR in (15)N-labeled 1,8-bis(dimethylamino)naphthalene (DMAN), which contains a well-known strong NHN H-bond. We complexed DMAN with three different counteranions to examine the effects of the chemical environment on the H-bond lengths and chemical shifts. All three DMAN compounds exhibit significantly elongated N-H distances compared to the covalent bond length, and the (1)H(N) chemical shifts are larger than ∼17 ppm, consistent with strong NHN H-bonds in the DMAN cation. However, the (15)N and (1)H chemical shifts and the precise N-H distances differ among the three compounds, and the (15)N chemical shifts show opposite dependences on the proton localization from the general trend in organic compounds, indicating the significant effects of the counteranions on the electronic structure of the H-bond. These data provide useful NMR benchmarks for strong H-bonds and caution against the sole reliance on chemical shifts for identifying strong H-bonds in proteins since neighboring side chains can exert influences on chemical shifts similar to those of the bulky organic anions in DMAN. Instead, N-H bond lengths should be measured, in conjunction with chemical shifts, as a more fundamental parameter of H-bond strength.

Hydrogen bond distinction and activation upon catalytic etherification of hydroxyl compounds

The slope (A) of the linear correlation between ln δand 1/Tcan distinguish intra- and intermolecular H-bonds and predict their reactivities.

NMR of glycans: shedding new light on old problems

The diversity in molecular arrangements and dynamics displayed by glycans renders traditional NMR strategies, employed for proteins and nucleic acids, insufficient. Because of the unique properties of glycans, structural studies often require the adoption of a different repertoire of tailor-made experiments and protocols. We present an account of recent developments in NMR techniques that will deepen our understanding of structure-function relations in glycans. We open with a survey and comparison of methods utilized to determine the structure of proteins, nucleic acids and carbohydrates. Next, we discuss the structural information obtained from traditional NMR techniques like chemical shifts, NOEs/ROEs, and coupling-constants, along with the limitations imposed by the unique intrinsic characteristics of glycan structure on these approaches: flexibility, range of conformers, signal overlap, and non-first-order scalar (strong) coupling. Novel experiments taking advantage of isotopic labeling are presented as an option for overcoming spectral overlap and raising sensitivity. Computational tools used to explore conformational averaging in conjunction with NMR parameters are described. In addition, recent developments in hydroxyl detection and hydrogen bond detection in protonated solvents, in contrast to traditional sample preparations in D2O for carbohydrates, further increase the tools available for both structure information and chemical shift assignments. We also include previously unpublished data in this context. Accurate determination of couplings in carbohydrates has been historically challenging due to the common presence of strong-couplings. We present new strategies proposed for dealing with their influence on NMR signals. We close with a discussion of residual dipolar couplings (RDCs) and the advantages of using (13)C isotope labeling that allows gathering one-bond (13)C-(13)C couplings with a recently improved constant-time COSY technique, in addition to the commonly measured (1)H-(13)C RDCs.

Are there really low-barrier hydrogen bonds in proteins? The case of photoactive yellow protein

For a long time, low-barrier hydrogen bonds (LBHBs) have been proposed to exist in many enzymes and to play an important role in their catalytic function, but the proof of their existence has been elusive. The transient formation of an LBHB in a protein system has been detected for the first time using neutron diffraction techniques on a photoactive yellow protein (PYP) crystal in a study published in 2009 (Yamaguchi, S.; et al. Proc. Natl. Acad. Sci. U.S.A. 2009, 106, 440-444). However, very recent theoretical studies based on electronic structure calculations and NMR resonance experiments on PYP in solution (Saito, K.; et al. Proc. Natl. Acad. Sci. U.S.A. 2012, 109, 167-172) strongly indicate that there is not such an LBHB. By means of electronic structure calculations combined with the solution of the nuclear Schrödinger equation, we analyze here under which conditions an LBHB can exist in PYP, thus leading to a more reasonable and conciliating understanding of the above-mentioned studies.

Geometry and Cooperativity Effects in Adenosine−Carboxylic Acid Complexes

NMR experiments and theoretical investigations were performed on hydrogen bonded complexes of specifically 1- and 7-15N-labeled adenine nucleosides with carboxylic acids. By employing a freonic solvent of CDClF2 and CDF3, NMR spectra were acquired at temperatures as low as 123 K, where the regime of slow hydrogen bond exchange is reached and several higher-order complexes were found to coexist in solution. Unlike acetic acid, chloroacetic acid forms Watson-Crick complexes with the proton largely displaced from oxygen to the nitrogen acceptor in an ion pairing structure. Calculated geometries and chemical shifts of the proton in the hydrogen bridge favorably agree with experimentally determined values if vibrational averaging and solvent effects are taken into account. The results indicate that binding a second acidic ligand at the adenine Hoogsteen site in a ternary complex weakens the hydrogen bond to the Watson-Crick bound carboxylic acid. However, substituting a second adenine nucleobase for a carboxylic acid in the trimolecular complex leads to cooperative binding at Watson-Crick and Hoogsteen faces of adenosine.

Determination of the Temperature Dependence of the H−D Spin−Spin Coupling Constant and the Isotope Effect on the Proton Chemical Shift for the Compressed Dihydride Complex [Cp*Ir(P−P)H2]2+

Complex [CpIr(dmpm)H(2)](2+) (dmpm = bis(dimethylphosphino)methane) has been reported to display temperature-dependent spin-spin coupling constant ((1)J(HD)) and isotope effect on the (1)H NMR chemical shift (Deltadelta). A combined electronic structure density functional theory + quantum nuclear dynamics study is used to determine from first-principles the unusual temperature dependence of the spin-spin coupling constant. It is found that the potential energy surface describing the motion of the Ir-H(2) unit has a deeper minimum in the dihydride region and is characterized by important anharmonicities. These anomalies affect the nature of the vibrational states of the unit and are the main reason for the unusual temperature dependence of (1)J(HD) and Deltadelta. These results suggest experimental tests to identify compressed dihydride transition metal complexes.

Synthesis and Properties of Compressed Dihydride Complexes of Iridium: Theoretical and Spectroscopic Investigations

Reaction of [Cp*Ir(P-P)Cl][B(C6F5)4] (P-P = bisdimethydiphosphinomethane (dmpm), bisdiphenyldiphosphinomethane (dppm)) with [Et3Si][B(C6F5)4] in methylene chloride under 1 atm of hydrogen gas affords the dicationic compressed dihydride complexes [Cp*Ir(P-P)H2][B(C6F5)4]2. These dicationic complexes are highly acidic and are very readily deprotonated to the corresponding monohydride cations. When the preparative reaction is carried out under HD gas, the hydride resonance exhibits JHD = 7-9 Hz, depending upon the temperature of observation, with higher values of JHD observed at higher temperatures. A thermally labile rhodium analogue, [CpRh(dmpm)(H2)][B(C6F5)4]2, was prepared similarly. A sample prepared with HD gas gave JHD = 31 Hz and J(HRh) = 31 Hz, allowing the Rh complex to be identified as a dihydrogen complex. Quantum dynamics calculations on a density functional theory (DFT) potential energy surface have been used to explore the structure of the Ir complexes, with particular emphasis on the nature of the potential energy surface governing the interaction between the two hydride ligands and the Ir center.

Isotopomeric Polymorphism

The complex formed between 4-methylpyridine and pentachlorophenol (4MPPCP) crystallizes in a triclinic space group. If the same complex is synthesized from deuterated pentachlorophenol (4MPPCP-d1), it crystallizes in an entirely different monoclinic polymorph. Solid-state NMR of samples synthesized with a full range of deuteration levels, crystallized from solution or the melt, and in the presence or absence of seeds confirms that the isotopomers indeed have different thermodynamically stable crystal structures. The phenomenon is apparently due to very different hydrogen bond strengths between the two polymorphs.

Density Functional Calculation of the 2D Potential Surface and Deuterium Isotope Effect on 13C Chemical Shifts in Picolinic Acid N-Oxide. Comparison with Experiment

2D free energy surfaces V = V(rOH, rO...O) for the intramolecular H-bond in the title compound were calculated by the DFT method and used in the calculation of primary and secondary chemical shifts of the compound dissolved in chloroform and acetonitrile. Solvent effects were accounted for by the SCRF/PCM method. The corresponding two-dimensional chemical shift surfaces with included solvent reaction field were obtained using the Continuous Set of Gauge Transformations approach at the B3LYP/6-311+G(2d,2p) level of theory. The chemical shifts were estimated as quantum averages along the two internal coordinates in the hydrogen bond and along several vibrational levels according to the Boltzmann distribution at room temperature. Fairly good agreement between the experimental and calculated isotope effects was obtained. 1D and 2D NMR spectra of solutions of picolinic acid N-oxide and its deuterated analogue were recorded and assigned.

Theoretical study of the symmetry of the (OH...O)- hydrogen bonds in vinyl alcohol-vinyl alcoholate systems

The interactions between substituted vinyl alcohols and vinyl alcoholates (X = NH(2), H, F, Cl, CN) are studied at the B3LYP/6-311++G(d,p) level of theory. In a first step, the conformation of the monomers is investigated and the proton affinities (PA(A(-))) of the enolates are calculated. The enols and enolates are held together by strong (OH...O)(-) hydrogen bonds, the hydrogen bond energies ranging from 19.1 to 34.6 kcal mol(-1). The optimized O...O distances are between 2.414 and 2.549 A and the corresponding OH distances from 1.134 and 1.023 A. The other geometry parameters such as C[double bond]C or CO distances also indicate that, in the minimum energy configuration, the hydrogen bonds are characterized by a double well potential. The Mulliken charges on the different atoms of the proton donors and proton acceptors and the frequencies of the nu(OH) stretching vibrations agree with this statement. All the data indicate that the hydrogen bonds are the strongest in the homomolecular complexes. The transition state for hydrogen transfer is located with the transition barrier estimated to be about zero. Upon addition of the zero-point vibration energies to the total potential energy, the barrier vanishes. This is a characteristic feature of low-barrier hydrogen bonds (LBHBs). The hydrogen bond energies are correlated to the difference 1.5 PA(AH) - PA(A(-)). The correlation predicts different energies for homomolecular hydrogen bonds, in agreement with the theoretical calculations. Our results suggest that a PA (or pK(a)) match is not a necessary condition for forming LBHBs in agreement with recent data on the intramolecular hydrogen bond in the enol form of benzoylacetone (J. Am. Chem. Soc. 1998, 120, 12117).

Quantum chemical calculations on structural models of the catalytic site of chymotrypsin: comparison of calculated results with experimental data from NMR spectroscopy

Hybrid density functional quantum mechanical calculations were used to study the strength of the hydrogen bond between His(57) N(delta)(1) and Asp(102) O(delta)(1) in chymotrypsin and how it changes along the reaction coordinate. Comparison of experimental shifts with the results of chemical shift calculations on a variety of small molecules, including species containing very strong hydrogen bonds, has validated the overall approach and provided the means for calibrating and correcting the calculated values. Models of the active site of chymotrypsin in its resting state and tetrahedral intermediate state were derived from high-resolution X-ray structures. The distance between His(57) N(delta)(1) and Asp(102) O(delta)(1) in each model was varied between 2.77 A (weak hydrogen bond) and 2.50 A (extremely strong hydrogen bond), and the one-dimensional potential energy surface of the hydrogen-bonded proton (or deuteron/triton) was determined. The zero-point energy, probability distribution, and chemical shift were determined for each distance. Calculated values for NMR chemical shifts, NMR chemical shift differences between (1)H and (3)H, and (2)H/(1)H fractionation factors were compared with published experimental values. Energies provided by the calculations indicated that the hydrogen bond between His(57) N(delta)(1) and Asp(102) O(delta)(1) in the chymotrypsin active site increases in strength by 11 kcal mol(-)(1) in going from the resting state of the enzyme to the tetrahedral intermediate state. This result confirms the hypothesis that the strengthened hydrogen bond plays an important role in lowering the energy of the transition state and, hence, in the catalytic efficiency of the enzyme. Models of the transition state that best fit the experimental data are consistent with a "strong" hydrogen bond between His(57) N(delta)(1) and Asp(102) O(delta)(1) but apparently not a "low-barrier" or "very strong" hydrogen bond.

The (1)H NMR chemical shift for the hydroxy proton of 4-(dimethylamino)-2'-hydroxychalcone in chloroform: a theoretical approach to its inverse dependence on the temperature

[reaction: see text] The inverse dependence of the chemical shift on the temperature experimentally found for the phenolic proton of 4-(dimethylamino)-2'-hydroxychalcone (DMAHC) is theoretically studied. As the temperature decreases, the solvent dielectric constant epsilon increases and the zwitterionic resonance form is more stabilized. Electronic calculations at the DFT level of theory were performed by immersing the solute DMAHC in chloroform cavities of different epsilon values. The values of the calculated chemical shifts for DMAHC as a function of epsilon show that the growing contribution of the zwitterionic structure justifies the experimental results.

A QM/MM study of the racemization of vinylglycolate catalyzed by mandelate racemase enzyme

The experimentally postulated mechanism for the interconversion between (S)-vinylglycolate and (R)-vinylglycolate catalyzed by mandelate racemase enzyme consists of a two-step quite symmetric process through a dianionic enolic intermediate that is formed after the abstraction of the alpha-proton of vinylglycolate by a basic enzymatic residue and is then reprotonated by another residue. The challenging problem behind this reaction is how the enzyme manages to stabilize such an intermediate, that is, how it lowers enough the high pK(a) of the alpha-proton for the reaction to take place. The QM/MM simulations performed in this paper indicate that catalysis is based on the stabilization of the negative charge developed on the substrate along the reaction. We have identified three different reaction mechanisms starting from different quasi-degenerate structures of the substrate-enzyme complex. In two of them the stabilizing role is done by means of a catalytic proton transfer that avoids the formation of a dianionic intermediate, and they involve six steps instead of the two experimentally proposed. On the contrary, the third mechanism passes through a dianionic species stabilized by the concerted approach of a protonated enzymatic residue during the proton abstraction. The potential energy barriers theoretically found along these mechanisms are qualitatively in good agreement with the experimental free energy barriers determined for racemization of vinylglycolate and mandelate. The theoretical study of the effect of the mutation of Glu317 by Gln317 in the kinetics of the reaction reveals the important role in the catalysis of the hydrogen bond formed by Glu317 in the native enzyme, as only one of the mechanisms, the slower one, is able to produce the racemization in the active site of the mutant. However, we have found that this hydrogen bond is not an LBHB within our model.

VT 1H NMR Investigations of Resonance-Assisted Intramolecular Hydrogen Bonding in 4-(Dimethylamino)-2‘-hydroxychalcone

[reaction: see text] Resonance-assisted intramolecular hydrogen bonding in both polar aprotic and nonpolar solutions of 4-(dimethylamino)-2'-hydroxychalcone (DMAHC) has been investigated by variable-temperature proton NMR spectroscopy. In both nonpolar and polar solvents, the signal for the phenolic hydrogen moves downfield as the temperature is lowered. In each solvent system studied, a linear relationship between chemical shift and temperature was observed.