Quantum chemical topological analysis of hydrogen bonding in HX…HX and CH3X…HX dimers (X = Br, Cl, F)

We present a systematic investigation of the nature and strength of the hydrogen bonding in HX···HX and CH3X…HX (X = Br, Cl and F) dimers using ab initio MP2/aug-cc-pVTZ calculations in the framework of the quantum theory of atoms in molecules (QTAIM) and electron localisation functions (ELFs) methods. The electron density of the complexes has been characterised, and the hydrogen bonding energy, as well as the QTAIM and ELF parameters, is consistent, providing deep insight into the origin of the hydrogen bonding in these complexes. It was found that in both linear and angular HX…HX and CH3X…HX dimers, F atoms form stronger HB than Br and Cl, but they need short (∼2 Å) X…HX contacts.

Controlling the formation of ionic-liquid-based aqueous biphasic systems by changing the hydrogen-bonding ability of polyethylene glycol end groups

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Interplay of pi–pi stacking and hydrogen bonding in conjugated supramolecular systems studied by solid-state NMR

In this work supramolecular organic systems based on rigid pi-conjugated building blocks and flexible side chains were studied via solid-state NMR spectroscopy. Specifically, these studies focussed on phenylene ethynylene based macrocycles, polymer systems including polythiophenes, and rod-coil copolymers of oligo(p-benzamide) and poly(ethylene glycol). All systems were studied in terms of the local order and mobility. The central topic of this dissertation was to elucidate the role of the flexible side chains in interplay of different non-covalent interactions, like pi-pi-stacking and hydrogen bonding.Combining the results of this work, it can be concluded that the ratio of the rigid block and the attached alkyl side chains can be crucial for the design of an ordered pi-conjugated supramolecular system. Through alkyl side chains, it is also possible to introduce liquid-crystalline phases in the system, which can foster the local order of the system. Moreover in the studied system longer, unbranched alkyl side chains are better suited to stabilize the corresponding aggregation than shorter, branched ones.The combination of non-covalent interactions such as pi-pi-stacking and hydrogen bonding play an important role for structure formation. However, the effect of pi-pi-stacking interaction is much weaker than the effect of hydrogen bonding and is only observed in systems with a suitable local order. Hence, they are often not strong enough to control the local order. In contrast, hydrogen bonds predominantly influence the structural organization and packing. In comparison the size of the alkyl side chains is only of minor importance. The suppression of certain hydrogen bonds can lead to completely different structures and can induce a specific aggregation behavior. Thus, for the design of a supramolecular ordered system the presence of hydrogen bonding efficiently stabilizes the corresponding structure, but the ratio of hydrogen bond forming groups should be kept low to be able to influence the structure selectively.

Hydrogen bonding of nucleotide base pairs: Application of the PM3 method

The ability of the pm3 semiempirical quantum mechanical method to reproduce hydrogen bonding in nucleotide base pairs was assessed. Results of pm3 calculations on the nucleotides 2′-deoxyadenosine 5′-monophosphate (pdA), 2′-deoxyguanosine 5′-monophosphate (pdG), 2′-deoxycytidine 5′-monophosphate (pdC), and 2′-deoxythymidine 5′-monophosphate (pdT) and the base pairs pdA–pdT, pdG–pdC, and pdG(syn)–pdC are presented and discussed. The pm3 method is the first of the parameterized nddo quantum mechanical models with any ability to reproduce hydrogen bonding between nucleotide base pairs. Intermolecular hydrogen bond lengths between nucleotides displaying Watson–Crick base pairing are 0.1–0.2 Å less than experimental results. Nucleotide bond distances, bond angles, and torsion angles about the glycosyl bond (χ), the C4′C5′ bond (γ), and the C5′O5′ bond (β) agree with experimental results. There are many possible conformations of nucleotides. pm3 calculations reveal that many of the most stable conformations are stabilized by intramolecular CHO hydrogen bonds. These interactions disrupt the usual sugar puckering. The stacking interactions of a dT–pdA duplex are examined at different levels of gradient optimization. The intramolecular hydrogen bonds found in the nucleotide base pairs disappear in the duplex, as a result of the additional constraints on the phosphate group when part of a DNA backbone. Sugar puckering is reproduced by the pm3 method for the four bases in the dT–pdA duplex. pm3 underestimates the attractive stacking interactions of base pairs in a B-DNA helical conformation. The performance of the pm3 method implemented in SPARTAN is contrasted with that implemented in MOPAC. At present, accurate ab initio calculations are too timeconsuming to be of practical use, and molecular mechanics methods cannot be used to determine quantum mechanical properties such as reaction-path calculations, transition-state structures, and activation energies. The pm3 method should be used with extreme caution for examination of small DNA systems. Future parameterizations of semiempirical methods should incorporate base stacking interactions into the parameterization data set to enhance the ability of these methods.

Ability of the PM3 Quantum Mechanical Method to Model Intermolecular Hydrogen Bonding between Neutral Molecules

The PM3 semiempirical quantum-mechanical method was found to systematically describe intermolecular hydrogen bonding in small polar molecules. PM3 shows charge transfer from the donor to acceptor molecules on the order of 0.02-0.06 units of charge when strong hydrogen bonds are formed. The PM3 method is predictive; calculated hydrogen bond energies with an absolute magnitude greater than 2 kcal mol-' suggest that the global minimum is a hydrogen bonded complex; absolute energies less than 2 kcal mol-' imply that other van der Waals complexes are more stable. The geometries of the PM3 hydrogen bonded complexes agree with high-resolution spectroscopic observations, gas electron diffraction data, and high-level ab initio calculations. The main limitations in the PM3 method are the underestimation of hydrogen bond lengths by 0.1-0.2 for some systems and the underestimation of reliable experimental hydrogen bond energies by approximately 1-2 kcal mol-l. The PM3 method predicts that ammonia is a good hydrogen bond acceptor and a poor hydrogen donor when interacting with neutral molecules. Electronegativity differences between F, N, and 0 predict that donor strength follows the order F > 0 > N and acceptor strength follows the order N > 0 > F. In the calculations presented in this article, the PM3 method mirrors these electronegativity differences, predicting the F-H- - -N bond to be the strongest and the N-H- - -F bond the weakest. It appears that the PM3 Hamiltonian is able to model hydrogen bonding because of the reduction of two-center repulsive forces brought about by the parameterization of the Gaussian core-core interactions. The ability of the PM3 method to model intermolecular hydrogen bonding means reasonably accurate quantum-mechanical calculations can be applied to small biologic systems.

N-H center dot center dot center dot pi hydrogen-bonding and large-amplitude tipping vibrations in jet-cooled pyrrole-benzene

The N-H center dot center dot center dot pi hydrogen bond is an important intermolecular interaction in many biological systems. We have investigated the infrared (IR) and ultraviolet (UV) spectra of the supersonic-jet cooled complex of pyrrole with benzene and benzene-d(6) (Pyr center dot Bz, Pyr center dot Bz-d(6)). DFT-D density functional, SCS-MP2 and SCS-CC2 calculations predict a T-shaped and (almost) C(s) symmetric structure with an N-H center dot center dot center dot pi hydrogen bond to the benzene ring. The pyrrole is tipped by omega(S(0)) = +/- 13 degrees relative to the surface normal of Bz. The N center dot center dot center dot ring distance is 3.13 angstrom. In the S(1) excited state, SCS-CC2 calculations predict an increased tipping angle omega(S(1)) = +/- 21 degrees. The IR depletion spectra support the T-shaped geometry: The NH stretch is redshifted by -59 cm(-1), relative to the "free" NH stretch of pyrrole at 3531 cm(-1), indicating a moderately strong N-H center dot center dot center dot pi interaction. The interaction is weaker than in the (Pyr)(2) dimer, where the NH donor shift is -87 cm(-1) [Dauster et al., Phys. Chem. Chem. Phys., 2008, 10, 2827]. The IR C-H stretch frequencies and intensities of the Bz subunit are very similar to those of the acceptor in the (Bz)(2) dimer, confirming that Bz acts as the acceptor. While the S(1) <- S(0) electronic origin of Bz is forbidden and is not observable in the gas-phase, the UV spectrum of Pyr center dot Bz in the same region exhibits a weak 0(0)(0) band that is red-shifted by 58 cm(-1) relative to that of Bz (38 086 cm(-1)). The origin appears due to symmetry-breaking of the p-electron system of Bz by the asymmetric pyrrole NH center dot center dot center dot pi hydrogen bond. This contrasts with (Bz)(2), which does not exhibit a 0(0)(0) band. The Bz moiety in Pyr center dot Bz exhibits a 6a(0)(1) band at 0(0)(0) + 518 cm(-1) that is about 20x more intense than the origin band. The symmetry breaking by the NH center dot center dot center dot pi hydrogen bond splits the degeneracy of the v(6)(e(2g)) vibration, giving rise to 6a' and 6b' sub-bands that are spaced by similar to 6 cm(-1). Both the 0(0)(0) and 6(0)(1) bands of Pyr center dot Bz carry a progression in the low-frequency (10 cm(-1)) excited-state tipping vibration omega', in agreement with the change of the omega tipping angle predicted by SCS-MP2 and SCS-CC2 calculations.

Crystal chemistry and hydrogen bonding of rustumite Ca10(Si2O7)2(SiO4)(OH)2Cl2 with variable OH, Cl, F

Three samples of the skarn mineral rustumite Ca10(Si2O7)2(SiO4)(OH)2Cl2, space group C2/c, a ≈7.6, b ≈ 18.5, c ≈ 15.5 Å, β ≈ 104°, with variable OH, Cl, F content were investigated by electron microprobe, single-crystal X-ray structure refinements, and Raman spectroscopy. “Rust1LCl” is a low chlorine rustumite Ca10(Si2O7)2(SiO4)(OH1.88F0.12)(Cl1.28,OH0.72) from skarns associated with the Rize batholith near Ikizedere, Turkey. “Rust2F” is a F-bearing rustumite Ca10(Si2O7)2(SiO4)(OH1.13F0.87) (Cl1 96OH0.04) from xenoliths in ignimbrites of the Upper Chegem Caldera, Northern Caucasus, Russia. “Rust3LClF” represents a low-Cl, F-bearing rustumite Ca10(Si2O7)2(SiO4)0.87(H4O4)0.13(OH1.01F0.99) (Cl1.00 OH1.00) from altered merwinite skarns of the Birkhin massif, Baikal Lake area, Eastern Siberia, Russia. Rustumite from Birkhin massif is characterized by a significant hydrogarnet-like or fluorine substitution at the apices of the orthosilicate group, leading to specific atomic displacements. The crystal structures including hydrogen positions have been refined from single-crystal X-ray data to R1 = 0.0205 (Rust1_LCl), R1 = 0.0295 (Rust2_F), and R1 = 0.0243 (Rust3_LCl_F), respectively. Depletion in Cl and replacement by OH is associated with smaller unit-cell dimensions. The substitution of OH by F leads to shorter hydrogen bonds O-H⋯F instead of O-H⋯OH. Raman spectra for all samples have been measured and confirm slight strengthening of the hydrogen bonds with uptake of F.This study discusses the complex crystal chemistry of the skarn mineral rustumite and may provide a wider understanding of the chemical reactions related to contact metamorphism of limestones.

Three strategically placed hydrogen-bonding residues convert a proton pump into a sensory receptor.

In haloarchaea, light-driven ion transporters have been modified by evolution to produce sensory receptors that relay light signals to transducer proteins controlling motility behavior. The proton pump bacteriorhodopsin and the phototaxis receptor sensory rhodopsin II (SRII) differ by 74% of their residues, with nearly all conserved residues within the photoreactive retinal-binding pocket in the membrane-embedded center of the proteins. Here, we show that three residues in bacteriorhodopsin replaced by the corresponding residues in SRII enable bacteriorhodopsin to efficiently relay the retinal photoisomerization signal to the SRII integral membrane transducer (HtrII) and induce robust phototaxis responses. A single replacement (Ala-215-Thr), bridging the retinal and the membrane-embedded surface, confers weak phototaxis signaling activity, and the additional two (surface substitutions Pro-200-Thr and Val-210-Tyr), expected to align bacteriorhodopsin and HtrII in similar juxtaposition as SRII and HtrII, greatly enhance the signaling. In SRII, the three residues form a chain of hydrogen bonds from the retinal's photoisomerized C(13)=C(14) double bond to residues in the membrane-embedded alpha-helices of HtrII. The results suggest a chemical mechanism for signaling that entails initial storage of energy of photoisomerization in SRII's hydrogen bond between Tyr-174, which is in contact with the retinal, and Thr-204, which borders residues on the SRII surface in contact with HtrII, followed by transfer of this chemical energy to drive structural transitions in the transducer helices. The results demonstrate that evolution accomplished an elegant but simple conversion: The essential differences between transport and signaling proteins in the rhodopsin family are far less than previously imagined.