Structural criteria for the rational design of selective ligands: convergent hydrogen bonding sites for the nitrate anion

A large number of crystal structures are analyzed to characterize the structural aspects of hydrogen bonding interactions with the NO(3)(-) anion. Further insight is provided by the use of electronic structure calculations to determine stable geometries and interaction energies for NO(3)(-) complexes with several simple hydrogen bond donor groups, including water, methanol, N-methylform-amide, and methane. The results establish the existence of a clear set of structural criteria for the rational design of molecular receptors that complex the NO(3)(-) anion through hydrogen bonding interactions.     

Computational analysis of the potential energy surfaces of glycerol in the gas and aqueous phases: effects of level of theory, basis set, and solvation on strongly intramolecularly hydrogen-bonded systems.

The 126 possible conformations of 1,2,3-propanetriol (glycerol) have been studied by ab initio molecular orbital and density functional theory calculations in the gas and aqueous phases at multiple levels of theory and basis sets. The partial potential energy surface for glycerol as well as an analysis of the conformational properties and hydrogen-bonding trends in both phases have been obtained. In the gas phase at the G2(MP2) and CBS-QB3 levels of theory, the important, low-energy conformers are structures 100 and 95. In the aqueous phase at the SM5.42/HF/6-31G* level of theory, the lowest energy conformers are structures 95 and 46. Boltzmann distributions have been determined from these high-level calculations, and good agreement is observed when these distributions are compared to the available experimental data. These calculations indicate that the enthalpic and entropic contributions to the Gibbs free energy are important for an accurate determination of the conformational and energetic preferences of glycerol. Different levels of theory and basis sets were used in order to understand the effects of nonbonded interactions (i.e., intramolecular hydrogen bonding). The efficiency of basis set and level of theory in dealing with the issue of intramolecular hydrogen bonding and reproducing the correct energetic and geometrical trends is discussed, especially with relevance to practical computational methods for larger polyhydroxylated compounds, such as oligosaccharides.     

4-Methylpseudoproline derived from α-methylserine - synthesis and conformational studies

This paper presents the synthesis and solution conformational studies of the tripeptides Fmoc-Ala-(R)-(αMe)Ser(Ψ(H,H)Pro)-Ala-OBu(t) (6a) and Fmoc-Ala-(S)-(αMe)Ser(Ψ(H,H)Pro)-Ala-OBu(t) (6b). Additionally, the X-ray structure of 6a is given. NMR analysis corroborated by theoretical calculations (XPLOR) shows that in both peptides the amide bond between pseudoproline and the preceding amino acid is in the trans conformation. The same amide bond geometry was observed in the crystal state of 6a. The latter is additionally influenced by the presence of two symmetrically independent molecules in an asymmetric unit. Both molecules adopt a conformation which resembles β-turn type II, stabilized by hydrogen bonding. The conformational preferences and prolyl cis-trans isomerization of Ac-(αMe)Ser(Ψ(H,H)Pro)-NHMe (7) were explored at the IEFPCM/B3LYP/6-31+G(d) level of theory in vacuum, water and chloroform. It has been shown that the trans isomer predominates in water solutions and the cis isomer is preferred in chloroform. The conformation of 7 is down-puckered independently of the geometry of the amide bonds, with lower puckering in the transition state of the cis-trans isomerization.     

Ab initio prediction of the vibrational spectrum of the hydrogen bonded complex O2N...HONO2

The changes in the structural parameters and vibrational characteristics (vibrational frequencies, infrared intensities and Raman activities) arising from the hydrogen bonding between NO(2) and HONO(2) have been studied employing ab initio 6-31G(d, p)/UHF and 6-31+G(d, p)/UHF, and B3LYP/6-31G(d, p) calculations. The charge rearrangement upon hydrogen bonding have been, estimated using the Mulliken population analyses. It was established that the complexation between NO(2) and HONO(2) leads to changes in the structural parameters and the vibrational characteristics of the monomers. The most sensitive to the hydrogen bond formation are the vibrational characteristics of the normal modes of the monomer bonds participating in the hydrogen bonding. The predicted shifts in the vibrational frequencies by ab initio and B3LYP/6-31G(d, p) calculations are in very good agreement with the experimentally observed, which is an evidence for the reliance of the studied structure.     

Self‐Complementary Dimers of Oxalamide‐Functionalized Resorcinarene Tetrabenzoxazines

Self‐complementarity is a useful concept in supramolecular chemistry, molecular biology and polymeric systems. Two resorcinarene tetrabenzoxazines decorated with four oxalamide groups were synthesized and characterized. The oxalamide groups possessed self‐complementary hydrogen bonding sites between the carbonyls and amide groups. The self‐complementary nature of the oxalamide groups resulted in self‐included dimeric assemblies. The hydrogen bonding interactions within the tetrabenzoxazines gave rise to the formation of dimers, which were confirmed by single‐crystal X‐ray diffractions analysis and supported by NMR spectroscopy and mass spectrometry. The self‐included dimers were connected by numerous and strong intermolecular N?H???O and C?H???O hydrogen bonds supplemented with C?H???π interactions, forming one‐dimensional polymers, which were then further linked into three‐dimensional networks.     

Modeling the interplay of inter- and intramolecular hydrogen bonding in conformational polymorphs

The predicted stability differences of the conformational polymorphs of oxalyl dihydrazide and ortho-acetamidobenzamide are unrealistically large when the modeling of intermolecular energies is solely based on the isolated-molecule charge density, neglecting charge density polarization. Ab initio calculated crystal electron densities showed qualitative differences depending on the spatial arrangement of molecules in the lattice with the greatest variations observed for polymorphs that differ in the extent of inter- and intramolecular hydrogen bonding. We show that accounting for induction dramatically alters the calculated stability order of the polymorphs and reduces their predicted stability differences to be in better agreement with experiment. Given the challenges in modeling conformational polymorphs with marked differences in hydrogen bonding geometries, we performed an extensive periodic density functional study with a range of exchange-correlation functionals using both atomic and plane wave basis sets. Although such electronic structure methods model the electrostatic and polarization contributions well, the underestimation of dispersion interactions by current exchange-correlation functionals limits their applicability. The use of an empirical dispersion-corrected density functional method consistently reduces the structural deviations between the experimental and energy minimized crystal structures and achieves plausible stability differences. Thus, we have established which types of models may give worthwhile relative energies for crystal structures and other condensed phases of flexible molecules with intra- and intermolecular hydrogen bonding capabilities, advancing the possibility of simulation studies on polymorphic pharmaceuticals.     

L-arginine trifluoroacetate salt bridges in its solid state compound: the low-temperature three dimensional structural determination of L-arginine bis(trifluoroacetate) crystal and its vibrational spectral analysis

Structural varieties of L-arginine trifluoroacetate (abbreviated as LATF) and L-arginine bis(trifluoroacetate), LABTF, in the solid state compounds were observed and analyzed by the nuclear magnetic resonance (NMR) spectroscopy. The guanidinium-carboxylate interaction plays an important role involving in the crystal structure construction. Conformational changes of L-Arg(+) and L-Arg(2+) cations result from the intrinsic structural difference by hydrogen bonding and electrostatic interactions. The low-temperature structure of its crystalline salt, L-arginine bis(trifluoroacetate), was determined to describe the hydrogen bonding interactions. In comparison with the crystal structure at room temperature, the low-temperature L-Arg(2+) cations present tiny conformational difference and the rotational disorder of CF(3) group disappears. FT-IR and Raman spectra were investigated and hydrogen bonding interactions were analyzed on the basis of its vibrational spectra. Results indicate that this type interaction is greatly contributive to the structural features and vibrational spectral properties.     

Conformational behaviour, hydrogen bond competition and intramolecular dynamics in vanillin derivatives: acetovanillone and 6-hydroxy-3-methoxyacetophenone

The conformational landscape of the structural isomers acetovanillone (apocynin, AV) and 6-hydroxy-3-methoxyacetophenone (HMAP) has been investigated in a supersonic jet using Fourier transform microwave spectroscopy. Two conformers have been detected in the jet-cooled expansion for each molecule (s-cis and s-trans in AV; s-trans and a-trans for HMAP), differing in the relative orientation of the acetyl and methoxy groups. Both molecules are stabilized by O-H···O or O-H···O=C hydroxyl intramolecular hydrogen bonds, either constraining the local conformations of the methoxy group in AV, or that of the acetyl group in HMAP. Internal rotation splittings have been observed in both conformers of each molecule, originated by the acetyl group, that yield information on the influence of the intramolecular hydrogen bonds on the methyl torsion. The similar internal rotation barriers in both molecules (6.6 and 7.4 kJ mol(-1) in AV; 7.3 and 7.0 kJ mol(-1) in HMAP) suggest that the acetyl torsion is only slightly affected by intramolecular hydrogen bonding. The absence of torsional tunnellings due to the methoxy group indicates torsional barriers above 10.2 and 8.9 kJ mol(-1) for AV conformers, 10.1 and 10.4 kJ mol(-1) for HMAP. Conformational ratios and relative free energies have been estimated from relative intensity measurements of the spectral lines. Ab initio (MP2) and density functional calculations using the recent M05-2X empirical functional have been used to aid the experimental work in describing the structures, internal rotation barriers and isomerization potentials.     

Hydrogen bonding and solvatochromatic shift of the lowest 1(n, pi*) excitation of s-tetrazine in its hydrated clusters and dilute solutions

Hydrogen bonding involving azine and its derivatives such as nucleic bases is very important for understanding the structure and function of biological systems. In this work, we have investigated the hydrogen bonding structures of the hydrated cluster and dilute aqueous solution of s-tetrazine using computer simulation techniques, and evaluated the absorption and fluorescence shifts of the lowest 1(n, pi*) excitation of s-tetrazine solution using our solvent shift method. For the s-tetrazine-water cluster, a linear orthodox hydrogen bond arrangement is predicted in both ground and excited states with small structural and energetic differences, and a bifurcated hydrogen bond isomerization is anticipated. Further, ab initio calculations have verified these conformations. For the s-tetrazine-water solution, a mixture of two hydrogen bonding arrangements is found to be in both ground and excited states, resulting in small magnitudes of absorption and fluorescence solvent shifts. This finalizes our series investigation of hydrogen bonding and solvent shifts of dilute azines in water.     

Hydrogen bonding topology influences gelating properties of malonamides

Preparation, structural characterisation and topology of hydrogen bonding networks of bis(phenylglycinol)malonamide, as well as its Cα mono- and dialkyl-substituted derivatives are described. Their hydrogen bonding motifs are described in view of their gelling properties. Topology of hydrogen bonding typical of malonamide gelators is compared with those of well-examined oxalamide gelators.     

The first oxalamide‐bridged ferrocene: Facile synthesis,preliminary conformational analysis and biological evaluation

tert ‐Butyl‐1′‐methoxycarbonyl‐1‐ferrocenecarbamate ( 1 ) was Boc‐deprotected to give free amine which underwent oxalyl chloride‐mediated dimerization. The structure of the so‐obtained oxalamide‐bridged ferrocene 2 was elucidated using infrared and NMR (¹H, ¹³C, COSY, NOESY, HSQC, HMBC) spectroscopies, crystal structure analysis, and electrospray ionization and high‐resolution mass spectrometry. The preliminary conformational analysis in solution suggested the intramolecular engagement of oxalamide protons, while single‐crystal analysis revealed an intermolecular hydrogen‐bonding pattern. Also, the effect of oxalamide‐bridged ferrocene 2 on cell viability of three human cell lines (HEK293T, HeLa and HepG2) was tested. In vitro screening revealed proliferative as well as cytotoxic effects of the tested compound in the applied concentration range (1–350 μM) on HEK293T and HepG2 cells. Stimulatory effect on cell growth was the most pronounced for normal HEK293T cells, while the highest cytotoxic effect was observed towards HeLa tumour cells and it was dose‐dependent. The observed dual biological activity of 2 implies its potential application in drug development.     

Conformational control and photoenolization of pyridine-3-carboxaldehydes in the solid state: stabilization of photoenols via hydrogen bonding and electronic control

We have investigated the solid-state photobehavior of a broad set of pyridine-3-carboxaldehydes 1-5. The introduction of a heteroatom into mesitaldehydes as in aldehydes 1 raises the question of conformational preference in the solid state. The preferred conformations have been unequivocally established from X-ray crystal structure analyses of two of the aldehydes, 1c and 2c; it is shown that intramolecular hydrogen bonding could be utilized to achieve conformational control. In contrast to mesitaldehydes, which undergo efficient photocyclization to benzocyclobutenols in the solid state, the heteroatom analogues 1b and 1c exhibit a perceptible color change (from colorless to pale yellow for 1b and yellow-orange for 1c) upon UV irradiation; the color attributed to (E)-enols is persistent for several hours. Continued irradiation leads to an intractable polymeric material. The AM1 calculations, which have been reliably applied to the thermal cyclization of xylylenols to benzocyclobutenols, reveal that the (E)-enols of 1 are more stable than those of the mesitaldehydes relative to their corresponding benzocyclobutenols. The stabilization is interpreted as arising from the possibility of engaging the heteroatom in resonance delocalization. That the contribution from such a role of the nitrogen atom is so pronounced is elegantly demonstrated by forming the fluoroborate salts; 1a-HBF(4) and 1b-HBF(4) readily exhibit highly red-shifted absorption upon exposure to UV radiation as a result of stabilization of the photoenols. Notably, such a remarkable stabilization via electronic control of the photoenols is unprecedented. All of the 2-methoxy- and 2-chloro-substituted aldehydes 2-5 exhibit photochromism. Ab initio calculations show that the methoxy group in aldehydes 2 and 3 stabilizes the (E)-enols via O[bond]H...O hydrogen bonding as compared to those of 1 by 5-6 kcal/mol relative to their corresponding benzocyclobutenols. Thus, the presence of methoxy and halo groups at position 2 serves not only to direct the formyl oxygen toward the methyl group for H-abstraction but also to stabilize the (E)-enols.     

Secondary structural preferences of 2,2-disubstituted pyrrolidine-4-carboxylic acid oligomers: beta-peptide foldamers that cannot form internal hydrogen bonds

We examine a new class of beta-peptides, 2,2-disubstituted pyrrolidine-4-carboxylic acid oligomers, and show that they manifest discrete conformational preferences despite the impossibility of internal hydrogen bonding. Numerous beta-peptide families have been described that display specific secondary structural preferences, but all of the conformations characterized in detail so far have contained internal hydrogen bonds. Internal hydrogen bonding is observed within the most common secondary structures of conventional peptides as well. Identifying foldamers in which shape control is independent of hydrogen bonding is significant in two ways. At a fundamental level, foldamers in this small but growing class are interesting because their shapes are controlled by distinctive networks of noncovalent forces. At a practical level, non-hydrogen bonded foldamers may be useful in biomedical applications because the low intrinsic polarity of their backbones may promote bioavailability.     

Probing the configurational space of a metalloprotein core: an ab initio molecular dynamics study of Duo Ferro 1 binuclear Zn cofactor

We present three theoretical models of various degree of completeness to explore the chemical phase space available to the Glu4His2Zn2 cofactor found in the four-helix bundle of de novo designed metalloprotein Duo Ferro 1. We have found that the planewave DFT geometry optimization of 94-atom Model I, which contains both the protein scaffold constraints as well as the second shell hydrogen bonding network, reproduces the crystal structure bonding with remarkable accuracy (0.34 A). Surprisingly, the geometry optimization of 66-atom Model II (lacking the second shell hydrogen bonding) and 48-atom Model III (being also free of the protein scaffold constraints) still result in the fidelity with the crystallographic structure (RMSDs 0.29 and 0.34 A, respectively). To examine whether these structures are close to the global minimum as well as to investigate various conformational transitions to which the di-Zn cofactor may be susceptible to, we have carried out a 10 ps Car-Parrinello Molecular Dynamics (CPMD) simulation of Model III. We suggest that weak hydrogen bonds between imidazole hydrogens and carboxylate oxygens modulate the dynamical behavior of the system. One part of the molecule was found to be rigid due to the particular H(imidazole)-O(carboxylate) interaction restricting both the motion of the imidazole ring as well as the terminal carboxylate conformational mobility. The second half of the system was very flexible demonstrating a coupling of a transient formation of H(imidazole)-O(carboxylate) bonds with the spinning of the imidazole ring and syn-anti isomerization of the terminal carboxylate group. In addition, two low-energy snapshots from the 10 ps CPMD run were quenched, and their geometries were optimized, leading to two new isomers 48 kJ/mol lower in energy than the one associated with the crystal structure. We suggest that periodic quenching of the CPMD simulation snapshots of a minimalist model may be used as an efficient method to generate a large number of competitive local minima, which may be consequently pruned by imposing the protein scaffold constraints as well as further tuned by the second shell hydrogen bonding network.     

Structure and vibrational spectrum of the hydrogen-bonded system between 4-tert-butylphenol and N-bases: Ab initio and DFT studies

The structural and vibrational characteristics of the hydrogen-bonded system between 1,5,7-triazabicyclo [4.4.0] dec-5-ene (TBD) and 4-tert-butylphenol have been investigated employing ab initio and DFT calculations at different basis sets. The calculations show that the optimized structure of the studied system is cyclic. The corrected values of the dissociation energy for the hydrogen-bonded complex have been calculated in order to estimate its stability. The influence of the hydrogen bonding on the properties of the monomers (TBD and 4-tert-butylphenol) has been investigated. The hydrogen bonding between TBD and 4-tert-butylphenol leads to changes in the structural (bond lengths and angles) and vibrational (vibrational frequencies and infrared intensities) characteristics of the monomers. It was established that the TBD molecule is considerably deformed upon hydrogen bonding, while the deformation of the 4-t-BuPhOH is smaller. In agreement with the experiment, the calculations show that the stretching O-H vibration from 4-tert-butylphenol is shifted to lower frequency upon hydrogen bonding. The predicted frequency shift Deltanu(O-H) (-338cm(-1)) is in very good agreement with the experimentally observed (-351cm(-1)). In the same time the IR intensity of the nu(O-H) increases dramatically in the hydrogen-bonded system.     

Protonation of the proximal histidine ligand in heme peroxidases

The heme peroxidases have a histidine group as the axial ligand of iron. This ligand forms a hydrogen bond to an aspartate carboxylate group by the other nitrogen atom in the side chain. The aspartate is not present in the globins and it has been suggested that it gives an imidazolate character to the histidine ligand. Quantum chemical calculations have indicated that the properties of the heme site strongly depend on the position of the proton in this hydrogen bond. Therefore, we have studied the location of this proton in all intermediates in the reaction mechanism, using a set of different quantum mechanical and combined experimental and computational methods. Quantum refinements of a crystal structure of the resting FeIII state in yeast cytochrome c peroxidase show that the geometric differences of the two states are so small that it cannot be unambiguously decided where the proton is in the crystal structure. Vacuum calculations indicate that the position of the proton is sensitive to the surroundings and to the side chains of the porphyrin ring. Combined quantum and molecular mechanics (QM/MM) calculations indicate that the proton prefers to reside on the His ligand in all states in the reaction mechanism of the peroxidases. QM/MM free energy perturbations confirm these results, but reduce the energy difference between the two states to 12-44 kJ/mol.     

Conformational Analysis of Free and Bound Retinoic Acid

The conformational profiles of unbound all-trans and 9-cis retinoic acid (RA) have been determined using classical and quantum mechanical calculations. Sixty-six all-trans-RA (ATRA) and forty-eight 9-cis-RA energy minimum conformers were identified via HF/6-31G* geometry optimizations in vacuo. Their relative conformational energies were estimated utilizing the M06, M06-2x and MP2 methods combined with the 6-311+G(d,p), aug-cc-pVDZ and aug-cc-pVTZ basis sets, as well as complete basis set MP2 extrapolations using the latter two basis sets. Single-point energy calculations performed with the M06-2x density functional were found to yield similar results to MP2/CBS for the low-energy retinoic acid conformations. Not unexpectedly, the conformational propensities of retinoic acid were governed by the orientation and arrangement of the torsion angles associated with the polyene tail. We also used previously reported QM/MM X-ray refinement results on four ATRA-protein crystal structures plus one newly refined 9-cis-RA complex (PDB ID 1XDK) in order to investigate the conformational preferences of bound retinoic acid. In the re-refined RA conformers the conjugated double bonds are nearly coplanar, which is consistent with the global minimum identified by the Omega/QM method rather than the corresponding crystallographically determined conformations given in the PDB. Consequently, a 91.3% average reduction of the local strain energy in the gas phase, as well as 92.1% in PCM solvent, was observed using the QM/MM refined structures versus the PDB deposited RA conformations. These results thus demonstrate that our QM/MM X-ray refinement approach can significantly enhance the quality of X-ray crystal structures refined by conventional refinement protocols, thereby providing reliable drug-target structural information for use in structure-based drug discovery applications.     

Solid state NMR spectroscopy as a precise tool for assigning the tautomeric form and proton position in the intramolecular bridges of o-hydroxy Schiff bases

Two analogous Schiff bases, (S,E)-2-((1-hydroxy-3-methyl-1,1-diphenylbutan-2-ylimino)methyl)phenol (1) and (S,Z)-2-hydroxy-6-((1-hydroxy-3-methyl-1,1-diphenylbutan-2-ylamino)methylene)cyclohexa-2,4-dienone (2), exist in the solid state as phenol-imine and keto-amine tautomers, respectively. Their crystal structures were solved using the X-ray diffraction method. Sample 1 forms orthorhombic crystals of space group P2(1)2(1)2(1), while 2 forms monoclinic crystals of space group P2(1). In each sample, one molecule is in the asymmetric unit of the crystal structure. One-dimensional and two-dimensional solid state NMR techniques were used for structure assignment and for inspection of the (13)C and (15)N δ(ii) of the chemical shift tensor (CST) values. NMR study indicates that the span (Ω = δ(11)-δ(33)) and the skew (κ = 3(δ(22)-δ(iso)/Ω) are extremely sensitive to change in the tautomeric form of the Schiff bases. Theoretical calculations of NMR shielding parameters for 1 and 2 and a model compound with reduced aliphatic residue were performed using the GIAO method with B3LYP functional and 6-311++g(d,p) basis sets. From comparative analysis of the experimental and theoretical parameters, it was concluded that the position of hydrogen in the intramolecular bridge has tremendous influence on (13)C and (15)N CST parameters. Inspection of Ω and κ parameters allowed for the establishment of the nature of the hydrogen bonding and the assignment of the equilibrium proton position in the intramolecular bridges in the solid state.     

A conformational analysis of 3′-azido-3′-deoxythymidine

An extensive conformational analysis of 3′-azido-3′-deoxythymidine (AZT) was performed at the semiempirical AM1 level with full relaxation of all geometric parameters and careful consideration of furan puckering and the rotational states of the thymine—furan, furan—azide, furan—methylene, and methylene—hydroxyl bonds. The search located 70 conformers, 21 of which have relative energies within 2.5 kcal/mol of the global minimum. Several geometric features, including various forms of hydrogen bonding, within this selected lowenergy subset were examined in terms of their relative contributions to the conformational states of AZT. Hydrogen bonding of thymine's position 2 carbonyl oxygen atom to the hydroxymethyl group (O₂? ;HO), which until recently has not been mentioned in the literature, is observed in a few low-energy AM1 conformations; however, this form is less favored at the AM1 level than the usually depicted modes involving the thymine moiety with the oxygen atoms of the hydroxyl and furan groups (H₆? ;OH and H₆? O_fur, as observed in the two crystallographically independent structures), as well as that involving the hydroxyl hydrogen and furan oxygen atoms (OH? O_fur, which also has not been mentioned for AZT in the literature until recently). The AM1-optimized geometries agree more closely with nuclear magnetic resonance data than with crystallographic structures and bear little resemblance to molecular mechanics results. The present study shows no evidence of a single dominant conformation or single structural parameter that determines AZT's conformational states. In contrast to our previous analogous study of cGMP, this computational study of AZT does not show strong evidence of a syn conformation with hydrogen bonding involving the base.