Competition between hydrogen bonds and halogen bonds in complexes of formamidine and hypohalous acids

Quantum chemical calculations have been per-formed for the complexes of formamidine (FA) and hypohalous acid (HOX, X = F, Cl, Br, I) to study their structures, properties, and competition of hydrogen bonds with halogen bonds. Two types of complexes are formed mainly through a hydrogen bond and a halogen bond, respectively, and the cyclic structure is more stable. For the F, Cl, and Br complexes, the hydrogen-bonded one is more stable than the halogen-bonded one, while the halogen-bonded structure is favorable for the I complexes. The associated H-O and X-O bonds are elongated and exhibit a red shift, whereas the distant ones are contracted and display a blue shift. The strength of hydrogen and halogen bonds is affected by F and Li substitutents and it was found that the latter tends to smooth differences in the strength of both types of interactions. The structures, properties, and interaction nature in these complexes have been understood with natural bond orbital (NBO) and atoms in molecules (AIM) theories.     

Effect of superalkali substituents on the strengths and properties of hydrogen and halogen bonds

Quantum chemical calculations have been performed for the complexes Li₃OCCX–Y (X?=?Cl, Br, H; Y?=?NH₃, H₂O, H₂S) and Li₃OCN–X′Y′ (X′Y′?=?ClF, BrCl, BrF, HF) to study the role of superalkalis in hydrogen and halogen bonds. The results show that the presence of an Li₃O cluster in a Lewis acid weakens its acidity, while its presence in a Lewis base enhances its basicity. Furthermore, the latter effect is more prominent than the former one, and the presence of an Na₃O cluster causes an even greater effect than Li₃O. The strengths of hydrogen and halogen bonds were analyzed using molecular electrostatic potentials. The contributions of superalkalis to the strength of hydrogen and halogen bonds were elucidated by analyzing differences in electron density.     

Cooperativity between fluorine-centered halogen bonds: investigation of substituent effects

This article analyzes the substitution effects on cooperativity between fluorin-centered halogen bonds in NCF?·?·?·?NCF?·?·?·?NCX and CNF?·?·?·?CNF?·?·?·?CNX complexes, where X?=?H, F, Cl, CN, OH, and NH₂. These effects are investigated theoretically in terms of geometric and energetic features of the complexes, which are computed by ab initio methods. The topological analysis, based on the quantum theory of atoms in molecules (QTAIM), is used to characterize the interactions and analyze their enhancement with varying electron density at bond critical points. It is found that the complexes with electron-donating groups exhibit a strong cooperativity, while a much weaker cooperativity occurs in the NCF?·?·?·?NCF?·?·?·?NCCN and CNF?·?·?·?CNF?·?·?·?CNCN trimers. An excellent correlation is found between the cooperative energy in the ternary complexes and the calculated three-body interaction energies. The energy decomposition analysis (EDA) indicates that the electrostatic and dispersion effects play a main role in the cooperativity of fluorine-centered halogen bonding.

The interplay of hydrogen bonds and halogen bonds in the structure of NH-pyrazoles bearing C-aryl and C-halogen substituents

The behavior in solution and in the solid state of 3(5)-phenyl-1H-pyrazole (7), 3(5)-phenyl-4-chloro-1H-pyrazole (6), 3(5)-phenyl-4-bromo-1H-pyrazole (1), and 3(5)-p-chlorophenyl-4-bromo-1H-pyrazole (8) is discussed in relation to their 3-phenyl (a)/5-phenyl (b) annular tautomerism. Two new X-ray structures are reported: a new polymorph of 1 and the structure of 6. The new polymorph is a 3-phenyl-1H-pyrazole 1a′ trimer while the new structure is a 5-phenyl-1H-pyrazole 6b trimer. The combined use of NMR at low temperature and DFT calculations allows to discuss the tautomerism of the first three pyrazoles and to predict that the fourth one should be a tetramer formed by both tautomers, 8a and 8b.     

Topological analysis of aromatic halogen/hydrogen bonds by electron charge density and electrostatic potentials

In this work, the intermolecular distribution of the electronic charge density in the aromatic hydrogen/halogen bonds is studied within the framework of the atoms in molecules (AIM) theory and the molecular electrostatic potentials (MEP) analysis. The study is carried out in nine complexes formed between benzene and simple lineal molecules, where hydrogen, fluorine and chlorine atoms act as bridge atoms. All the results are obtained at MP2 level theory using cc-pVTZ basis set. Attention is focused on topological features observed at the intermolecular region such as bond, ring and cage critical points of the electron density, as well as the bond path, the gradient of the density maps, molecular graphs and interatomic surfaces. The strength of the interaction increases in the following order: F⋅⋅⋅π < Cl⋅⋅⋅π < H⋅⋅⋅π. Our results show that the fluorine atom has the capability to interact with the π−cloud to form an aromatic halogen bond, as long as the donor group is highly electron withdrawing. The Laplacian topology allows us to state that the halogen atoms can act as nucleophiles as well as electrophiles, showing clearly their dual character.     

Uniformity,ideality, and hydrogen bonds in transmembrane alpha-helices

Protein environments substantially influence the balance of molecular interactions that generate structural stability. Transmembrane helices exist in the relatively uniform low dielectric interstices of the lipid bilayer, largely devoid of water and with a very hydrophobic distribution of amino acid residues. Here, through an analysis of bacteriorhodopsin crystal structures and the transmembrane helix structure from M2 protein of influenza A, some helices are shown to be exceptionally uniform in hydrogen bond geometry, peptide plane tilt angle, and backbone torsion angles. Evidence from both the x-ray crystal structures and solid-state NMR structure suggests that the intramolecular backbone hydrogen bonds are shorter than their counterparts in water-soluble proteins. Moreover, the geometry is consistent with a dominance of electrostatic versus covalent contributions to these bonds. A comparison of structure as a function of resolution shows that as the structures become better characterized the helices become much more uniform, suggesting that there is a possibility that many more uniform helices will be observed, even among the moderate resolution membrane protein structures that are currently in the Protein Data Bank that do not show such features.     

Geometric criteria of hydrogen bonds in proteins and identification of "bifurcated" hydrogen bonds

Empirical criteria for identification of hydrogen bonds were analyzed to produce a set of geometrically consistent criteria. For a data set of 30 structures, application of a set of purely geometrical criteria, along with exclusion of abnormal backbone conformations, also excluded a common interaction of Ser/Thr side chains with Asp/Glu side chains ([ST]/[DE] pairs). These interactions were termed "bifurcated hydrogen bonds", which implies delocalization of a positively charged hydrogen of hydroxyl between the two acceptor atoms of the carboxylic group. These "bifurcated" interactions are among the most common packing patterns for [ST]/[DE] pairs of side chains. Therefore, the identification of hydrogen bonds cannot be based on geometrical criteria only and requires introduction of some physico-chemical criteria.     

Short hydrogen bonds in proteins

Short hydrogen bonds are present in many chemical and biological systems. It is well known that these short hydrogen bonds are found in the active site of enzymes and aid enzyme catalysis. This study aims to systematically characterize all short hydrogen bonds from a nonredundant dataset of protein structures. The study has revealed that short hydrogen bonds are commonly found in proteins and are widely present in different regions of the protein chain, such as the backbone or side chain, and in different secondary structural regions such as helices, strands and turns. The frequency of occurrence of donors and acceptors from the charged side chains as well as from the neutral backbone atoms is equally high. This suggests that short hydrogen bonds in proteins occur either due to increased strength or due to geometrical constraints and this has been illustrated from several examples.     

Modelling hydrogen bonds in NN-dimethylformamide

N, N-dimethylformamide (DMF) is a ‘universal’ solvent with the simplest amide structure. DMF has different interactions with many polymers and biomolecules. It is therefore necessary to study systematically the interactions in DMF itself first. In this study, both FT-IR and two molecular theoretical methods (MP2 and DFT/B3LYP) were used to study various hydrogen bonding interactions in DMF molecules based on its weak H-bonding donors CH/CH₃ and strong H-bonding acceptor C = O. The possible H-bonding donors and acceptors in DMF molecules were first analysed followed by modelling the effect of different structural environments on vC = O bands in infrared spectra. Finally, H-bonding properties including distance, angles and the energy as well as the probability of H-bonding patterns were obtained. The results showed that there exist five possible different weak types of H-bonding dimers; among them, three dimers consist of a pair of weak H-bonds, whereas two other dimers have two pairs of H-bonds, leading to 14 (including eight different) H-bonds. Two types of dimers were dominant, whereas three others can be omitted.     

Influence of hydrogen bonds on edge-to-face interactions between pyridine molecules

Edge-to-face interactions between two pyridine molecules and the influence of simultaneous hydrogen bonding of one or both of the pyridines to water on those interactions were studied by analyzing data from ab initio calculations. The results show that the edge-to-face interactions of pyridine dimers that are hydrogen bonded to water are generally stronger than those of non-H-bonded pyridine dimers, especially when the donor pyridine forms a hydrogen bond. The binding energy of the most stable edge-to-face interacting H-bonded pyridine dimer is ?5.05 kcal/mol, while that for the most stable edge-to-face interacting non-H-bonded pyridine dimer is ?3.64 kcal/mol. The interaction energy data obtained in this study cannot be explained solely by the differences in electrostatic potential between pyridine and the pyridine–water dimer. However, the calculated cooperative effect can be predicted using electrostatic potential maps.     

Structurally symmetrical,easily polarizable hydrogen bonds between side chains in proteins and proton-conducting mechanisms. III

(L -Cys)_n, (L -Lys)_n, and (L -Glu)_n were studied by ir spectroscopy in terms of their degree of deprotonation or protonation. It is shown that structurally symmetrical, easily polarizable SH ?S^? ? ^?S ?HS, N⁺H ?N ? N ?H⁺N, and OH ?O^? ? ^?O ?HO hydrogen bonds are formed between the side chains. The different wave number distributions of the ir continua caused by these hydrogen bonds show that the barrier in the double-minimum proton potential decreases in the series of these hydrogen bonds. The stability of these hydrogen bonds against hydration increases in this series. The OH ?O^? ? ^?O ?HO bonds are not broken by small amounts of water. With (L -Cys)_n the formation of easily polarizable hydrogen bonds and a β-structure–coil transition are strongly interdependent. As a result of this coupling effect, the β-structure–coil transition becomes cooperative. With (L -Glu)_n, the formation of the polarizable hydrogen bonds and the observed conformational change are independent processes. The (L -Glu)_n conformation changes from α-helix to coil only if more than 80% of the residues are deprotonated. Finally, on the basis of the various types of easily polarizable hydrogen bonds, charge shifts in active centers of enzymes and the proton-conducting mechanism through hydrophobic regions of biological membranes are discussed.     

Modeling hydrogen bonds in three dimensions

The effective potential between two hydrogen bonded atoms is calculated on the basis of the Lippencott-Schroeder bent bond model, taken to be a typical model interaction. We differ from other calculations in that the minimum energy configuration for the proton is treated adiabatically, its position being recomputed at each value of the larger atoms separation. We find the typical hard core to have been a consequence of an artificial restriction of the proton to a fixed angle with the larger atom axis, basically a one-dimensional assumption. Free to move in three dimensions, the proton is squeezed off the axis as the separation narrows, and the hard core feature is gone. Depending on the degree of bond bending, the anharmonicity of the bond may be diminished, eliminated, or even reversed.     

Proton transfer and polarizability of hydrogen bonds formed between cysteine and lysine residues

(L -Cys)_n + N-base systems and (L -Cys)_n + (L -Lys)_n systems were studied by ir spectroscopy. It is shown that in the water-free systems, SH ?N ? S^? ?H⁺N hydrogen bonds are formed. With the (L -Cys)_n + N-base systems, both proton-limiting structures in the SH ?N ? S^? ?H⁺N bonds have equal weight when the pK_a of the protonated N-base is 2 pK_a units larger than that of (L -Cys)_n. The same is true with the water-free (L -Cys)_n + (L -Lys)_n system. Thus, with regard to the type of proton potentials present, these hydrogen bonds are proton-transfer hydrogen bonds showing very large proton polarizabilities. This is confirmed by the occurrence of continua in the ir spectra. Small amounts of water open these hydrogen bonds and increase the transfer of the proton to (L -Lys)_n. In the (L -Lys)_n + N-base systems, with increasing proton transfer the backbone of (L -Cys)_n changes from antiparallel β-structure to coil. In (L -Cys)_n + (L -Lys)_n, the conformation is determined by the (L -Lys)_n conformation and changes depending on the chain length of (L -Lys)_n. Finally, the reactivity increase in the active center of fatty acid synthetase, which should be caused by the shift of a proton, is discussed on the basis of the great proton polarizability of the cysteine–lysine hydrogen bonds.     

Co-crystallization between a thymine and a metal complex connected by triple hydrogen bonds

The synthesis and crystal structure of co-crystal between bis(dithiobiureto)platinum(II) ([Pt(dtb)₂]) with thymine is reported. The crystalline structure of [Pt(dtb)₂](thymine)₂ shows two-dimensional arrays created by hydrogen-bonding interactions. One mononuclear complex and two thymine molecules form a building unit connected by triple hydrogen bond employing ADA-DAD arrangements (A, hydrogen bond acceptor; D, hydrogen bond donor). The building unit is linked to adjacent units via additional hydrogen bonds to form planar sheet.     

Microsomal oxidation of bromo-, chloro- and fluorobiphenyls

1. The metabolism of 2-, 3-, 4-bromo-, 2-, 4-chloro-, and 2-fluorobiphenyl by hepatic microsomes isolated from control and Aroclor 1254-treated rats and pigeons was studied.2. Meta and para as well as dihydroxylated metabolites were detected, but para hydroxylation was the preferred route of metabolism with all of the substrates used.3. The overall rates of hydroxylation were greater with hepatic microsomes from rats than from pigeons.4. Treatment with Aroclor 1254, a potent inducer of hepatic monooygenases, resulted in increased rates of metabolism and in the enhanced formation of diol metabolites. Metabolism of halobiphenyls by induced P450 isoenzymes altered the regioselective hydroxylation pathways.5. Ortho- and meta halosubstituted biphenyls were less rapidly metabolised when compared with paru substituted isomers.     

Counting peptide-water hydrogen bonds in unfolded proteins

It is often assumed that the peptide backbone forms a substantial number of additional hydrogen bonds when a protein unfolds. We challenge that assumption in this article. Early surveys of hydrogen bonding in proteins of known structure typically found that most, but not all, backbone polar groups are satisfied, either by intramolecular partners or by water. When the protein is folded, these groups form approximately two hydrogen bonds per peptide unit, one donor or acceptor for each carbonyl oxygen or amide hydrogen, respectively. But when unfolded, the backbone chain is often believed to form three hydrogen bonds per peptide unit, one partner for each oxygen lone pair or amide hydrogen. This assumption is based on the properties of small model compounds, like N-methylacetamide, or simply accepted as self-evident fact. If valid, a chain of N residues would have approximately 2N backbone hydrogen bonds when folded but 3N backbone hydrogen bonds when unfolded, a sufficient difference to overshadow any uncertainties involved in calculating these per-residue averages. Here, we use exhaustive conformational sampling to monitor the number of H-bonds in a statistically adequate population of blocked polyalanyl-six-mers as the solvent quality ranges from good to poor. Solvent quality is represented by a scalar parameter used to Boltzmann-weight the population energy. Recent experimental studies show that a repeating (Gly-Ser) polypeptide undergoes a denaturant-induced expansion accompanied by breaking intramolecular peptide H-bonds. Results from our simulations augment this experimental finding by showing that the number of H-bonds is approximately conserved during such expansion⇋compaction transitions.